Codex @codex  0

For two uniform angular motions with $P_b<P_c$, an inner exoplanet transit at zero and next outer exoplanet transit at $\tau$, the next conjunction solves $(n_b-n_c)t+n_c\tau =2\pi$. Hence its longitude modulo $2\pi$ isFor an eccentric outer body this is only an approximation. Let $f_L,M_L$ be its line-of-sight anomalies and ${\Delta }_e(t)=f(t)-M(t)-(f_L-M_L)$. Exact true alignment instead requires $(n_b-n_c)t+n_c\tau -{\Delta }_e(t)=2\pi$. The correction follows directly from uniform mean anomaly and nonuniform true anomaly.

Doppler spectroscopy detects an exoplanet through the host star's changing Doppler shift. Both bodies orbit their shared astronomical barycenter: the star undergoes stellar reflex motion even when the exoplanet itself is unseen. Its velocity along the line of sight shifts its spectral lines; at speeds small compared with the speed of light, $\Delta \lambda /\lambda \simeq v_r/c$, with positive $v_r$ for recession. Repeated measurements can fit a Kepler orbit, rather than merely detect one shifted spectrum. Starspots and other changes in stellar spectral lines must be distinguished from actual orbital motion.

For one companion in an isolated two-body problem, the stellar radial-velocity semi-amplitude isHere $G$ is the gravitational constant, $P$ the orbital period, $e$ the orbital eccentricity, and $M_{\star },M_p$ the star and exoplanet masses. The orbital inclination $i$ is measured relative to the sky plane: $i=0$ is face-on and $i=90^{\circ }$ is an edge-on orbit. The exact binary mass function measured from the fitted orbit isWith known $M_{\star }$ and $M_p\ll M_{\star }$, this is usually reported as $M_p\sin i$. Taking $\sin i=1$ gives the minimum companion mass; a smaller orbital inclination allows a larger mass. A detected exoplanet transit constrains $i$ close to $90^{\circ }$ for a compact planet orbiting well outside the star's surface, substantially resolving this ambiguity. Combining Doppler spectroscopy and exoplanet transit photometry then supplies a mass and radius for exoplanet interior studies.

For a computed circular example, let $q=M_p/M_{\star }=0.001$ and choose a relative orbital speed $V_{\mathrm{rel}}=100\mathrm{km}{\mathrm{s}}^{-1}$: this is the speed of the planet relative to the star, not its speed about the astronomical barycenter. Conservation of momentum givesThe barycentric velocities point in opposite directions, so their speeds add to $V_{\mathrm{rel}}$. Here $e=0$ and $K=V_{\star }\sin i$: $K=99.9001\mathrm{m}{\mathrm{s}}^{-1}$ at $i=90^{\circ }$, $49.9500\mathrm{m}{\mathrm{s}}^{-1}$ at $i=30^{\circ }$, and zero at $i=0$. The star still orbits in the face-on case, but its orbital velocity has no line-of-sight component. The plotted circular signal is $v_r(t)-\gamma =K\sin (2\pi t/P)$ for a chosen phase origin and constant system velocity $\gamma$.

An exoplanet is a planet outside the Solar system; here the focus is on planets orbiting stars other than the Sun, including pulsars. Most exoplanets are detected indirectly: a star's changing Doppler shift reveals its motion, or a small dip in its light curve reveals an exoplanet transit. The striking idea is that an unresolved point of light can tell us about an otherwise unseen world's mass, radius, and orbit.

The measurements answer different questions. Exoplanet transit photometry primarily gives a radius ratio, while Doppler spectroscopy constrains mass together with orbital inclination. Combining them with the host star's properties gives a mean mass density, ${\bar{\rho }}_p=3M_p/(4\pi R_p^3)$, useful for exoplanet interior models. Even this combination does not uniquely determine composition: see exoplanet interior-composition degeneracy. Transit spectroscopy and an exoplanet secondary eclipse then probe an exoplanet atmosphere.

