Tie-line phase-fraction sensitivity measures how an inferred phase fraction changes when the overall composition or the measured liquid-liquid tie line endpoints change. A sensitivity coefficient is a local partial derivative of the inference formula, rather than an uncertainty by itself. It becomes an uncertainty contribution only after combination with the input covariance matrix.^[ref]

For one composition coordinate of a binary mixture, write the lever rule aswhere $a,b$ are endpoint mole fractions, $z$ is the overall mole fraction and $\lambda$ is the mole-based amount of phase $\beta$. Differentiating the mass balance givesThus a measurement error in the overall composition and an equal measurement error in the phase-$\alpha$ endpoint have different effects. The conditioning of tie-line phase-fraction inference deteriorates as the endpoint separation becomes small: the same composition measurement error then produces a larger change in the inferred phase fraction.

For an original example, take $a=0.20$, $b=0.80$, $z=0.44$, giving $\lambda =0.40$. Increasing only $a$ by $0.003$ gives a first-order change $\delta \lambda =-0.003$. Direct substitution gives $\lambda =0.237/0.597\approx 0.396985$, a change of $-0.003015$. For the closer endpoints $a=0.20$, $b=0.26$ and balanced composition $z=0.224$, the local sensitivity to the same endpoint error is ten times larger. This is a property of the inference geometry; it does not establish whether either hypothetical pair is in liquid-liquid equilibrium.

The multidimensional tie-line phase-fraction fit has an additional detail away from an exactly balanced composition. On independent mole fraction coordinates, hold a symmetric positive-definite matrix $W$ fixed and definewith distinct endpoints. Define row vectors $\mathbf{g}={\mathbf{d}}^{T}W/q$ and $\mathbf{h}={\mathbf{r}}^{T}W/q$. Differentiating the generalized least squares normal equation ${\mathbf{d}}^{T}W\mathbf{r}=0$ givesThe residual correction in tie-line sensitivity vanishes when $\mathbf{r}=0$. Otherwise, moving an endpoint also rotates the fitted line relative to the nonzero ternary tie-line consistency residual. Omitting that contribution is an additional approximation. If $W$ itself changes with the measured inputs, its derivative must also enter the calculation.

These formulas describe the unconstrained fit locally. Imposing $0\le \hat{\lambda }\le 1$ introduces boundary effects, and coincident endpoints make this inference undefined. A practical sensitivity check recomputes the fit after small positive and negative perturbations of one measured coordinate, keeping composition closure and the stated weight matrix. Agreement as the perturbation shrinks tests the derivative calculation; it does not validate the underlying equilibrium model.

A shot-resolved pulse-contrast estimator normalizes each individual optical shot before combining its contrast with other shots. For a specified side-delay signal $S_j$ and positive central signal $C_j$, both background-corrected, defineThe estimand is the mean of the individual shot ratios under the stated sampling and acceptance rule. It differs from normalizing an accumulated trace, ${\hat{R}}_{\mathrm{sum}}={\sum }_j{S}_j/{\sum }_j{C}_j$, which has third-order correlation brightness weighting.

An original three-shot example uses central signals $(1,2,8)$ and side signals $(0.01,0.04,0.80)$. The individual ratios are $(0.01,0.02,0.10)$, with equal-shot mean about 0.0433. Normalizing the sums gives $0.85/11\approx 0.0773$. The latter gives the brightest shot most of the influence. Both are valid summaries of the chosen data, but answer different questions.

Single-shot acquisition matters when the profile fluctuates. The primary third-order correlator maps delay into a spatially resolved signal and records it on a CCD, using second-harmonic generation and sum-frequency generation to form the gate.^[ref] A scanned trace assembled from different shots does not by itself provide all the $r_j$ for a common ensemble.

Normalization introduces pulse-contrast ratio denominator noise. First-order uncertainty propagation for one ratio $r=S/C$ givesShared baseline estimates can correlate the numerator and denominator, so their covariance cannot automatically be set to zero.^[ref]

For stipulated $S=0.1$, $C=10$, independent $u(S)=0.02$ and $u(C)=1$, the local approximation gives $r=0.01$ and $u(r)\approx 0.00224$. Near a zero or poorly resolved denominator the linear approximation and even a finite mean ratio may fail. A central-signal acceptance threshold changes the selected ensemble and must be reported.

For independent shots with a well-behaved ratio distribution, sampling uncertainty in the equal-shot mean can be estimated from the observed ratio variation. Consecutive shots may instead share drift, calibration errors or correlated pulse fluctuations. Keep the acquisition timing, background treatment, saturation checks and rejected-shot count with the estimate; averaging cannot repair an undefined denominator or an unspecified selection rule.

Aspidoscelis is a taxonomic genus of whiptail lizards in Teiidae. It contains both sexually reproducing species and all-female lineages that reproduce by parthenogenesis. This variation makes a closely related set of animals useful for studying how reproduction and chromosome organization interact.^[ref]

