Two body non-leptonic decays of

Heavy flavor mesons serve as crucial probes of the Quark-Gluon Plasma (QGP) in heavy-ion collisions at RHIC and LHC facilities. Precise theoretical predictions of their decay properties in vacuum are essential baseline measurements for understanding medium modifications in hot QCD matter. Using the factorization approach, this work presents a comprehensive study of two-body nonleptonic decays of B, and Bs mesons. We calculated branching fractions for a few decay channels, employing relativistic quark model form factors and mass values of parent particles with Hydrogen-like and Gaussian-like wavefunctions as input. Our results demonstrate good agreement with PDG data and existing theoretical predictions for most channels, validating the factorization framework's effectiveness. The Gaussian wavefunction approach shows accuracy, with branching fraction predictions aligning well with experimental values like PDG. Discrepancies in some channels reflect known limitations of the factorization approximation and highlight the importance of final-state interactions. Our results provide important theoretical benchmarks for collision data and guiding future advancements in computational modelling within the field of QCD matter.