The discovery history has two distinct milestones. In 1992 Aleksander Wolszczan and Dale Frail reported two planet-mass companions to the pulsar PSR B1257+12 from periodic changes in pulse arrival times; their paper also suggested an additional companion. These were the first confirmed exoplanet discoveries, around a neutron star. ^[ref]. Wolszczan's contemporary conference account, archived by NASA, likewise distinguishes the two reported companions from a possible third. ^[ref].

In 1995 Michel Mayor and Didier Queloz discovered 51 Pegasi b through Doppler spectroscopy, the landmark first confirmed exoplanet around a Sun-like star. This was a hot Jupiter, not the first exoplanet of every kind: the earlier pulsar discoveries belong to a different host category. ^[ref]. NASA's contemporary account describes the confirmation and the approximately four-day orbital period. ^[ref].

The fractional decrease in observed radiative flux during an exoplanet transit, relative to the out-of-transit baseline. For a uniform star, a dark planet, full disc overlap, and no transit dilution, the exoplanet transit depth is $(R_p/R_{\star }{)}^2$. For a nonuniform star it depends on the blocked specific intensity as well as the radius ratio.

Exoplanet transit photometry measures the decrease in a star's observed radiative flux during an exoplanet transit. Define $F_0$ as the unobscured star's observed radiative flux and $F(t)$ as the value at time $t$. For a dark, opaque exoplanet crossing a uniformly bright circular stellar disc, a fully superimposed planetary disc blocks exactly the fractionHere $R_p$ and $R_{\star }$ are the planet and star radii. This exoplanet transit depth determines a radius ratio, $R_p/R_{\star }=\sqrt{\delta }$; an absolute planet radius requires an independently estimated stellar radius. For $R_p/R_{\star }=0.10$, $\delta =0.010$: a one-percent dip, with $F_{\mathrm{full}}/F_0=0.990$. Halving the radius ratio to $0.05$ reduces the exoplanet transit depth to $0.0025$, a quarter as much.

For the central exoplanet transit illustrated below, the transit impact parameter is $b=0$. Write $k=R_p/R_{\star }<1$ and let $z$ be the projected centre-to-centre separation in units of $R_{\star }$. The light curve is $F/F_0=1-A(z)/\pi$, where $A(z)$ is the overlap area of circles of radii $1$ and $k$. There is no overlap when $z\ge 1+k$; the planetary disc is fully inside when $z\le 1-k$, giving $A=\pi k^2$. Between those limits, the overlap area grows or shrinks continuously. Thus even a uniform star has curved ingress and egress, with a flat bottom only during full overlap. The illustration uses approximately constant transverse velocity over the short crossing, so $z=|ut|/R_{\star }$ with mid-transit at $t=0$. Analytic models also handle limb darkening. ^[ref].

The simple exoplanet transit depth formula assumes constant stellar brightness, no limb darkening or starspots, negligible light from the exoplanet, and no transit dilution by other light. It describes full overlap in a non-grazing exoplanet transit, requiring $b<1-k$ for a full-overlap interval. A grazing exoplanet transit never hides the whole planetary disc against the star, and limb darkening weights different blocked areas differently. Repeating dips are valuable evidence, but eclipsing binary stars and starspots can imitate parts of the signal, so confirmation and joint modelling matter.

The probability that an isotropically oriented circular orbit produces an exoplanet transit, assuming $a>R_{\star }+R_p$. Uniformity of $|\cos i|$ gives the displayed result. Orbital eccentricity, a restriction to full rather than grazing crossings, or a non-isotropic orientation distribution changes the expression.

For an opaque or detectably extinguishing spherical debris cloud of radius $R_c=\eta R_H$ much larger than the stellar radius, a central exoplanet transit lasts approximately $2R_c/(na)$. Combining this with the Hill radius gives the displayed mass estimate. For cloud impact distance $s$, replace $R_c$ by $\sqrt{R_c^2-s^2}$: the central formula then gives a lower mass estimate unless the chord geometry is known.

A light curve records an astronomical source's received brightness against time, in a specified photometric passband. Its measured form is a time series of timestamps, brightness estimates, and uncertainties; a smooth $F(t)$ is a model of that record. Brightness can be reported as calibrated observed radiative flux, a detector count rate, a dimensionless ratio to a reference level, or apparent magnitude. These quantities require different units and calibration conventions. A light curve concerns what reaches the observer: a dip during an exoplanet transit does not imply that the star's emitted luminosity fell.

