The v' (0-1)- v" (0-80) vibrational bands of b 1 Σ g , 0 + + → X 3 Σ g , 1 (M s = ± 1) - and b 1 Σ g , 0 + + →A3Π u,Ω=1 systems are investigated, which are in the region of 1281-336371 nm (7803-29.73cm−1) and 1616-311056nm (6190-32.15 cm−1), respectively. We finally predict that the v' (0-1)- v" (0-12) vibronic bands of b - X 3 Σ g , 1 - transition could be observed experimentally. [Display omitted] • The ground X 3 Σ g - , singlet excited b 1 Σ g + , a 1Δ g , 11Π g and triplet excited A3Π u , 13Π g states are calculated utilizing high-level ab initio multi-reference configuration interaction (MRCI) method taking into account core-valence correlation (CV), the Davidson correction (+Q) and spin–orbit coupling (SOC). • Intensity of infrared band transitions from the singlet excited b 1 Σ g + and a 1Δ g states to the triplet ground X 3 Σ g - state of the Pb 2 molecule has been studied accounting the spin (S) and orbital (L) angular momentum. • The vibronic progressions (0, 1)-(0 ∼ 12) for the b 1 Σ g , 0 + + - X 3 Σ g , 1 - transition and the (0, 1)−(0 ∼ 11) series for the a 1 Δ g , 2 - X 3 Σ g , 1 - intercombination could be observed experimentally. • The electro-quadrupole a 1 Δ g , 2 - X 3 Σ g , 0 + - transition can be greatly enhanced, since it "borrows intensity" from the a 1 Δ g , 2 - b 1 Σ g , 0 + + quadrupole transition. Intensity of infrared band transitions from the singlet excited b 1 Σ g + and a 1Δ g states to the triplet ground X 3 Σ g - state of the Pb 2 molecule has been studied utilizing high-level ab initio multi-reference configuration interaction (MRCI) method taking into account core-valence correlation (CV), the Davidson correction (+Q) and spin–orbit coupling (SOC) effect. Intensity of the magnetic dipole transition b 1 Σ g , 0 + + - X 3 Σ g , 1 - is calculated accounting the spin (S) and orbital (L) angular momentum. The former is determined by the zero-field splitting of the ground X 3 Σ g - multiplet and the SOC-induced mixing coefficient between b 1 Σ g , 0 + + and X 3 Σ g , 0 + - states; the latter is determined by the admixtures of the intermediate 3Π g and 1Π g states in the framework of SOC perturbation theory. Furthermore, the intensity of magnetic dipole a 1 Δ g , 2 - X 3 Σ g , 1 - transition only needs to take into account the orbital angular momenta interactions with magnetic wave, which arise from the a 1Δ g -1Π g and X 3 Σ g - -3Π g magnetic transitions. The weak a, b - X 3 Σ g , 1 - vibronic bands induced by magnetic dipole transition moments are easily overlapped by the stronger a − A3Π u,1,2 and b − A3Π u,1 vibronic bands of electric-dipole nature; therefore, the intensity of related a, b - A bands is also calculated. At the same time, it is concluded that the vibronic progressions (0, 1)-(0 ∼ 12) for the b - X 3 Σ g , 1 - transition and the (0, 1)−(0 ∼ 11) series for the a - X 3 Σ g , 1 - intercombination could be observed experimentally. The electro-quadrupole a 1 Δ g , 2 - X 3 Σ g , 0 + - transition can be greatly enhanced, since it "borrows intensity" from the a 1 Δ g , 2 - b 1 Σ g , 0 + + quadrupole transition. Thus, we calculate, predict, and interpret series of the weak forbidden bands in the Pb 2 infrared spectroscopy, which will extend our understanding of SOC effects in the IV series dimers of periodic system. [ABSTRACT FROM AUTHOR]