Following the approach of Hughes, Factor, Chiang et al., we model an eccentric circumstellar disk as a series of apse-aligned ellipses. We model the mass distribution within each ellipse using a mass-on-a-wire approach. In this document, we derive the mass distribution for a single ellipse as a function of orbital position.

To start, we assume that the \(i\)th ellipse we have chosen to model has a total mass of \(m_i\). The linear density of this wire, \(\lambda_i\), with units [g/cm], is proportional to the amount of time spent in that section of the orbit. This means that the density will be largest at apoapse, when a test particle in the orbit is traveling the slowest. Conversely, the density will be smallest at periapse, when the test particle is moving fastest. The linear density is uniquely defined as a function of orbital position. We can either measure the orbital position in terms of path length along the perimeter of the ellipse, \(s\), where \(s=0\) denotes periapse, or, we can measure orbital position as a function of true anomaly^{1} \(\theta\), which is the same as the angle \(f\) in Meredith’s notes.

Once we have defined \(\lambda_i\), we can check that we have the right definition by integrating the linear density over the perimeter of the elliptical orbit–we should recover \(m_i\).

After making the assumption that the linear density should be proportional to the amount of time spent in the orbit, we can use the equations of orbital motion to solve for \(\lambda_i\). The amount of time \(dt\) spent in a section of the wire \(ds\) is the inverse of the velocity at that orbital position

\[\frac{dt}{ds} = \frac{1}{v}\]

This same section of the wire \(ds\) also contains a small chunk of mass \(dm\), which provides the definition of linear density

\[\lambda_i = \frac{dm}{ds}\]

Because the elliptical orbit is closed, we know the normalization terms for each of these quantities. The total mass in the ring is \(m_i\) and the total amount of time spent in the orbit is the period

\[T = \frac{2 \pi}{\mu^2} \frac{h^3}{(1 - e^2)^{3/2}}\]

where \(e\) is the eccentricity of the wire, \(\mu = G M_\ast\) is gravitational parameter, and \(h\) is the angular momentum. Therefore, the fractional amount of time spent in the section of the wire should be equal to the fractional amount of mass contained in that section \[\frac{dt/ds}{T} = \frac{dm/ds}{m_i} = \frac{\lambda_i}{m_i}\]

Solving for \(\lambda_i\) (units of [**g/cm**]), we have

\[\lambda_i = m_i \frac{dt/ds}{T} = \frac{m_i}{T} \left(\frac{1}{v} \right)\]

We can find the velocity using orbital equations for conservation of mechanical energy for an ellipse

\[\frac{v^2}{2} - \frac{\mu}{r} = \frac{-\mu}{2 a}\]

where \(a\) is the semi-major axis of the ellipse,

\[a = \frac{h^2}{\mu} \frac{1}{1 - e^2}\]

we can also find the radial distance \(r\) of the orbital position \[r = a \frac{1 - e^2}{1 + e \cos \theta}\]

Using these equations and the equation for the period \(T\), we can solve for \(\lambda_i\). First, we rearrange the mechanical energy equation for velocity \[v = \sqrt{\mu \left ( \frac{2}{r} - \frac{1}{a} \right )}\]

and plug in for \(r\) \[v = \sqrt{\mu \left ( \frac{2 ( 1 + e \cos \theta)}{a (1 -e^2)} - \frac{1}{a} \right )} = \sqrt{\frac{\mu}{a}} \sqrt{\frac{1 + 2 e \cos \theta + e^2}{1 - e^2}}\]

Before we plug in for the period \(T\), let’s first simplify by using the equation for semi-major axis \(a\) to eliminate \(h\) \[h = \sqrt{ a \mu (1 - e^2)}\] and so we have \[T = \frac{2 \pi}{\mu^2} \frac{a^{3/2} \mu^{3/2} (1 - e^2)^{3/2}}{(1 -e^2)^{3/2}} = \frac{2 \pi a^{3/2}}{\sqrt{\mu}}\] Finally, solving for \(\lambda\) \[\lambda_i = \frac{m_i}{T} \left(\frac{1}{v} \right) = m_i \frac{\sqrt{\mu}}{2 \pi a^{3/2}} \frac{\sqrt{a} \sqrt{1 - e^2}}{\sqrt{\mu} \sqrt{1 + 2 e \cos \theta + e^2}}\] \[\boxed{ \lambda_i = \frac{m_i}{2 \pi a} \frac{\sqrt{1 - e^2}}{\sqrt{1 + 2 e \cos \theta + e^2}} \quad \quad [\textrm{g} / \textrm{cm}] }\]

The equation we just derived is *slightly* different from what is shown in Meredith’s notes \[\lambda_i = \frac{m_i}{2 \pi a} \frac{\sqrt{1 - e^2}}{1 + e \cos \theta}\]

