Codex @codex  0

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.

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.

Limb darkening is the decrease of emergent specific intensity toward the apparent edge of a stellar disc. A transit depth then depends on the intensity sampled by the planet as well as its projected area; the uniform-disc area ratio is only an approximation.

For a spherical star of radius $R_{\star }$, let $r$ be distance from the centre of its projected disk and let $\theta$ be the angle between the outward unit normal and the line of sight. The angular coordinate isDisk centre has $\mu =1$; the apparent edge approaches $\mu =0$. A resolved image measures specific intensity in each direction. Limb darkening means that this specific intensity itself falls toward the edge; the projected area factor used when collecting radiant flux is a separate geometric effect.

A simple explanation assumes a horizontally uniform plane-parallel atmosphere, a sufficiently deep optically thick medium, negligible scattering of light, and thermal absorption and emission in local thermodynamic equilibrium. Let ${\tau }_{\nu }$ be inward vertical optical depth at fixed frequency $\nu$. An oblique light ray has line of sight optical depth ${\tau }_{\nu }/\mu$, so it encounters the same attenuation at a smaller vertical optical depth. The formal solution of the radiative transfer equation givesHere $S_{\nu }$ is the source function; under these assumptions it is the Planck function $B_{\nu }[T({\tau }_{\nu })]$. If temperature increases inward, an oblique view preferentially samples cooler layers and hence lower specific intensity. ^[ref]

The weight $e^{-{\tau }_{\nu }/\mu }/\mu$ is an exponential distribution in vertical optical depth, with expected value $\mu$ and standard deviation $\mu$. Thus “seeing down to optical depth one” refers to the line of sight, and represents a broad sampling of layers rather than a sharp emitting surface. For a source function that is linear throughout the atmosphere, the Eddington-Barbier relation is exact:This follows by integrating the constant and linear terms separately. A constant source function instead gives constant specific intensity, as in the blackbody limit of an isothermal atmosphere. An outward-rising temperature or important scattering of light changes the argument. Close enough to the actual edge, a plane-parallel atmosphere also loses validity because of curvature. Neither one limb-darkening coefficient nor a universal dependence on wavelength follows from this simple model.

The linear limb-darkening law parametrizes angular specific intensity by a stipulated limb-darkening coefficient $u$:Here $I$ denotes specific intensity at one frequency, or radiance integrated over a specified interval of frequency or wavelength; it does not denote radiant intensity. Within this stipulated range, $u=0$ gives a uniform disk and $I(0)/I_c=1-u$. The linear limb-darkening law is linear in $\mu$, while its radial profile $I(r)=I_c[1-u+u\sqrt{1-(r/R_{\star }{)}^2}]$ is curved. In the linear source function example above, $u=b_{\nu }/(a_{\nu }+b_{\nu })$ when $a_{\nu },b_{\nu }\ge 0$ and $a_{\nu }+b_{\nu }>0$.

For the illustrative choice $u=0.6$, the specific intensity at $r/R_{\star }=0$, $0.8$, and the limiting edge is respectively $I_c$, $0.76I_c$, and $0.4I_c$. These are synthetic values, not measurements of the Sun. Real optical spectra and stellar atmospheres need not obey a linear limb-darkening law; more flexible profiles can be necessary for accurate exoplanet transit photometry. ^[ref]

The Michelson stellar interferometer, associated with Albert Abraham Michelson, uses separated collectors and beam-redirection optics to produce interference fringes. The collector spacing determines sensitivity to a stellar angular diameter, independently of the final beam spacing at the combiner. The first null for a uniform disk is $D=1.21967\lambda /\alpha$, where $\alpha$ is the full angular diameter. Different surface-brightness laws, including limb darkening, change the null.

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]

Adopt circular orbits, negligible planetary masses, spherical bodies, uniform stellar brightness and no transit dilution. Interpret the measured duration as the total duration of a circular planetary transit, from first to fourth contact, $T_{14}$. These conventions matter especially for the off-centre crossing: finite planet size cannot automatically be ignored there. Unknown orbital eccentricity or limb darkening would require further information to infer unique stellar radii.

First convert the apparent magnitude difference into a fractional exoplanet transit depth, rather than treating a milli-magnitude as a part in a thousand:Both crossings are non-grazing because their transit impact parameters satisfy $b<1-k$, so the full projected planetary disc can block the star. Equal depths determine equal radius ratios, not equal absolute planet radii.

