Buckling analysis of nanobeams is investigated using nonlocal continuum beam models of the different classical beam theories namely as Euler-Bernoulli beam theory (EBT), Timoshenko beam theory (TBT), and Levinson beam theory (LBT). To this end, Eringen’s equations of nonlocal elasticity are incorporated into the classical beam theories for buckling of nanobeams with rectangular cross-section. In contrast to the classical theories, the nonlocal elastic beam models developed here have the capability to predict critical buckling loads that allowing for the inclusion of size effects. The values of critical buckling loads corresponding to four commonly used boundary conditions are obtained using state-space method. The results are presented for different geometric parameters, boundary conditions, and values of nonlocal parameter to show the effects of each of them in detail. Then the results are fitted with those of molecular dynamics simulations through a nonlinear least square fitting procedure to find the appropriate values of nonlocal parameter for the buckling analysis of nanobeams relevant to each type of nonlocal beam model and boundary conditions.analysis.