Data for figures in paper "The velocity field of our Milky Way outer stellar halo based on DESI DR2"

This repository provides all data that can be used to reproduce figures in following paper based on DESI Milky Way Survey Data Release 2: https://ui.adsabs.harvard.edu/abs/2026arXiv260401628L/abstract ================================================================================= Description of contents (which can also be found in a readme.txt file after decompressing the tar file): figure1.csvIt stores the MCMC chain for making the contour plot (figure 1 in the paper). Meaning of columns are stored in the same order as the x-axis quantities of figure 1 in the paper.:#1: <V_r> [km/s]#2: <V_\Phi> [km/s]#3: <V_\theta> [km/s]#4: \sigma_r [km/s]#5: \sigma_\phi [km/s]#6: \sigma_\theta [km/s]#7: \mu_[Fe/H]#8: \sigma_[Fe/H]#9: log likelihood values------------------------------------------------------------------------------figure2.npyThis dataset can be loaded with:import numpy as npt=np.load('figure2.npy',allow_pickle=True) It contains the following data:r_GC - distance to Galactic center in unit of kpc<vr>, <v_phi>, <v_theta> - three mean velocity components in unit of km/s<sigma_r>, <sigma_phi>, <sigma_theta> - velocity dispersions for three velocity components in unit of km/s<beta> - velocity anisotropy<beta_with_sgr> - velocity anisotropy with Sagittarius included-------------------------------------------------------------------------------figure3.csv MCMC chains to make the contour plot (figure 3) of the paper. The meaning of columns are stored in the same order as the x-axis quantities of figure 3 in the paper. The last column is the log likelihood values. --------------------------------------------------------------------------------figure4.npyThis dataset can be loaded with:import numpy as npt=np.load('figure4.npy',allow_pickle=True) It contains the following data:r_GC - distance to Galactic center [kpc]P_MR - fraction of metal-rich population, which is the curve in the upper panel of figure 4miu_fe/h_MR - mean Fe/H of metal-rich populationmiu_fe/h_MP - mean Fe/H of metal-poor populationsigma_fe/h_MR - std in Fe/H of metal-rich populationsigma_fe/h_MP - std in Fe/H of metal-poor population ------------------------------------------------------------------------------figure5.npyThis dataset can be loaded with:import numpy as npt=np.load('figure5.npy',allow_pickle=True) It contains the following data:r_GC - distance to Galactic center [kpc]<v_r_MR>, <v_phi_MP>, <v_theta_MR> - three mean velocity components of the metal-rich population [km/s]<vr_MP>, <v_phi_MP>, <v_theta_MP> - three mean velocity components of the metal-poor population [km/s]sigma_r_MR,sigma_phi_MR,sigma_theta_MR - velocity dispersions for the three components of the metal-rich population [km/s]sigma_r_MP,sigma_phi_MP,sigma_theta_MP - velocity dispersions for the three components of the metal-poor population [km/s]beta_MP,beta_MR - velocity anisotropy of the metal-rich and metal-poor populations------------------------------------------------------------------------------figure6.npyThis dataset can be loaded with:import numpy as npt=np.load('figure6.npy',allow_pickle=True) It contains the following data:r_GC - distance to Galactic center [kpc]<vr_north> - mean radial-velocity profile of the northern stellar halo<vr_south> - mean radial-velocity profile of the southern stellar halo ----------------------------------------------------------------------------------figure7.csv MCMC chains to make the contour plot (figure 7) of the paper. The meaning of columns are stored in the same order as the x-axis quantities of figure 7 in the paper. The last column is the log likelihood values. -----------------------------------------------------------------------------------figure8.npy This dataset can be loaded with:import numpy as npt=np.load('figure8.npy',allow_pickle=True) It stores the healpix values to make the observational mean radial velocity field in Mollweide projection, with nside=6 (left panel of figure 8). ----------------------------------------------------------------------------------figure9.npy This dataset can be loaded with:import numpy as npt=np.load('figure9.npy',allow_pickle=True) It contains the following data:r_GC - distance to Galactic center [kpc]miu_compr - cotracting motion [km/s] miu_dipole - mean polar velocity [km/s]l_apex, b_apex - apex direction of the dipole velocity field in Galactic coordinates