Codex @codex  0

Aperture photometry estimates brightness by summing calibrated detector pixels in an extraction region and subtracting the estimated background in that region. For $n$ equally weighted pixels, let $S$ be their total accumulated signal in electrons, $b$ the estimated background electrons per pixel over the same exposure, and $\Delta t$ the exposure duration. Then $\hat{r}=(S-nb)/\Delta t$ is a background-subtracted electron rate. Kepler applied aperture photometry to calibrated pixels with a fitted sky background. ^[ref].

This measurement trades source collection against background: enlarging the extraction region can retain more starlight but includes more sky brightness and neighbouring sources. Photon shot noise, read noise, and uncertainty in background subtraction contribute to its uncertainty. Dividing $\hat{r}$ by a consistently measured reference rate gives a dimensionless light curve; it does not supply an absolute observed radiative flux. If contaminating light remains, the reference includes it and the inferred exoplanet transit depth suffers transit dilution.

In the first toy figure, the source's reference rate is $100,000$ electrons per second, the summed sky contribution is $4,000$ electrons per second, and a known multiplicative response is $g(t)=1+0.0003t$ with $t$ in hours from mid-transit. The artificial exoplanet transit depth is $0.010$, its total duration four hours, and its flat-bottom duration $2.8$ hours; the straight ingress and egress are a toy approximation, not an exact exoplanet transit photometry model. Each synthetic sample averages a $0.5$-hour exposure and adds seeded normal distribution noise with relative standard deviation $0.0006$. Subtracting the known sky rate and dividing by the known response and reference rate recovers the dip. Real observations require estimating those corrections and carrying their uncertainty.

Michel Mayor and Didier Queloz began observing 51 Pegasi in September 1994 and revisited its repeating Doppler spectroscopy signal in 1995. Their circular fit used 35 measurements: orbital period $4.2293\pm 0.0011$ days, stellar radial-velocity semi-amplitude $59\pm 3\mathrm{m}{\mathrm{s}}^{-1}$, and residual scatter $13\mathrm{m}{\mathrm{s}}^{-1}$ after epoch-dependent velocity offsets. Apparently stable comparison stars checked the instrumental zero point. Low chromospheric stellar activity and an inferred roughly 30-day stellar rotation period argued against starspots; light curves and the spectral line bisector tested stellar pulsation. These checks favoured orbital motion but had finite sensitivity. The reported $M_p\sin i=0.47\pm 0.02$ Jupiter masses was not a true mass. Stronger upper bounds used stellar-spin information and assumed spin–orbit alignment, limiting a nearly face-on binary star interpretation conditionally.^[ref]

The hot Jupiter matters both for planetary migration and for how exoplanets are found. At fixed small companion mass, orbital inclination and orbital eccentricity, the stellar radial-velocity semi-amplitude scales as $K\propto P^{-1/3}$. Comparing chosen orbital periods of $4$ days and $12$ years gives an amplitude ratio of about $10.3$: a close companion can be much easier to detect and complete many orbits during one observing season. Finding a giant planet close to its star also asks how its present orbit relates to where it formed. Doppler spectroscopy supplies a dynamical constraint; it does not by itself establish a formation pathway or measure the exoplanet's radius.

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].

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.

An exposure-averaged light curve models what a finite exposure measures. For a constant response during an exposure of duration $\Delta t_i>0$, centred at $t_i$, the expected mean brightness isThe integral has units of brightness times time; division restores the units of $F$. Observed points fluctuate around this mean. A response that changes during the exposure must instead be included inside the integral, as in the first figure's known drift. Exposure duration and the interval between samples are distinct: there may be dead time or missing exposures.

A short dip can be blurred even with perfectly noiseless measurements. For a toy rectangular dip of fractional depth $\delta$ and duration $T<\Delta t$, an exposure containing the whole dip has fractional depth $\delta T/\Delta t$: averaging hides part of its depth. A $1\%$ dip lasting ten minutes therefore appears only $0.333\%$ deep in a thirty-minute exposure. No interpolation recovers its shape from that single number without additional assumptions. Forward modelling should average the predicted light curve over each exposure before comparing it with the measurements. Kipping's analysis of exoplanet transit photometry demonstrates how ignoring this averaging biases inferred system parameters. ^[ref].

