Testing the Prediction of Fuzzy Dark Matter Theory in the Milky Way Center

The fuzzy dark matter model (FDM; also known as quantum wave dark matter model) argues that light bosons with a mass of are a possible candidate for dark matter in the universe. One of the most important predictions of FDM is the formation of a soliton core instead of a density cusp at the center of galaxies. If FDM is the correct theory of dark matter, then the predicted soliton core can help form the Central Molecular Zone (CMZ) in the Milky Way. We present high-resolution hydrodynamical simulations of gas flow patterns to constrain the properties of the soliton core based on a realistic Milky Way potential. We find that a dense center is required to form a reasonable CMZ. The size and kinematics of the CMZ offer a relatively strong constraint on the inner enclosed mass profile of the Galaxy. If a soliton core is not considered, a compact nuclear bulge alone with a radially varying mass-to-light ratio can match the observed size and kinematics of the CMZ. A soliton core model with a mass of and a core radius of , together with a less massive nuclear bulge with a constant mass-to-light ratio, also agrees nicely with the current data. Such an FDM soliton core corresponds to a boson mass of , which could be further constrained by the improved determination of the mass-to-light ratio in the Galactic center.


Introduction
Although the current cold dark matter (CDM) model successfully explains many problems on the large-scale structures of the universe (e.g., Frenk & White 2012;Bennett et al. 2013), some of its predictions on the galactic scales are still in tension with modern observations (see the review by Bullock & Boylan-Kolchin 2017). One of the famous discrepancies is the so-called "cusp-core" problem: dissipationless CDM simulations predicted a universal Navarro-Frenk-White (NFW) density profile for bound dark matter halos, with a cuspy density profile at the central part (Navarro et al. 1997), while observations favor a flat cored profile in low surface brightness galaxies (LSB) and dwarf galaxies (de Blok 2010). A shallower dark matter core might be created by considering the various baryonic effects in simulations (e.g., Oñorbe et al. 2015;Read et al. 2019). However, it is still not clear how these "subgrid physics" alter the structure formation in the lowest mass galaxies (e.g., Sánchez et al. 2012;Weinberg et al. 2015).
On the other side, alternative dark matter models are being developed to solve the problems in a more self-consistent way. In recent years, the Fuzzy Dark Matter (FDM) model, also known as the "quantum wave dark matter" model, has been gaining more attention. In this scenario, the dark matter is composed of very light bosons with a particle mass of~-10 eV 22 (e.g., Hu et al. 2000;Marsh & Silk 2014;Schive et al. 2014a). Such bosons have a characteristic wavelength of mv 0.2 kpc assuming a velocity of -100 km s 1 , which helps to suppress the formation of small-scale structures. One of the most important predictions of FDM is the presence of a stable ground-state soliton core due to the Bose-Einstein condensate at the central part of each dark matter halo. The core radius is usually comparable to the characteristic wavelength, and the halo transitions to an NFW profile within a few core radii according to recent numerical simulations (e.g., Schive et al. 2014b). This may provide a plausible solution to the "cusp-core" problem in the CDM model. One can easily deduce that a lower mass halo would host a larger soliton core due to the particles inside it having a smaller velocity and a larger characteristic wavelength. Therefore, previous works (e.g., Chen et al. 2017) preferred to use dwarf galaxies to constrain the soliton core mass and size, and a boson mass of ( -)´-1 2 10 eV 22 is favored in such studies. Naturally, one may wonder whether the soliton core exists in our Milky Way and if so, how we can detect its effects on the stars and gas. The Milky Way should be a unique laboratory for testing the existence of such a soliton core due to the unprecedented data quality we could achieve compared to external galaxies. However, one should be careful when applying the halo-core scaling relation found in previous FDM simulations (e.g., Schive et al. 2014b) to Milky Way-like galaxies to derive the soliton core properties. One reason is that previous FDM simulations focus mainly on the halo mass range of - M 10 10 halos, such as the Milky Way with a virial halo mass of  »´M 1.5 10 12 (Portail et al. 2017;Posti & Helmi 2019). Another reason is that the halo-core relation could be interpreted as the specific kinetic energy of the soliton core equals that of the halo. However, the specific kinetic energy of the soliton core can be strongly affected by the baryons in the Galactic center, and therefore, the halo-core scaling relation may need to be adjusted. (Bar et al. 2019, see also Section 2.1.3).
Nevertheless, according to the scaling relation in Schive et al. (2014b), the expected soliton mass and core radius are  M 1.44 10 9 and 0.16 kpc for the Milky Way, respectively, by assuming a boson mass of´-1.0 10 eV 22 . It is also important to note that these two numbers should be interpreted as the maximum soliton core mass and the minimum soliton core radius, due to the fact that the halo of the Milky Way may still be relaxing (Hui et al. 2017).
Observations on stellar and gas kinematics in the Galactic center have revealed some evidence for a central component with a mass of  M 10 9 , although it is not clear whether this component is composed of baryons or dark matter (see the discussions in Section V.B of Bar et al. 2018). Portail et al. (2017) and De Martino et al. (2018) included such a component in their dynamical models to explain the high velocity dispersion of the stars in the Galactic bulge. Moreover, the CMZ (Bally et al. 1988), which is thought to be a molecular gas ring or disk orbiting the Galactic center, has a size of »100 pc (Galactic longitude | |   l 1 . 5) with a rotation velocity of » -100 km s 1 (e.g., Molinari et al. 2011;Kruijssen et al. 2015;Henshaw et al. 2016). This implies an enclosed mass of  »´M 2.3 10 8 inside »100 pc by assuming circular motions. The formation of the CMZ is related to the bar-induced gas inflow due to the large-scale bisymmetric gravitational torques (e.g., Sormani et al. 2015a;Li et al. 2016;Ridley et al. 2017). In addition, the CMZ needs a dense center (~1% disk mass) to support its backbone x 2 orbits (Regan & Teuben 2003;Li et al. 2017;Sormani et al. 2018a). Previous studies usually included a very compact stellar nuclear bulge to form the CMZ (e.g., Baba et al. 2010;Li et al. 2016;Shin et al. 2017;Sormani et al. 2018a;Armillotta et al. 2019), but the physical origin of this component is still not clear. The mass and size of the nuclear bulge are also not well constrained, due both to dust extinction on the Galactic plane and to the uncertainties on the mass-to-light ratio of the nuclear bulge (see Section 2.1.2). This gives us some freedom to add a soliton core to our previous Milky Way gas-dynamical model (Li et al. 2016) to see whether we could still reproduce the observed CMZ shape and kinematics, thus putting important constraints on the FDM theory.
The purpose of the present paper is to use high-resolution hydrodynamical simulations to determine the allowed mass and size of the soliton core in the Milky Way. The paper is organized as follows: in Section 2, we describe our model setup and numerical methods; in Section 3, we present our simulation results; and in Section 4, we discuss the possibility of the existence of such a soliton core together with its additional dynamical effects. The summary is presented in Section 5. The Galactic potential used in this work is based on the bestfit gas-dynamical model of Li et al. (2016), where they used the features in the (l v , ) diagram to constrain the gas flow patterns in the Milky Way on large scales. The (l v , ) diagram shows the intensity of cold gas emission as a function of Galactic longitude l and line-of-sight velocity v (e.g., Dame et al. 2001). The gaseous features viewed from face-on appear as highdensity ridges in the (l v , ) diagram, and these features are formed mainly due to the non-axisymmetric structures such as the Galactic stellar bar and spiral arms.