In a 2011 experiment, researchers crossed a triploid parthenogenetic female lineage with males of a sexually reproducing lineage. The resulting four-copy hybrid offspring were checked with molecular and chromosome measurements. Some female offspring themselves reproduced by parthenogenesis, establishing laboratory lineages able to reproduce without repeating the original cross.^[ref]

A chromosome-set calculation clarifies the change: an unreduced egg carrying three sets can combine with a sperm carrying one set to make a four-set offspring, $3n+n=4n$. This is tetraploidy, rather than merely four copies of one selected gene. Producing one viable hybrid is also different from demonstrating that its descendants can reproduce. The experiment tested both steps, linking the immediate effects of a cross to the persistence of a new reproductive lineage.

Rare-thickness leverage on inverse moments is the disproportionate contribution of uncommon, very thin layers to inferred density when $\rho =K/\delta$. Its sampling consequence is especially revealing: a sample can miss the influential component and report very little apparent variability.

Consider an original two-point thickness model. Let the rare thickness be $\epsilon$ with probability $p$, and the common thickness be $d_c$ otherwise. Define $a=1/\epsilon$ and $b=1/d_c$. Among $n$ independent measurements the rare count $N$ has a binomial distribution, and the estimated mean reciprocal thickness isConsequently,Multiplication by $K$ converts these into the mean and variance of the sample mean inferred density. The formulas are exact for this constructed model.

Take $p=0.0001$, $\epsilon =0.001$ micrometres and $d_c=10.001$ micrometres. The true mean reciprocal thickness is approximately $0.19998\mu {\mathrm{m}}^{-1}$. With $n=1000$, zero rare observations has probability about 0.905 and gives $\hat{m}\approx 0.09999$. Exactly one rare observation instead gives $\hat{m}\approx 1.09989$. The sampling distribution is highly asymmetric despite the large-looking sample count. When every sampled value is the same common thickness, the ordinary sample variance is even zero; it has not revealed the missing component.

Rare-event discovery versus inverse-moment precision distinguishes two experimental objectives. To have at least 95% probability of seeing one rare event requiresor 29,956 observations for this $p$. Estimating the rare fraction to relative standard uncertainty about 10% instead needs roughly one million independent observations, since $u(\hat{p})/p=\sqrt{(1-p)/(np)}$. Seeing one event is not precise quantification of its contribution.

A distributional assumption therefore needs physical support near zero thickness, not just a fit to the sample centre. NIST's gamma distribution and lognormal distribution entries specify distinct lower-tail models; choosing either is a modeling decision.^[ref] Dimensional resolution and independent geometric bounds can be more informative than extending a simulation that never challenges its assumed tail.

Reference-state transfer in overlap normalization uses a known optical state to estimate a state-dependent instrumental offset, then transfers that correction to the specimen measurement. Its central assumption is that the relevant state-gain ratio remains applicable between the reference and specimen. Repeatability of the reference alone does not establish transfer.

Let $u=\log \rho$, $h=\log q>0$ andSuppose specimen and reference observations satisfyWith known reference state $u_R$, the corrected specimen log-ratio isThe residual is exactly the failure of the instrumental offset to transfer. Acquisition order can address time drift, but cannot by itself prove ${\delta }_S={\delta }_R$.

An original matched-state example has $u_S=u_R=0$ and $q=2$. If the specimen's state-gain ratio differs from the reference's by $e^{0.02}$, the corrected ratio still contains 0.02. Ignoring that residual gives an inferred radius ratio approximately 0.9917 instead of one. Repeating both measurements can reduce independent scatter while preserving that systematic offset.

Reference uncertainty must also propagate. Near a matched reference, $D_h^{\prime }(0)=-4\tanh h$. For $q=2$, its magnitude is 2.4. A standard uncertainty 0.01 in the reference log-radius therefore contributes 0.024 to the corrected log-ratio before combining other errors. If the same reference estimate corrects many specimens, this uncertainty is shared; it does not become independent for each specimen. The GUM explicitly retains covariance from common references.^[ref]

Operating-point transfer of an overlap reference requires checking the variables capable of changing the state response: optical power, detector bias, electronic range, frequency and illumination. A useful proposed test brackets specimen readings with reference readings and repeats the reference at several operating powers. A trend then tests a transfer assumption; an unchanged reference at only one power leaves power-dependent transfer untested.

NIST's receiver apparatus controls spatial mode and polarization with a polarizer and single-mode-fibre coupling, while separately calibrating the electrical response.^[ref] A reference must preserve the relevant optical and electrical conditions, or the analysis must model their differences. The transfer calculation above is an original diagnostic construction, not a claimed NIST protocol.

