<p>Before I tackle this problem, I will assume that the size of the planet is negligible compared to the size of the star, so I won't take into account partial transits.  I will also assume circular orbits.  I will also assume uniform random probability distribution of orbital revolution axes on the unit sphere.</p>
<p>Suppose the orbital revolution axis of the planet is pointed directly at us, the observer.  Then, we would be perfectly perpendicular to the orbital plane and there would be no transit.  Instead, from our perspective, the planet's orbit would be perfectly circular at a distance of <span class="math-container">$4d$</span> from the star.</p>
<p>Only when the orbital revolution axis (along r in the below image) was some angle <span class="math-container">$\theta$</span> from our line of site (along h), would the planetary orbit cross the star from our perspective. The revolution axes where no transit exists could be projected onto a unit sphere as a <a href="https://en.wikipedia.org/wiki/Spherical_cap" rel="nofollow noreferrer">spherical cap</a>, with area <span class="math-container">$A = 2\pi (1-\cos{\theta})$</span>. The area of a unit hemisphere is <span class="math-container">$A_h=2\pi$</span>, so the fractional area of no transits (which is the probability of no transits) is simply <span class="math-container">$P(\sim t)=A/A_h=1-\cos{\theta}$</span>.
<a href="https://en.wikipedia.org/wiki/Spherical_cap#/media/File:Spherical_cap_diagram.tiff" rel="nofollow noreferrer">Image credit wikipedia</a></p>
<p><a href="https://i.sstatic.net/NAv8K.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/NAv8K.png" alt="enter image description here" /></a></p>
<p>The angle <span class="math-container">$\theta$</span> is simply the complementary angle of the maximum transit angle, so  <span class="math-container">$\theta=\cos^{-1}{(1/8)}$</span>. So <span class="math-container">$P(\sim t)=1-1/8$</span>, and the probability of a transit is then <span class="math-container">$P(t) = 1- P(\sim t)=\frac{1}{8}$</span>.</p>
<p><strong>Note:</strong> The OP claims the answer should instead be <span class="math-container">$P(t) = \frac{1}{8\pi}$</span>. Please let me know if you find my missing <span class="math-container">$\pi$</span>!</p>