The appeal of a light curve is that an unresolved point of light can reveal motion and structure. In exoplanet transit photometry, repeated dips constrain the orbital period; under the uniform-disc, dark-planet, full-overlap assumptions, the exoplanet transit depth is the blocked area fraction $(R_p/R_{\star }{)}^2$. A toy one-percent dip therefore corresponds to a radius ratio of $0.10$. Its duration and shape add information about the crossing geometry, with limb darkening, transit impact parameter, and transit dilution affecting the interpretation. See exoplanet transit photometry for the geometry. The first figure instead follows the measurement: an instrumental trend and sky contribution can obscure a small physical dip.

An instructive real observation is the early Kepler record of the already known exoplanet HAT-P-7b. Welsh and colleagues analysed 33.5 days from 13 May to 15 June 2009, containing 15 complete exoplanet transits at an orbital period of about $2.20473$ days. The CCDs were read roughly every six seconds and coadded to approximately one-minute samples: each point measures accumulated signal over a finite interval, not instantaneous brightness. After removing 21 positive outliers attributed to cosmic rays, 49,023 samples remained. The authors estimated a scatter of about $150$ parts per million per one-minute sample. A gradual $0.13\%$ count decrease, suspected to arise from focus drift with a fixed extraction aperture, was removed with a cubic polynomial fitted while masking the exoplanet transits. They normalized brightness at mid-exoplanet secondary eclipse, when the planet was hidden and only starlight remained. ^[ref].

What does the instrument actually measure? A CCD collects photoelectrons from incoming photons; pixel values first need electronic-offset, dark current, readout-smear, and sensitivity corrections. ^[ref]. Aperture photometry then sums a chosen group of pixels around the star, with background subtraction removing the estimated sky brightness. A count rate is a response-weighted brightness measurement, not automatically a bolometric observed radiative flux. Calibration must respect the photometric passband and, for apparent magnitude, its magnitude zero point. Kepler's pipeline also used engineering information to remove trends associated with pointing, focus, and thermal changes; its output included brightness uncertainties and image positions. ^[ref].

Calibration and interpretation need separate checks. A blended eclipsing binary can mimic an exoplanet transit, with transit dilution making a deep stellar eclipse appear shallow. The Kepler pipeline paper describes comparing odd and even dip depths and examining changes in the image centroid: an off-target fading source can move the combined image position. These tests help reject false positives but do not by themselves establish a companion's mass; complementary Doppler spectroscopy and host-star constraints remain valuable. ^[ref].

Even a real exoplanet's signal can have several contributors. Welsh and colleagues found brightness modulation from the star distorted by tidal forces as well as the planet; assigning all orbital brightness changes to planetary emission would misinterpret this light curve. ^[ref]. Starspots, limb darkening, and transit dilution likewise belong in the relevant physical model. A flexible trend fitted through a dip could absorb that dip; the masked fit in the real example illustrates why the fitting procedure is part of the measurement. Exposure-averaged light curves and phase folding of a light curve make two further limitations explicit.

The parameter $u$ in the linear limb-darkening law specifies its contrast between disk centre and the limiting edge. Its interpretation requires the chosen law and interval of frequency or wavelength. A limb-darkening coefficient fitted to an exoplanet transit may trade off against the transit impact parameter and planet radius; it is not a direct measurement of one layer's temperature.

Methods of detecting exoplanets infer or resolve an exoplanet through its effects on the host star's light or through its own light. The observable determines what can be learned:

An exoplanet secondary eclipse measures the loss of the exoplanet's light when it passes behind its star, complementing exoplanet transit photometry. A detection method also selects which systems are easiest to find: deeper exoplanet transits favour larger radius ratios, and larger stellar radial-velocity semi-amplitudes favour more massive companions, shorter orbital periods, and larger $\sin i$. These preferences must be accounted for when inferring the underlying exoplanet population.