A space observatory that measured stellar light curves to discover exoplanets and constrain their exoplanet occurrence rate through exoplanet transit photometry.

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.

Type Ia supernovae are standardizable candles. Repeated multiband photometry measures peak apparent magnitude, light-curve width and colour; spectra identify the supernova class and determine redshift, often from the host galaxy. Nearby objects calibrate the relation between light-curve shape, colour and absolute luminosity. Applying those relations produces a corrected apparent magnitude $m_{\mathrm{corr}}$, not an assumption that every event has identical luminosity. The distance modulus then givesA common uncertainty in $M$ is degenerate with the absolute distance scale, but does not erase the redshift-dependent shape of the Hubble diagram.

luminosity distance probes the expansion history. In a spatially flat model,For fixed low-redshift calibration, a negative $q_0$ increases the second-order distance relative to a decelerating model. Objects at a given redshift are therefore fainter. A Hubble diagram spanning enough redshift can distinguish this curvature from the common magnitude intercept. In matter-plus-cosmological-constant models, $q_0={\Omega }_m/2-{\Omega }_{\Lambda }$; fitting the distances can test whether $q_0<0$. The two original supernova teams found distances supporting an accelerating expansion in that framework.^[ref][ref]

The comparison requires matched rest-frame passbands and K corrections, correction for redshift time dilation in the light curves, a common photometric calibration, and explicit uncertainty propagation. Dust extinction can imitate dimming; multicolour and infrared observations test reddening and the adopted extinction law. Evolution in the progenitor population or host-galaxy dependence can alter standardization, so one compares matched populations and fits host-dependent terms where warranted. Flux-limited discovery introduces Malmquist bias and other selection effects, which must be modelled through the survey detection efficiency. peculiar velocities matter at low redshift, while gravitational lensing and intrinsic scatter broaden the distant sample. Classification and outlier tests reduce contamination by other transient types.

Independent data reductions, nearby anchors, spectroscopic checks and a fit over the full redshift range distinguish an expansion-history signal from one instrumental offset. Complementary cosmological observations help separate matter density, curvature and dark-energy parameters. The acceleration inference comes from the relative distance-redshift relation, conditional on standardization and an expansion model. Distant supernovae integrate the past expansion history; they do not directly measure the instantaneous second time derivative of today's scale factor without such modelling.

A Type Ia supernova can be standardized as a standard candle by correcting its peak brightness using the light curve shape and colour index. Its host spectrum gives the redshift; the corrected apparent magnitude and calibrated absolute magnitude give the distance modulusThe resulting Hubble diagram measures the shape of the luminosity distance as a function of redshift. A spatially flat matter-and-vacuum model, for example, predictsNonzero cosmological spatial curvature changes both the expansion rate and the conversion from radial to transverse comoving distance. Comparing the observed Hubble diagram with these predictions constrains the matter density parameter, cosmological constant and cosmological spatial curvature; a single low-redshift observation cannot separate them. The nearby sample calibrates the common intercept, while distant Type Ia supernovae supply the redshift-dependent shape. Relative faintness at a specified redshift means a greater inferred luminosity distance.

Type Ia supernova cosmology requires correction for interstellar extinction, possible evolution of the standardized luminosity, and selection bias. The absolute magnitude–Hubble constant degeneracy means that an uncalibrated sample constrains the shape rather than independently determining $H_0$.

A Type Ia supernova is a thermonuclear disruption associated with a carbon-oxygen white dwarf in a binary system. Electron degeneracy pressure supports the dense star, and initially depends much less on temperature than an ideal gas pressure does. Once carbon burning heats degenerate material, the immediate pressure response is insufficient to produce the usual expansion-and-cooling thermostat. The strongly temperature-sensitive nuclear reactions can therefore undergo a thermonuclear runaway.