CC BY-SA 4.0
2022-01-17T19:29:46.100Z
"The angle θ is simply the complementary angle of the maximum transit angle,..." is difficult to follow. The "maximum transit", if it is defined as the transit with longest time, would it not be with θ =90°? in your convention for θ. In this case the non-transit cap is defined by that max θ that results in the minimum transit time (the planet just touches the star disk w/o entering it). Did I miss something?
Ng Ph
CC BY-SA 4.0
2022-01-18T16:34:04.673Z
@NgPh Sorry I was late with this reply.  In the figure in your answer, my $\theta$ is the same as $i$, Since cosine is the Adjacent over Hypotenuse of a right triangle, $\cos{\theta} = 1/8$.  The maximum transit angle is the (complementary) unlabeled interior angle.
Connor Garcia
CC BY-SA 4.0
2022-01-13T21:26:17.097Z
2022-01-16T15:00:19.303Z
2022-01-16T15:00:19.303Z
Connor Garcia
<p>This is an amusing <strong>(and tough)</strong> math quiz.</p>
<p>We are confronted with 3 numerical results and various mathematical reasonings (of which one is not given):</p>
<ul>
<li>
<ol>
<li>the OP's calculation gives <span class="math-container">$Pt=0.8/4\pi$</span> (about 6.4%).</li>
</ol>
</li>
<li>
<ol start="2">
<li>the OP’s announced correct result <span class="math-container">$Pt=1/8\pi$</span> (about 4%), w/o demonstration or source.</li>
</ol>
</li>
<li>
<ol start="3">
<li>@Connor Garcia calculation gives <span class="math-container">$Pt=1/8$</span> (12.5%).</li>
</ol>
</li>
<li>
<ol start="4">
<li>The tutorial linked by <a href="http://www.iiassvietri.it/down/ases_2015/Lectures_notes/Tr_1.pdf" rel="nofollow noreferrer">@ProfRob, Slide 7</a>, gives <span class="math-container">$Pt=R/a$</span>=(d/2)/(4d)=1/8</li>
</ol>
</li>
</ul>
<p><strong>Definition of Pt</strong></p>
<p>Given an Observer(O) at a far distance from a Star (S) of radius R, having a much smaller  Exoplanet orbiting at distance <span class="math-container">$a$</span>, the probability that the Observer can detect the Exoplanet by the transit method is denoted Pt. Note that the duration of transit, the luminosity of the Star, the observation time and the sensitivity of the detector are not considered. It is purely a geometrical exercise.</p>
<p>For example, when the orbital plane contains the line between the Observer and the Star (O-S), the transit time is maximum. We can call these orbits « maximum transit orbits » (for that Observer position). When the orbital plane is such that the exoplanet just touches the disk of the Star, viewed from the Observer, the transit time is minimal (near zero). We call these « minimum transit orbits ». Both types of orbits count equally as detectable in the calculation of Pt. When O-S coincides with the axis of rotation, the Observer will never see a transit.</p>
<p><strong>The reasonings</strong></p>
<p>The reasoning in (3) computes (1-Pt) as the ratio of the areas of two spherical caps: the &quot;non-transit cap&quot; and the hemisphere cap facing the observer, both on a unit sphere.</p>
<p>The reasoning in (1) computes Pt as the ratio of arc lengths: the arc of maximum transit and the half circumference, both of a circle of radius a (=4d).</p>
<p>The reasoning in (4) computes Pt as the probability Prob{abs(cos(i)) &lt; <span class="math-container">$R/a$</span>}, <span class="math-container">$i$</span> being the angle that the orbital plane makes with the <strong>perpendicular</strong> to O-S (see Fig).</p>
<p><a href="https://i.sstatic.net/Jb70w.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/Jb70w.png" alt="enter image description here" /></a>
<a href="http://www.iiassvietri.it/down/ases_2015/Lectures_notes/Tr_1.pdf" rel="nofollow noreferrer">source</a></p>
<p><strong>Walking through the reasoning in (3)</strong></p>
<p>Take a perfect sphere and make the following « experiment » (I owe this intuitive and nice reasoning to @ConnorGacia):</p>
<ul>
<li><p>pick a random point on the sphere,</p>
</li>
<li><p>paint in blue a cap of half-cone angle <span class="math-container">$\theta$</span> centered at the chosen point.
<a href="https://i.sstatic.net/W1BX9.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/W1BX9.png" alt="enter image description here" /></a></p>
</li>
<li><p>Paint also the « opposite » cap (180° symmetric).</p>
</li>
<li><p>Let the sphere roll on a <strong>bumpy</strong> surface so that, when it comes to rest, <strong>any point on the sphere has equal opportunity</strong> of being the rest point of contact with the surface. The chance that the sphere comes at rest at a point inside any of the two blue areas is equal to the ratio of one cap area to the area the sphere, multiplied by 2. This is the same as taking area of one cap and divide it by the area of the hemisphere.
When the sphere has a radius r=1, the area of each blue cap is <span class="math-container">$2\pi(1-cos \theta)$</span> and the hemisphere has an area equals to <span class="math-container">$2\pi$</span>. Hence the probability of not coming to rest in the blue area is <span class="math-container">$cos  \theta$</span>.</p>
</li>
</ul>
<p>The connection to our Pt is as follows :</p>
<p>Call the angle between the line O-S and the axis of orbit rotation, <span class="math-container">$\delta$</span>. Call <span class="math-container">$\theta$</span> that value of <span class="math-container">$\delta$</span> resulting in minimum transit orbits (in the direction of the given Observer). If we increase <span class="math-container">$\delta$</span> from 0° (the orbit is orthogonal to the O-S), then the orbit becomes a transiting one (from the Observer’s point of view), when and only when <span class="math-container">$\delta$</span> crosses <span class="math-container">$\theta$</span>.</p>
<p>For any given R and a, it turns out that the value of <span class="math-container">$\theta$</span> is related only to R and a in the following way <span class="math-container">$cos \theta = sin (90°- \theta) = R/a$</span>.</p>
<p>For a given orbit, consider a (any) sphere centered at the star and draw on it the two blue caps, each cap is centered on the orbit axis of rotation. Choosing a random direction for O-S is akin to rolling that sphere on a bumpy surface. The Observer can detect the planet if the O-S does not intersect any of the 2 blue caps.</p>
<p><strong>Dissecting and comparing (1) and (4)</strong></p>
<p>Both (1) and (4) reasonings are based on a reduction to a « 2D » geometry. Both make use, more or less, of the same angles.</p>
<p>In (1), it is reasoned that, when we increase the inclination angle (called <span class="math-container">$\theta/2$</span> in the figure), from 0° (maximum transit), the orbit leaves the transit zone when <span class="math-container">$sin(inclination)=R/a$</span>.</p>
<p>In (4), it is reasoned that when we reduce the angle <span class="math-container">$i$</span> (see Fig) from 90° (maximum transit), the orbit leaves the transit zone when <span class="math-container">$i$</span> reaches the value that satisfies <span class="math-container">$cos i =R/a$</span>.</p>
<p>They are both correct up to this point, but diverges on a (very subtle) mathematical point. In (1), it is assumed obvious that <strong>the inclination angle is uniformly distributed</strong>, whereas in (4) it is <strong><span class="math-container">$\cos i$</span> that is uniformly distributed</strong>. We can suspect that the calculation of (4) is correct because the result of (4) matches that of (3). But why is it that saying that the inclination is uniformly distributed is flawed ?</p>
<p><strong>Another ball rolling experiment</strong></p>
<p>I will resort again to the intuitive (and nice) « experiment » with the rolling sphere above (again, I am indebted to @ConnorGarcia for this). This time we do not paint the sphere. Instead, once an orbit is chosen, we identify the Equator and the Poles and we draw the « latitude» markers, from Equator to Poles. We roll the ball again on our bumpy surface and when it comes to rest, we note the latitude of the resting point. This latitude is the inclination angle of the orbit with respect to the (random) observation line O-S, which is the outcome of our experiment run. If we repeat several runs of this, we can observe that we get more inclination below 45° than inclination above 45°. Another test is to count how many times we get a rest at, or close to any Pole, and how many time we get, at or close to the Equator. We can observe that we get more often closer to the Equator than to the Pole. This shows that the random inclination angle <strong>is not</strong> uniformly distributed.</p>
<p>The reasoning in (1) however is based on this false evidence and therefore its result is incorrect.</p>
<p>If you like this answer, please upvote @Connor Garcia too.</p>