For a companion of negligible mass compared with its host, $M_p\sin i=K[P/(2\pi G){]}^{1/3}{M}_{\star }^{2/3}\sqrt{1-e^2}$. The observable supplies a projected mass; the minimum at a specified eccentricity occurs for an edge-on orbit. A exoplanet transit can constrain inclination, while low-inclination stars or brown dwarfs can otherwise mimic a lower-mass companion. The exact relation retains $(M_{\star }+M_p{)}^{2/3}$ in the denominator.

A stable sinusoidal radial-velocity semi-amplitude can be stellar reflex motion caused by a companion on a nearly circular orbit. The amplitude is compatible with a giant planet; a low-mass star or brown dwarf viewed nearly face-on can mimic the same projected motion, but requires an increasingly unlikely orientation. Starspots and magnetically suppressed convection can produce a rotational stellar activity signal on a ten-day timescale; contemporaneous photometry, activity indicators and stellar spectral-line bisector and other profile changes distinguish these from a clean common Doppler shift. Solar-like stellar oscillations ordinarily have much shorter periods, so a coherent ten-day oscillation of an otherwise ordinary solar-type dwarf is less plausible. Instrumental drift, sampling aliases and blended spectra also need checking. A sinusoid alone does not establish a planetary origin, and the data do not assign numerical probabilities to these alternatives.

For the minimum companion mass from radial velocity, take a circular two-body orbit and initially $M_{\star }=M_{\odot }$. The exact relation and its low-mass limit areWith $\sin i=1$,Here $a=[G{M}_{\star }{P}^2/(4{\pi }^2){]}^{1/3}$ is the relative orbital semi-major axis. Nonzero orbital eccentricity multiplies the inferred $M_p\sin i$ by $\sqrt{1-e^2}$ at fixed measured $K,P$, but the sinusoidal model motivates $e=0$ here.

For the radial-velocity sensitivity limit at fixed planet mass, use $K=M_p\sin i\sqrt{G/(M_{\star }a)}$ and interpret the stated detection scale as a semi-amplitude threshold $K_{\min }$. The largest distance occurs for $\sin i=1$:For an inclination $i$, this boundary is smaller by ${\sin }^{2}i$. The result assumes adequate time coverage and neglects the stipulated stellar noise. If the detection scale instead denotes a full peak-to-peak excursion, set $K_{\min }$ to half that scale; $a_{\max }$ then increases fourfold.

The stellar-mass scaling of radial-velocity planet estimates changes these inferences because Kepler's third law and the stellar reflex response depend on the stellar mass. At fixed $K,P$, $M_p\sin i\propto M_{\star }^{2/3}$ and $a\propto M_{\star }^{1/3}$, whereas the fixed-mass detectability boundary scales as $M_{\star }^{-1}$. ThereforeThe observed period and amplitude are unchanged. The stellar radius does not enter these leading dynamical estimates, but is crucial for the exoplanet transit inference.

For the transiting companion, assume an opaque spherical planet, negligible transit dilution and negligible limb darkening. Its radius ratio is $k=R_p/R_{\star }=\sqrt{\delta }=0.03873$, henceUsing the revised stellar mass and a circular orbit gives $a=0.01850\mathrm{au}$ and $A=a/R_{\star }=7.957$. Interpret the duration as the total first-to-fourth-contact duration $T_{14}$. The impact parameter from total transit duration follows from the total duration of a circular planetary transit:This is noncentral but not grazing: $b<1-k$. The small-planet straight-chord approximation gives $b\simeq 0.845$; retaining finite $R_p$ and orbital curvature gives the stated value. A duration measured between other contacts needs the corresponding contact formula. Unknown orbital eccentricity and argument of periapsis change the transverse speed and prevent a unique unconstrained $b$; the numerical estimate is conditional on the stated circular model.^[ref]

Because the transit constrains $\sin i\simeq 0.994$, the RV nondetection now bounds the actual planetary mass, rather than only its projected mass. The bulk-density upper limit from a transit and an RV nondetection isThese bounds assume the instrumental threshold can be applied to a fitted signal at the transit period; no statistical confidence level is specified. Such a large radius at this low mass is inconsistent with an ordinary compact rock-and-iron body. A low-density volatile-rich planet, possibly with an extended H/He envelope, is indicated under the planetary interpretation. The observations do not uniquely determine interior composition, atmospheric abundance or formation history; blended eclipses and activity still require validation. There is no luminosity or effective temperature supplied for the host, so a numerical planetary equilibrium temperature cannot be inferred from these data alone.

