Dataset of the energies and multipole radiative transition properties for Be I

We present a comprehensive dataset of energies and multipole radiative transition properties for neutral beryllium (Be I), calculated using the relativistic configuration-interaction plus many-body perturbation theory (RCI+MBPT) method. The dataset includes energies for the 93 lowest-lying levels arising from the 2sns n=2-8, 2snp n=2-8), 2snd n=3-7, 2snf n=4-7, and 2p2 electronic configurations, along with reduced matrix elements (RMEs), line strengths, transition rates, and oscillator strengths for 464 electric dipole (E1), 402 electric quadrupole (E2), 37 magnetic dipole (M1), and 546 magnetic quadrupole (M2) transitions. Calculations were performed in both length and velocity gauges for E1 transitions to assess uncertainties, with accuracies classified per NIST guidelines. Comparisons with NIST Atomic Spectra Database and prior theoretical/experimental values show excellent agreement. Many higher-level energies and multipole transitions are reported for the first time. This dataset supports applications in astrophysics, atomic clocks, and precision spectroscopy. All data are provided in tabular form in the manuscript and supplementary information.