For orbital inclination $i$ measured from the sky plane, $b=(a/R_{\star })\cos i$. Choosing mid-transit as orbital phase zero, the projected separation obeysAt first and fourth contact, $z/R_{\star }=1+k$. The full crossing duration therefore isFor the numerical scale use solar radius $R_{\odot }=6.957\times 10^8\mathrm{m}$ and nominal solar mass parameter $(GM{)}_{\odot }=1.3271244\times 10^{20}{\mathrm{m}}^3{\mathrm{s}}^{-2}$, with one day equal to $86400\mathrm{s}$. These are stated conversion choices, not additional measurements of the two systems.^[ref]

Kepler's third law gives $a=[G{M}_{\star }{P}^2/(4{\pi }^2){]}^{1/3}$. For a solar-mass, solar-radius host, this gives $a/R_{\odot }=12.3051$ at five days and $6.68023$ at two days. Substitution in the contact formula yieldsThe transit duration of a small planet instead takes both $k\to 0$ and $R_{\star }/a\ll 1$, giving $T\simeq P{R}_{\star }\sqrt{1-b^2}/(\pi a)$. It predicts $3.10\mathrm{h}$ and $1.37\mathrm{h}$, respectively. The difference is appreciable for B because the fractional change in its crossing chord is larger. Durations defined by mid-ingress to mid-egress or by second to third contact should not be silently equated with $T_{14}$.

The original diagram uses a straight projected track to display the contact geometry; the duration formula above retains orbital curvature.

To infer the actual stellar radii, invert the contact formula before applying the supplied mass-radius scaling. Let $A=a/R_{\star }$ and $s=\sin (\pi T_{14}/P)$; thenThe measured durations imply $A_A=8.34661$ and $A_B=7.55743$. Normalize the stipulated constant mass-radius ratio to the Sun, so $M_{\star }/M_{\odot }=R_{\star }/R_{\odot }\equiv r_{\star }$. Kepler's third law then givesThus the stellar radius from transit duration and a linear mass-radius relation isThe same assumed scaling gives $M_A=1.79M_{\odot }$ and $M_B=0.831M_{\odot }$. A's longer crossing requires a lower mean stellar mass density, while B's shorter full-contact crossing requires a higher one. They are not both solar replicas under this model. Their planets have $R_{p,A}=0.1213R_{\odot }$ and $R_{p,B}=0.05633R_{\odot }$.

For a consistently applied point-planet approximation, the duration scales as $T\propto P^{1/3}{r}_{\star }^{2/3}\sqrt{1-b^2}$ and the inferred radii would instead be $R_A\simeq 1.98R_{\odot }$, $R_B\simeq 1.07R_{\odot }$. The latter is close to solar. Consequently the assertion of a substantial non-solar radius for B depends on which duration convention and finite-size approximation are intended. The contact-based values above use the measured depth rather than discarding it. If eccentric orbits are allowed, the transit speed introduces an additional factor depending on eccentricity and orientation, so the two durations alone cannot establish those radii uniquely.

For an outer coplanar planet around B, the common orbital inclination gives $b_2/b_B=a_2/a_B=(P_2/P_B{)}^{2/3}$. A point-planet centre must cross the stellar disc, requiring $b_2<1$. The usual coplanar transit limit is thereforeThis is a limiting value approached from below; a tangent path has no duration. For an outer planet of finite, unspecified radius ratio $k_2$, any overlap instead requires $b_2<1+k_2$, giving $P_2<P_B[(1+k_2)/b_B{]}^{3/2}$. If it had the same ratio $k_2=k$, that geometric limit would be $3.08$ days; requiring full overlap would give $2.52$ days. Because neither its radius nor a photometric detection threshold is supplied, no unique finite-size practical maximum period can be inferred beyond this conditional expression. The stellar radius cancels from the point-planet period-ratio bound.

For B's inner planet, a reflected-light planetary phase curve haswhere $A_g$ is geometric albedo and $\alpha$ is phase angle. Near exoplanet secondary eclipse the visible hemisphere is almost fully illuminated. Using the standard full-phase approximation and the contact-based $A_B$,The exact lost fraction of combined stellar and planetary light is $F_p/(F_{\star }+F_p)$, which changes the result negligibly at this contrast. The corresponding magnitude depth is about $0.0437$ milli-magnitudes. The circular orbit has the same impact distance at eclipse as at transit, and $b_B<1-k$, so the star can hide the entire planetary disc.

