Arxiv:2010.10337V3 [Astro-Ph.GA] 19 Jan 2021 Galaxies Should Have the Steepest, Most Concentrated Pro- the Term Ψdm

Detection of a universal core-halo transition in dwarf galaxies as predicted by Bose-Einstein dark matter

Alvaro Pozo1,2, Tom Broadhurst1,2,3, Ivan de Martino4, Tzihong Chiueh9,10, George F. Smoot2,5,6,7,8, Silvia Bonoli,2,3, Raul Angulo2,3. 1Department of Theoretical Physics, University of the Basque Country UPV/EHU, E-48080 Bilbao, Spain; email:[email protected]; [email protected]; 2DIPC, Basque Country UPV/EHU, E-48080 San Sebastian, Spain 3Ikerbasque, Basque Foundation for Science, E-48011 Bilbao, Spain 4Universidad de Salamanca,Facultad de Ciencias.,F´ısica Teórica, Salamanca, Plaza de la Merced s/n. 37008, Spain; [email protected] 5Institute for Advanced Study and Department of Physics, IAS TT & WF Chao Foundation Professor, Hong Kong University of Science and Technology, Hong Kong 6Energetic Cosmos Laboratory, Nazarbayev University, Nursultan, Kazakhstan 7Physics Department, University of California at Berkeley CA 94720 Emeritus 8Paris Centre for Cosmological Physics, APC, AstroParticule et Cosmologie, UniversitéParis Diderot, CNRS/IN2P3, CEA/lrfu,UniversitéSorbonne Paris Cité, 10, rue Alice Domon et Leonie Duquet, 75205 Paris CEDEX 13, France Emeritus 9Department of Physics, National Taiwan University, Taipei 10617, Taiwan 10National Center for Theoretical Sciences, National Taiwan University, Taipei 10617, Taiwan (Dated: January 20, 2021) The presence of large dark matter cores in dwarf galaxies has long been puzzling and many are now known to be surrounded by an extensive halo of stars. Distinctive core-halo structure is characteristic of dark matter as a Bose Einstein condensate, ψDM, with a dense, soliton core predicted in every galaxy, representing the ground state, surrounded by a large, tenuous halo of excited density waves. A marked density transition is predicted between the core and the halo set by the de Broglie wavelength, as the soliton core is a prominent standing wave that is denser by over an order of magnitude than the surrounding halo. Here we identify this predicted behavior in the stellar profiles of the well known ”isolated” dwarfs that lie outside the Milky Way, each with a clear density transition at ' 1.0 kpc, implying a very light −22 boson, mψ ' 10 eV. The classical dwarf galaxies orbiting within the Milky Way also show this predicted core-halo structure but with larger density transitions of over two orders of magnitude, that we show implies tidal stripping of dwarf galaxies by the Milky way, as the tenuous halo is more easily stripped than the stable soliton core. We conclude that dark matter as a light boson explains the observed family of classical dwarf profiles with tidal stripping included, in contrast to the standard heavy particle interpretation where low mass galaxies should be concentrated and core-less, quite unlike the core-halo structure observed.

The well resolved classical dwarf galaxies orbiting matter cores has been uncovered in the first simula- the Milky Way are typically concluded to have sizeable tions of dark matter (DM) in a Bose-Einstein conden- dark matter cores of radius ' 0.5 − 1.0 kpc (Walker sate state (Schive et al. 2014, Schive et al. 2014b, & Peñarrubia2011, Amorisco & Evans 2012), or rela- Schwabe et al. 2016) revealing unanticipated, rich wave- tively shallow central density profiles s (e.g. Hayashi like structure on the de Broglie scale, described by a et al. 2020), much flatter than predicted by standard coupled Schrodinger-Poisson equation for the mean field ΛCDM simulations of heavy particle DM, where dwarf behaviour that evolves simply under self-gravity, hence arXiv:2010.10337v3 [astro-ph.GA] 19 Jan 2021 galaxies should have the steepest, most concentrated pro- the term ψDM. A prominent standing wave is found to files. Proposals invoking repeated gas outflows to gen- comprise the ground state formed at the center of every erate cores are not reinforced by the highest resolution galaxy, surrounded by an extensive halo of slowly vary- simulations where little physical effect on the dark matter ing density waves that fully modulate the halo density profile is found to be induced by repeated gas loss (Bose on the de Broglie scale by self-interference (Schive et al. et al. 2019). This tension is exacerbated by discover- 2014, Schive et al. 2014b, Mocz et al. 2017, Veltmaat ies of extremely low surface brightness “ghostly” dwarfs, et al. 2019 , Hui et al. 2020). Condensates are inher- Crater II, And XIX and Antlia II, that have even wider ently non-relativistic and hence ψDM is viable as dark cores than the classical dwarfs, of up to 3 kpc with lower matter, behaving as CDM on large scale exceeding the velocity dispersions (Torrealba et al. 2019, Collins et al de Broglie wavelength, demonstrated by the pioneering 2019). simulations of Schive et al. (2014). This ψDM interpre- A physically very different explanation for wide dark tation makes several unique predictions, including a core ii

