Using the complete active space self-consistent field (CASSCF) method followed by the internally contracted multi-reference configuration interaction (MRCI) approach in combination with the correlation-consistent basis sets, this paper studies the potential energy curves of X -S states of BeBr molecule and the corresponding 30 states for the first time. All the -S states correlate to the first two dissociation channels, Be( 1 S g ) + Br( 2 P u ) and Be( 3 P u ) + Br( 2 P u ), of BeBr molecule. Of these -S states, the 3 2 and 2 4 are found to be repulsive without the spin-orbit coupling, whereas 1 4 , 2 4 , 3 2 and 2 4 + are found to be repulsive with the spin-orbit coupling included. A 2 and 2 2 + possess the double well whether the spin-orbit coupling effect is included or not. Only 1 4 + , 1 4 − , 1 2 and 2 2 are found to be the inverted -S states. The spin-orbit coupling is accounted for by the state interaction approach with Breit-Pauli Hamiltonian using the all-electron cc-pCVTZ basis set. The potential energy curves determined by the internally contracted MRCI method are corrected for size-extensivity errors by means of the Davidson correction. Core-valence correlation correction is calculated with a cc-pCVTZ basis set. Scalar relativistic correction is included using the third-order Douglas-Kroll Hamiltonian approximation at the level of cc-pVTZ basis set. The spectroscopic parameters of all the -S and bound states are evaluated. The spectroscopic parameters are compared with those reported in the literature. Fair agreement is found between the present results and available measurements. In particular, the energy splitting of 204.43 cm −1 in the A 2 -S state agrees well with the measurements of 201 cm −1 . Analyses demonstrate that the spectroscopic parameters reported here can be expected to be reliably predicted ones.