A periodic RC pulse-train baseline is the repeatable pre-pulse capacitor voltage reached after startup. It need not be zero even when the input low level is zero. With high duration $T$, low duration $D$, period $P=T+D$, and $\alpha =e^{-T/\tau }$, $\beta =e^{-D/\tau }$, the recurrence in residual baseline in RC pulse-width discrimination has contraction factor $q=\alpha \beta =e^{-P/\tau }$. Its fixed baseline and peak areThis is a deterministic periodic state, not a random noise floor. The startup error is reduced by a specified fraction $\eta$ after $n\ge (\tau /P)\log (1/\eta )$ cycles, rounded upward.

The mean supplies an independent check. Integrating $\tau \dot{v}+v=u$ over a full period makes the derivative term vanish, so $⟨v⟩=⟨u⟩=(T/P)V_0$. The mean depends on duty cycle while the baseline and ripple depend on period relative to the time constant. For $P\ll \tau$ at fixed $r=T/P$, expansion givesThus a rapid train approaches a level near its duty-cycle average, rather than isolated relaxed pulses.

For an original fast-train example, choose $V_0=1\mathrm{V}$, $\tau =1\mathrm{ms}$, $P=0.20\mathrm{ms}$ and $r=0.30$. The exact baseline is about $0.27929\mathrm{V}$, the peak $0.32127\mathrm{V}$, and the mean exactly $0.30\mathrm{V}$. The startup deviation contracts by $e^{-0.2}$ each cycle; reducing it below one percent of its initial value takes 24 cycles, or $4.8\mathrm{ms}$.

A pulse-generator/oscilloscope test should therefore distinguish the first pulses from the repeatable train, and record minima, maxima and mean rather than only input width. The natural-response measurement in the primary RC laboratory checks the time constant used here.^[ref] Threshold placement relative to the repeatable minimum and maximum determines acceptance and missing comparator reset in an RC pulse train.

Missing comparator reset in an RC pulse train occurs when the filtered output never spends a finite interval below the reset threshold between pulses. In the periodic first-order model, with $\alpha =e^{-T/\tau }$ and $\beta =e^{-D/\tau }$, the baseline is the minimum output and equalsFor a nonhysteretic ideal electronic comparator with $a=V_{\mathrm{th}}/V_0$, $V_b\ge V_{\mathrm{th}}$ leaves no finite low interval. At equality the waveform merely touches the threshold. This is different from a pulse whose peak never reaches threshold and therefore gives no finite high interval.

Solving the strict reset condition givesSeparate rising and falling transitions also require the peak to exceed the threshold. Both inequalities matter, so increasing the quiet interval is not automatically sufficient to produce an accepted output pulse from a weak input.

For an original three-regime example, take $\tau =1\mu \mathrm{s}$, $V_0=1\mathrm{V}$, $T=0.50\mu \mathrm{s}$ and $V_{\mathrm{th}}=0.40\mathrm{V}$. The reset boundary is $D\approx 0.46386\mu \mathrm{s}$. At $D=0.20\mu \mathrm{s}$, the periodic baseline/peak are about $0.63992/0.78160\mathrm{V}$: after startup the electronic comparator does not reset. At $D=1.0\mu \mathrm{s}$, they are $0.18632/0.50648\mathrm{V}$, allowing both transitions. At $D=10\mu \mathrm{s}$, the peak is only about $0.39348\mathrm{V}$, so no finite high interval is produced. The same input height and high duration can therefore give three different output patterns as repetition changes.

The threshold-reset margin for a periodic RC pulse train must also accommodate uncertainty and switching dynamics. With a Schmitt trigger, the lower switching threshold controls reset while the upper one controls acceptance. The capacitor-level inequalities alone do not specify propagation delay or a real electronic comparator's minimum switching dwell time.

A test records the source, capacitor output and electronic comparator output together while sweeping $D$. A steady high output with a pulsing input should be compared with the measured capacitor minimum before being attributed to a source failure. The RC decay model can be checked independently using the primary laboratory's natural-response method.^[ref]

The quiet interval required for an RC baseline timing budget controls how much a residual capacitor baseline may advance the next threshold crossing. In residual baseline in RC pulse-width discrimination, a positive baseline $V_b$ changes the relaxed crossing byFor a permitted advance $\delta >0$, requiring $|\Delta t|\le \delta$ is equivalent to $V_b/V_0\le ϵ$, where $ϵ=1-e^{-\delta /\tau }$. To retain a fresh rising crossing, also require $ϵ<a=V_{\mathrm{th}}/V_0$, equivalently $\delta <-\tau \log (1-a)$.

After a previous peak $V_p$, a zero-input low interval $D$ gives $V_b=V_p{e}^{-D/\tau }$. Hence a measured upper bound on that peak gives the sufficient intervalUsing only the model bound $V_p\le V_0$ recovers the conservative expression in the header, $D\ge -\tau \log (1-e^{-\delta /\tau })$. The peak bound and initial-state assumptions are part of the guarantee. The next pulse must still last long enough to produce the rising crossing; this timing bound alone does not guarantee acceptance.

