The Blanco DECam Bulge Survey (BDBS) VIII. Chemo-kinematics in the southern Galactic bulge from 2.3 million red clump stars with Gaia DR3 proper motions

Context. The inner Galaxy is a complex environment, and the relative contributions of di fferent formation scenarios to its observed morphology and stellar properties are still debated. The di fferent components are expected to have di fferent spatial, kinematic, and metallicity distributions, and a combination of photometric, spectroscopic, and astrometric large-scale surveys is needed to study the formation and evolution of the Galactic bulge. Aims. The Blanco DECam Bulge Survey (BDBS) provides near-ultraviolet to near-infrared photometry for approximately 250 million unique stars over more than 200 square degrees of the southern Galactic bulge. By combining BDBS photometry with the latest Gaia astrometry, we aim to characterize the chemodynamics of red clump stars across the BDBS footprint using an unprecedented sample size and sky coverage. Methods. Our field of view of interest is |l| <= 10 degrees, -10 degrees <= b -3 degrees. We constructed a sample of approximately 2 :3 million red clump giants in the bulge with photometric metallicities, BDBS photometric distances, and proper motions. Photometric metallicities are derived from a (u- i)(0) versus [Fe=H] relation; astrometry, including precise proper motions, is from the third data release (DR3) of the ESA satellite Gaia. We studied the kinematics of the red clump stars as a function of sky position and metallicity by investigating proper-motion rotation curves, velocity dispersions, and proper-motion correlations across the southern Galactic bulge. Results. By binning our sample into eight metallicity bins in the range of 1.5 dex < [Fe/H] < +1 dex, we find that metal-poor red clump stars exhibit lower rotation amplitudes, at similar to 29 km s 1 kpc 1. The peak of the angular velocity is similar to 39 km s 1 kpc 1 for [Fe/H] similar to 0.2 dex, exhibiting declining rotation at higher [Fe/H]. The velocity dispersion is higher for metal-poor stars, while metalrich stars show a steeper gradient with Galactic latitude, with a maximum dispersion at low latitudes along the bulge minor axis. Only metal-rich stars ([Fe/H] >= 0.5 dex) show clear signatures of the bar in their kinematics, while the metal-poor population exhibits isotropic motions with an axisymmetric pattern around Galactic longitude l = 0. Conclusions. This work describes the largest sample of bulge stars with distance, metallicity, and astrometry reported to date, and shows clear kinematic di fferences with metallicity. The global kinematics over the bulge agrees with earlier studies. However, we see striking changes with increasing metallicity, and, for the first time, kinematic di fferences for stars with [Fe/H] > -0.5, suggesting that the bar itself may have kinematics that depends on metallicity.