CC BY-SA 4.0
2022-01-18T16:38:34.037Z
Related: https://astronomy.stackexchange.com/questions/35471/why-are-asteroids-with-zero-orbital-inclination-rare/40518#40518
Connor Garcia
CC BY-SA 4.0
2022-01-18T16:53:24.453Z
I think a key takeaway from this answer is that the probability distribution of orbital inclinations should not be treated as uniform random from 0 to 180.  Instead, one should consider the distribution of orbital axes on the unit sphere.
Connor Garcia
CC BY-SA 4.0
2022-01-18T18:30:04.397Z
@Connor Garcia, neatly formulated. Nevertheless, as much as I am convinced I have identified **where** (1) is wrong, I am still unable to find an intuitive way to answer the "why this is wrong?". Especially when the 2 functions $1-cos \theta$ and $2\theta/\pi$ agree at the two extremes of $\theta=0$ and $\theta=\pi/2$. Need more coffee !
Ng Ph
CC BY-SA 4.0
2022-01-18T22:22:43.560Z
As an analogy, imagine if one was to roll a globe (spherical earth model) down a bumpy hill.  It would be much more likely that it came to rest on an absolute latitude less than 45 degrees rather than greater.  This is because equi-latitude strips cover more of the globe closer to the equator.  But if you throw a random dart at a Mercator Projection map of the Earth, the opposite is true!  Which technique would you want to use to ensure that equal sized areas have equal probability of a sample landing in them?
Connor Garcia
CC BY-SA 4.0
2022-01-19T12:02:26.827Z
@Connor Garcia, I think you are taking me closer to a "AHA! moment". I may be able to connect the dots to the related question (your 1st comment post). But I haven't fully got the connection with this problem, yet. Please, do not give another hint before a couple of days. I like math quiz!
Ng Ph
CC BY-SA 4.0
2022-01-20T21:18:09.273Z
@Connor Garcia, AHA! Please check my edits.
Ng Ph
CC BY-SA 4.0
2022-01-17T22:12:32.067Z
2022-01-20T21:17:19.857Z
2022-01-20T21:17:19.857Z
Ng Ph
<p>I'm trying to calculate the probability of transit of an exoplanet in orbit around a star of diameter <span class="math-container">$d$</span>, at a distance of <span class="math-container">$4d$</span> from its centre. Here's what I did:</p>
<p><a href="https://i.sstatic.net/N97n5.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/N97n5m.png" alt="enter image description here" /></a></p>
<p><a href="https://i.sstatic.net/O7SsB.png" rel="nofollow noreferrer"><img src="https://i.sstatic.net/O7SsBm.png" alt="enter image description here" /></a></p>
<p>In this figure, we are only concerned about the right hemisphere. The idea I have used is that for a transit to be visible for an observer situated directly to the right, the orbit must extend to within <span class="math-container">$l$</span> of the right extreme of this sphere.</p>
<p>To calculate <span class="math-container">$l$</span>, I applied basic trigonometry:</p>
<p><span class="math-container">$$ \frac{\theta}{2}=tan^{-1}\frac{1}{8},$$</span></p>
<p>and <span class="math-container">$\theta$</span> came out to be roughly <span class="math-container">$0.2$</span> radians.</p>
<p><span class="math-container">$$l = 4d * \theta \approx 0.8 \ d.$$</span></p>
<p>Now the probability is simply given by <span class="math-container">$\large\frac{0.8 \ d}{4 \pi d}$</span> and this comes out to be about <span class="math-container">$\large\frac{1}{5\pi}$</span>.</p>
<p>Where have I gone wrong, because the answer to this problem is supposed to be <span class="math-container">$\large\frac{1}{8\pi}$</span>?</p>
<p>Even after approximating <span class="math-container">$l$</span> to be nearly <span class="math-container">$d$</span>, the answer I get is <span class="math-container">$\large\frac{1}{4\pi}$</span> and not <span class="math-container">$\large\frac{1}{8\pi}$</span>.</p>