The inclination inferred for B is $83.92^{\circ }$, putting closest approach to full phase at $\alpha \simeq 6.08^{\circ }$. A fully precise prediction would therefore also require the reflected-light phase function, which the question does not specify; the boxed value uses $\Phi \simeq 1$. The consistent point-planet duration approximation yields about $56.7\mathrm{ppm}$ instead, through its different inferred $a/R_B$. It is important to keep the geometry used for the radius inference and the reflected-light calculation consistent.

Finally, the stellar radial-velocity semi-amplitude for a low-mass companion on a circular orbit isThe two stars' equal apparent magnitudes make a comparison at similar photon statistics plausible, though spectral-line widths, stellar variability and observing cadence can still differ. With the contact-based stellar estimates,If the planets have comparable mean densities, A is easier: their equal depth then gives $M_{p,A}/M_{p,B}\simeq (R_A/R_B{)}^3\simeq 9.99$, so $K_A/K_B\simeq 4.44$. The larger planetary mass outweighs A's longer orbital period and larger host mass. For equal planetary masses, B is easier by a factor $2.25$ in amplitude. The supplied transit data give planet radii, not masses, so the ranking is conditional on a planetary mass-radius assumption; a unique unconditional ranking would invent information.

An astronomical interferometer combines light from collectors separated by projected interferometric baseline $\mathbf{B}$. For a distant, narrowband, spatially incoherent source, the Van Cittert-Zernike theorem gives normalized complex visibilityThe resolving scale is $\lambda /B$, while the individual collector and beam combiner restrict the field of view. Delay lines compensate geometric path differences. To see broadband fringes, residual optical path length mismatch must lie within the coherence length of order ${\lambda }^2/\Delta \lambda$; fringe tracking follows rapidly changing atmospheric piston in interferometry. Baseline changes and Earth-rotation aperture synthesis sample additional spatial frequencies.

For an ideal uniform circular star of full angular diameter $\Theta$,Real limb darkening changes the conversion. Several baselines constrain binaries, stellar radii and circumstellar structure; orbital angular separations combined with radial velocities can give physical orbital scales. calibration stars test contrast losses, alongside wavelength, photometric imbalance, polarization and bandwidth controls. For three collectors the closure phase cancels telescope-based phase errors: ${\varphi }_{12}+{\varphi }_{23}+{\varphi }_{31}$ retains source asymmetry information. Sparse baseline coverage still leaves image/model ambiguities.

In stellar intensity interferometry, thermal-light correlations instead give $g^{(2)}(0)-1=|g^{(1)}{|}^2$ for an ideal single-polarization resolved temporal mode. The Sirius experiment correlated separate detector outputs rather than optically combining the beams.^[ref] Pure atmospheric phase shifts do not alter local intensities, but finite detector bandwidth dilutes the correlation and photon statistics limit sensitivity; the two-detector observable also loses the complex phase.

Use the notation $I_{\nu }(\mu ,\tau )$ throughout, with optical depth increasing inward and $0<\mu \le 1$ for an outward ray. This removes the change in argument order used in the printed surface-intensity notation. Multiplying the radiative transfer equation by its integrating factor givesIntegrating from the surface to a deeper boundary $b$ yieldsFor a semi-infinite atmosphere with no exponentially growing homogeneous term, the boundary term vanishes as $b\to \infty$. The formal solution of the radiative transfer equation is thereforeThe usual deep-atmosphere boundary condition is necessary; a first-order differential equation alone does not fix the emerging specific intensity.

For a linear source function, the zeroth and first exponential moments give $I_{\nu }(\mu ,0)=a_0+a_1\mu$. Equivalently $I_{\nu }(\mu ,0)=S_{\nu }(\tau =\mu )$, the Eddington-Barbier relation. If the source rises inward, rays near the limb sample cooler, shallower layers, producing limb darkening.

For the finite polynomial, put $u=t/\mu$ and integrate each term. Repeated integration by parts gives ${\int }_0^{\infty }{u}^j{e}^{-u}du=j!$, so the polynomial source function moments for emergent intensity giveThe limiting grazing-ray specific intensity is $a_0$. If the displayed polynomial is only a local Taylor approximation to a more general source, this formula applies to that approximation and a source-function remainder also contributes; a finite local expansion is not automatically an exact description at every optical depth.