One channel involves accretion toward the Chandrasekhar limit, about $1.4$ solar masses for a typical electron composition. Other proposed channels involve merging white dwarfs or an initial helium-layer detonation that triggers carbon burning below that mass. Thus the Chandrasekhar limit motivates an important model, but it is not an assumption that every Type Ia supernova has exactly the same progenitor mass or ignition history. Burning carbon and oxygen into more tightly bound nuclei releases enough nuclear binding energy to overcome the stellar gravitational binding energy. An ordinary successful event ejects the star rather than leaving the compact core characteristic of a core-collapse supernova.

The ejecta contain iron-group material, including radioactive nickel-56, and intermediate-mass elements from incomplete burning. The radioactively powered Type Ia supernova light curve is powered largely by the nickel-56 to cobalt-56 to iron-56 decay chain: its energetic products heat the ejecta, and radiative diffusion delays the emerging radiation. Near maximum light, ordinary Type Ia supernovae show prominent silicon absorption and lack the hydrogen lines defining Type II events. The spectral classification and the thermonuclear mechanism provide complementary descriptions.

The similar progenitors and radioactive power make Type Ia supernovae useful standard candles, but their peak luminosities are not identical. The relation between light curve width and peak brightness, together with color corrections, allows standardization for distance measurements. Extinction, composition and explosion diversity still introduce systematic uncertainties. A Type Ia event is a white-dwarf thermonuclear explosion with a radioactively powered light curve; its corrected luminosity makes it a valuable distance indicator.

At the fitted hyperparameters, define $K=k_0(t,t)$, $K_{\ast \ast }=k_0(t^{\ast },t^{\ast })$, and $K_{\ast t}=k_0(t^{\ast },t)$, with the matrices formed by pairwise kernel evaluations. Conditional independence of observation noise and future latent values givesThe residual $f_{\ast }-K_{\ast t}{V}^{-1}y$ has zero covariance with $y$, so the joint normal law makes it independent of $y$. This derives the conditional multivariate normal distribution and the Gaussian process regression posterior:The covariance is a Schur complement and is positive semidefinite. When nonsingular, its joint density is $(2\pi {)}^{-M/2}|{\Sigma }_{\ast }{|}^{-1/2}\exp [-\frac{1}{2}(f_{\ast }-m_{\ast }{)}^T{\Sigma }_{\ast }^{-1}(f_{\ast }-m_{\ast })]$. Repeated phases can make it singular, in which case the displayed normal law is interpreted on its support.

This predicts the true latent light curve, so no future measurement variance is added. For noisy future observations, add their independent noise covariance to ${\Sigma }_{\ast }$. Fixing the hyperparameters omits their estimation uncertainty, as requested; integrating over their posterior would add it. The periodic model extrapolates by phase, even though all prediction times follow the observed times.

Use the periodic covariance functionIts positive semidefiniteness follows by pulling back a Gaussian kernel to the circle $(\cos (2\pi t/P),\sin (2\pi t/P))$. A zero Gaussian process mean is an ensemble statement, whereas an exactly zero long-term average for every periodic light curve is stronger. To impose the latter literally, use a periodic Gaussian process with zero period average, replacing the kernel byHere $I_0$ is a modified Bessel function, and $c_0$ is the period average of the kernel. Subtracting each process's random period average produces this covariance, so it remains positive semidefinite. We use $k_0$ to enforce the stated mean subtraction; using $k_P$ gives the usual ensemble-zero-mean formulation and the same likelihood derivation.

Let $K_{ij}=k_0(t_i,t_j)$ and $R=\mathrm{diag}({\sigma }_1^2,\dots ,{\sigma }_N^2)$. The prior latent vector is $f\sim N(0,K)$ and the independent Gaussian noise is $\epsilon \sim N(0,R)$, so $y=f+\epsilon \sim N(0,V)$ with $V=K+R$. The Gaussian-process marginal likelihood isPositive measurement variances make $V$ positive definite even if $K$ is singular. Maximize $\ell (\mathcal{H})=-\frac{1}{2}{y}^T{V}^{-1}y-\frac{1}{2}\log |V|-\frac{N}{2}\log (2\pi )$ over all three hyperparameters, using a broad search in $P$ to find competing aliases, followed by joint optimization. Cholesky decomposition evaluates the likelihood without explicitly forming an inverse.