CC BY-SA 4.0
2022-01-13T08:17:23.500Z
Please edit the image to make it smaller, I don't know how to do that.
Ambica Govind
CC BY-SA 4.0
2022-01-13T08:31:31.767Z
Your two images seem to disagree about the definition of d. In the upper it is l = d. In the lower it is $l = 4d\cdot \theta / (2\pi) \approx 2d$. Is the lower sketch is supposed to be perpendicular to the top view?
planetmaker
CC BY-SA 4.0
2022-01-13T09:10:25.800Z
@planetmaker Thanks for pointing it out. I have edited it.
Ambica Govind
CC BY-SA 4.0
2022-01-13T10:05:09.607Z
@AmbicaGovind I've made some edits, using a trick to display the images smaller and adjusted the MathJax formatting. Please have a look and feel free to edit further.
uhoh
CC BY-SA 4.0
2022-01-13T10:43:31.300Z
I haven't throught this through. However there is a probability for a central passage to happen right now (orbit inclination of 0°), but the probability of this has to be multiplied with the probability of the orbits inclination being such that it transits at all due to inclination (I recon that could be $\theta$ in the lower sketch). I don't see this being taken care of in the calculation.
planetmaker
CC BY-SA 4.0
2022-01-13T15:35:13.047Z
I can't see what the diagrams are meant to represent. You are trying to work out the probability that the star is obscured by a disc of radius $d/2$ when viewed from effectively infinitely far away. This involves tangent lines coming from infinity which graze the planet and the limb of the star. There are lots of internet resources that derive this e.g. http://www.iiassvietri.it/down/ases_2015/Lectures_notes/Tr_1.pdf
ProfRob
CC BY-SA 4.0
2022-01-17T19:59:21.577Z
I think you meant "an exoplanet in orbit around a **star** (not planet).
Ng Ph
<exoplanet><mathematics><planetary-transits><planetary-science>
Probability of transit of an exoplanet in front of its star
Ambica Govind
# Probability of transit of an exoplanet in front of its star

(Asked by: username_0 on 2022-01-13)