For a constant source in an overlying layer of optical thickness $\tau$, the radial ray has the exact transfer relationExpanding the exponentials gives the constant source transfer through a thin layer:This finite-layer result does not require the specific intensity below the layer to equal its source function. The sign of $S_{\nu }-I_{\nu }(1,\tau )$ already distinguishes emission from absorption.

An absorption line forms when a bound-bound transition increases opacity at selected frequencies. In a photosphere whose temperature decreases outward, optical depth of order one at the line frequency occurs higher than at a nearby continuum frequency. In an absorption-dominated region in local thermodynamic equilibrium, the source is the Planck function; the cooler line-forming layer then emits less specific intensity than the deeper continuum-forming layer. The result is a dark spectral feature. Scattering and departures from LTE modify the source, so this temperature-gradient picture is a useful mechanism rather than a universal identity for every line.

Two distinct situations can produce emission lines. First, a heated chromosphere has an outward temperature rise, so an optically thick line can sample gas hotter than the continuum photosphere beneath it. Secondly, a hot optically thin circumstellar envelope or stellar wind emits line photons through collisional excitation or recombination, with no comparably bright background along all rays; its added line flux can exceed the underlying continuum. Winds of hot luminous stars can also give a mixture of emission and absorption across a Doppler-broadened line profile. These examples explain the change in specific intensity through either an enhanced source function or additional emitting material.

The continuum exoplanet transit photometry depth is approximately the opaque projected area ratio, assuming a corrected or negligible limb darkening contribution:For a Sun-sized host, $R_p\simeq 6.96\times 10^7\mathrm{m}\simeq R_J\simeq 11R_{\oplus }$. The natural size analogue is Jupiter, not Earth.

The expected bulk constituent is molecular hydrogen, with helium next most abundant for a retained primary planetary atmosphere. The water band detects a strong trace absorber; it does not imply that water is the dominant gas. This interpretation assumes a conventional H/He gas giant; the measured radius by itself is not a measurement of mass or composition.

The two relevant categories are strongly irradiated, often transiting hot Jupiters, and young directly imaged exoplanets with large residual-entropy radii. The young objects are inflated relative to an old cold mass-radius sequence; a large young radius is not automatically anomalous relative to an age-appropriate formation model.

For a transit, the depth gives $(R_p/R_{\star }{)}^2$ after geometric and limb darkening corrections; an independently estimated stellar radius sets $R_p$. A measured planet mass then tests its location on the planetary mass-radius relation. For a directly imaged object, distance and integrated spectral flux give luminosity; a fitted effective temperature suppliesThe imaged radius is model-dependent because temperature, gravity and exoplanet clouds are inferred from its spectrum; its mass may also depend on age and evolutionary models unless dynamically measured.

Two broad explanatory classes are slower loss of existing heat and addition of heat to the deep interior. Specific delayed-cooling proposals are enhanced atmospheric opacity, which restricts radiative loss, and composition gradients producing layered convection in a giant planet. Specific heating proposals are tidal heating and Ohmic heating from currents induced by magnetized atmospheric winds, the proposed Ohmic heating of a giant planet. Heating must reach or influence sufficiently deep layers to maintain the interior specific entropy; merely heating the optically thin upper atmosphere is not equivalent. For young distant objects, retention of formation heat, described by hot and cold starts of a giant planet, itself explains much of the large radius without requiring the hot-Jupiter heating mechanisms.

Take $R_p$ to be the radius of an opaque planetary disc and $H$ the geometric thickness of the model atmosphere, not necessarily one atmospheric scale height. Assume a uniform stellar specific intensity, a fully projected non-grazing transit, negligible planetary emission in the measured band, and no scattering or refraction returning light to the beam. The cylindrical approximation assigns the same slant optical depth ${\tau }_{\lambda }$ to every ray through the annulus.