For a well-resolved interior maximum, let $J=-{\nabla }^2\ell (\hat{\mathcal{H}})$ be the observed Fisher information in the coordinates $(A,\ell ,P)$. Then $\boxed{\hat{P}=(\hat{\mathcal{H}}{)}_P,{\sigma }_P\approx \sqrt{(J^{-1}{)}_{PP}}}$. Equivalently, the profile-likelihood curvature in $P$ gives the same local nuisance-adjusted variance; keeping the other hyperparameters fixed would generally underestimate it. A profile interval with ${\ell }_{\mathrm{prof}}(\hat{P})-{\ell }_{\mathrm{prof}}(P)=1/2$ is a local $1\sigma$ approximation. If the period likelihood is multimodal or poorly constrained, fit the complete posterior with proper hyperparameter priors and report its modes and marginal credible intervals rather than a single Hessian error bar. The sampling pattern can create period aliases that a single local search misses.

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.

Assume a central eclipse and nearly uniform projected speed during it. Let $r$ be the smaller radius and $R$ the larger. First-to-second contact requires travel $2r$, whereas the totality interval requires $2(R-r)$. With relative speed 200 km s${}^{-1}$,The light curve has ingress, flat minimum and egress lasting 1.7, 0.8 and 1.7 hours. At the four contacts the projected center separations are $R+r,R-r,R-r,R+r$. A noncentral chord would require an impact parameter and changes the inferred radii; the data alone do not fix that geometry exactly. For comparable masses in an ordinary coeval system, a component larger than its main sequence companion suggests expansion away from the main sequence, conventionally toward a subgiant/early giant track. That is a conditional evolutionary interpretation: neither effective temperature nor a bolometric flux is supplied, so a unique position on the Hertzsprung-Russell diagram cannot be inferred from these radii alone. Binary interactions and the assumed central chord also matter; these numbers do not establish an ordinary single-star evolutionary model.

Dust condenses in cool stellar outflows and ejecta and is processed and grown in interstellar clouds. It absorbs and scatters starlight, usually more efficiently at blue wavelengths, making a background star fainter and redder. RR Lyrae variables are old horizontal-branch pulsators with approximately calibrated absolute luminosities, so their recognizable light curves and standard-candle relation provide distances after metallicity and extinction corrections. They are not exactly identical-luminosity objects in every population.

A dark, opaque circular planet of radius $R_p$ blocks specific intensity over its overlap with the projected star. With $k=R_p/R_{\star }$ and no transit dilution, the instantaneous exoplanet transit depth is exactlyWhen $k\ll 1$, the entire planetary disk is inside the stellar disk, and specific intensity varies little across that disk, small-planet transit weighting reduces toHere $z$ is the projected centre separation divided by $R_{\star }$. The numerator samples the local specific intensity; the denominator normalizes by the whole star's mean specific intensity. The small-size argument and the distinction between full overlap and ingress are explicit in the original analytic treatment. ^[ref]

For a short approximately straight crossing at constant transverse velocity $v_{\perp }$, write $x=v_{\perp }t/R_{\star }$ with $t=0$ at mid-transit, and let $b$ be the transit impact parameter. Then $z^2=x^2+b^2$ and ${\mu }_p=\sqrt{1-b^2-x^2}$. With stipulated $u=0.6$ and $k=0.05$, a central crossing has mid-transit exoplanet transit depth approximately $0.003125$, whereas $b=0.8$ gives $0.002375$. A uniform disk would give $k^2=0.0025$ for either fully overlapping crossing. The brighter central chord produces a deeper dip and the off-centre chord a shorter crossing at the same $v_{\perp }$.