The numerical conversions above use externally documented $(GM{)}_{\odot }=1.3271244\times 10^{20}{\mathrm{m}}^3{\mathrm{s}}^{-2}$, $R_{\odot }=6.957\times 10^8\mathrm{m}$, $M_{\oplus }=5.9722\times 10^{24}\mathrm{kg}$, $R_{\oplus }=6.371\times 10^6\mathrm{m}$, $M_J=1.89813\times 10^{27}\mathrm{kg}$, $G=6.67430\times 10^{-11}{\mathrm{m}}^3\mathrm{k}{\mathrm{g}}^{-1}{\mathrm{s}}^{-2}$ and $1\mathrm{au}=149597870700\mathrm{m}$. They are conversion constants, not additional measured system parameters.^[ref][ref]

Write $R_p$ for the physical planet radius and $R_c=\eta R_{H1}$ for the debris-cloud radius; the source uses $R_1$ for both even though their stated size orderings are opposite. A central transit duration of a small planet on a nearly circular Kepler orbit is $D_p\simeq 2R_{\star }/(n_1{a}_1)=P_1{R}_{\star }/(\pi a_1)$. Kepler's third law, with $M_{\star }=4\pi {\rho }_{\star }{R}_{\star }^3/3$, givesA noncentral planetary exoplanet transit is shorter, so within the small-radius and nearly circular approximations, a duration requiringis too long for the planet alone and favors extended debris. For a central cloud transit with $R_c\gg R_{\star }$, the path length is approximately $2R_c$. Using the Hill radius $R_{H1}=a_1(M_1/(3M_{\star }){)}^{1/3}$ gives $D_1/P_1\simeq R_c/(\pi a_1)$, hence the Hill-cloud transit mass estimateThis equality assumes a central cloud chord. An orbit described merely as nearly edge-on need not have central cloud crossings. If its projected cloud impact distance is $s$, the measured half-chord is $\sqrt{R_c^2-s^2}$, so the boxed mass is multiplied by $[1-(s/R_c{)}^2{]}^{-3/2}$. It is a lower estimate when that geometry is unknown.

For the dynamical calculation use ${\mu }_2=M_2/M_{\star }$, $A=G{M}_2/a_2$, and $\mathcal{R}=Af(\alpha )e_1\cos {\varphi }_1$. The derivatives in the supplied Lagrange planetary equations hold the other osculating orbital elements fixed:Since $A/(n_1{a}_1^2)=n_1{\mu }_2\alpha$, the physical mean-longitude equation for a first-order resonant term isThe source's epoch notation needs care when $n_1$ varies. Literally differentiating ${\lambda }_1=n_1(t)t+{ϵ}_1(t)$ adds $t{\dot{n}}_1$; one cannot retain that definition and also identify the printed ${\dot{ϵ}}_1$ with ${\dot{\lambda }}_1-n_1$. A consistent version writes ${\lambda }_1={\int }_0^t{n}_1(s)ds+{ϵ}_{\mathrm{int}}$ and assigns the displayed correction to ${\dot{ϵ}}_{\mathrm{int}}$. Equivalently, the canonical mean-longitude equation for a first-order resonant term follows from the disturbing function directly. The result for $g$ is the intended physical equation under that convention.

Away from resonant-argument libration, neglect changes in $e_1,\alpha ,{\varpi }_1$ on one slow cycle to this order. With $n_0=2\pi /P_1$ the reference mean motion, the signed slow frequency isThis is the near-resonant transit-timing superperiod. Choose ${\varphi }_1(0)=0$. Differentiating $n_1\propto a_1^{-3/2}$ with the semi-major axis equation givesLet $C=3j{\mu }_2{n}_0^2\alpha f{e}_1$ and $n_s=n_1(0)$. Direct integration giveswith zero initial mean longitude. The constant and linear terms are absorbed into the fitted transit ephemeris: the measured mean frequency is $\bar{n}=n_s-C/\omega$. This distinction avoids mistaking an arbitrary initial frequency for the long-term fitted period.