FIG. 1. Isolated Dwarf Galaxies: The star count based profiles of the well-studied isolated dwarf galaxies in the local group, lying beyond the Milky Way. Each of these galaxies has an extended halo of stars stretching to ' 2 kpc and beyond. These are compared with our predictions for light boson dark matter, ψDM, in green, where the core has a soliton profile that is close to Gaussian in form, surrounded by halo of excited states that average azimuthally to an approximately NFW-like profile, beyond the soliton radius. Our best fit MCMC based profiles are found to fit well, including the characteristic density drop of about a factor of ' 30 between the prominent core and tenuous halo at a radius ' 1 kpc (indicated by vertical orange lines), see Table I. The best fit NFW profile to the halo is shown as a thin grey line, above which the core region is seen to be a prominent overdensity in each case. profile with the unique soliton profile (see methods) and kpc, corresponding to the lower mass limit for galaxies that the soliton radius is set by the de Broglie wavelength permitted by ψDM, with the above boson mass (Broad- and hence the soliton radius, rsol, is larger in lower mass hurst et al. 2020). This momentum dependence for the galaxies, mgal, of lower momentum, following this scal- soliton width is also found to explain the excessive star −1/3 and gas motions measured near the center of the Milky ing relation: rsol ∝ m found by Schive et al. (2014b) gal Way (Portail et al. 2017) that imply a soliton core of and verified in independent simulations, (Schwabe et al. 9 2016, Niemeyer 2020). about 150pc and 10 M (de Martino et al. 2018), which falls on the predicted soliton-halo mass scaling relation. The predicted soliton profile has been shown to match An intermediate sized soliton plus halo is derived for the in detail the cored mass profile of the Fornax dwarf prototypical “ultra diffuse galaxy” DF44 (van Dokkum for which the dynamical data extends beyond the core et al. 2019) from its extended flat velocity dispersion of radius allowing a determination of the boson mass of −22 −22 ' 33 km/s that is also consistent with mψ ' 10 eV mψ ' 10 eV (Schive et al. 2014). This is also sup- (Pozo et al. 2020). ported by Jeans analyses of other classical dwarfs (Chen et al. 2017). The momentum dependence of the soliton Here we explore a related prediction of ψDM that a radius has also been tested with the ghostly dwarf Antlia characteristic drop in density should be found at the II, uncovered by GAIA, which has a larger extent and transition radius between the core and the halo. The lower mass than the classical dwarfs, requiring a wider soliton forms a prominent core that contrasts by over an soliton radius of ' 3 kpc than for the classical dwarfs, ' 1 order of magnitude in density at the core radius with iii

FIG. 2. Stellar profiles of classical dwarfs orbiting the Milky Way and also the four known ”isolated” dwarfs that lie beyond in the Local Group (all listed in Table I), rescaled by their transition radius, revealing the profiles have a common core-halo form that is more pronounced for for the orbiting dwarfs, indicated by redder colours, for which the halo density is relatively compared to the ”isolated” dwarfs shown in blue colour. The right hand panel compares the change in stellar density, or ”gap”, between the core and halo versus pericenter radius from GAIA proper motions, showing the isolated dwarfs beyond the Milky Way (blue points) have generally smaller transition amplitudes compared to most of the orbiting dwarfs (red points). Sextans is the only obvious exception, shown in green, as discussed in the text in the context of tidal stripping. the surrounding halo. The halo comprises de Broglie Milky Way where significant tidal stripping is evident in scale waves with an approximately Navarro-Frenk-White some cases (Ibata et al. 2001). (NFW) form when averaged azimuthally, reflecting the We begin by examining the profiles of the four well cold, non-relativistic nature of the condensate on scales studied isolated dwarf galaxies in the Local Group, Ce- exceeding the de Broglie wavelength (Schive et al. 2016). tus, Tucana, Aquarius and Leo A, which have small ve- Even in projection this transition is predicted to be dis- locity dispersions ' 10 km/s and old spheroidal stellar tinct because the soliton core is close to Gaussian in form, populations, like the common classical dwarf spheroidal so its sharp 3D boundary is preserved in 2D at the same galaxies orbiting the Milky Way and Andromeda. How- radius. It is important to appreciate this prominent core ever, these dwarfs are understood not to have suffered is quite unlike the behaviour of smooth cores employed tidal stripping or ram pressure gas stripping by the main modelling dark halos, where the core is continuous in jor Local Group members. This simplicity has motivated density with the halo. In contrast, the ψDM core is a sta- deep imaging of these isolated galaxies and they are close ble standing wave that is a gravitationally self-reinforcing enough to detect individual stars with the aim of estab- (Schive et al. 2014b) with a pronounced overdensity of lishing intrinsic dwarf galaxies properties using individ- > 30 above the halo predicted for low mass galaxies of ually detected stars which allows the entire profile to be 10 ' 10 M (see Figures 1 & 2 of Schive et al. 2014b) and traced free of surface brightness limitations. These star is shown to be stable to tidal stripping in a new simula- counts have established that the stellar profiles of each tion of a dwarf galaxy orbiting the Milky Way (Schive et of these four galaxies extends to 2 kpc or more in radius al. 2020) where the density profile of the soliton self re- (see Figure 1), somewhat surprisingly for low mass dwarf inforces changing relatively slowly compared to the outer galaxies (Gregory et al. 2019, McConnachie et al. 2006, halo that is readily stripped. It is this distinctive profile McConnachie et al. 2006b) prediction we aim to test here by comparing the profiles Our MCMC based profile fits to these isolated galaxy of ”isolated” dwarf galaxies lying beyond the Milky Way profiles are shown in Figure 1, with the best parameters and also with all the classical dwarfs orbiting within the and uncertainties listed in Table 1 and in Figure 6-9 in the iv