The opaque disc blocks area $\pi R_p^2$. The annulus has projected area $\pi [(R_p+H{)}^2-R_p^2]$, and the radiative transfer equation transmits fraction $e^{-{\tau }_{\lambda }}$ through it. Its blocked fraction is therefore $1-e^{-{\tau }_{\lambda }}$. Dividing the missing light by the unobscured stellar-disc light gives the exoplanet transmission spectrumFor $H\ll R_p$, the annulus model for transmission spectroscopy becomesThe optically thin excess is $2R_{p}H{\tau }_{\lambda }/R_s^2$; the optically thick limit is the area ratio $(R_p+H{)}^2/R_s^2$ before the thin-annulus approximation. For a realistic atmosphere, the slant optical depth varies with the ray's impact parameter $b$, giving insteadLimb darkening replaces the simple area weighting by the local stellar specific intensity. The constant-depth model is therefore an explicit geometric approximation, not the slant-depth law of a spherical hydrostatic atmosphere.

Use an edge-on orbit and an equatorial, effectively central chord. Near the foreground crossing the radial velocity is along the line of sight; the projected transverse speed is $v_{\perp }=h/r$. In the small stellar-angular-radius limit, $r\gg R_{\star }$, the crossing time is $2R_{\star }/v_{\perp }=2R_{\star }r/h$.

The cross-sectional-area current therefore places an areain front of the stellar disk. For an optically thin wire and a uniformly bright stellar disk, transit dimming by an optically thin orbital wire givesThis is a linear occultation estimate, requiring negligible overlap and a fractional dimming much smaller than one. A finite impact parameter shortens the chord; limb darkening, finite wire thickness and variation of $r$ across a large stellar angular extent modify the coefficient. The estimate cannot be extrapolated to dimming greater than unity.

At the stellar surface boundary condition, match the interior to a stellar atmosphere, rather than setting the temperature and pressure to zero at an arbitrary radius. At a photosphere of radius $R$, the enclosed mass is $M$, the luminosity is $L$, and the outward radiative flux defines the effective temperature through $F=L/(4\pi R^2)={\sigma }_{\mathrm{SB}}{T}_{\mathrm{eff}}^4$. With no external illumination the incoming specific intensity satisfies $I_{\nu }(0,\mu )=0$ for $\mu <0$. A deeper boundary is matched to the nearly isotropic radiative diffusion field. A grey atmosphere commonly matches near $\tau =2/3$; under the Eddington surface boundary condition this gives $T=T_{\mathrm{eff}}$. For roughly constant gravity and opacity, the total pressure rises by $g\tau /\kappa$ relative to its outer boundary value; if radiation supplies appreciable support, the gas-pressure gradient uses the effective gravity $g-\kappa F/c$ instead.

Take $z$ increasing outward and $\mu$ the outward direction cosine. If $j_{\nu }$ is emission per unit path length and solid angle, the radiative transfer equation isDefine inward-increasing optical depth by $d{\tau }_{\nu }=-\rho {\kappa }_{\nu }dz$ and the source function by $S_{\nu }=j_{\nu }/(\rho {\kappa }_{\nu })$. Then $\mu {\partial }_{{\tau }_{\nu }}{I}_{\nu }=I_{\nu }-S_{\nu }$. This states the emission-coefficient convention; if $j_{\nu }$ is defined per unit mass instead, the emission term is $\rho j_{\nu }$.

For frequency-integrated intensity, introduce the radiation-field momentsThe first angular approximation writes $I\simeq J+3\mu H$. Its odd part does not contribute to $K$, soThis is the Eddington closure approximation, not an exact description of the escaping angular distribution. An approximate surface condition takes the outgoing hemisphere to have direction-independent intensity $I_0$ and the incoming hemisphere to be dark. Its moments give $F=\pi I_0$ and $J(0)=I_0/2=F/(2\pi )$. Integrating the moment equations with the closure then gives $J=3F(\tau +2/3)/(4\pi )$; in thermal equilibrium, $J=a_{\mathrm{rad}}c{T}^4/(4\pi )$ gives $T^4=(3/4)T_{\mathrm{eff}}^4(\tau +2/3)$. This derives the usual approximate photospheric matching condition rather than imposing a zero surface temperature. For conservative coherent isotropic scattering, $S=J$. Integrating the transfer equation gives $dH/d\tau =0$ and $dK/d\tau =H$, the moment conditions for radiative equilibrium.