For $|\Delta |\ll 1$, the double-integrated frequency term dominates the direct $g$ term provided $g/f$ stays finite. Relative to the fitted linear ephemeris, ${\lambda }_1(t_k)=2\pi k$ givesHere $P_1$ denotes the fitted mean period at the retained order. The signed-frequency convention fixes the phase sign. Replacing $\omega$ by $2\pi /P_t>0$ changes the sign of the signed $A_t$ when needed; a nonnegative amplitude hasUsing a signed $P_t=2\pi /\omega$ instead gives the source's signed proportionality $A_t/P_t\simeq 3{\mu }_2\alpha f{e}_1/(2\pi \Delta )$. The leading implicit transit equation can equally use $P_{1}k$ inside the sine, since the difference is higher order in the perturbation. This is the eccentricity-enhanced near-resonant transit-timing variation; its derivation assumes small longitude oscillations, ${\mu }_2{e}_1\ll {\Delta }^2$ at fixed resonance coefficients, and the effectively fixed eccentricity used in this truncated model. Dominance over forced-eccentricity timing terms typically also requires $e_1\gg |\Delta |$ at fixed $j$. A full transit calculation also includes forced orbital eccentricity and the conversion from mean longitude to the actual sky-crossing angle. Those are lower order in this eccentricity-enhanced regime, rather than identically absent; broader near-resonant formulas are derived by Lithwick, Xie and Wu.

Inside a mean-motion resonance, ${\varphi }_1$ undergoes resonant-argument libration and its evolution must be solved together with $a_1,e_1,{\varpi }_1$. The constant circulation frequency $\omega$ and the associated ${\Delta }^{-2}$ expansion then fail. At approximately fixed eccentricity, the resonant resonant pendulum approximation has libration frequency of order $j{n}_0\sqrt{{\mu }_2\alpha |f|e_1}$, up to numerical factors. The transit-timing variation follows that libration period and depends on its libration amplitude of a resonant argument. Exactly at a resonant equilibrium there need be no libration signal; near a separatrix the period grows, and eccentricity dynamics can modify the simple estimate. Resonant amplitudes remain controlled by the bounded resonant motion rather than diverging as the circulating detuning tends to zero.

For constant mean motions, consecutive mean conjunctions are separated by $P_{\mathrm{syn}}=P_b{P}_c/(P_c-P_b)$ and advance in longitude by $\Delta \Lambda =2\pi /(P_c/P_b-1)$. To follow one of the $q$ branches, compare conjunctions $q$ steps apart. At exact resonance that advance is $2\pi p$; subtract it to measure the slow drift. Hence conjunction-pattern precession near a mean-motion resonance isThe printed expression is its leading small-offset term:Positive $\delta$ gives retrograde drift and negative $\delta$ prograde drift. The tracked-branch rotation time is about $q{P}_b/(p|\delta |)$; the full set of $q$ interchangeable directions repeats after a rotation $2\pi /q$, giving the usual mean-pattern super-period $P_b(1+\delta )/(p|\delta |)$. The slow, continuously tracked branch interpretation assumes the resonance neighbourhood $|\delta |(p+q)\ll q$; merely $|\delta |\ll 1$ need not give a slow pattern for an arbitrarily large integer ratio $p/q$. At $\delta =0$ there is no offset-driven drift.

The allowed directions can consequently move through the observer's line of sight. Triple exoplanet transits require the $bc$ and $cd$ patterns, together with their timing phases, to overlap that line of sight within the finite exoplanet transit windows; differing drift rates can open and close such windows. A small offset neither guarantees a triple event nor prohibits it permanently.

Moreover, in a true mean-motion resonance, bounded ${\varphi }_{bc}$ can coexist with a period offset because $(p+q)n_c-p{n}_b=q{\dot{\varpi }}_c+{\dot{\varphi }}_{bc}$. The mean conjunction pattern then follows apsidal precession plus the resonant-argument libration correction. For appreciable orbital eccentricity, the true-conjunction geometry has the periodic corrections identified above; the drift calculation is a mean-phase result.

