The third (PNOF3), fourth (PNOF4) and fifth (PNOF5) versions of the Piris natural orbital functional were used to characterize the beryllium dimer. The results obtained were compared to those gained afforded by CASSCF and CASPT2 as well as experimental data. The equilibrium distances (R (e)), dissociation energies (D (e)), effective bond orders (EBOs), and rovibrational levels were calculated. PNOF3, PNOF4, and CASPT2 predicted a bonded Be(2) molecule, while PNOF5 and CASSCF did not, which demonstrates the importance of the dynamical electron correlation. We observed that PNOF3 yields the most accurate equilibrium distances, while PNOF4 most accurately calculates the rovibrational levels. However, both of these functionals overestimate dissociation energies. Both PNOF3 and PNOF4 predict EBOs that agree with that obtained using CASPT2.