For an original timing budget, let $\tau =2\mu \mathrm{s}$, $\delta =20\mathrm{ns}$ and $a=0.4$. Then $ϵ\approx 0.0099502<a$, and the conservative interval is $D\ge 9.2203\mu \mathrm{s}$. If a validated peak bound is only $V_p\le 0.60V_0$, it reduces to $8.1987\mu \mathrm{s}$. Simply calling either interval several time constants long would hide the actual timing specification.

For $\delta \ll \tau$, expansion gives $D\approx \tau \log (\tau /\delta )+\delta /2$ under the full-height peak bound. Tightening a small timing allowance by a factor of two therefore adds roughly $\tau \log 2$ of quiet time. If the allowed advance is as large as the whole relaxed crossing time, a separate reset condition is still required; the logarithmic timing budget alone no longer ensures the initial output is below threshold.

A bench test measures the peak and the baseline while sweeping the low interval, keeping the threshold and pulse height fixed. MIT's RC laboratory supplies the natural-decay measurement used to check the underlying time constant.^[ref] Probe capacitance and source resistance belong in that measured time constant.

A correlated-channel update of pressure-map information measures what a new pressure channel adds after accounting for noise that it shares with the old channels. The rank-one update of pressure-map information for independent noise is a special case. An extra channel may be a precise-looking duplicate of an existing noisy signal and contribute almost no new information.

Write the old and new observations as $y=Hp+e$ and $z=a^{\mathsf{T}}p+\eta$. Suppose their known joint error covariance matrix is positive definite,Subtract the part predicted from old noise by definingDirect multiplication gives $\mathrm{Cov}(e,\nu )=0$ and ${\tau }^2=\mathrm{Var}\nu =s-r^{\mathsf{T}}{R}^{-1}r>0$, a Schur complement. Thus the invertible residual transformation separates the weighted residual objective into its old term and one scalar innovation term:For jointly normal random variables, uncorrelatedness here also gives independence, and the matrix is the parameter Fisher information. For more general zero-mean errors with this fixed covariance, the displayed algebra still describes the generalized least squares information matrix. This is a partial whitening transformation; the primary covariance interface documents whitening and its identity-covariance check.^[ref]

An original scalar example observes the same pressure twice, with unit error variances and error correlation $\rho \in (-1,1)$. Then $H=a=R=s=1$, $r=\rho$, soAt $\rho =0.8$, the true estimator variance is $0.9$, whereas an independent-noise calculation would claim $0.5$. At $\rho =-0.8$, the variance is $0.1$: opposite noise contributions can cancel. Correlation can reduce or increase the benefit relative to the independent case; it does not make the correctly computed added information negative.

The effective sensitivity is $b$, not merely the original $a$. If $b=0$ and ${\tau }^2>0$, the new residual carries no additional pressure information. A zero-innovation-variance pressure channel requires separate treatment: division by ${\tau }^2=0$ is invalid, and the limiting relation may be a redundant duplicate or an exact constraint, depending on $b$.

A proposed array characterization measures simultaneous repeated readings under fixed known loads, estimates the complete cross-channel covariance, and tests reconstruction on withheld repetitions. A common reference offset can survive indefinitely despite more channels. The derivation assumes known, unchanged sensitivities and covariance; estimating them from few samples adds uncertainty not included by this update.

The transverse-noise limit of a rotating rank-one covariance model describes how an ordinary two-dimensional noise distribution concentrates onto a line as one physical noise scale tends to zero. For orthonormal vectors ${\mathbf{u}}_{\theta },{\mathbf{v}}_{\theta }$, setA zero-mean Gaussian observation has densityThe pseudoinverse fitting with a parameter-dependent covariance range article explains the range restriction at the exactly singular endpoint.

For a fixed point off the line, ${\mathbf{v}}_{\theta }^T\mathbf{y}\ne 0$, the exponential drives this density to zero as $ϵ\downarrow 0$. On the line the two-dimensional density grows like $1/ϵ$. There is no finite pointwise two-dimensional density limit: probability concentrates onto a one-dimensional Gaussian measure. A density on that line uses a different reference measure. Simply substituting zero for $ϵ$ in the full-rank density formula is invalid.

For an original angular example, let an observed vector have length $R=3$ along orientation ${\theta }_0$, with $\sigma =1$ and $ϵ=0.06$. At angular displacement $\delta =\theta -{\theta }_0$,At $\delta =0.02$ radians the density ratio is about $0.608$. The angular resolution set by transverse spectral noise has characteristic small-angle scale $ϵ/R\approx 0.02$ radians, or $1.15^{\circ }$, when the transverse noise is much smaller than the longitudinal noise.

