Herein, we investigated the thermodynamic stability and opto-electronic properties of a newly BN-doped azulene. The gas-phase formation enthalpies of 11 BN-doped azulene were calculated by the atomization energy method using three computational models (CBS-APNO, CBS-QB3, and G3MP2). The results suggest that AZ-1N9B exhibits the highest stability among the studied isomers. On the other hand, AZ-1B9N and AZ-9B10N display nearly equal stability with relative energies of 19.36 and 19.82 kcal/mol at CBS-QB3, respectively. These two isomers are considered the least stable among the investigated compounds. The frontier molecular orbitals (FMO), ionization energies (IE), and electron affinities (EA) of these isomers were discussed. Additionally, the electronic absorption spectra of the BN-doped azulenes were computed using the TD-B3LYP/6–31 + G(d,p) and TD-CAM-B3LYP level of theories, which using a long-range corrected hybrid functional in acetone. The computational results obtained in this research are align closely with the existing literature, thereby reinforcing the credibility and reliability of our findings.