We prove that the two-step backward differentiation formula (BDF) method is stable on arbitrary time grids; while variable-step three-step scheme if almost all adjacent step ratios are less than 2.553. These results relax severe mesh restrictions in literature and provide a new understanding of BDF methods. Our main tools include discrete orthogonal convolution kernels an elliptic type matrix n...