We compute the electric dipole transitions χ bJ (1P ) → γΥ(1S), with J = 0, 1, 2, and h b (1P ) → γη b (1S) in a model-independent way. We use potential non-relativistic QCD (pNRQCD) at weak coupling with either the Coulomb potential or the complete static potential incorporated in the leading order Hamiltonian. In the last case, the perturbative series shows very mild scale dependence and a good convergence pattern, allowing predictions for all the transition widths. Assuming ΛQCD mv 2 , the precision that we reach is k 3where kγ is the photon energy, m is the mass of the heavy quark and v its relative velocity. Our results are: Γ(χ b0 (1P ) → γΥ(1S)) = 28 +2 −2 keV, Γ(χ b1 (1P ) → γΥ(1S)) = 37 +2 −2 keV, Γ(χ b2 (1P ) → γΥ(1S)) = 45 +3 −3 keV and Γ(h b (1P ) → γη b (1S)) = 63 +6 −6 keV.