Comment on: Discontinuous codimension-two bifurcation in a Vlasov equation (arXiv:2212.01250)

We comment on the recent work by Yamaguchi and Barr\'e [Phys. Rev. E 107, 054203 (2023)], which uses linear stability analysis of the Vlasov equation to characterize phase transitions in a generalized Hamiltonian Mean Field (gHMF) model. By performing extensive molecular dynamics simulations with $N=10^8$ particles, we demonstrate that the bifurcation analysis of the initial stationary distribution is insufficient to predict either the location or the nature of the phase transition to a quasi-stationary state (qSS). Specifically, we show that for bimodal momentum distributions, the instability threshold identified by the authors does not correspond to a ferromagnetic transition; instead, the system remains in a paramagnetic state characterized by magnetization oscillations with a zero time-average. We find that the true paramagnetic-ferromagnetic transition is discontinuous (first-order) and occurs at significantly larger coupling strengths, characterized by a clear coexistence of states. These results indicate that linear bifurcation and symmetry-breaking phase transitions are distinct phenomena in long-range interacting systems, and that the former lacks the predictive power to describe the long-time fate of the system.