The energy levels and electric dipole () matrix elements of the ground state and low-lying excited states of Be atoms are calculated using the relativistic configuration interaction plus core polarization (RCICP) method. The static and dynamic , magnetic dipole () and electric quadrupole () polarizabilities as well as the hyperpolarizabilities of the and states are determined. Two magic wavelengths, 300.03 and 252.28 nm, of clock transition are found. Then, the multipolar and nonlinear Stark shifts of the clock transition at the magic wavelength are discussed in detail. We find that when the laser intensity I is in the range of 14.3–15.9 kW cm−2 and the detuning δ (the frequency detuning of the lattice laser frequency relative to the magic frequency) is in the range of 40.7–40.9 MHz, the fractional Stark shifts of the clock transition are less than 1.0 × 10−18. While, when I is in the range of 15.01–15.46 kW cm−2 and δ is in the range of 40.73–40.76 MHz, the fractional Stark shifts are lower than 1.0 × 10−19.