FIG. 3. The left hand panel compares the dwarf profiles scaled by their individual best fitting core radius. The profiles are very similar in the core region, differing in the extended halo region where the density variation is 2 orders of magnitude relative to the core, with the isolated galaxies shown with bluer colours having denser halos than the orbiting dwarfs, shown with redder colours. The right hand panel shows the projected simulation profiles for ψDM by Schive, Chiueh & Broadhurst (2020), where little dependence on the degree of tidal stripping is predicted for the core region, in contrast to the halo region bracketed in grey. The coloured curves indicate the duration of stripping, with the timescale indicated by the legend in the inset of the left hand panel for the 3D profiles, with matched colour coding, spanning 4Gyrs, that is very similar to the range of the data. Note, the predicted halo slope is relatively shallow and independent of the degree of stripping and agrees well with the mean galaxy profile of the isolated dwarfs (blue data points) and with the orbiting dwarf galaxies (red data points) with a larger core-halo transition for the orbiting dwarfs (red data points). Standard Plummer profiles are also shown (see legend), that fit approximately the core region but fall well short at large radius. The NFW profile is also shown to match the halo region but falls well below the prominent core.

Supplementary information. The soliton profile for the core halo structure extending beyond ' 1.0 kpc, with core is fully characterised by its core radius free param- very similar halo slopes that are relatively shallow. In eters, rc and for the stellar halo fit an NFW profile with Figure 3 we separately scale the radial profiles by their scale radius rs∗ and normalization ρo∗. Finally, the tran- individually fitted core radius which reveals clearly the sition radius, rt, between the soliton and NFW profile distinctive core-halo profile with on average a relatively is also allowed to vary, within the prior range indicated larger density transition for the orbiting dwarfs (red data, by the ψDM simulations (see Methods section). The fit Figure 3b) compared to the mean isolated dwarf profile results given in Table 1 are scaled to a fixed boson mass (blue data, Figure3 b). Figure 3a reveals the measured that we set to the canonical value indicated by dwarf dy- profiles are very similar in form within the core region, −22 namics, mψ = 1.5 × 10 eV since it is degenerate with below the mean transition radius of 0.75kpc, whereas at halo mass (see methods section). larger radius the star density profile varies relative to the core by two orders of magnitude, but with a similar We now compare these isolated dwarfs to the well stud- slope. This behaviour is very much like the prediction un- ied, classical dwarf spheroidal galaxies orbiting within the der tidal stripping, shown in the right hand panel, where Milky Way for which it has long been claimed that halos the soliton profile remains relatively unchanged whereas of stars are a common feature, originally termed “ex- the halo declines in amplitude by 2 orders of magnitude tra tidal stars, for convenience” (Irwin & Hatzidimitriou with increasing stripping, while maintaining a predicted 1995). Such halos are now very clear in the latest deep NFW-like form (see Figure 3b) in the halo region. imaging, as shown in Figure 2a, and there is a continuous range of profiles, most of which have a clearly extended We define a transition gap; ∆C−H = log ρC /ρH , where v

