Charles Beck edited PH1 stable?.tex  about 9 years ago

Commit id: 95349209308a6e43450cec5f3e41a28fcf2d4581

deletions | additions      

       

\item use Holliday and REsnick's simulation results to compare PH1b's inclination and eccentricity with that of the expected stability zones.  \end{itemize}  \subsection{Semimajor Axis Inner Region}  The planet's semimajor axis is approximately 0.634 AU. The eccentricity is $e \approx 0.0539$, and the mass ratio is $\mu=m_2 /(m_1 +m_2 )=0.408/(1.528+.408)=.211$. The binary semimajor axis is $a_b =0.1744$. So, using Holman & Wiegert's best fit equation: $a_c =[0.464+(-0.380)\mu +(-.631)e+(.586)\mu e +.15e^2 +(-0.198)\mu e^2 ]a_b $ we can find the critical semimajor axis,   \begin{equation}  a_c =[0.464+(-0.380).211 +(-.631).0539+(.586)(.211) .0539 +.15(.0539^2) +(-0.198)(.211) .0539^2 ].1744=0.0622  \end{equation}.  So, $a \stackrel{?}{>} a_c$,  \begin{equation}  0.634 \stackrel{\checkmark}{>} 0.0622  \end{equation}!  This means for an inner region orbit the planet would be stable. however, the planet is orbiting both of the binary stars which means it is an outer orbit so we will next determine its stability with outer region calculations.