The origin of very metal-poor (VMP; Fe/H -2. 0) stars on planar orbits has been the subject of great attention since their first discovery. However, prior to the release of the Gaia BP/RP (XP) spectra, and large photometric samples such as SkyMapper, SAGES, J-PLUS and S-PLUS, most studies have been limited due to their small sample sizes or strong selection effects. Here, we cross-match photometric metallicities derived from Gaia XP synthetic photometry and geometric distances from Bailer-Jones et al. , and select 12, 000 VMP stars (1604 dwarfs and 10, 396 giants) with available high-quality astrometry. After calculating dynamical parameter estimates using AGAMA, we employ the non-negative matrix factorization technique to the v_ϕ distribution across bins in Z ₌₀ₗ (the maximum height above or below the Galactic plane during the stellar orbit). We find three primary populations of the selected VMP stars: halo, disk system, and the Gaia Sausage/Enceladus (GSE) structure. The fraction of disk-like stars decreases with increasing Z ₌₀ₗ (as expected), although it is still 20\% for stars with Z ₌₀ₗ 3 kpc. Similar results emerge from the application of the Hayden criterion, which separates stellar populations on the basis of their orbital inclination angles relative to the Galactic plane. We argue that such high fractions of disk-like stars indicate that they are an independent component, rather than originating solely from Galactic building blocks or heating by minor mergers. We suggest that most of these VMP stars are members of the hypothesized ``primordial" disk.
Xu et al. (Sun,) studied this question.