stripped, and so it is interesting that several claims have been made that Sextans lacks any sign of tidal interaction from its smooth morphology and internal dynamics, concluding that Sextans is relatively newly bound to the Milky Way (Roderick et al. 2016 , Okamoto et al. 2017). This conclusion is also supported by the relatively large +11 pericenter distance of Sextans, 71−12 kpc and near circu- lar orbit that is now well established with GAIA proper motions(Fritz et al. 2018). An independently estimated boson mass for Sextens was derived in our earlier Jeans +0.41 −22 analysis,mψ = 1.23−0.33 × 10 eV, using the velocity dispersion profile (Broadhurst et al. 2020), under the implicit assumption of no significant stripping which is consistent with our conclusion here based on the shape of the stellar profile shape. In the case of the isolated dwarf Tucana, sufficient stellar velocity dispersion data is available to allow an independent dynamical consistency check via Jeans analysis and this we show in Figure 4 (described in the Supple- mentary information), finding consistency with our best fit model profile model, for a mildly negative velocity FIG. 4. The predicted velocity dispersion profile for Tucana anisotropy, like the well studied dwarf spheroidal galax- using our best core-halo profile calculated using the MCMC ies (Schive et al.2014 , Chen et al. 2017). It will be analysis of Figures 1 and 7, and listed in Table 1. Notice that important to compare the core radius defined by the star Eq. (B2) has been used instead of the usual Plummer model to determine the stellar density profile, assuming the stars profile to that of the mass profile inferred dynamically, as approximately trace the DM mass profile. these may differ somewhat depending on the cooling of the gas relative to the DM, particularly for the old star population that typically dominates with a larger radial extent and higher velocity dispersion than younger pop- ρ is the asymptotic central core stellar density. The C ulations that are often present with lower velocity dis- value of the gap is seen in Figure 2b to be generally persions. The effect of tidal stripping should be equally higher for the orbiting dwarfs than the isolated dwarfs, effective for DM and the stars as tracers of the DM, mo- with ∆ > 2.0. The largest apparent core-halo tran- C−H tivating deeper imaging of the many other less studied sition is found for Leo I, with ∆ ' 3.5 (Figure 2b) C−H dwarfs associated with Andromeda for a fuller study in- and so it is interesting that although Leo I is relatively cluding tidal effects based on accurate stellar profiles. distant now at 250 kpc, its orbital pericenter is now es- Lensing will complement this dynamical work within the tablished by the GAIA satellite to be only ' 50 kpc, halo where full density modulation from self interference moving on an eccentric orbit (see table 1) where tidal is predicted to measurably perturb image fluxes and po- stripping is enhanced during close approaches (Sohn et sitions near the Einstein ring(Chan et al. 2020, Hui et al. 2018). Tidal effects are also claimed for Ursa Minor al. 2020). Oscillation modes of the central soliton in- (Martinez-Delgado et al. 2001), for Sculptor (Walker & duced by the de Broglie scale fluctuations of the halo can Peñarrubia2011, Amorisco & Evans 2012) and for Ca- cause dynamical heating of stars (Church et al. 2019) rina (Battaglia et al. 2012b). In the context of ψDM, the and random motions of star clusters within the soliton tidal effects are quite different from the CDM based ex- that random walk to the core radius where they may be pectations, as the tenuous halo is found to be efficiently disrupted tidally (Schive et al. 2020), marking the core stripped in ψDM simulations (Du et al. 2018, Schive radius(de Boer et al. 2013). etal 2020) to produce a relatively large break at the core radius (see Figure 3b) with a resilient core due to the sta- bility of the soliton, that is hard to disrupt tidally (Du et al. 2018) and which may account for the morphological symmetry of Draco (Segall et al. 2007) and the marked drop in velocity dispersion beyond the core (Wilkinson et al. 2004), warranting further examination with a tailored ψDM simulation. The Sextans dwarf galaxy is an exception as its stellar profile is like the isolated dwarfs, with an extensive core- halo structure and a relatively small transition gap of ∆C−H ' 1.3 (Figure 2b). This suggests, in the context of ψDM, that Sextans has not been significantly tidal vi

