This paper comprehensively explores bottomonia mass spectra and their decay properties by solving the non-relativistic Schrodinger wave equation numerically with an 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. Using the Van Royen—Weisskopf formula, we investigate leptonic decay constants, di-leptonic, di-gamma, tri-gamma, di-gluon decay widths and incorporate first-order radiative corrections. We also computed radiative transition widths, which provide better insight into the non-perturbative aspects of QCD. The states ϒ(10750), ϒ(10860), ϒ(11020) are analyzed as pure states 33 D 1, 53 S 1 and 43 D 1, respectively. The present mass spectroscopy and decay properties results align with available experimental values and other theoretical predictions. Our results may provide better insight into upcoming experimental information in the near future.