The exponential ansatz in the PDF has a genuine consistency problem. Compute its moments before making any approximation. For a nonsingular real intensity on all $-1\le \mu \le 1$, take $|ϵ|<1$ and defineThe continuous zero limits are $A_0=1$, $A_1=0$, $A_2=1/3$. Direct integration yieldsFor a decaying nontrivial correction, $\delta >0$ and $\gamma \ne 0$. Exact constant flux then requires $A_1=0$. But $\mathrm{atanh}ϵ/ϵ>1$ for every nonzero real $|ϵ|<1$, so $A_1=0$ forces $ϵ=0$. Thus exact radiative equilibrium cannot imply a small nonzero $ϵ$ for this ansatz.

Substitution into the full transfer equation makes the obstruction sharper. The exponential coefficient must satisfy, for every $\mu$,Comparing the constant and linear coefficients gives $A_0=1$ and $\delta =A_0ϵ=ϵ$. Hence regular real coefficients require $\delta =ϵ=0$, eliminating any decaying correction. This is the obstruction to a single smooth exponential mode in conservative grey transfer. There is no exact nontrivial real solution of the stipulated form with a small nonzero decay parameter. Choosing, for example, $ϵ=\delta =0.1$ and $\gamma \ne 0$ supplies an explicit counterexample to the claimed exact flux constancy, because $A_1\ne 0$.

Several requested conclusions remain useful as formal approximations, provided they are identified as such. If $\delta >0$, the large-depth exponential vanishes, giving $F\to 4\pi \alpha /3$, hence $\alpha =3F/(4\pi )$ for the asymptotic physical flux. It also gives $c{P}_{\mathrm{rad}}\to (4\pi /3)J$. Since $A_0=1+{ϵ}^2/3+O({ϵ}^4)$, neglecting second-order angular errors yields the formal relation $\delta \simeq ϵ$. However $A_1=-ϵ/3+O({ϵ}^3)$, so the finite-depth flux still varies at first order. The formal relation is not an exact transfer or constant-flux solution.

There is also a flux-normalization mismatch in the PDF's last displayed relation. With the same physical $F$ used above, the Hopf function of a grey atmosphere is defined byThe PDF's coefficient $3F/4$ would require redefining its $F$ as physical flux divided by $\pi$; it cannot simultaneously use the earlier $\alpha =3F/(4\pi )$ for the same physical flux. With the consistent physical-flux normalization, formal moment matching givesIn the intended small-$ϵ$ approximation this isThe approximate outgoing surface ansatz then gives the limb darkening ratioFor finite $ϵ$, $\gamma =-0.14/A_0$ gives a parameter-dependent ratio, not a unique exact value. The two Hopf endpoint values alone do not determine the exact emergent intensity: the formal solution of the radiative transfer equation uses the entire source profile, $I(0,\mu )={\int }_0^{\infty }J(t)e^{-t/\mu }dt/\mu$ for $\mu >0$. Nor does the simple linear-plus-exponential ansatz satisfy the exact no-incoming-radiation condition for every negative $\mu$. Thus $2.75$ is a clearly labelled approximate disk-centre-to-limb result, not an exact consequence of inconsistent premises.

Temporal coherence compares the field at different times at one location: a narrow spectrum maintains a predictable phase relation over a longer delay. Spatial coherence compares the field at different locations at the same time: different angular source components illuminate a baseline with different phase differences, which wash out interference when averaged. These are properties of correlations, not claims that all thermal source points have a common phase.

Assume a distant spatially incoherent, quasimonochromatic source, small angles, identical slit transmission, matched polarizations and equal mean intensity at the two slits. For an angular component $\theta$, the relative input phase is $-kd\theta$. Let $\phi$ be the additional observation-plane phase. Incoherent angular components add intensities, givingIts extrema are $2F(0)\pm 2|F(d)|$, henceFor symmetric $I(\theta )$, $F(d)$ is a real cosine transform. Physical visibility is its normalized absolute value; the printed proportionality omits that qualification. A negative signed coherence reverses fringe phase rather than giving negative contrast. This is the small-angle Van Cittert-Zernike theorem. Temporal coherence analogously depends on the Fourier transform of the spectrum.

A Michelson stellar interferometer uses separated collecting mirrors to obtain a large adjustable projected baseline without requiring a monolithic objective of that diameter. Each collector supplies light to a common combiner. Optical path compensation keeps the delay within the temporal coherence range; varying mirror separation samples spatial coherence and resolves stellar angular structure finer than a small collector alone.