References 12. de Martino, I., Broadhurst, T., Henry Tye, S. -H.; Chiueh, T. & Schive, H.-Y. Dynamical Evidence of 9 1. Amorisco, N. C. & Evans, N. W. Dark matter cores a Solitonic Core of 10 M . Phys. Dark Universe and cusps: the case of multiple stellar populations 28, 100503 (2018). in dwarf spheroidals. Mon. Not. R. Astron. Soc. 419, 184-196. (2012). 13. de Martino, I., Broadhurst, T., Tye, S. -H. H.,; Chiueh, T., Schive, H.-Y. & Lazkoz, R. Recognising 2. Battaglia, G., Rejkuba, M., Tolstoy, E., Irwin, M. Axionic Dark Matter by Compton and de-Broglie J.& Beccari, G. A wide-area view of the Phoenix Scale Modulation of Pulsar Timing. Phys. Rev. dwarf galaxy from Very Large Telescope/FORS Lett. 119, 221103 (2017). imaging. Mon. Not. R. Astron. Soc. 424, 1113- 14. Du, X., Schwabe, B.,; Niemeyer, J. C.& Bürger,D. 1131. (2012). Ghostly galaxies as solitons of Bose-Einstein dark 3. Battaglia, G., Irwin, M., Tolstoy, E., de Boer, T. & matter. Phys. Rev. Lett. 97, id.063507. (2018). Mateo, M. The extensive age gradient of the Carina 15. Frinchaboy, P. M., Majewski, S. R., Muñoz,R. R., dwarf galaxy Astrophys. J. 761, L31. (2012). Law, D. R.,Lokas, E. L., Kunkel, W. E., Patterson, R. J.& Johnston, K. V. A 2MASS All-sky View of 4. Binney, J. &; Tremaine, S. Galactic Dynamics: Sec- the Sagittarius Dwarf Galaxy. VII. Kinematics of ond Edition Princeton Univ. (2008). the Main Body of the Sagittarius dSph. Astrophys. 5. Broadhurst, T., de Martino, I., Luu, H. N., Smoot, J. 756, 19. (2012). G. F.& Tye, S. -H. H. Ghostly galaxies as solitons 16. Fritz, T. K., Battaglia, G., Pawlowski, M. S., Kalli- of Bose-Einstein dark matter. Phys. Rev. D. 101, vayalil, N., van der Marel, R., Sohn, S. T., Brook, id.083012. (2020). C. & Besla, G. Gaia DR2 proper motions of dwarf galaxies within 420 kpc. Orbits, Milky Way mass, 6. Bose, S., Frenk, C. S., Jenkins, A., Fattahi, A., tidal influences, planar alignments, and group in- Gómez, F. A., Grand, R. J. J., Marinacci, F., fall. Astron. Astrophys. 619, A103 (2018). Navarro, J. F., Oman, K. A., Pakmor, R., Schaye, J., Simpson, C. M.& Springel, V. No cores in dark 17. Gregory, A. L., Collins, M. L. M., Read, J. I., Irwin, matter-dominated dwarf galaxies with bursty star M. J., Ibata, R. A., Martin, N. F., McConnachie, formation histories. Mon. Not. R. Astron. Soc. A. W.& Weisz, D. R. Kinematics of the Tucana 486, 4790-4804 (2019). dwarf galaxy: an unusually dense dwarf in the Lo- cal Group. Mon. Not. R. Astron. Soc. 485, 7. Chan, J. H. H., Schive, H.-Y., Wong, S.-K., Chiueh, 2010-2025 (2019). T.& Broadhurst, T. Multiple Images and Flux Ra- tio Anomaly of Fuzzy Gravitational Lenses. ArXiv 18. Hayashi, K.,; Chiba, M. & Ishiyama, T. Diver- 2002.10473,(2020). sity of dark matter density profiles in the Galactic dwarf spheroidal satellites. Astrophys. J. 904, 19. 8. Chen, S.-R., Schive, H.-Y.& Chiueh, T. Jeans anal- (2020). ysis for dwarf spheroidal galaxies in wave dark matter. Mon. Not. R. Astron. Soc. 468, 1338-1348 19. Hu, W., Barkana, R. & Gruzinov, A. Fuzzy Cold (2017). Dark Matter: The Wave Properties of Ultralight Particles. Phys. Rev. Lett. 85, 1158-1161 (2000). 9. Church, B. V., Mocz, P. & Ostriker, J. P. Heating of Milky Way disc stars by dark matter fluctua- 20. Hui, L., Joyce, A., Landry, M. J.& Li, X. Vortices tions in cold dark matter and fuzzy dark matter and waves in light dark matter Arxiv 2007.13780, paradigms. Mon. Not. R. Astron. Soc. 485, (2020). 2861-2876. (2019). 21. Irwin, M.& Hatzidimitriou, D. Structural parameters for the Galactic dwarf spheroidals. Mon. Not. 10. Collins, M. L. M., Tollerud, E. J., Rich, R. M., R. Astron. Soc. 277, 1354-1378. (1995). Ibata, R. A., Martin, N. F., Chapman, S. C., Gilbert, K. M.,& Preston, J. A detailed study of 22. Ibata, R., Lewis, G. F., Irwin, M., Totten, E. & Andromeda XIX, an extreme local analogue of ul- Quinn, T. Great Circle Tidal Streams: Evidence for tra diffuse galaxies. Mon. Not. R. Astron. Soc. a Nearly Spherical Massive Dark Halo around the 491, 3496-3514 (2019). Milky Way. Astrophys. J. 551, 294-311. (2001). 11. de Boer, T. J. L., Tolstoy, E., Saha, A. & Olszewski, 23. Kang, H. D. & Ricotti, M. Ghostly haloes in dwarf E. W. A new study of stellar substructures in the galaxies: constraints on the star formation effi- Fornax dwarf spheroidal galaxy. Astron. Astro- ciency before reionization. Mon. Not. R. Astron. phys. 551, 7 (2013). Soc. 488, 2673-2688 (2019). vii

