Spatial (in the RA and Dec plane referred to the system center, left panel) and kinematic (in the RV-PM plane referred to the systemic velocities, right panel) distributions of the 39 stars that are likely members of Liller 1, color-coded according to their metal content: red (metal-poor), green (metal-intermediate), and blue (metal-rich). The 1, 2, and 3 σ contour levels from the GMM analysis are also plotted.