SciPy's pinned implementation treats singular covariance using its rank and a separate numerical support check, rather than only discarding zero eigenvalues in a quadratic.^[ref] Numerical support tolerance and actual transverse measurement noise are different quantities. Choose the latter from the detector and processing model; a tiny arbitrary regularization can otherwise imply unjustified angular precision.

A photoelastic branch-slip error is an incorrect integer-cycle choice during photoelastic phase unwrapping. If the measured wrapped phase retardation is ${\varphi }_k$ and the true unwrapped value is ${\delta }_k={\varphi }_k+2\pi m_k$, replacing the correct integer $m_k$ by $m_k+j$ changes the recovered phase by $2\pi j$. A small error in one branch decision can therefore become a full-cycle offset in many later observations.

For an original noiseless counterexample, suppose a known starting phase is $2.2\pi$ and the next true phase is $3.4\pi$. Their wrapped values in $(-\pi ,\pi ]$ are $0.2\pi$ and $-0.6\pi$. The nearest-increment rule chooses $-0.8\pi$ instead of the true $1.2\pi$, producing $1.4\pi$ as the second recovered phase. The recovered photoelastic fringe order is $0.7$ rather than $1.7$. Every later small increment can be locally consistent while retaining that wrong cycle.

In the uniform-thickness stress-optical model, $N=t|{\sigma }_1-{\sigma }_2|/f_{\sigma }$, where $f_{\sigma }$ is the material fringe value. While staying on a known nonnegative stress-difference branch, one cycle changes the inferred stress difference by $f_{\sigma }/t$. Stipulated $f_{\sigma }=10\mathrm{N}\mathrm{m}{\mathrm{m}}^{-1}$ and thickness $t=5\mathrm{mm}$ make the example's one-cycle error $2\mathrm{MPa}$. The material coefficient depends on the optical conditions, including wavelength; the primary photoelastic derivation explains that dependence and the analyzer's periodic intensity response.^[ref]

A sufficient local noise margin iswhen each phase error is bounded in magnitude by $ϵ$. The triangle inequality bounds the error in the phase increment by $2ϵ$, leaving it inside the nearest-branch interval. Equality gives no strict guarantee. Reducing the loading increment helps only if the assumed phase recovery and error bound are valid.

The Illinois polariscope provides controlled loading, a load readout and camera observation, which support a proposed check using smaller repeated load increments.^[ref] Agreement between a coarse sequence and a finer sequence tests branch stability over those measurements. It cannot establish an unknown common starting order. A lone crossed-polarizer intensity also does not supply a unique signed phase: the additional phase-recovery information assumed by the algorithm must be measured or independently justified.

The noise-whitened angle between spectral signal and nuisance shapes measures their separation in units of measurement noise rather than raw plotting units. For nonzero shapes $\mathbf{s},\mathbf{x}$ and a positive-definite covariance $C$, choose a square whitening transformation $W$ with $WC{W}^T=I$. The angle isIt is the ordinary Euclidean angle after whitening. Both parallel and antiparallel shapes are indistinguishable up to an unknown nuisance amplitude; the precision penalty depends on ${\sin }^2\varphi$, as in noise penalty from nearly parallel signal and reference shapes.

For an original example, take raw shapes $\mathbf{s}=(1,1{)}^T$, $\mathbf{x}=(1,-1{)}^T$ and $C=\mathrm{diag}(1,100)$. The plotted shapes are orthogonal, but whitening gives $(1,0.1{)}^T$ and $(1,-0.1{)}^T$. Thus $\cos \varphi =0.99/1.01\approx 0.98020$ and $\varphi \approx 11.42^{\circ }$. The second channel, which separates their signs, is noisy. Rejecting the nuisance multiplies the optimum variance by $1/{\sin }^2\varphi =25.5025$, or standard uncertainty by $5.05$.

The invariance of noise-whitened spectral angles to the whitening factor follows from $W^{T}W=C^{-1}$. Different valid square whitening factors can rotate the whitened coordinates without changing the angle. Likewise an invertible change of measurement units transforms both shapes and covariance, leaving the inverse-covariance expression unchanged. Rescaling only the displayed curves while leaving their noise model fixed is a different operation.

The covariance must describe the actual measured channels and shared errors, using the full propagation rule.^[ref] If $C$ is singular, first restrict to its physically allowed range; an unrestricted inverse whitening factor does not exist.

The primary stimulated-Raman experiment described systematic spectral noise and subtraction of a pure-solvent reference.^[ref] The angle above assesses which desired patterns remain distinguishable after such nuisance rejection under a stated covariance. It does not establish that an estimated reference shape or noise model is itself accurate.

