Electronic states of the PbSi molecule up to 4 eV have been studied by carrying out ab initio based MRDCI calculations which include relativistic effective core potentials (RECPs) of both the atoms. The use of semicore RECPs of Pb produces better dissociation limits than the full-core one. However, the (3)P(0)-(3)P(1) splitting due to Pb is underestimated by about 4000 cm(-1). At least 25 bound electronic states of the Λ-S symmetry are predicted for PbSi. The computed zero-field-splitting in the ground state is about 544 cm(-1). A strong spin-orbit mixing changes the nature of the potential energy curves of many Ω states. The overall splitting among the spin components of A(3)Π is computed to be 4067 cm(-1). However, the largest spin-orbit splitting is reported for the (3)Δ state. A number of spin-allowed and spin-forbidden transitions are predicted. The partial radiative lifetime for the A(3)Π-X(3)Σ(-) transition is of the order of milliseconds. The computed bond energy in the ground state is 1.68 eV, considering the spin-orbit coupling. The vertical ionization energy for the ionization to the X(4)Σ(-) ground state of PbSi(+) is about 6.93 eV computed at the same level of calculations.