Our calculations were conducted within density functional theory (DFT) and density functional perturbation theory (DFPT) using norm-conserving pseudo-potential and the local density approximation. The elastic constants of Zn1−xBexOZn1−xBexO were calculated, C11C11, C33C33 and C44C44 increase with the increase of Be content, whereas the C12C12 shows a non-monotonic variation and C13C13 decreases when Be concentration increases. The values of bulk modulus BB, Young’s modulus EE and shear modulus GG increase with the increase of Be content. Poisson’s ratio σσ decreases with increased Be concentration. The ductility decreases with increasing Be concentration and the compressibility for Zn1−xBexOZn1−xBexO along cc-axis is smaller than along aa-axis. Phonon dispersion curves show that Zn1−xBexOZn1−xBexO is dynamically stable (no soft modes). Quantities such as refractive index, Born effective charge, dielectric constants and optical phonon frequencies were calculated as a function of the Be molar fraction xx. The agreement between the present results and the known data that are available only for ZnO and BeO is generally satisfactory. Our results for Zn1−xBexOZn1−xBexO(0<x<1)(0<x<1) are predictions.