The bottomonium mass spectra are computed in the framework of potential non-relativistic quantum chromodynamics. The potential consists of a static term incorporating the Coulomb plus confinement part along with a correction term added non-perturbatively from pNRQCD, which is classified in powers of the inverse of the heavy quark mass O(1/m). The masses of excited bottomonia are calculated by perturbatively adding spin–hyperfine, spin–orbit and tensor components of one-gluon exchange interactions in powers of O(1/m
2). Calculated masses are found to be consistent with other theoretical studies and experimental data. The Regge trajectories of the calculated mass spectra are also constructed. The values of the wave functions are extracted and employed to calculate the electromagnetic transition widths and γγ, e
+
e
−, light hadron and γγγ decay widths of several states at various leading orders, within the non-relativistic QCD formalism. Some of the experimentally reported states of bottomonium family like ϒ(10 860), ϒ(11 020) and X(10 610) are identified as mixed S–D wave and P wave states.