First adopt the uniform-angular-motion approximation implicit in the requested formula. Set the line of sight to longitude zero, the $b$ exoplanet transit to $t=0$, and the next $c$ exoplanet transit to $t=\tau$, with $0\le \tau <P_c$. Write $n_b=2\pi /P_b$, $n_c=2\pi /P_c$, with $P_c>P_b$. The next conjunction satisfiesThe common longitude is $n_{b}t$ modulo $2\pi$. Subtracting one $2\pi$ gives the conjunction longitude from a transit time lagSubstitution of the period ratio yieldsIt is exact in this uniform-angle model; no first-order expansion in $\delta$ is necessary here. For other choices of which exoplanet transit is used, choose the appropriate whole-turn branch.

For the finite-eccentricity orbit specified earlier, this is an approximation. Let $f_L$ and $M_L$ be the true anomaly and mean anomaly at the line of sight. The eccentric planet has $M_c(t)=M_L+n_c(t-\tau )$ and true angular advance $n_c(t-\tau )+{\Delta }_e(t)$, whereTrue conjunction therefore requires $(n_b-n_c)t+n_c\tau -{\Delta }_e(t)=2\pi$. The PDF omits this term, which is generally $O(e_c)$.

For a concrete counterexample take $p=q=1$, $\delta =0$, $e_c=1/5$, $\tau =P_b/4$, and a line of sight along $c$'s periapsis. The uniform model predicts $t=7P_b/4$ and longitude $3\pi /2$. At that time $M_c=3\pi /2$, but Kepler's equation gives $f_c\ne 3\pi /2$, so the planets are not truly aligned. Accurate finite-eccentricity conjunctions must be found from the corrected equation.

A simultaneous transit in a resonant chain requires both conjunction patterns to admit the observer's longitude at the same time. In the leading mean-conjunction geometry put $x=\Lambda -{\varpi }_c$. The necessary phase compatibility at symmetric centres isLet $d=\gcd (q,n)$. Eliminating $x$, and using that $q\ell -nk$ ranges over multiples of $d$, gives the exact modular compatibility criterionHere $\gcd$ denotes the greatest common divisor. For the orders $q,n\in \{1,2,3\}$ above, the two phase sets intersect only whenThe common directions are the two quadratures relative to periapsis. In reduced resonances this requires $p$ and $m$ odd; otherwise the apparent second-order ratios reduce to different resonance orders. There is no additional restriction on those odd values from angular compatibility alone. A shared direction and suitable temporal phase can be chosen, and the rational period ratios allow its recurrence.

This is a potential low-amplitude configuration, not a guarantee of an observed triple exoplanet transit. The observer must lie in the common orbital plane and near an allowed direction. Finite stellar radii, nonzero libration amplitudes of a resonant argument, apsidal motion and asymmetric centres broaden or change the possibilities. At finite $e_c$, true simultaneous alignment obeys the more precise equations $(p+q)y-px={\varphi }_{bc}$ and $(m+n)x-my={\varphi }_{cd}$, with $y=M_c(x)$. Thus the order-only conclusion is not a universal necessary condition for every eccentric resonant chain.

Define $s_i=\sin [2\pi (t_i-T)/P+\pi ]$ and use the independent quoted Gaussian errors $\sigma =10\mathrm{m}{\mathrm{s}}^{-1}$, with period and exoplanet transit epoch fixed. The model is linear in $K,V_0$ despite its sinusoidal time dependence. Minimize ${\chi }^2={\sum }_i(V_i-K{s}_i-V_0{)}^2/{\sigma }^2$. With $N=5$ the normal equations areUsing $\sum s_i=-0.741472$ and $\sum s_i^2=2.824633$, the solution isThe residuals in observation order are $(3.67,11.01,2.74,-7.90,-9.51)\mathrm{m}{\mathrm{s}}^{-1}$. Therefore ${\chi }_{\min }^2=2.950$ for $5-2=3$ degrees of freedom, giving reduced ${\chi }^2=0.983$. The upper-tail probability is about $0.40$: the fit is consistent with the stated errors, although five points do not establish a unique physical explanation.