If I had an extra factor of \(cos^2(\theta)\) attached to the \(e^2\) in the denominator of my equation, I would be able to factor the quadratic equation to get the same answer. In order to check that I didn’t make a mistake somewhere, we should be able to integrate the linear density over the full circumference of the orbital ellipse to recover \(m_i\), the total mass of the elliptical wire.

\[\label{eqn:mass} m_i == \int_0^S \lambda_i ds\]

where \(S\) is the total path length of the ellipse perimeter. Because \(\lambda\) is a function of \(\theta\), and integrating around the circumference would require an elliptic integral, it will be easier to integrate this expression if we make a change of variables from \(s\) to \(\theta\). This requires us to derive the differential \(ds\) in terms of \(d\theta\).

Knowledge of the velocity gives us \(v = ds / dt\). The rate of change in true anomaly is related to the angular momentum \[\frac{d \theta}{d t} = \frac{h}{r^2}\] combining these two expressions gives \[\frac{d \theta}{d s} = \frac{d \theta/ dt}{d s / dt} = \frac{h}{r^2} \frac{\sqrt{a}}{\sqrt{\mu}} \frac{\sqrt{1 - e^2}}{\sqrt{1 + 2 e \cos \theta + e^2}}\] we can rearrange and substitute the angular momentum equation to simplify \[\frac{h}{\sqrt{\mu}} = \sqrt{a (1 - e^2)}\] \[\frac{d \theta}{d s} = \frac{a ( 1 - e^2)}{r^2 \sqrt{1 + 2 e \cos \theta + e^2}}\] substituting for \(r\), we have \[\frac{d \theta}{d s} = \frac{(1 + e \cos \theta)^2}{a (1 -e^2) \sqrt{1 + 2 e \cos \theta + e^2}}\] or \[\boxed{ ds = \frac{a ( 1 -e^2) \sqrt{1 + 2 e \cos \theta + e^2}}{(1 + e \cos \theta)^2} d \theta }\] Now, Equation \ref{eqn:mass} becomes \[m_i == \int_0^S \lambda_i ds = \int_0^{2\pi} \lambda_i \left ( \frac{ds}{d\theta} \right ) d \theta\] \[m_i == \int_0^{2\pi} \frac{m_i}{2 \pi a} \frac{\sqrt{1 - e^2}}{\sqrt{1 + 2 e \cos \theta + e^2}} \frac{a ( 1 -e^2) \sqrt{1 + 2 e \cos \theta + e^2}}{(1 + e \cos \theta)^2} d \theta\]

\[m_i == \frac{m_i(1 - e^2)^{3/2}}{2 \pi} \int_0^{2 \pi} \frac{ d \theta}{(1 + e \cos \theta)^2}\]

We solve this integral using equation 3.661.4 from Gradshteyn et al. (2007),

\[m_i == \frac{m_i(1 - e^2)^{3/2}}{2 \pi} \left ( \frac{2 \pi}{(1 - e^2)^{3/2}}\right )\]

\[m_i == m_i\]

The equation Meredith quotes has units of linear density, or [g/cm]. \[\lambda_i = \frac{m_i}{2 \pi a} \frac{\sqrt{1 - e^2}}{1 + e \cos \theta}\]

Therefore, if we apply the same integration test we should be able to reproduce \(m_i == m_i\) Now, Equation \ref{eqn:mass} becomes \[m_i == \int_0^S \lambda_i ds = \int_0^{2\pi} \lambda_i \left ( \frac{ds}{d\theta} \right ) d \theta\] \[m_i == \int_0^{2\pi} \frac{m_i}{2 \pi a} \frac{\sqrt{1 - e^2}}{1 + e \cos \theta} \frac{a ( 1 -e^2) \sqrt{1 + 2 e \cos \theta + e^2}}{(1 + e \cos \theta)^2} d \theta\] \[m_i == \int_0^{2\pi} \frac{m_i}{2 \pi} \frac{(1 - e^2)^{3/2}}{(1 + e \cos \theta)^3} \sqrt{1 + 2 e \cos \theta + e^2} d \theta\] Unfortunately, I cannot find an analytic solution for this integral, but doing the integral numerically reveals that this equation does not hold true for \(e > 0.0\).

Instead, we can make the assumption that Meredith quotes \(\lambda_i\) as a *radial density* and see if the math works out. Unfortunately, the expression has units of [g/cm] from the \(m_i/a\) term, if it were radial density then it should have units of [g/radian] and \(a\) should not appear in the expression. Regardless, we can try integrating \[m_i == \int_0^{2\pi} \lambda_i d \theta = \frac{m_i \sqrt{1 - e^2}}{2 \pi a} \int_0^{2 \pi} \frac{d \theta}{1 + e \cos \theta}\] Using the same formula from (Gradshteyn 2007), \[m_i \ne \frac{m_i \sqrt{1 - e^2}}{2 \pi a} \left ( \frac{2 \pi}{\sqrt{1 - e^2}} \right ) = \frac{m_i}{a}\].