Minimum-noise weights under baseline moment constraints choose a linear pulse estimator with the smallest variance among those that preserve a specified signal and remove specified backgrounds. Write the reading vector as $\mathbf{x}=\theta \mathbf{s}+H\mathbf{\beta }+\mathbf{ϵ}$, where $\mathbf{s}$ is a known signal template, the columns of $H$ contain baseline shapes, and the error covariance matrix $C$ is positive definite. For $\hat{\theta }={\mathbf{a}}^T\mathbf{x}$, requireColumns $1,t,t^2,\dots$ in $H$ implement the moment conditions for baseline-drift cancellation. The noise cost of polynomial baseline cancellation depends on which feasible weights are chosen.

If the constraint columns are independent, Lagrange multipliers giveFor any feasible perturbation $\mathbf{v}$, $F^T\mathbf{v}=0$ and ${\mathbf{v}}^{T}C\mathbf{a}=0$, so the variance increases by ${\mathbf{v}}^{T}C\mathbf{v}$. This proves the minimum and connects the estimator to generalized least squares. The covariance propagation is the full quadratic form, not a sum that discards correlations.^[ref]

For an original one-bin pulse example, use five equally spaced times $-2h,-h,0,h,2h$, a signal present only at the centre, and independent errors of variance ${\sigma }^2$. To remove just constant and linear backgrounds, the optimal dimensionless weights are $(-1/4,-1/4,1,-1/4,-1/4)$. Their variance is $5{\sigma }^2/4$, lower than $3{\sigma }^2/2$ for the centre and its two nearest references alone. Multiplying the amplitude estimate by a bin duration $\Delta t$ gives energy variance $5\Delta t^2{\sigma }^2/4$. Extra uncontaminated references can reduce noise when the cancelled background family stays the same.

The inverse requires genuinely independent constraints. Constraint redundancy in baseline-noise minimization calls for removing redundant columns and checking feasibility; a signal template inside the nuisance span cannot be preserved while that span is rejected. NIST's before/after pulse baselines provide the physical measurement context.^[ref] The five-point design is an original idealization and assumes no pulse signal in its reference samples.

Integer frame-count choice under a Raman saturation limit chooses the number of pump-on/off exposure pairs as well as their durations. More pairs incur more readout overhead; too few pairs can saturate before using the available wall-clock budget. Under the model of maximum-exposure constraint in paired Raman acquisition, for positive readout time $d$ and frame cap $c$, each candidate integer $n$ with $2nd<W$ has usable live timeIts unused budget is $W-2nd-T_n\ge 0$. This wall-clock slack in saturation-limited Raman frame allocation is allowed: failure to consume the entire budget does not make a shorter acquisition physically impossible.

For each candidate, allocate $T_{\mathrm{on}}+T_{\mathrm{off}}=T_n$ with each total at most $nc$. Clip the unconstrained on total $T_n\sqrt{\alpha }/(\sqrt{\alpha }+\sqrt{\beta })$ to $[\max (0,T_n-nc),\min (T_n,nc)]$, then evaluate $V_n=\alpha /T_{\mathrm{on}}+\beta /T_{\mathrm{off}}$. This solves the continuous subproblem; comparing the finitely many integer candidates solves the stated frame-count problem. Positive counts and no other costs make additional feasible live time beneficial.

In an original example, use $W=0.60\mathrm{s}$, $d=0.010\mathrm{s}$, $c=0.097\mathrm{s}$, $\alpha =2\times 10^{-6}\mathrm{s}$ and $\beta =0.5\times 10^{-6}\mathrm{s}$. The relevant candidates areFour pairs are globally optimal in this model: one or two pairs use both states at their caps and improve as $n$ increases; three is evaluated explicitly; for every $n\ge 4$ the unconstrained ratio is feasible and $V_n=4.5\times 10^{-6}/(0.60-0.02n)$ increases with $n$. At four pairs the on/off totals are about $0.3467$ and $0.1733\mathrm{s}$.

The primary Raman apparatus used 97 ms exposures to avoid saturation and 10 ms shuttered readout, motivating these stipulated example timings.^[ref] Actual frame choice can also be constrained by read noise, drift and minimum switching cadence. Such costs must be included before treating this shot-noise-only optimum as an experimental operating prescription.

A quadratic Gaussian overlap correction models an estimated logarithmic variance as a nonnegative quadratic function of Gaussian calibration residuals. It gives an exact mean for a corrected overlap while retaining the dependence between the reading and the variance estimate. The model is useful for testing uncertainty propagation software and understanding which shared residuals matter.

Let $Z\sim N(0,I_n)$, $d=a^{\mathsf{T}}Z$ and $\hat{v}=c+Z^{\mathsf{T}}BZ$, with $B$ a symmetric positive semidefinite matrix and $c\ge 0$. For $\tilde{\gamma }=\gamma \exp (d-\hat{v}/2)$, the joint exponential moment in overlap correction givesThe eigenvalues of $B$ describe how strongly each calibration direction enters the variance estimate, while the components of $a$ describe how those same directions enter the reading. Eigenmode alignment in overlap variance correction therefore affects the bias even when the variance estimate's marginal distribution is unchanged.

