A method is presented allowing accurate computations of the weak dipole transition moments in HD. As an illustration the method is applied to the calculation of the Σu components of the wave functions. The Σu contributions to the transition moments are computed and upper bounds for the energy corrections due to Σg–Σu coupling are given. It is argued that the j dependence of the total transition moment is of the order of 10−6 D.