24. Mart´ınez-Delgado,D., Alonso-Garc´ıa,J., Aparicio, 36. Schive, H.-Y., Chiueh, T., Broadhurst, T.&; A.& Gómez-Flechoso, M. A. A Tidal Extension in Huang, K.-W. Contrasting Galaxy Formation from the Ursa Minor Dwarf Spheroidal Galaxy. Astro- Quantum Wave Dark Matter, ψDM, with ΛCDM, phys. J. 549, L63-L66. (2001). using Planck and Hubble Data. Astrophys. J. 818, 14 . (2016). 25. McConnachie, A. W. & Irwin, M. J. Structural properties of the M31 dwarf spheroidal galaxies. 37. Schive, H.-Y., Chiueh, T. & Broadhurst, T. Soliton Mon. Not. R. Astron. Soc. 365, 1263-1276 (2006). Random Walk and the Cluster-Stripping Problem in Ultralight Dark Matter. Phys. Rev. Lett. 124, 26. McConnachie, A. W., Arimoto, N., Irwin, M. & 201301 (2020). Tolstoy, E. The stellar content of the isolated tran- 38. Schwabe, B., Niemeyer, J. C. & Engels, J. F. Simu- sition dwarf galaxy DDO210. Mon. Not. R. As- lations of solitonic core mergers in ultralight axion tron. Soc. 373, 715-728. (2006). dark matter cosmologies Phys. Rev. D. 94, 043513 (2016). 27. Mocz, P., Vogelsberger, M., Robles, V. H., Zavala, J., Boylan-Kolchin, M., Fialkov, A.& Hernquist, 39. Segall, M., Ibata, R. A., Irwin, M. J., Martin, N. L. Galaxy formation with BECDM - I. Turbulence F. & Chapman, S. Draco, a flawless dwarf galaxy*. and relaxation of idealized haloes. Mon. Not. R. Mon. Not. R. Astron. Soc. 375, 831-842. (2007). Astron. Soc. 471, 4559-4570 (2017). 40. Sohn, S.T., Majewski, S. R., Muñoz, R. R., Kunkel, W. E., Johnston, K. V., Ostheimer, J. 28. Navarro, J.F., Frenk, C. S. & White, Simon D. M. C., Guhathakurta, P., Patterson, R. J., Siegel, M. The Structure of Cold Dark Matter Halos. Astro- H.& Cooper, M. C. Exploring Halo Substructure phys. J. 462, 563 (1996). with Giant Stars. X. Extended Dark Matter or Tidal Disruption?: The Case for the Leo I Dwarf 29. Niemeyer, J. C. Small-scale structure of fuzzy and Spheroidal Galaxy. Astrophys. J. 663, 960-989. axion-like dark matter. P. in Particle & Nuclear (2007). Phys. 113, id. 103787. (2020). 41. Torrealba, G., Belokurov, V., Koposov, S. E., Li, 30. Okamoto, S., Arimoto, N., Tolstoy, E., Jablonka, T. S., Walker, M. G., Sanders, J. L., Geringer- P., Irwin, M. J., Komiyama, Y., Yamada, Y.& Sameth, A., Zucker, D. B., Kuehn, K., Evans, N. Onodera, M. Population gradient in the Sextans W.& Dehnen, W. The hidden giant: discovery of dSph: comprehensive mapping of a dwarf galaxy an enormous Galactic dwarf satellite in Gaia DR2. by Suprime-Cam. Mon. Not. R. Astron. Soc. Mon. Not. R. Astron. Soc. 488, 2743-2766 (2019). 467, 208-217. (2017). 42. van Dokkum, P., Wasserman, A., Danieli, S., Abra- 31. Portail, M., Wegg, C., Gerhard, O. & Ness, M. ham, R., Brodie, J., Conroy, C., Forbes, D. A., Chemodynamical modelling of the galactic bulge Martin, C., Matuszewski, M., Romanowsky, A. J.& and bar. Mon. Not. R. Astron. Soc. 470, 1233- Villaume, A. Spatially Resolved Stellar Kinematics 1252. (2017). of the Ultra-diffuse Galaxy Dragonfly 44. I. Obser- vations, Kinematics, and Cold Dark Matter Halo 32. Pozo, A., Broadhurst, T., de Martino, I., Nhan Fits. Astrophys. J. 880, 91 (2019). Luu, H., Smoot, G. F., Lim, J. & Neyrinck, M. 43. Veltmaat, J., Schwabe, B. & Niemeyer, J. C. Wave Dark Matter and Ultra Diffuse Galaxies. Baryon-driven growth of solitonic cores in fuzzy Arxiv 2003.08313, (2020). dark matter halos. Phys. Rev. Lett. 101, 083518 (2020). 33. Roderick, T. A., Jerjen, H., Da Costa, G. S.& Mackey, A. D. Structural analysis of the Sextans 44. Walker, M. G.& Peñarrubia, J. A Method for dwarf spheroidal galaxy. Mon. Not. R. Astron. Measuring (Slopes of) the Mass Profiles of Dwarf Soc. 460, 30-43. (2016). Spheroidal Galaxies. Astrophys. J. 742, 20 . (2011). 34. Schive, H.-Y., Chiueh, T. & Broadhurst, T. Cosmic 45. Widrow, L. M. & Kaiser, N. Using the Schroedinger structure as the quantum interference of a coherent equation to simulate collisionless matter. Astro- dark wave. Nature Phys 10, 496-499 (2014). phys. J. 416, L71-L74 (1993). 35. Schive, H.-Y., Liao, M.-H., Woo, T.-P., Wong, S.- 46. Wilkinson, M. I., Kleyna, J. T., Evans, N. W., K., Chiueh, T., Broadhurst, T. & Hwang, W. -Y. P. Gilmore, Ge. F., Irwin, M. J.& Grebel, E. K. Kine- Understanding the Core-Halo Relation of Quantum matically Cold Populations at Large Radii in the Wave Dark Matter from 3D Simulations. Phys. Draco and Ursa Minor Dwarf Spheroidal Galaxies. Rev. Lett. 113, 261302 (2014). Astrophys. J. 611, L21-L24. (2004). viii 47. Woo, T. -P. & Chiueh, T. High-resolution simu- Methods: lation on structure formation with extremely light bosonic dark matter. Astrophys. J. 697, 850-861 (2009). Appendix A: The Wave Dark Matter Halo

