We look at the mass spectra of the [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], and [Formula: see text] mesons using a relativistic square root potential. Before looking at the mass spectra, we have to figure out the model parameters, which are [Formula: see text][Formula: see text]GeV and [Formula: see text][Formula: see text]GeV. The calculated result of [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], [Formula: see text], findings of this study exhibit a notable concurrence with the experimental observations and pertinent theoretical projections. We estimate the decay constant, leptonic decay width, semileptonic decay width, and branching fractions of pseudoscalar and vector mesons, specifically B and D mesons, while keeping the model parameters unchanged. The pseudoscalar decay constants and partial decay widths of “B and D-mesons” reasonably agree with the theoretical predictions, lattice quantum chromodynamics (LQCD) calculations, and experimental data. Moreover, we have efficiently found the values for these mesons’ leptonic decay width and branching fraction, matching the experimental findings and theoretical forecasts. The calculated values of semileptonic decays are [Formula: see text], [Formula: see text], [Formula: see text] and [Formula: see text], [Formula: see text], [Formula: see text], the proximity of the observed results to experimental and certain theoretical models is evident.