  In this paper, we comprehensively explore bottomonia mass spectra and their
decay properties by solving the non-relativistic Schrodinger wave equation
numerically with approximate quark-antiquark potential form. We also
incorporate spin-dependent terms - spin-spin, spin-orbit, and tensor terms to
remove mass degeneracy and to obtain excited states ($nS, nP, nD, nF, n = 1, 2,
3, 4, 5$) mass spectra. By using Van Royen - Weisskopf formula, we investigate
leptonic decay constants, di-leptonic, di-gamma, tri-gamma, di-gluon decay
widths and also incorporate first-order radiative corrections. We also computed
radiative transition widths, which give a better insight into the
non-perturbative aspects of QCD. The present results for mass spectroscopy and
decay properties are in tune with available experimental values and other
theoretical predictions. Our results may provide better insight to upcoming
experimental information in the near future.