Ultralight bosons, such as Axions, explored in relation dark matter (Widrow & Kaiser 1993, Hu et al. 2000) in their simplest version, without self-interaction, the boson mass is the only free parameter, which if suffi- ciently light means the de-Broglie wavelength exceeds the mean free path set by the density of dark matter, so these bosons can satisfy the ground state condition for a Bose-Einstein condensate described by the coupled Schroedinger-Poisson equation, that in comoving coordi- nates reads: ∂ ∇2 i + − aV ψ = 0 , (A1) ∂τ 2 ∇2V = 4π(|ψ|2 − 1) . (A2) Here ψ is the wave function, V is the gravitation potential and a is the cosmological scale factor. The system is normalized to the time scale dτ = χ1/2a−2dt, and to the 1/4 1/2 3 2 scale length ξ = χ (mB/~) x, where χ = 2 H0 Ω0 where Ω0 is the current density parameter (Widrow & Kaiser 1993). Recently, it has proved possible with advanced GPU computing to make reliable, high dynamic range cosmological simulations that solve the above equations, (Schive et al. 2014, Schwabe et al 2016, Mocz et al. 2017, May , Springel 2021) that evolve to produce large scale structure indistinguishable from CDM, but with virial- ized halos characterized by a solitonic core in the ground state that naturally explains the dark matter dominated cores of dwarf spheroidal galaxies (Schive et al. 2014b). Another, important feature arising from simulations is that the central soliton is surrounded by an extended halo with a “granular” texture on the de-Broglie scale, due to interference of excited states, but which when azimuthally averaged follows closely the Navarro-Frank- White (NFW) density profile (Navarro et al. 1996, Woo et al. 2009, Schive et al. 2014, Schive et al. 2014b). The fitting formula for the density profile of the solitonic core in a ψDM halo is obtained from cosmological simulations (Schive et al. 2014, Schive et al. 2014b): −1 −23 −2 −4 1.9 a (mψ/10 eV) (rc/kpc) −3 ρc(r) ∼ −2 2 8 M pc . [1 + 9.1 × 10 (r/rc) ] (A3) where the values of the constants are: c1 = 1.9, c2 = −23 −2 10 , c3 = 9.1 × 10 ; mψ is the boson mass and rc is the solitonic core radius. The latter scales with the halo mass of the galaxy and the mass boson obeying to the following the scaling relation which has been calibrated with cosmological simulation (Schive et al. 2014b):

−22 −1/6 −1/3 10 1/2 ζ(z) Mh rc = 1.6 eV a 9 kpc mψ ζ(0) 10 M (A4) ix

Beyond the soliton, at radii larger than a transition observations, as presented in Figure 4: scale (rt), the simulations also reveal the halo is approx- 2 Z ∞ R2 σ2(r)ρ (r) imately NFW in form, presumably reflecting the non- σ2 (R) = 1 − β r ∗ rdr (B4) los 2 2 2 1/2 relativistic nature of condensates beyond the de Broglie Σ(R) R r (r − R ) scale, and therefore the total density profile can be writ- ten as: where Z ∞  2 2 −1/2 ρc(r) if r < rt, Σ(R) = 2 ρ∗(r)(r − R ) rdr . (B5) ρDM (r) = ρ0 (A5) R 2 otherwise,  r 1+ r rs rs and: Appendix C: MCMC Numerical Analysis More in detail, the scale radius of the solitonic solution, which represents the ground state of the Schrodinger- We have explored the full range of relevant parameter Poisson equation, is related to size to the halo through space with the Monte Carlo Markov Chain (MCMC) the uncertainty principle. From cosmological simula- technique based on the Metropolis-Hastings sampling tions, the latter is found to hold non-locally, relating a algorithm, to obtain the core radius that characterises local property with a global one (for more details we refer fully the soliton profile, rc, and the scale radius and to Schive et al. 2014b). the normalization of the stellar density profile,rs∗ and ρ0∗, and the also the transition radius rt between the core and halo profiles. We allow for an adaptive step Appendix B: Dynamical Model of Galaxies in Wave size in order to reach an acceptance rate between 20% Dark Matter Halo and 50%, and we ensure the convergence relying on the Gelman-Rubin criteria. Once the convergence criteria The classical dwarf galaxies are known to be domi- are satisfied the chains are combined to compute the nated by DM, and so the stars are treated as tracer par- total likelihood, together with the 1D marginalized ticles (Gregory et al. 2019, McConnachie et al. 2006, likelihood distribution with the corresponding the Kang et al. 2019 , McConnachie et al.2006b) moving in expectation value and variance. The results are shown the gravitational potential generated by DM halo density in Table 1 for the above free parameters and their distribution. uncertainties. We also show the covariances between the In this context, the corresponding velocity dispersion free paramaters in Figures 6-10 for the isolated galaxies profile can be predicted by solving the spherically sym- and for sextans, as the profiles of these galaxies may be metric Jeans equation: regarded as unaffected by tidal interaction, as described in the text, reflecting intrinsic properties in the context 2 2 d(ρ∗(r)σ (r)) GMDM (r) ρ∗(r)σ (r) of ψDM. r = −ρ (r) − 2β r , dr ∗ r2 r (B1) where MDM (r) is the mass DM halo obtained by in- tegrating the spherically symmetric density profile in Eq. (A5), β is the anisotropy parameter (see Binney & Tremaine 2008, Equation (4.61)), and ρ∗(r) is the stellar density distribution defined by the solitons wave dark matter imprint:  ρ1∗(r) if r < rt, ρ∗(r) = ρ0∗ (B2) 2 otherwise,  r 1+ r rs∗ rs∗ where