I'm trying to calculate the probability of transit of an exoplanet in orbit around a star of diameter $d$, at a distance of $4d$ from its centre. Here's what I did:

[![enter image description here](https://i.sstatic.net/N97n5m.png)](https://i.sstatic.net/N97n5.png)

[![enter image description here](https://i.sstatic.net/O7SsBm.png)](https://i.sstatic.net/O7SsB.png)

In this figure, we are only concerned about the right hemisphere. The idea I have used is that for a transit to be visible for an observer situated directly to the right, the orbit must extend to within $l$ of the right extreme of this sphere.

To calculate $l$, I applied basic trigonometry:

$$ \frac{\theta}{2}=tan^{-1}\frac{1}{8},$$

and $\theta$ came out to be roughly $0.2$ radians.

$$l = 4d \* \theta \approx 0.8 \ d.$$

Now the probability is simply given by $\large\frac{0.8 \ d}{4 \pi d}$ and this comes out to be about $\large\frac{1}{5\pi}$.

Where have I gone wrong, because the answer to this problem is supposed to be $\large\frac{1}{8\pi}$?

Even after approximating $l$ to be nearly $d$, the answer I get is $\large\frac{1}{4\pi}$ and not $\large\frac{1}{8\pi}$.

**Tags:** exoplanet, mathematics, planetary-transits, planetary-science

> **Comments:**
> **username_0**: Please edit the image to make it smaller, I don't know how to do that.  
> **username_1**: Your two images seem to disagree about the definition of d. In the upper it is l = d. In the lower it is $l = 4d\cdot \theta / (2\pi) \approx 2d$. Is the lower sketch is supposed to be perpendicular to the top view?  
> **username_0**: @username_1 Thanks for pointing it out. I have edited it.  
> **username_2**: @username_0 I've made some edits, using a trick to display the images smaller and adjusted the MathJax formatting. Please have a look and feel free to edit further.  
> **username_1**: I haven't throught this through. However there is a probability for a central passage to happen right now (orbit inclination of 0°), but the probability of this has to be multiplied with the probability of the orbits inclination being such that it transits at all due to inclination (I recon that could be $\theta$ in the lower sketch). I don't see this being taken care of in the calculation.  
> **username_3**: I can't see what the diagrams are meant to represent. You are trying to work out the probability that the star is obscured by a disc of radius $d/2$ when viewed from effectively infinitely far away. This involves tangent lines coming from infinity which graze the planet and the limb of the star. There are lots of internet resources that derive this e.g. http://www.iiassvietri.it/down/ases_2015/Lectures_notes/Tr_1.pdf  
> **username_4**: I think you meant "an exoplanet in orbit around a **star** (not planet).  

---

## 2 answers

### Answer by username_5 on 2022-01-13. Score: 2

Before I tackle this problem, I will assume that the size of the planet is negligible compared to the size of the star, so I won't take into account partial transits. I will also assume circular orbits. I will also assume uniform random probability distribution of orbital revolution axes on the unit sphere.

Suppose the orbital revolution axis of the planet is pointed directly at us, the observer. Then, we would be perfectly perpendicular to the orbital plane and there would be no transit. Instead, from our perspective, the planet's orbit would be perfectly circular at a distance of $4d$ from the star.

Only when the orbital revolution axis (along r in the below image) was some angle $\theta$ from our line of site (along h), would the planetary orbit cross the star from our perspective. The revolution axes where no transit exists could be projected onto a unit sphere as a [spherical cap](https://en.wikipedia.org/wiki/Spherical_cap), with area $A = 2\pi (1-\cos{\theta})$. The area of a unit hemisphere is $A\_h=2\pi$, so the fractional area of no transits (which is the probability of no transits) is simply $P(\sim t)=A/A\_h=1-\cos{\theta}$.
[Image credit wikipedia](https://en.wikipedia.org/wiki/Spherical_cap#/media/File:Spherical_cap_diagram.tiff)

[![enter image description here](https://i.sstatic.net/NAv8K.png)](https://i.sstatic.net/NAv8K.png)

The angle $\theta$ is simply the complementary angle of the maximum transit angle, so $\theta=\cos^{-1}{(1/8)}$. So $P(\sim t)=1-1/8$, and the probability of a transit is then $P(t) = 1- P(\sim t)=\frac{1}{8}$.

**Note:** The OP claims the answer should instead be $P(t) = \frac{1}{8\pi}$. Please let me know if you find my missing $\pi$!

> **Comments:**
> **username_4**: "The angle θ is simply the complementary angle of the maximum transit angle,..." is difficult to follow. The "maximum transit", if it is defined as the transit with longest time, would it not be with θ =90°? in your convention for θ. In this case the non-transit cap is defined by that max θ that results in the minimum transit time (the planet just touches the star disk w/o entering it). Did I miss something?  
> **username_5**: @username_4 Sorry I was late with this reply.  In the figure in your answer, my $\theta$ is the same as $i$, Since cosine is the Adjacent over Hypotenuse of a right triangle, $\cos{\theta} = 1/8$.  The maximum transit angle is the (complementary) unlabeled interior angle.  

---

### Answer by username_4 on 2022-01-17. Score: 1

This is an amusing **(and tough)** math quiz.

We are confronted with 3 numerical results and various mathematical reasonings (of which one is not given):

* 1. the OP's calculation gives $Pt=0.8/4\pi$ (about 6.4%).
* 2. the OP’s announced correct result $Pt=1/8\pi$ (about 4%), w/o demonstration or source.
* 3. @username_5 calculation gives $Pt=1/8$ (12.5%).
* 4. The tutorial linked by [@username_3, Slide 7](http://www.iiassvietri.it/down/ases_2015/Lectures_notes/Tr_1.pdf), gives $Pt=R/a$=(d/2)/(4d)=1/8

**Definition of Pt**

Given an Observer(O) at a far distance from a Star (S) of radius R, having a much smaller Exoplanet orbiting at distance $a$, the probability that the Observer can detect the Exoplanet by the transit method is denoted Pt. Note that the duration of transit, the luminosity of the Star, the observation time and the sensitivity of the detector are not considered. It is purely a geometrical exercise.

For example, when the orbital plane contains the line between the Observer and the Star (O-S), the transit time is maximum. We can call these orbits « maximum transit orbits » (for that Observer position). When the orbital plane is such that the exoplanet just touches the disk of the Star, viewed from the Observer, the transit time is minimal (near zero). We call these « minimum transit orbits ». Both types of orbits count equally as detectable in the calculation of Pt. When O-S coincides with the axis of rotation, the Observer will never see a transit.

**The reasonings**

The reasoning in (3) computes (1-Pt) as the ratio of the areas of two spherical caps: the "non-transit cap" and the hemisphere cap facing the observer, both on a unit sphere.

The reasoning in (1) computes Pt as the ratio of arc lengths: the arc of maximum transit and the half circumference, both of a circle of radius a (=4d).

The reasoning in (4) computes Pt as the probability Prob{abs(cos(i)) < $R/a$}, $i$ being the angle that the orbital plane makes with the **perpendicular** to O-S (see Fig).

[![enter image description here](https://i.sstatic.net/Jb70w.png)](https://i.sstatic.net/Jb70w.png)
[source](http://www.iiassvietri.it/down/ases_2015/Lectures_notes/Tr_1.pdf)

**Walking through the reasoning in (3)**

Take a perfect sphere and make the following « experiment » (I owe this intuitive and nice reasoning to @ConnorGacia):

* pick a random point on the sphere,
* paint in blue a cap of half-cone angle $\theta$ centered at the chosen point.
  [![enter image description here](https://i.sstatic.net/W1BX9.png)](https://i.sstatic.net/W1BX9.png)
* Paint also the « opposite » cap (180° symmetric).
* Let the sphere roll on a **bumpy** surface so that, when it comes to rest, **any point on the sphere has equal opportunity** of being the rest point of contact with the surface. The chance that the sphere comes at rest at a point inside any of the two blue areas is equal to the ratio of one cap area to the area the sphere, multiplied by 2. This is the same as taking area of one cap and divide it by the area of the hemisphere.
  When the sphere has a radius r=1, the area of each blue cap is $2\pi(1-cos \theta)$ and the hemisphere has an area equals to $2\pi$. Hence the probability of not coming to rest in the blue area is $cos \theta$.

The connection to our Pt is as follows :

Call the angle between the line O-S and the axis of orbit rotation, $\delta$. Call $\theta$ that value of $\delta$ resulting in minimum transit orbits (in the direction of the given Observer). If we increase $\delta$ from 0° (the orbit is orthogonal to the O-S), then the orbit becomes a transiting one (from the Observer’s point of view), when and only when $\delta$ crosses $\theta$.

For any given R and a, it turns out that the value of $\theta$ is related only to R and a in the following way $cos \theta = sin (90°- \theta) = R/a$.

For a given orbit, consider a (any) sphere centered at the star and draw on it the two blue caps, each cap is centered on the orbit axis of rotation. Choosing a random direction for O-S is akin to rolling that sphere on a bumpy surface. The Observer can detect the planet if the O-S does not intersect any of the 2 blue caps.

**Dissecting and comparing (1) and (4)**

Both (1) and (4) reasonings are based on a reduction to a « 2D » geometry. Both make use, more or less, of the same angles.

In (1), it is reasoned that, when we increase the inclination angle (called $\theta/2$ in the figure), from 0° (maximum transit), the orbit leaves the transit zone when $sin(inclination)=R/a$.

In (4), it is reasoned that when we reduce the angle $i$ (see Fig) from 90° (maximum transit), the orbit leaves the transit zone when $i$ reaches the value that satisfies $cos i =R/a$.

They are both correct up to this point, but diverges on a (very subtle) mathematical point. In (1), it is assumed obvious that **the inclination angle is uniformly distributed**, whereas in (4) it is **$\cos i$ that is uniformly distributed**. We can suspect that the calculation of (4) is correct because the result of (4) matches that of (3). But why is it that saying that the inclination is uniformly distributed is flawed ?

**Another ball rolling experiment**

I will resort again to the intuitive (and nice) « experiment » with the rolling sphere above (again, I am indebted to @username_5 for this). This time we do not paint the sphere. Instead, once an orbit is chosen, we identify the Equator and the Poles and we draw the « latitude» markers, from Equator to Poles. We roll the ball again on our bumpy surface and when it comes to rest, we note the latitude of the resting point. This latitude is the inclination angle of the orbit with respect to the (random) observation line O-S, which is the outcome of our experiment run. If we repeat several runs of this, we can observe that we get more inclination below 45° than inclination above 45°. Another test is to count how many times we get a rest at, or close to any Pole, and how many time we get, at or close to the Equator. We can observe that we get more often closer to the Equator than to the Pole. This shows that the random inclination angle **is not** uniformly distributed.

The reasoning in (1) however is based on this false evidence and therefore its result is incorrect.

If you like this answer, please upvote @username_5 too.

> **Comments:**
> **username_5**: Related: https://astronomy.stackexchange.com/questions/35471/why-are-asteroids-with-zero-orbital-inclination-rare/40518#40518  
> **username_5**: I think a key takeaway from this answer is that the probability distribution of orbital inclinations should not be treated as uniform random from 0 to 180.  Instead, one should consider the distribution of orbital axes on the unit sphere.  
> **username_4**: @username_5, neatly formulated. Nevertheless, as much as I am convinced I have identified **where** (1) is wrong, I am still unable to find an intuitive way to answer the "why this is wrong?". Especially when the 2 functions $1-cos \theta$ and $2\theta/\pi$ agree at the two extremes of $\theta=0$ and $\theta=\pi/2$. Need more coffee !  
> **username_5**: As an analogy, imagine if one was to roll a globe (spherical earth model) down a bumpy hill.  It would be much more likely that it came to rest on an absolute latitude less than 45 degrees rather than greater.  This is because equi-latitude strips cover more of the globe closer to the equator.  But if you throw a random dart at a Mercator Projection map of the Earth, the opposite is true!  Which technique would you want to use to ensure that equal sized areas have equal probability of a sample landing in them?  
> **username_4**: @username_5, I think you are taking me closer to a "AHA! moment". I may be able to connect the dots to the related question (your 1st comment post). But I haven't fully got the connection with this problem, yet. Please, do not give another hint before a couple of days. I like math quiz!  
> **username_4**: @username_5, AHA! Please check my edits.  

---