Model
The large-scale potential of the Milky Way has been studied extensively in recent works. Portail et al. (2015) constructed a made-to-measure Milky Way model that matches the 3D density of red clump giants in the bulge region (Wegg & Gerhard 2013) and the BRAVA kinematics (Kunder et al. 2012) quite well. The axial ratios of the bar are (10:6.3:2.6) with the semimajor axis of »5 kpc. The bar angle to the Sun-Galactic center line is ( )    27 2 . We adopt this model as the basis of the Galactic potential with a few minor modifications as in Li et al. (2016). The bar pattern speed is chosen to be , which places the corotation radius at = R 6.2 kpc. In this study, we remove the large-scale spiral potentials from Li et al. (2016) as they are only important for the outer gas features, which is not the main focus of this work.
On the other hand, the central mass distribution in our Galaxy is not well constrained. The stellar kinematics may only be reliable at higher latitudes (   b 2 ), due to the heavy dust contamination closer to the disk plane. The gas kinematics in the central disk region may not serve as a good indicator of mass either, due to its highly noncircular motions caused by the Galactic bar (e.g., Fux 1999;Sormani et al. 2015bSormani et al. , 2018b. We therefore modify the central potential of the best-fit model in Li et al. (2016) to test the possible effects caused by a soliton core. We include the following two components, namely the nuclear bulge and the soliton core, in the central part of the Milky Way potential. Launhardt et al. (2002) found hints of a stellar nuclear bulge component in the central~200 pc of the Milky Way using COBE near-infrared image. The authors parameterized the nuclear bulge with the following equation: , to fit the data. The authors claimed that the presence of two distinct components may not have a clear physical meaning, but could reflect the extinction effects not accounted for in their model. By assuming ¡ = 2, Launhardt et al. (2002) derived the total mass of the nuclear bulge to be ( )  ´M 1.4 0.6 10 9 . Following Launhardt et al. (2002), the nuclear bulge included in our simulations is also composed of two components with = R 120 pc nb,1 and = R 220 pc nb,2 . Their density profiles are described by Equation (1). We generally assume the mass-to-light ratio ϒ is the same for the two components. The total mass of the nuclear bulge is therefore As the nuclear bulge is a massive and compact object, its dynamical effects cannot be negligible. Other studies have also revealed the presence of the nuclear bulge (or at least a massive and dense central object), including the young stellar nuclear disk found in APOGEE by Schönrich et al. (2015), and the increase of stellar velocity dispersion caused by a massive central object in Portail et al. (2017) and De Martino et al. (2018). Li et al. (2016) also included the nuclear bulge to generate a reasonable CMZ in their best-fit gas-dynamical models. However, the exact mass of the nuclear bulge is still not well determined. This is partly because the COBE photometry used to derive the density profile of the nuclear bulge is heavily affected by the dust extinction. Another reason is that the assumed ¡ = 2 in the near-infrared band is uncomfortably large (e.g., Conroy et al. 2009;Conroy & Gunn 2010). We therefore keep ϒ as a free parameter in this study. , and r c is the core radius defined as the radius enclosing 25% of the total soliton mass M s (r c is also roughly the radius where density decreases to half of the central peak value). However, as pointed out in Section II.C from Hui et al. (2017), these two numbers should be interpreted only as the upper and lower limits, i.e., the maximum soliton mass and the minimum soliton core radius as the Milky Way halo may have not been fully relaxed and the soliton core may be developing. We therefore define . Note that this mass is very close to the mass of the nuclear bulge with ¡ = 2.

Nuclear Bulge
The mass M s and core radius r c of self-gravitating soliton cores are inversely related, because the characteristic de Broglie wavelength cannot exceed the virial radius of the system (Hui et al. 2017; also see Equation (5) in Section 3.5). This relation needs to be modified for the Milky Way, due to the presence of the nuclear bulge whose mass and size are similar to the expected soliton core predicted by the halo-core relation. Recently, Bar et al. (2019) studied how the baryonic nuclear bulge would affect the soliton core size in the Milky Way. They found that the soliton core radius r c could be compressed to half of its original value in the potential of a nuclear bulge with a mass profile similar to Launhardt et al. (2002). Inspired by their work, we adopt a modified soliton core density profile based on the results in Schive et al. (2014b) as where M s and r c are independent. We could then use Equation (2) to make soliton cores with the same mass but different core radius. The corresponding gravitational acceleration of such a soliton core is given by where 1 . Assuming M s and r c are independent implies the mass profile of the nuclear bulge varies accordingly. This parameterization is more flexible to obtain the required mass profile of the soliton core, while the method in Bar et al. (2019) is a more self-consistent way of determining the soliton core properties by solving the Schrödinger-Poisson equation. In this study, we first adopt M s and r c as two free soliton parameters. We then examine whether the preferred soliton core properties (i.e., certain combinations of M s and r c ) are allowed by the nuclear bulge potential we adopt in Section 3.5.
If M s and r c are chosen to be the same as the predictions of the halo-core relation (Schive et al. 2014b), one gets a selfgravitating soliton core. If M s is the same but a smaller r c is adopted, one gets a compressed soliton core under an external baryonic potential (i.e., the nuclear bulge). In Figure 1. we show the rotation curves of a self-gravitating soliton core and a compressed soliton core with the same total mass obtained using Equation (  , the same as in Bar et al. (2019). The blue solid line shows the rotation curve of a self-gravitating soliton core with a core radius = r 0.17 kpc c , while the red dotted-dashed line shows a soliton core that is compressed by a factor of 2. The compressed soliton core has the same total mass but a smaller core radius = r 0.09 kpc c . The green dotted-dashed curve shows the results from Bar et al. (2019), who obtained the compressed soliton core profile by directly solving the Schrödinger-Poisson equation in an external nuclear bulge potential. Their results are in good agreement with Equation (3).

Numerical Scheme
We study how a thin gas disk evolves under a realistic barred Milky Way potential described in the previous sections. The barred potential is rigidly rotating about the Galactic center with a fixed pattern speed W = -- . The gas is assumed to be inviscid, and we use the latest public version of the grid-based MHD code Athena++ (J. Stone et al. 2020, in preparation) to solve the Euler equations. Athena++ is a complete rewrite of the previous Athena code with improved performance and scalability. We employ a uniform Cartesian grid with 4096×4096 cells for 2D simulations covering a square box with size of = L 14 kpc, thus the grid spacing is D = D = x y 3.4 pc. We also test 3D runs where the grid iś2 048 2048 21 with a vertical extent of 200 pc. We adopt the Roe Riemann solver, piecewise linear reconstructions, and outflow boundary conditions at the domain boundaries. High spatial resolution is necessary to capture the turbulent gas motions in the central region of galaxies, as demonstrated in previous works (e.g., Li et al. 2015;Sormani et al. 2015a;Few et al. 2016). All our models are run for a period of 0.5 Gyr. The 2D simulations take about ∼15-20 hr using 64 Intel cores, and the 3D runs take about ∼5 days using 256 Intel cores. In general, Athena++ is ∼5 times faster than the previous version Athena (Gardiner & Stone 2005;Stone et al. 2008;Stone & Gardiner 2009).
We set up the initial rotating gaseous disk with an exponential surface density profile The gas disk then has a surface density around pc 2 at solar neighborhood to match the observed value roughly (Bovy & Rix 2013). The initial rotation velocity of gas is set to balance the azimuthally averaged gravitational force. We linearly ramp up the bisymmetric bar potential over the first bar rotation period = T 186 Myr bar to avoid transients. We compute the gas flow in the inertial frame, while Li et al. (2016) computed it in the bar corotation frame.
We adopt an isothermal equation of state assuming the specific gas internal energy is constant through an energy balance between radiative heating and cooling processes, same as previous studies (e.g., Kim et al. 2012aKim et al. , 2012bSormani et al. 2015c). The effective isothermal sound speed c s describes the level of velocity dispersion between molecular clouds and thus does not stand for the microscopic gas temperature. We test two different c s in our models with low ( -10 km s 1 ) and high ( -20 km s 1 ) cases, in the range of the observed ISM velocity dispersion (e.g., Walter et al. 2008;Leroy et al. 2009). We ignore gas self-gravity, magnetic fields, star formation and stellar feedback, and other additional physics in our models for simplicity, but prefer to study these effects in a follow-up paper. We discuss how these unaccounted physics will affect our results in Section 4.4.

General Evolution
We first report the simulation of the gas evolution in the potential described in Section 2.1 with a ¡ = 2 nuclear bulge, but without a soliton core. We explore the case with a relatively low sound speed = c 10 km s s 1 , which is the same as most of the previous studies (e.g., Fux 1999;Rodriguez-Fernandez & Combes 2008;Ridley et al. 2017). In Figure 2, we show the gas surface density in the simulation at = t 380 Myr when the bar has finished about two rotation periods. Strong shocks develop at the leading side of the bar as the flow is supersonic, which are indicated by the high-density gas ridges in Figure 2. Gas loses energy and angular momentum when encountering the shocks, then flows inwards and accumulates at the center. The CMZ can be clearly seen as the high-density ring-like structure at » R 300 pc in the center. The four spiral arms around the bar region are known as the near and far 3 kpc arms and the molecular ring in the (l v , ) diagram, which are related to the 4:1 resonance of the rotating bar (Sormani et al. 2015b). The flow pattern then becomes quasi-steady until the end of the simulation.
In general, this model is very similar to the best-fit model of Li et al. (2016), as we adopt the same Galactic potential but with a higher resolution. In the following sections, we vary the mass-to-light ratio ϒ of the nuclear bulge, the mass of the soliton core M s , and the radius of the soliton core r c to determine the possible mass and size of the soliton core allowed in the Milky Way, thus putting constraints on FDM theory.

Tests with Only the Nuclear Bulge
As we have shown in Section 2.1.2, an assumed a mass-tolight ratio ¡ = 2 gives a nuclear bulge mass of  M 1.4 10 9 , which is very close to the maximum mass of the soliton core . However, the total mass of the central components is constrained to be »1.5  M 10 9 (e.g., Launhardt et al. 2002;Portail et al. 2017;De Martino et al. 2018). It is therefore not reasonable that the Milky Way hosts both a maximum mass soliton core and a ¡ = 2 nuclear bulge at the same time. The mass degeneracy of a soliton core and a stellar nuclear bulge in the Milky Way center was also discussed in Bar et al. (2018). As the ϒ of the nuclear bulge may have some uncertainties, we first examine its effects on gas features when a possible soliton core is not included. . The bar lies horizontally along the x-axis with a semimajor length around »5 kpc, and the corotation radius of the bar is placed at = R 6.2 kpc. The CMZ can be seen as the high-density gas ring around » R 300 pc. The central potential in this model includes a nuclear bulge with ¡ = 2 but without a soliton core. The Sun is at ( ) ( ) = -x y , 7.4 kpc, 3.8 kpc in this plot. The Sun-Galactic center line (l = 0) is indicated by the black dashed line.
In Figure 3, we plot the gas surface density (top panels), the rotation curves in the central 500 pc with the position of the CMZ marked (insets), and the corresponding (l v , ) diagram (bottom panels) of the gas in the central 500 pc of our simulations with different ϒ. The colors in the (l v , ) diagram indicate gas density, which is binned according to Galactic longitude and line-of-sight velocity. We only select the dense gas ( C II mainly traces the hot ionized gas, while NH 3 originates from cold and dense molecular gas. It seems that gas with different physical conditions in the CMZ forms a consistent shape in the (l, v) space.
To form a CMZ that has similar size and kinematics to observations, the model needs to satisfy (at least) two conditions: (1) the rotation curve in the model has a circular velocity of » -100 km s 1 at » R 100 pc; (2) the inflowing gas piles up on roughly circular orbits with » R 100 pc. In Figure 3, we show that the size of the CMZ shrinks with lower ϒ, due to a decreasing central mass (e.g., Athanassoula 1992; Regan & Teuben 2003;Li et al. 2015Li et al. , 2017. For ¡ = 0.5, a central mass of  M 3.6 10 8 is not sufficiently massive to produce a nearly circular CMZ with high gas density. The model with ¡ = 2 roughly satisfies condition (1) On the other hand, the model with ¡ = 1 roughly satisfies condition (2), but the rotation velocity is lower, and thus it forms a shallower slope in the (l v , ) diagram compared to observations. These results imply that varying the ϒ of the nuclear bulge alone with the density profile given by Equation (1) may not produce a reasonable CMZ that is consistent with observations in both size and kinematics.
We slightly fine-tune the model of the nuclear bulge above to see whether we could better reproduce the observed CMZ kinematics. We find that a less massive but more concentrated nuclear bugle is needed to satisfy conditions (1) and (2) simultaneously. Figure 3 ). This results in a nuclear bulge with a mass of  M 5.1 10 8 . The gas kinematics in Figure 3(d) agree with observations fairly well. We regard this model, which does not include a soliton core, as one of our best-fit models, although the nuclear bulge density profile in this model is slightly different from that in Launhardt et al. (2002). If the luminosity profile in Launhardt et al. (2002) is correct, then this best-fit model implies that the mass-to-light ratio of the nuclear bulge is radially varying.

Tests with Possible Soliton Cores
In Section 2.1.3, we argue that the soliton core (if it exists) could be significantly compressed by the presence of a nuclear bulge. The compressed soliton core would result in a steeper rotation curve at the inner region, which helps generate a steeper slope in the (l, v) space for the CMZ. In the following sections, we adopt the nuclear bulge profile as determined in  Kruijssen et al. (2015). The blue shaded region shows the [C II ] observations in CMZ from Langer et al. (2017). The bottom part of panel (c) is almost empty because there is little dense gas in this model that satisfies our CMZ criteria. Note that the Galactic longitude of 2°corresponds to a radius of »300 pc. Launhardt et al. (2002), where the mass-to-light ratio is the same for the two components. As shown in Section 3.2, such a nuclear bulge is not concentrated enough, thus leaving some room for a compressed soliton core. Figure 3(b) suggests a ¡ = 1 nuclear bulge already generates a slightly larger CMZ compared to observations. Adding a compact soliton core to such a nuclear bulge would only make the CMZ even larger. We thus adopt ¡ = 0.5 for the nuclear bulge so that we can produce a CMZ smaller than 200 pc by including a suitable soliton core. . The blue line at the inset of each panel denotes the rotation curve, and the vertical red line is the averaged CMZ radius, same as in Figure 3.
It is clear that the location of the CMZ is roughly the turnover radius of the rotation curve (i.e., where the rotation curve turns flat). For the detailed formation mechanism, we refer the reader to Li et al. (2015). The effects of a different soliton core mass M s and core radius r c are similar to those of a varying ϒ of the nuclear bulge. A larger M s leads to a larger CMZ, while a smaller r c produces a more circular CMZ together with a steeper (l, v) feature. Combining Figures 4 and 5, we conclude that cases (b), (e), (f), (i), and (j) can roughly satisfy conditions (1) and (2) mentioned in Section 3.2 simultaneously. The CMZ formed in these five simulations better reproduce the observed gas kinematics together with a size of »200 pc. These five models are highlighted with a red box in the figures and are regarded as our best-fit models that contain a soliton core.
The soliton cores with larger r c form more elliptical CMZs compared to models with smaller r c . The orientation of these more elliptical CMZs oscillates with time, which is not shown here but can be inferred in panels (g) and (k) in Figure 4. This is because that gas tends not to follow periodic orbits when the pressure force is comparable to the gravitational force. The oscillating CMZ is more or less similar to the open stream model proposed by Kruijssen et al. (2015), who argued that a precessing orbit instead of a closed orbit is a more physical solution in a Galactic potential. In our models, when these elliptical CMZs have a specific angle with respect to the Sun, the resulting (l, v) plot agrees better with observations (e.g., panel g), while other angles are worse (e.g., panel k). This can be understood through the following arguments: the Sun is at ( ) ( ) = -x y , 7.4 kpc, 3.8 kpc in our models (also see the dashed Sun-Galactic center line in Figure 6), if the line of sight is tangent to the elliptical CMZ near its pericenter (e.g., panel g), we would expect a steeper slope in the (l, v) space due to a higher line-of-sight velocity and a narrower extent compared to the case where the tangent point is near the apocenter (e.g., panel k). It is possible that the shape of the CMZ changes with time in reality, but it then implies we are observing the CMZ at a special time. If this is true, the constraints on the soliton core radius r c could be loosened to a slightly larger value.

Varying Sound Speed
Apart from the gravitational constraints, the gas in the CMZ is also quite turbulent. ALMA observations of SiO revealed that the velocity dispersion of gas can be as large as -20 km s 1 (Kauffmann 2016). The NH 3 data from the HOPS survey (Longmore et al. 2017) indicate a velocity dispersion of  -17.4 4.8 km s 1 . As described in Section 2.2, the effective isothermal sound speed c s used in our simulations reflects mainly the velocity dispersion of gas clouds. We therefore further test a higher value = c 20 km s s 1 for the preferred models highlighted by the red box in Figure 5 to see how this would affect our results.
We find that a higher sound speed c s generally makes the CMZ smaller, but only slightly. This is because shocks with a higher c s develop closer to the bar major axis. The physical reason is that the gas is less supersonic and thus needs a steeper potential to be shocked (e.g., Kim et al. 2012b;Sormani et al. 2015a). We also find that the CMZ tends to be more elliptical in a higher c s case, similar to the model with a large core radius r c but a lower c s . A higher c s also helps to broaden the (l, v) features as a result of a higher velocity dispersion, which better matches the data. Figure 6 shows the model that matches most observed (l, v) features. The model hosts a nuclear bulge with ¡ = 0.5, together with a soliton core of It is interesting that the model shown in Figure 6 also qualitatively reproduces the asymmetry of the CMZ, i.e., it shows more emission on the positive longitude compared to the negative longitude in the (l, v) space. From the face-on view, we can see that the gas forms some feather-like features on the left side of the CMZ. This may be caused by the "wiggling" instability proposed by Wada & Koda (2004) then further studied by Kim et al. (2014) and Sormani et al. (2017). A higher c s enhances the instability, thus making the asymmetry more prominent than the lower c s case. Some of the "feathers" roughly reproduce the observed Sagittarius B2 cloud at »  l 1 . It is worth noting that the asymmetry is transient, with a typical timescale of a few tens Myr, similar to some of the more sophisticated simulations (e.g., Sormani et al. 2018b).
We thus conclude that a model with a light and compact soliton core (i.e., . Nevertheless, it is still worth exploring for a better match to the observed features in the future studies.

What Does a Light and Compact Soliton Core Imply?
We use . We use r c,0 to denote the core radius of a soliton without compression. Once m 22 is fixed, M smax and r c,0 are determined for a given halo mass. A larger m 22 corresponds to a less massive and more concentrated soliton core. However, due to the compression by the nuclear bulge, the real core radius r c should be smaller than the original r c,0 , as we have shown in Section 2.1.3. We therefore define a dimensionless compression parameter = C r r c c,0 to describe the level of compression. The . Note that M smax is very close to the mass of the nuclear bulge derived in Launhardt et al. (2002). The sound speed c s is -10 km s 1 . The inset in the bottom-right panel does not have a ring radius as the density of gas is not high enough, similar to Figure 3. The red box highlights the models that roughly match the observed CMZ size.
compressed core radius r c can be written as which is independent of the halo virial mass. Substituting = C r r c c,0 in Equation (5) recovers the original results in Schive et al. (2014b). Here C=1 corresponds to a selfgravitating soliton core, while for < C 1, the soliton core is compressed by baryons. Equations (4) and (5) suggest that for a given soliton mass and core radius ( ) M r , s c , there is a degeneracy between the boson mass m 22 and the compression parameter C. In other words, both C and m 22 can be varied accordingly to give the desired r c for a certain M s .
We first show that the compression parameter C in our best-fit models cannot be achieved by the background ¡ = 0.5 nuclear bulge if m 22 is around 1. The parameter set ( ) M r , s c of our best-fit models in Figures 4 and 5, i.e., panels (b), (e), (f), (i), and (j), can be interpreted as = C 0.16, 0.04, 0.08, 0.04, 0.02 with = m 1 22 , as suggested by Equation (5). We numerically solve the Schrödinger-Poisson equation based on the method in Bar et al. (2019) to self-consistently derive the soliton core radius r c in the ¡ = 0.5 nuclear bulge potential. We find r c should be around 160 pc, 184 pc, and 196 pc for a soliton core with mass of 1/2, 1/4, and M 1 8 smax , assuming = m 1

22
. The corresponding C is therefore 0. 50, 0.29, 0.29, 0.15, 0.15 for panels (b, e, f, i, j), about 3-7 times larger than the numbers quoted above. This . Considering all the models highlighted by the red box (i.e., b, e, f, i, and j) in Figures 4 and 5, we conclude that a boson mass m 22 in the range of 2-7 is more reasonable in the current configuration of the Galactic potential that contains both the soliton core and the nuclear bulge.
We need to emphasize that our upper limit for m 22 comes from the assumption that a minimum mass for the soliton core must exist, otherwise the CMZ will not form. However, as shown in Section 3.2, a nuclear-bulge-only model could also form a reasonable CMZ consistent with observed kinematics, as long as we assume the nuclear bulge to be more concentrated than that in Launhardt et al. (2002). Therefore, the upper limit we get (i.e., m 7

22
) is only a weak constraint-it could be larger than 7 if the nuclear bulge in reality is more compact in mass compared to the model we adopt in the current paper. On the other hand, the lower limit should be a stronger constraint (i.e., m 2 22 ), as we have already assumed a relatively less massive nuclear bulge compared to previous studies.

The Possibility of a Higher Boson Mass
As pointed out above, an m 22 range of 2-7 is preferred in this work. This boson mass range is marginally higher than the conventional value  Bar et al. 2018Bar et al. , 2019. The larger boson mass also helps alleviate the tension arising from the high-z luminosity function (e.g., Schive et al. 2016;Corasaniti et al. 2017; and Lyα forest (e.g., Iršič et al. 2017;Leong et al. 2019), where either a larger boson mass or an extreme axion model is typically required to produce a sufficient amount of small-scale structure. Self-consistent FDM hydrodynamical simulations will be essential for narrowing down the mass constraints further.

Central Mass Profile of the Milky Way
The CMZ size and kinematics constrain the total mass profile in the central~200 pc, but it is still not quite clear whether the central mass is contributed by a compact stellar nuclear bulge, or a less compact nuclear bulge together with a dark soliton core. This mass degeneracy between a nuclear bulge and a soliton core was first discussed in Bar et al. (2018), and a better determination on the mass-to-light ratio of the nuclear bulge (i.e., stellar mass) would be helpful to break this degeneracy. In Figure 7, we plot the enclosed mass profiles for two of our bestfit models with and without the soliton core. The red solid line in Figure 7 shows the enclosed mass profile from our nuclearbulge-only model (Figure 3(d)); the central mass in this model is contributed entirely by baryons. The blue solid line shows the enclosed mass profile from the model that contains both a 10 7 using the GALACTICNUCLEUS survey; their results are indicated by the purple triangles in this plot. We see that this stellar mass is quite consistent with our nuclear-bulge-only model without a soliton core. If their mass estimation is accurate, then there is little room left for the FDM soliton core in the Milky Way center. Nevertheless, we still need further observations to provide tighter constraints on the baryon mass profile to better clarify this problem.

Other Dynamical Effects of the Soliton Core
We find that a soliton core with a mass of  »´M 3.6 10 8 and a nuclear bulge with ¡ = 0.5 would generate a CMZ that matches the observations relatively well. The total mass of these two components is therefore  »´M 0.7 10 9 . However, the velocity dispersion revealed by Gaia and VVV surveys at | |   l 10 requires a dense center of  »´M 1.5 10 9 as found by De Martino et al. (2018) and Portail et al. (2017), which is about two times larger than the value preferred here. We emphasize that the mass obtained in this paper is constrained mainly by the CMZ region. The CMZ size and kinematics offers no strong constraints on the mass profile outside  R 200 pc (| |   l 1 . 5). Therefore, this discrepancy could be explained by additional mass outside the CMZ radius which could still result in a large velocity dispersion. Indeed, the central velocity dispersion of the stars drops inside  R 200 pc (De Martino et al. 2018), implying a toroidal mass distribution of a less massive center and a more massive periphery.
Interestingly, the soliton core could form such a toroidal shape, due to the baryonic potential as reported in Bar et al. (2019). We emphasize that a toroidal soliton core is just one possible solution. More detailed modeling is still needed to predict the exact soliton core properties in the Milky Way.

Different Physical Conditions of the Gas inside the CMZ
The gas in the CMZ is multiphase (e.g., Langer et al. 2017), but it seems that both the cold and hot gas form a consistent feature in the (l, v) space for which we are trying to match using our isothermal simulations. The isothermal equation of state is a simple assumption, and it forms a slightly smaller CMZ compared to the simulations with a complicated chemistry cooling network, as shown in Sormani et al. (2018b). The cooling mechanism also helps to break the high-density gas in the CMZ into small clumps seen in the NH 3 data, while in our simulation, the gas is more smoothly distributed in the CMZ. The cooling timescale for the ISM is much shorter than the dynamical time, thus we expect the kinematics of the CMZ would not be significantly affected by considering the chemical evolution. The mass of the gas in the CMZ is around  M 3.0 10 7 (Launhardt et al. 2002), which is also small compared to the mass of the soliton core and the nuclear bulge used in this paper.
In general, including self-gravity and cooling would make the CMZ larger (e.g., Kim et al. 2012a;Sormani et al. 2018b), while including magnetic field, star formation, and stellar feedback would make it smaller (e.g., van de Ven & Chang 2009;Kim & Stone 2012). These physics basically influence condition (2) mentioned in Section 3.2, i.e., where the inflowing gas piles up. Condition (1) is still determined mainly by the gravitational potential. A recent work by Armillotta et al. (2019) studied the star formation in the CMZ with a detailed modeling of gas physics and the stellar evolution. In their models, a CMZ of »200 pc forms, very similar to the isothermal case with the same Galactic potential adopted in Ridley et al. (2017). We therefore regard our isothermal simulations as a first-order approximation to the more complicated simulations. We will explore further the effects of these baryonic physics on the constraints of the soliton core in follow-up studies.

The 3D Structure of the CMZ
We have verified that 3D simulations produce very similar gas flow patterns to those in the 2D case. Gas cannot go very far above the plane because the deep gravitational potential caused by the soliton core and the nuclear bulge. The edge-on ¥ shape of the CMZ discovered by Molinari et al. (2011) can also be seen in our 3D simulations. However, the ¥ shape in the simulations exists only in the first~200 Myr then gradually decays, as also reported in the smoothed particle hydrodynamics simulations of Shin et al. (2017). Because it is a transient feature, we prefer not to use it as an indicator to constrain the possible existence of the soliton core. The observed ¥ shape of the CMZ might be caused by a recent inflow, as suggested in Sormani & Barnes (2019).

Summary
We have performed high-resolution hydrodynamical simulations of bar-driven gas flows in a realistic Milky Way potential, considering both the effects of a baryonic nuclear bulge and a Figure 7. Enclosed mass in spheres of radius R for the inner 300 pc of the Galaxy. The thick blue line is the result from Figure 14 in Launhardt et al. (2002). The red solid line shows the mass distribution in Figure 3(d), which is our best-fit model that includes only the nuclear bulge. The blue solid line shows that used in Figures 4(f) and 5(f), which is one of our best-fit models that contains both the nuclear bulge and the soliton core. The blue dashed line plots the baryonic mass in this model. The purple triangle is the stellar mass (  ´M 6.5 0.4 10 7 ) obtained by Nogueras-Lara et al. (2019) at 45 pc. Note that this mass is not the enclosed mass inside a = r 45 pc sphere, but the total mass inside a = R 45 pc and = z 20 pc cylinder. We also include a  M 2.6 10 6 black hole and a  M 3.0 10 7 nuclear stellar cluster in the models to be consistent with Launhardt et al. (2002).
possible dark soliton core. The observed size and kinematics of the CMZ could be reproduced by including a compact central mass component, whose mass profile is relatively well constrained by our models. If a soliton core is not considered, a more compact nuclear bulge than usually assumed could match the observed size and kinematics of the CMZ. Including a moderate soliton core together with a less massive nuclear bulge also nicely agrees with observations. An effective sound speed of gas around -20 km s 1 can help further improve the match. The preferred soliton core could be achieved if the boson mass is larger than the conventional´-1.0 10 eV 22 by a factor of 2-7. Such a boson mass range is also broadly consistent with many other astrophysical constraints. This mass range can be further narrowed down by the improved determination of the mass-to-light ratio of the nuclear bulge. The 3D ¥ shape structure of the CMZ probably does not offer a tighter constraint on the soliton core because it is likely a transient feature that occurs only at the beginning of the inflow process. This may imply that the current CMZ may have just experienced a recent inflow, which also helps to explain the low star formation rate in the CMZ caused by a large gas velocity dispersion.