ρ0∗ −3 ρ1∗(r) = −2 2 8 N∗kpc (B3) [1 + 9.1 × 10 (r/rc) ]

Here, rs∗ is the 3D scale radius of the stellar halo corresponding to ρ0∗ the central stellar density and the transition radius, rt, is the point where the soliton structure ends and the halo begins at the juncture of the core and halo proﬁles. Finally, the predicted velocity dispersion proﬁle can be projected along the line of sight to compared with the x

Galaxy rc rt rs∗ Gap(Log) Distance Pericenter (kpc) (kpc) (kpc) ∆C−H (kpc) (kpc) +0.01 +0.07 +0.82 Aquarius 0.35−0.01 1.25−0.06 1.05−0.64 1.78 1071 - +0.02 +0.08 +0.14 Cetus 0.36−0.02 0.87−0.07 0.24−0.06 1.34 775 - +0.01 +0.06 +0.90 Tucana 0.25−0.01 0.78−0.06 1.05−0.57 1.5 887 - +0.02 +0.09 +0.49 Leo A 0.43−0.02 1.51−0.08 0.73−0.47 1.76 800 - +0.01 +0.05 +0.51 +11 Sextans 0.48−0.01 1.31−0.06 1.61−0.49 1.36 90 71−12 +0.05 +0.01 +0.54 +126 Phoenix 0.28−0.06 1.27−0.01 1.10−0.55 3.44 415 263−219 +0.01 +0.08 +0.78 +80 Leo I 0.24−0.01 1.30−0.08 1.75−0.96 3.43 250 45−34 +0.01 +0.02 +0.09 +12 Draco 0.17−0.01 0.56−0.02 0.1−0.05 2.32 80 28−7 +0.01 +0.04 +0.51 +21 Carina 0.21−0.01 0.81−0.04 1.17−0.61 2.26 101 60−16 +0.01 +0.07 +0.25 +15 Sculptor 0.20−0.01 0.72−0.07 0.12−0.09 2.18 80 51−10 +0.01 +0.05 +0.90 +8 Ursa Minor 0.28−0.01 0.96−0.04 0.52−0.40 2.14 66 29−6 +0.01 +0.03 +0.42 +121 Leo II 0.14−0.01 0.60−0.03 1.34−0.55 2.45 210 45−30

TABLE I. Observations and ψDM profile fits. Column 1: Dwarf galaxy colour coded as in figure 2b, Column 2: Core radius rc, Column 3: Transition point rt, Column 4: Stellar scale radius rs∗, Column 5: Gap ∆C−H , Column 6: Distances from Milky Way center & Column 7: Pericenter determined from GAIA (Fritz et al. 2018). xi

FIG. 5. Stellar density profiles of all dwarf galaxies listed on table I. The green profile shows the 95% uncertainty on the fitted profile obtained from the MCMC simulation. with the model transition radius rt and uncertainty marked as the vertical orange bar, separating the core and halo regions. References for the data are: Tucana (Gregory et al. 2019), Cetus(McConnachie et al. 2006),Leo A(Kang et al. 2019), Aquarius(McConnachie et al. 2006b),Carina(Frinchaboy et al. 2012),Ursa Minor(Martinez- Delgado et al. 2001), Leo I(Sohn et al. 2007), Phoenix(Battaglia et al. 2012), Sculptor (Frinchaboy et al. 2012) and Draco (Wilkinson et al. 2004)). xii

FIG. 6. Aquarius: correlated distributions of free parameters from MCMC simulation. The core radius and transition radius is well deﬁned here, despite the ﬂat input priors, indicating a reliable result. xiii

FIG. 7. Cetus: correlated distributions of the free parameters. As can be seen the core radius and transition radius are well deﬁned despite the ﬂat input priors, indicating a reliable result. xiv

FIG. 8. Tucana’s correlated distributions of free parameters. As can be seen the core radius and transition radius are well deﬁned despite the ﬂat input priors, indicating a reliable result with a well constrained core and transition radius. xv

FIG. 9. Leo A’s correlated distributions of free parameters. As can be seen the core radius and transition radius are well deﬁned despite the ﬂat input priors, indicating a reliable result. xvi

FIG. 10. Sextans: correlated distributions of free parameters. As can be seen the core radius and transition radius are well deﬁned despite the ﬂat input priors, indicating a reliable result.