Проблема с астрофизическими единицами

надеюсь, у вас все хорошо.

Я не очень хорошо знаком AstroPy, и у меня возникла проблема с вычислением преобразования между единицами углового момента .

Если я определяю абс. значение углового момента

$$\|\vec{J_2}\| \propto \|\vec{r}\times \vec{p}\| \overset{!}{\propto} a_2 m_2 \left(\frac{G M_s}{a_2}\right)^{1/2} = m_2 (G M_s a_2)^{1/2}$$

от движения тела$m_2$взаимодействие с$M_s$, предполагая кеплерово движение, я напишу:

затем общий пресс. анг. мама.$\|J_1\|$, скажем

$$\|J_1\| \equiv \mathrm{\|((1m) \times (1 \frac{kg \, m}{s}))\| = 1 \frac{kg\, m}{s}}$$ J1 = u.m * u.kg * u.m * (1/u.s)должны иметь эквивалентные размеры по сравнению с$J_2$. Поэтому я должен быть в состоянии проверить, используя

J2.decompose()или J2.to(J1), но все, что я получаю, этоUnitConversionError: 'AU(1/2) earthMass G(1/2) solMass(1/2)' and 'kg m2 / s' (angular momentum) are not convertible

J1 = angular_momentum(1* u.Msun, 3*u.Mearth, 0.1*u.au)

J2 = u.m * u.kg * u.m * (1/u.s)

J1.to(J2)
---> UnitConversionError

или

angular_momentum(1*u.Msun, 3*u.Mearth, 0.1*u.au).decompose() #why?

$$3.0900743 \times 10^{43} \mathrm{\frac{kg^2 \,m^{1/2}}{A^{1/2}\,s}} \qquad \mathrm{(output)}$$

это похоже * u.auне работает, как ожидалось. В любом случае, передача юнитов, например, u.mне решает проблему.

Вопрос: кто-нибудь знает, что здесь происходит?

Изменить: это функция, которую я использую для вычисления$J_2$def angular_momentum (мс, мп, а):

    #from Lp = r x p ~~ a * (mp (G ms / a)^1/2) keplerian rot ~~ mp (g ms a)^1/2
    
    Lp = mp * (u.G  *  ms  *  a)**(1./2.)   #This MUST! have L units
    return Lp