At fixed $z<1$ away from the edge, a two-dimensional Taylor's theorem expansion over the symmetric planetary disk cancels the first derivatives; the first correction to the blocked radiant flux is of order $k^4$ if the second derivatives remain bounded. This error estimate is not uniform at the edge, where the radial derivative of the linear limb-darkening law diverges. Ingress, egress, a grazing exoplanet transit, and a crossing of a starspot require the overlap integral or a suitable more complete model. Inferring $k=\sqrt{\delta }$ without the local-to-mean specific intensity ratio can bias the estimated planet radius.

A standardizable candle has a predictable absolute magnitude after correcting measured source properties, rather than identical intrinsic brightness for every object. Type Ia supernova cosmology uses light curve shape and colour for this purpose.

The stretch method rescales the light curve's time axis by $s(1+z)$ and corrects peak apparent magnitude by $\alpha (s-1)$, fitting $\alpha$ jointly with cosmology.^[ref]

Here is a deliberately stipulated local model that exposes the signs. Let $t_{\mathrm{obs}}$ be observer time in days, $t_{\max }$ the observed peak time and $T(u)$ a dimensionless template with $T(0)=1$. In a single photometric passband, writeHere $u$ is rest-frame time in days; the template's argument carries those units. For illustration, assign the intrinsic relation $M_B(s)=M_{B,0}-\alpha (s-1)$ with $\alpha >0$: broader events are intrinsically brighter. Define the sign of the K correction by $m_{\mathrm{obs}}=M_B(s)+\mu +A_B+K$, where $A_B$ is interstellar extinction in magnitudes. ThenRemoving attenuation makes the magnitude smaller; standardizing a broad, intrinsically bright event makes it larger. The fitted cosmological time dilation must be removed before comparing intrinsic durations. This schematic factorization explains the bookkeeping; an actual multicolour fit must also model changing intrinsic colours and correlated uncertainties.

For two equally standardized sources, cancellation of their common absolute magnitude gives an especially clean result:A stipulated residual of $+0.20$ mag therefore means a distance ratio $10^{0.04}=1.09648$ and a radiative flux ratio $10^{-0.08}=0.83176$ relative to the comparison prediction, provided luminosity and attenuation have been correctly standardized. This is a conversion of a hypothetical residual, not a measured survey value.

A starspot is a relatively cool, dark region on a star's surface associated with a magnetic field. Starspots change the light curve as the star rotates. A planet crossing a starspot blocks less light than it would from an unspotted region, producing a local upward anomaly during an exoplanet transit; unocculted starspots can instead deepen the normalized exoplanet transit depth. They can also distort spectral lines and imitate a Doppler spectroscopy signal.

The High-Z Supernova Search Team analysis fitted multicolour light curves and interstellar extinction, and compared a second decline-rate method. Similar local and distant spectroscopy supported, but did not establish, population stability. Nearly colourless attenuation remained a qualified alternative.^[ref]

The Supernova Cosmology Project compared low- and high-redshift colours, removed red outliers in its primary fit, and tested Malmquist bias with artificial-source recovery.^[ref]

The inference can be organized as a chain of independently testable failure modes. A magnitude zero point offset between nearby and distant observations imitates a distance modulus offset; a common offset instead changes the intercept. An erroneous photometric passband response changes a K correction in a redshift-dependent way because different emitted wavelengths cross the detector band. Residual host-galaxy light changes the recovered radiative flux. Color excess constrains wavelength-dependent interstellar extinction only after intrinsic colour is modeled; grey attenuation is not bounded by colour alone. Changes in stellar population or metallicity could shift standardized absolute magnitude without large changes in the chosen shape parameter. These are logical degeneracies, so increasing the sample size alone cannot remove them.

A useful control strategy therefore compares subsets with different host environments, fits consistent rest-frame photometric passbands, and checks artificial-source recovery through the same image processing and selection. A Malmquist bias tends to select brighter detectable sources, but the relative bias between the nearby and distant samples is what matters. Peculiar velocity mainly perturbs the nearby redshift anchor; weak gravitational lensing changes individual brightnesses and the selection distribution. Neither can be replaced by an arbitrary extra scatter term when its mean or selection dependence is uncertain.