For known errors the covariance matrix is ${\sigma }^2(X^{T}X{)}^{-1}$. Equivalently, profiling out the fitted offset givesThis is the $\Delta {\chi }^2=1$ one-parameter uncertainty, not the curvature with $V_0$ incorrectly held fixed. Propagating it through the supplied mass relation, and treating its coefficient and $P$ as fixed,The symbol $\pm nP$ in the exoplanet transit epoch denotes periodic copies of the same exoplanet transit phase, not an uncertainty distribution to propagate. No independent error for the empirical mass coefficient was supplied.

Additional independent observations retain positive weight even with twice the error: each has one quarter of the previous per-observation weight. Their information adds to the existing information matrix, so the formal mass uncertainty decreases for informative additional orbital phases and cannot increase under this fixed-error model. It can remain unchanged in a degenerate phase configuration that supplies no further amplitude information after fitting the offset. The conclusion assumes the same correct model and independent errors, not a new fitted error rescaling or unmodelled systematic effects.

The intact source figure gives orbital period $P=6.0\mathrm{yr}$, stellar radial velocity semiamplitude $K=30.0\mathrm{m}{\mathrm{s}}^{-1}$, and orbital eccentricity $e=0$. A planetary exoplanet transit implies a nearly edge-on orbit, so $\sin i\simeq 1$. Adopt a solar mass for the Solar-type star and initially neglect the companion's mass relative to it.

For a circular solution of the two-body problem, the star's distance from the centre of mass is $a_{\ast }=a{M}_p/(M_{\ast }+M_p)$, and the measured semiamplitude is $K=2\pi a_{\ast }\sin i/P$. Kepler's third law givesHenceThe inferred binary mass ratio is only $0.00183$, validating the approximation. Since $e=0$, the star–planet separation is constant and equals the orbital semimajor axis. Without the exoplanet transit, Doppler spectroscopy would generally give only a minimum companion mass through $M_p\sin i$.

Transit photometry. A planet crossing the stellar disc produces a periodic reduction in flux, approximately $\delta =(R_p/R_{\ast }{)}^2$ for a small, dark planet away from the limb. The orbital period and stellar radius allow an inferred planetary radius, but not its mass by this method alone. The geometric geometric transit probability is roughly $(R_{\ast }+R_p)/a$ for a circular orbit, so the method strongly favours close orbits and large planets. Small signals compete with instrumental noise and stellar activity; eclipsing stellar binaries and blended sources can mimic a planetary exoplanet transit, and long orbital periods require long observing baselines.

Doppler spectroscopy. Periodic shifts of stellar spectral lines reveal the stellar radial velocity caused by motion around the centre of mass. For a light companion, $K\propto M_p\sin i{P}^{-1/3}{M}_{\ast }^{-2/3}/\sqrt{1-e^2}$, favouring massive, short-period companions and low-mass host stars. Unknown orbital inclination gives a minimum mass; stellar activity, line broadening, and instrumental calibration limit precision. Long orbital periods require years of measurements.

Direct imaging. Spatially resolved reflected light or thermal emission separates a planet from its host. The large star–planet brightness contrast and small angular separation demand starlight suppression and precise correction of optical aberrations. The method favours wide orbits around nearby stars, and young, hot, massive planets with stronger thermal emission. Translating brightness to mass depends on planetary cooling models and age; it is difficult to image an Earth analogue close to a bright star.

A exoplanet transit occurs when a planet crosses the projected stellar disk and blocks some stellar flux. For a circular orbit with isotropically distributed angular-momentum directions, $\cos i$ is uniform on $[0,1]$. A transit requires $a|\cos i|<R_{\star }+R_p$, soFor a dark planet, a uniformly bright stellar disk and negligible dilution, the transit exoplanet transit depth is $\boxed{\delta \simeq (R_p/R_{\star }{)}^2}$. For $R_{\star }/a\ll 1$, a small planet and impact parameter $b=a\cos i/R_{\star }$,The central-transit result sets $b=0$. A first-to-fourth-contact duration retaining planet size replaces $\sqrt{1-b^2}$ by $\sqrt{(1+R_p/R_{\star }{)}^2-b^2}$ and includes the exact projected-speed geometry when needed. Eccentricity, limb darkening and grazing geometry modify these expressions.