An original eigenvector calculation writes $B=Q\mathrm{diag}({\lambda }_k)Q^{\mathsf{T}}$ and $\alpha =Q^{\mathsf{T}}a$. Taking logarithms givesThis form avoids taking a determinant and a potentially large exponential separately. It also displays which calibration directions overlap the readout direction. Unbiased variance estimation requires $c+{\sum }_k{\lambda }_k=\Vert a{\Vert }^2$; it does not force the displayed logarithm to vanish.

For a stipulated two-direction example, take $a=(1,0{)}^{\mathsf{T}}$, $c=0.9$ and one nonzero eigenvalue 0.1. If that eigenvalue lies along the reading direction, $B=\mathrm{diag}(0.1,0)$ and the mean ratio is about 0.9578. Rotating it into the independent direction, $B=\mathrm{diag}(0,0.1)$, gives about 1.00235. Both estimates have the form $0.9+0.1Z_1^2$ in distribution, and both have mean one. Only their joint relationship to $d$ changes.

A small-$B$ expansion, with $c=\Vert a{\Vert }^2-\mathrm{tr}B\ge 0$, makes the distinction sharper:The shared-direction term can appear at first order; the independent-direction contribution begins differently. All examples are probability calculations, not measured receiver performance.

In an actual analysis, build $a$ and $B$ from the stated measurement model and reused calibration records. JCGM 101 distinguishes input-distribution assignment from propagation and provides multivariate Gaussian inputs as one particular joint model.^[ref] Agreement of a Monte Carlo method with the exact formula checks implementation conditional on those inputs. It does not establish that Gaussian residuals or a quadratic variance estimator describe the physical experiment.

A floating-pinion ellipsoidal uncertainty bound gives the smallest two-contact compression margin when the uncertain parameter vector lies in a specified ellipsoid. It retains shared parameter effects, as does the floating-pinion affine uncertainty bound, but uses a different admissible-set shape.

Let approach and imbalance be affine functionsand describe the set by $u=u_0+Lz$ with $\Vert z{\Vert }_2\le 1$. Put $\bar{\delta }={\delta }_0+a^{\mathsf{T}}{u}_0$ and $\bar{b}=b_0+c^{\mathsf{T}}{u}_0$. The floating-pinion robust contact margin then satisfies the exact formulaTo prove it, write $\delta -|b|=\min (\delta -b,\delta +b)$. For either linear coefficient vector $q$, the Cauchy-Schwarz inequality givesIf $L^{\mathsf{T}}q\ne 0$, equality is attained at $z=-L^{\mathsf{T}}q/\Vert L^{\mathsf{T}}q{\Vert }_2$. Thus the result supplies both a margin and an attaining direction. The image representation also handles rank-deficient sets without inverting a shape matrix.

For an original example in micrometres, take $u_0=0$, $L$ the identity andOver $u_1^2+u_2^2\le 1$, the two branch minima are $4.5-\sqrt{0.1}$ and $5.5-\sqrt{0.1}$, so the margin is about $4.1838\mu \mathrm{m}$. The first branch attains it at $u=(-1,-3)/\sqrt{10}$. The encompassing box $|u_1|,|u_2|\le 1$ instead gives $4.1\mu \mathrm{m}$ because it admits combinations excluded by the disk.

A covariance-derived ellipse does not automatically certify an absolute engineering bound. Its probability content requires distributional assumptions and an appropriate radius; a hard guarantee requires a defensible admissible set. NASA's face-gear mounting and slow-roll load-sharing observations motivate the physical contact question, not the particular coefficients or set used here.^[ref] The two-contact force law must also remain valid; contact opening requires the floating-pinion contact-loss check.

Resolution-limit conditioning in finite-input detector calibration is the sensitivity of a small inferred detector width to uncertainty in two nearly equal measured widths. In finite-input correction in pulse-detector calibration, the nonnegative normalized convolution model givesThe signed relative sensitivities areThey diverge as the measured output width ${\sigma }_q$ approaches the known input width ${\sigma }_p$. Subtracting two well-measured squared widths can therefore be a poorly conditioned inference.

For an original example, ${\sigma }_q=1.00\mathrm{ns}$ and ${\sigma }_p=0.99\mathrm{ns}$ give $D=0.0199\mathrm{n}{\mathrm{s}}^2$ and ${\sigma }_h\approx 0.1411\mathrm{ns}$. The two relative sensitivity coefficients are about $50.25$ and $-49.25$. A small independent fractional width error can produce a much larger fractional response-width error; a common multiplicative timebase error instead partly cancels because the coefficients sum to one. Correlations must therefore be retained.^[ref]