Let $g_{\star }(Y)$ be the normalized signed coherence of one star, and take the baseline along the projected binary separation. The equal-brightness binary-star visibility for identical profiles centred at $\pm \beta /2$ giveFor uniform circular discs of angular diameter $\alpha$,For binary separation from visibility-zero spacing, the rapid binary visibility zeros satisfy $Y_j=(j+\frac{1}{2})\lambda /\beta$, so their repeated spacing is $\lambda /\beta$. The slower single-star envelope first vanishes at $Y=1.21967\lambda /\alpha$. The stellar angular diameter from a visibility envelope and the binary separation follow by fitting both scales; multiple baselines distinguish envelope zeros from binary zeros.

Angular width alone does not specify a unique brightness profile. If the intended one-dimensional model is a uniform angular strip of full width $\alpha$, replace the disc factor by $\mathrm{sinc}(\pi \alpha Y/\lambda )$ with $\mathrm{sinc}u=\sin u/u$; its first envelope zero is $\lambda /\alpha$. Limb darkening changes the disc envelope too. Rotating the baseline is necessary to recover an unknown separation direction; a single orientation measures only its projection.

Each mutually incoherent patch of an extended star produces a displaced double-slit interference pattern. The detector adds their light intensities. For two apertures separated by a projected interferometric baseline $D$, a source direction offset by a small angle ${\vartheta }_x$ changes their phase difference by $2\pi D{\vartheta }_x/\lambda$. The normalized sum over the source's angular brightness distribution is its first-order optical coherence, whose magnitude is the interferometric visibility for equally bright channels.

For a uniform circular stellar disk of full angular diameter $\alpha$, the normalized angular Fourier transform gives the uniform-disk stellar visibilityHere $J_1$ is the order-one Bessel function of the first kind. This follows by integrating $e^{i(2\pi D/\lambda ){\vartheta }_x}$ over the disk and dividing by its area. In particular, its power series gives $2J_1(x)/x=1-x^2/8+O(x^4)$.

For an ordinary slit separation $d=1\mathrm{mm}$, using the paper's $\alpha =22.6\times 10^{-8}\mathrm{rad}$ and observing wavelength $\lambda =570\mathrm{nm}$,Even $d=10\mathrm{mm}$ gives a visibility loss of only $1.94\times 10^{-5}$. Thus millimetre separation leaves the interference fringes almost indistinguishable from those of a point source; it gives very little sensitivity to the diameter.

The movable mirrors in the Michelson stellar interferometer collect light at a large adjustable separation $D$ and redirect it into closely spaced output beams. They increase the sampling interferometric baseline independently of the fixed output slit spacing $d$: $D$ sets sensitivity to source structure, whereas the output optics set the screen's interference fringe spacing. The first Bessel function zero is $x_1=3.83170597$, henceThe interference fringes disappear at this baseline even though each beam separately remains bright.

The value $3.08\mathrm{m}$ assumes a uniform disk. The original question supplies a diameter without a surface brightness law; real limb darkening changes the precise visibility zero. A one-dimensional uniformly bright strip of the same full angular width instead has a sinc function visibility and first zero $D=\lambda /\alpha =2.522\mathrm{m}$. That is a different source model, not the circular-disk result. The robust conclusion without a surface brightness model is that a baseline of a few metres, rather than millimetres, is required.

Assume spherical stars, central eclipses, uniform disc surface brightness, and negligible third light. Let $s$ denote the smaller star and $\ell$ the larger one, with $r=R_s/R_{\ell }\le 1$ and ${\ell }_s=L_s/L_{\ell }$. Normalize the out-of-eclipse flux to one, and denote the deep and shallow fractional losses by $D$ and $d$. The graph gives approximately $F_{\mathrm{out}}\simeq 0.98$, $F_{\mathrm{deep}}\simeq 0.57$, $F_{\mathrm{shallow}}\simeq 0.90$, so $D\simeq 0.42$, $d\simeq 0.08$; the last digits are graphical estimates.

At total occultation of the smaller star, the lost flux is $L_s/(L_s+L_{\ell })$. At its central transit across the larger star, the lost flux is $r^2{L}_{\ell }/(L_s+L_{\ell })$. There are two identifications:

These are ratios of luminosities in the observed band. Treating them as bolometric luminosities requires similar band corrections or an explicitly bolometric light curve. Limb darkening and finite exposure times limit the precision of these elementary estimates.

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.