In this talk, we focus on the fast computation of stationary states of a class of phase field crystal models using the Bregmanian proximal gradient (BPG) methods. In the full discretization formulation, it is proved that the sequence generated by the BPG method satisfies the energy dissipation property. Moreover, we extend the BPG to the infinite-dimensional case and analyze its convergence rate using the Lojasicwicz-Simon inequality. Finally, experiments on LB, LP and OK models validate the effectiveness of the proposed method.