Near the boundary, report uncertainty in the variance difference before transforming it. For a chosen interval $D\in [-0.010,0.050]\mathrm{n}{\mathrm{s}}^2$, imposing the physical constraint $D\ge 0$ gives a detector-width interval $[0,\sqrt{0.050}]\mathrm{ns}$, approximately $[0,0.2236]\mathrm{ns}$. A coverage-preserving detector-width bound from a variance interval uses this monotone square-root transformation on the physical intersection. It need not resemble a symmetric error bar around a point estimate. An interval wholly below zero contains no admissible width under this model; replacing it silently with a zero-width result hides the inconsistency.

NIST characterized its detector/cable/oscilloscope chain with an optical pulse much shorter than the receiver response.^[ref] That design avoids subtracting two nearly equal broad widths. The variance identity still requires nonnegative profiles and finite second moments; arbitrary full widths at half maximum and responses with negative ringing need a different model.

Capacitive current division in heterodyne receivers occurs when part of a photodiode's generated alternating current flows through shunt capacitance instead of reaching the measurement load. It changes both amplitude and phase. A DC arm-current measurement does not experience the same frequency-dependent division, so the electrical response must be included when checking the gain-matching condition for heterodyne normalization.

In an original lumped-element model, a generated peak-current phasor $I_d$ feeds a capacitance $C_d$ in parallel with a resistive load $R$. Current balance givesThe characteristic frequency $f_c=1/(2\pi R{C}_d)$ marks amplitude $1/\sqrt{2}$, phase $-\pi /4$ and delivered power one half of its low-frequency value for the same generated current. This is a current-divider model; a real receiver can have additional transit-time, amplifier and transmission-line responses.

For stipulated $R=50\Omega$ and $C_d=1$ pF, $f_c\approx 3.183$ GHz. At 1 GHz the amplitude ratio is approximately 0.9540; at 10 GHz it is approximately 0.3033. Those values describe the stated equivalent circuit, not a particular detector specification.

Recovering generated amplitude requires the correctionThe correction also carries calibration uncertainty. Its logarithmic sensitivity to capacitance isAt $x=10$, a small 5% capacitance uncertainty contributes approximately 4.95% relative amplitude uncertainty if other inputs are fixed. Dividing by an accurately known response cannot recover signal-to-noise lost before a fixed downstream noise source; calibration and noise are different issues.

Capacitance uncertainty in heterodyne current recovery should be evaluated at the actual detector bias and frequency model. A proposed diagnostic sweeps beat frequency while holding optical conditions stable, compares magnitude and phase, and separately characterizes load and network response. A fitted single pole is evidence about that circuit approximation, not proof that every observed roll-off is capacitive.

NIST's primary receiver apparatus uses a tunable two-laser beat and separately measures scattering parameters, reflection coefficients, power-sensor calibration and interconnecting-network corrections.^[ref] Those measurements help distinguish detector response from cable or mismatch effects. Retain the current-plane convention in heterodyne calibration when comparing the measured load current with generated photocurrent.

A shared-reference floor in averaged isotope stiffness ratios is the calibration uncertainty that remains when many sample measurements reuse one uncertain reference. In the first-order gain model of covariance of differential spectral-calibration gains, let the error in log stiffness-ratio estimate $i$ be $2(e_i-e_r)$. Assume independent sample errors with variance $s^2$ and one independent reference error with variance $s_r^2$. Averaging $m$ estimates givesThe sample contribution decreases with replication, while the shared reference contribution is unchanged. This is the covariance propagation rule applied to the actual reuse of the reference, rather than treating every ratio as independent.^[ref]

For an original example, choose sample calibration standard uncertainty $s=0.2\%$ and reference uncertainty $s_r=0.1\%$. One ratio has approximately $0.447\%$ standard uncertainty on the log-ratio scale. Sixteen sample measurements reduce it to about $0.224\%$, but the limiting floor is $2s_r=0.200\%$. Reaching within ten percent of that floor requiresso $m\ge 20$ for these chosen numbers. Replication beyond that point has progressively less effect unless the reference uncertainty is reduced.

If a new independent reference is measured for each sample, the model changes: both error contributions can average down, giving $4(s^2+s_r^2)/m$ under the corresponding independence assumptions. Alternatively, averaging several independent reference measurements reduces their random component, while a systematic component shared by all references still remains. A reference-versus-sample allocation in repeated isotope-ratio measurements must distinguish those resources and common influences.

The floor is conditional on the measurement model. Additive spectral calibration errors, reference drift, isotope-composition uncertainty and mode-assignment error are separate contributions. A useful acquisition protocol records which sample spectra share which reference observations; numerical averaging alone cannot remove an uncertainty component reused unchanged in every estimate.