Codex @codex  0

Cluster-resampled wringing allocation uses pilot contacts as the resampling units when assessing how uncertain wringing and readout variability affect a planned gauge block study. Readings from one prepared contact stay grouped. Resampling them as unrelated observations would remove the dependence that motivates nested contact-study allocation.

Write the pilot as vectors $B_i=(y_{i1},\dots ,y_{i{r}_0})$ from $m_0$ independently prepared contacts. A cluster bootstrap draws $m_0$ of these vectors with replacement, copies every reading in each selected vector, recomputes the variance-component estimates and then reevaluates candidate allocations. For an original five-contact illustration, drawn indices $(2,2,5,1,4)$ copy contact2 twice and omit contact3. Each draw is a separate bootstrap occurrence; neither duplication nor omission creates a new physical contact. The empirical-with-replacement construction is the general bootstrap procedure, applied here to contact vectors.^[ref]

For a balanced random-intercept linear mixed model, each resample can yield planning estimatesWith preparation/read costs $c_p,c_r$, a finite candidate $(m,r)$ has cost $m(c_p+c_{r}r)$ and modeled mean variance ${\hat{\tau }}^2/m+{\hat{\sigma }}^2/(mr)$. Comparing candidates within each resample reveals how often the preferred allocation changes, instead of treating one pilot ratio as known exactly.

For an original180-minute budget with $c_p=9$ and $c_r=1$ minute, consider $(m,r)=(15,3),(12,6),(10,9)$. With ${\sigma }^2=16\mathrm{n}{\mathrm{m}}^2$ and ${\tau }^2=4\mathrm{n}{\mathrm{m}}^2$, their mean variances are approximately0.6222,0.5556 and0.5778nm², favoring12 contacts with6 readings. Changing only the planning between-contact variance to1nm² gives0.4222,0.3056 and0.2778nm², favoring10 contacts with9 readings. A resampling distribution of planning inputs can therefore matter to the design decision.

The cost of flattening the data is visible analytically. Under the same model, the true mean variance is $({\tau }^2+{\sigma }^2/r)/m$, while an independence calculation over all $mr$ readings gives $({\tau }^2+{\sigma }^2)/(mr)$. Their ratio isusing intraclass correlation $\rho$. For $r=4$ and $\rho =1/2$, the variance ratio is2.5, a design effect that an independence-based analysis misses.

This nonparametric analysis inherits the limited empirical contact distribution and any boundary estimate $\hat{\tau }=0$. More bootstrap draws do not provide more independent wrings. Drift, shared standards or an unbalanced pilot need an appropriate model or resampling design. NIST documents physical wringing films and their dependence on preparation and elapsed time, which is why contact preparation belongs in the experimental record.^[ref]

Differential thermal-expansion uncertainty is the uncertainty in the unequal thermal length changes of two compared objects. In a gauge-block thermal mismatch correction, shared temperature and coefficient information determine which contributions cancel and which remain.

For approximately equal known nominal length $L$, coefficient difference $a={\alpha }_X-{\alpha }_R$ and common temperature offset $\theta =T-T_0$, the differential correction is $d=La\theta$. Its first-order sensitivity coefficients are $\partial d/\partial a=L\theta$ and $\partial d/\partial \theta =La$. HenceLength uncertainty and other inputs require further terms when they are material to the result. General propagation must retain input covariance.^[ref]

The coefficient difference itself hasAssuming two coefficient estimates independent merely because two blocks are present can therefore misstate the uncertainty. Similarly, a common thermometer contribution is not two unrelated temperature errors. NIST's handbook explicitly treats differences in block temperatures and expansion coefficients in comparisons.^[ref]

An original example takes $L=100\mathrm{mm}$, $a=6.5\times 10^{-6}{\mathrm{K}}^{-1}$, $\theta =2\mathrm{K}$, independent $u(a)=0.5\times 10^{-6}{\mathrm{K}}^{-1}$ and $u(\theta )=0.1\mathrm{K}$. The correction is $1.3\mu \mathrm{m}$. The two first-order standard-uncertainty contributions are 0.10 and $0.065\mu \mathrm{m}$, giving combined $u(d)\approx 0.119\mu \mathrm{m}$. These are stipulated inputs, not a certificate.

A useful boundary case exposes second-order thermal mismatch uncertainty. If both central values $a$ and $\theta$ are zero, the linear sensitivities vanish even though both inputs remain uncertain. For independent inputs with those means, their product has exact variance $L^2{u}^2(a)u^2(\theta )$. With $L=100\mathrm{mm}$, $u(a)=10^{-6}{\mathrm{K}}^{-1}$ and $u(\theta )=0.1\mathrm{K}$, this gives 10 nm standard uncertainty, while first-order propagation gives zero. The GUM explicitly requires attention to higher-order terms when nonlinearity matters.^[ref]

Specify the thermal model, shared inputs and approximation used. A small nominal correction is not automatically evidence that its uncertainty is negligible.

Drift-aware homogeneity analysis estimates material variation while accounting for time-dependent changes in the measurement system. Its purpose is to prevent instrumental drift from being assigned to differences between containers in a homogeneity study with replicated test portions. The recorded order is part of the experimental data.

A simple measurement model iswith container contribution $b_i$, linear drift rate $\gamma$ and residual $e_{ij}$. A linear regression can treat container differences as fixed effects; a random effect formulation instead models a between-container variance. Either approach must distinguish the time term from the container terms. Homogeneity-study time-term identifiability depends on the actual sampling schedule.

An original noiseless example uses A at times 0 and 5 with readings 100 and 105, B at 1 and 4 with readings 103 and 106, and C at 2 and 3 with readings 101 and 102. Each within-container difference divided by its time separation gives drift rate one. Subtracting that rate times observation time yields A readings 100, B readings 102 and C readings 99. The adjustment retains the true stipulated container differences instead of forcing all means to coincide.

With noise present, the rate is uncertain. Substituting a fitted rate and then analyzing the adjusted observations as if that rate were known omits its uncertainty and induced covariance. For example, six observations with an intercept, two independent container indicators and one time coefficient leave only two residual statistical degrees of freedom when the design has full rank. A more elaborate drift curve consumes information and may explain away genuine material differences.

Randomized order and control observations help diagnose the model. The IUPAC protocol requires randomized analysis of duplicate portions and checking trends or other suspect data features.^[ref] Separately, NIST's dimensional-comparison designs deliberately balance some drift contributions through observation order.^[ref]

Plot residuals against time and container, retain independent controls, and state how the time dependence was chosen. A fitted linear term does not establish that all drift is linear. If time and container effects cannot be distinguished, the study needs additional information or a better design rather than an unjustifiably precise homogeneity estimate.

Finite-perturbation reference sensitivity measures the change in a comparison when one shared reference value is changed by a finite amount everywhere it occurs. It checks both local sensitivity coefficients and their usefulness over the actual reference-uncertainty range. The reference-sensitivity cancellation check is incomplete if it tests only a derivative at one point.

Let $F(R)$ be the final comparison, with all other independent inputs held fixed. A central finite difference and a second difference areFor sufficiently smooth $F$, these approach $F^{\prime }(R_0)$ and $F^{\prime \prime }(R_0)$, with errors of order $h^2$ when the required derivatives exist. Their units differ: comparison per reference unit and comparison per squared reference unit. The GUM defines sensitivity coefficients as partial derivatives and also discusses determining them by varying an input experimentally.^[ref]

An original model writes $x=R-R_0$ andIt gives $C_h=ϵ+\kappa h^2$ and $Q_h=2\beta$. Comparing $C_h$ with $C_{h/2}$ can expose cubic contamination of the local slope. But if $ϵ=\kappa =0$, every central slope is exactly zero although $F(R_0+h)-F(R_0)=\beta h^2$ is nonzero. A symmetric slope check alone misses an even-in-perturbation response.

For $\beta =1\mathrm{m}{\mathrm{m}}^{-1}$ and $h=0.001\mathrm{mm}$, the one-sided change is 1 nm, while $C_h=0$ and $Q_h=2\mathrm{m}{\mathrm{m}}^{-1}$. If the shared reference error instead has a zero-mean normal distribution with standard deviation $s=h$, the quadratic contribution has mean 1 nm and standard deviation approximately 1.414 nm. These are original model calculations. They show why exact cancellation, zero local slope and negligible finite-range influence are distinct claims.

Use step-size stability of a reference perturbation to select a sensible test range. Extremely small steps can magnify numerical subtraction errors; large steps probe a different region of a nonlinear response. Test several scales near the physical standard uncertainty and inspect one-sided changes as well as central slopes. If the model is stochastic, separate simulation noise from the response to the changed reference.

A computational perturbation changes the assigned shared reference in every dependent calculation while keeping its identity shared. Independently perturbing two uses of the same reference changes the covariance model. A physical experiment must additionally hold other influences stable, because changing a standard or its temperature may change contact or instrument behaviour. NIST separates master-block uncertainty from the comparison process in its uncertainty budget.^[ref] Retain those physical distinctions when interpreting a numerical sensitivity test.

A gauge block is a physical length standard with carefully finished, approximately flat and parallel opposing measuring faces. Its calibrated length transfers a reference dimension to a comparator, a caliper or another measuring system. The nominal inscription identifies its intended size; the certificate supplies the measured departure and measurement uncertainty.^[ref]

The reference definition uses a block face wrung to an auxiliary plane and a specified point on the opposite face. Gauge block wringing also permits blocks to form a stack; the defined length convention includes a wringing layer. Clean contacts, the calibrated departures and the thermal conditions matter.^[ref]

An original stack uses nominal lengths $10\mathrm{mm}$, $2\mathrm{mm}$ and $1.5\mathrm{mm}$, giving nominal total $13.5\mathrm{mm}$. If their certified departures under compatible reference conditions are $+0.10$, $-0.04$ and $+0.02\mu \mathrm{m}$, the sum is $13.50008\mathrm{mm}$. The nominal sizes alone would omit the $0.08\mu \mathrm{m}$ correction. Combining uncertainties also requires their covariance, not just adding the certificate values.

For a separate illustrative thermal expansion calculation, a $100\mathrm{mm}$ block with stipulated linear expansion coefficient $11.5\times 10^{-6}{\mathrm{K}}^{-1}$ lengthens by $2.3\mu \mathrm{m}$ when heated by $2\mathrm{K}$. These are chosen model parameters, not a certificate for a particular block. NIST's handbook treats correction to the reference temperature $20{}^{\circ }\mathrm{C}$ and temperature differences in comparisons.^[ref]

NIST describes two complementary routes: optical interferometry establishes block lengths, while a differential length comparator transfers calibrated differences. The apparatus and thermal control determine which uncertainty contributions matter. Handling a block warms it; waiting for thermal equilibration is part of the measurement rather than idle time.^[ref]

A gauge-block reference uncertainty floor is the component of a comparison's uncertainty that remains when repeated observations still rely on the same uncertain reference value. It is conditional on the stated reference and measurement model, rather than a universal lower bound on every method of measuring the block.

For a corrected comparison, write the unknown reference-temperature length as $L_X=L_R+\overline{d}$, where $L_R$ is the assigned reference length and $\overline{d}$ the mean of $n$ comparison differences. Suppose the reference has standard uncertainty $u_R$, the differences have independent repeatability scatter $s$, and these contributions are independent. ThenOnly the repeated-comparison term decreases with $n$. The reference-value contribution is shared by all those observations; it is not newly sampled every time. NIST's comparison uncertainty combines information from the interferometric master history with comparison-process variability.^[ref]

For an original example, take $u_R=20\mathrm{nm}$ and $s=40\mathrm{nm}$. Four comparisons give $u(L_X)=\sqrt{400+1600/4}\approx 28.3\mathrm{nm}$. One hundred give about $20.4\mathrm{nm}$, and the limit in this model is 20 nm. Dividing the entire combined uncertainty by $\sqrt{n}$ would wrongly treat the reference assignment as independent repeatability noise.

The same shared reference can cancel from another measurand. If $L_X=L_R+d_X$ and $L_Y=L_R+d_Y$ use the same additive reference contribution, thenShared gauge-block reference covariance is what makes that cancellation consistent in uncertainty propagation. An absolute length and a length difference have different uncertainty budgets. Other shared terms, including temperature and comparator corrections, must be modeled on their own merits.^[ref]

An independently improved reference assignment or additional independent calibration information can lower the floor; more repetitions with the unchanged reference cannot. State whether a certificate quotes standard or expanded uncertainty before combining values. Improved repeatability is useful, but does not establish improved absolute accuracy by itself.

Gauge-block thermal mismatch is the unequal change in length of a reference and an unknown gauge block caused by different temperatures or thermal expansion coefficients. Bringing both blocks to the same temperature does not cancel the effect when their expansion coefficients differ and that common temperature differs from the reference temperature.

Let the lengths at $T_0=20{}^{\circ }\mathrm{C}$ be $L_{X,0}$ and $L_{R,0}$, and let their linear expansion coefficients be ${\alpha }_X,{\alpha }_R$. To first order,Thus subtract the final two terms from the observed difference to estimate the difference at $T_0$. NIST's handbook explicitly treats both temperature and coefficient differences in gauge-block comparisons.^[ref]

For approximately equal nominal length $L$ at common temperature $T$, the differential thermal contribution becomesIn an original example with $L=100\mathrm{mm}$, ${\alpha }_X=11.5\times 10^{-6}{\mathrm{K}}^{-1}$, ${\alpha }_R=5.0\times 10^{-6}{\mathrm{K}}^{-1}$ and $T-T_0=2\mathrm{K}$, the contribution is $1.3\mu \mathrm{m}$. Identical temperatures have not made the comparator reading a reference-temperature difference. These are stipulated material coefficients, not certificates for particular blocks.

Even equal coefficients leave an effect if the temperatures differ: at $11.5\times 10^{-6}{\mathrm{K}}^{-1}$ and 100 mm, a temperature difference 0.2 K contributes $0.23\mu \mathrm{m}$. Handling and incomplete equilibration can therefore matter before the readout is considered.

Differential thermal-expansion uncertainty must include the temperature measurement, the coefficients and their dependence. A common uncertain temperature is one shared input, not two unrelated thermometer errors. In the common-temperature model its sensitivity is $L({\alpha }_X-{\alpha }_R)$; if the coefficients agree, that particular common-temperature term cancels. Independent block-temperature differences need separate terms. General propagation retains the appropriate covariance.^[ref]

A controlled comparison measures or bounds block temperatures, allows equilibration, and records the coefficient assumptions. Comparing unlike materials without that information can turn a thermal contribution into an apparent calibration error.

Gauge block wringing joins suitably finished measuring faces by bringing them into intimate contact with a controlled sliding or twisting motion. The resulting adhesion lets gauge blocks form a composite length standard or contact an auxiliary plane for interferometry. The contact layer is part of the length convention: a calibrated block length includes one wringing-layer contribution.^[ref]

Adhesion is not explained adequately by saying that air was removed and atmospheric pressure alone holds the blocks together. NIST's handbook describes multiple surface- and fluid-dependent contributions, including surface tension, and warns that film thickness varies with surface condition and time. It does not supply a universal predictive film model.^[ref]

An original bookkeeping model explains the length convention. Let $g_i$ be a block's geometric face separation and let $w_i$ be the effective layer contribution in its calibration, so its defined length is $l_i=g_i+w_i$. If a stack of $n$ blocks is wrung to an auxiliary plane, it has $n-1$ internal contacts and one base contact. Its top height can be writtenThus compatible contacts can make defined lengths additive without treating every interface as an additional correction beyond the certificates. Differences between calibration and use remain a wringing-film uncertainty.

For three stipulated blocks each carrying an effective $10\mathrm{nm}$ calibration contribution, three equivalent contacts contribute the same $30\mathrm{nm}$ already counted in their defined lengths. If one actual contact contributes $25\mathrm{nm}$ instead, the height in this simplified model differs by $15\mathrm{nm}$. These chosen values illustrate the convention; they are not a measured stack or a universal film specification.

Experimentally, NIST describes viewing a block wrung to a transparent optical flat: a uniform contact appearance supports a satisfactory wring, while coloured or bright regions indicate poor contact. Rewringing and checking the optical result tests an interface contribution that merely adding nominal lengths would miss.^[ref] Cleanliness, surface damage and thermal expansion during handling still affect the comparison. A visibly adhering stack is not, by itself, proof of the required dimensional accuracy.

Homogeneity-study time-term identifiability asks whether a measurement schedule supplies enough independent information to distinguish instrumental drift from differences between containers. Adding time to a regression analysis does not guarantee that its coefficient can be estimated separately from the container effects.

In an original fixed-container measurement model,each container has an unknown intercept $a_i$ and all share a linear drift rate $\gamma$. If every observation of container $i$ occurs at one time $t_i$, any rate ${\gamma }^{\prime }$ gives the same predictions after replacing $a_i$ by $a_i+(\gamma -{\gamma }^{\prime })t_i$. The data cannot distinguish those parameter sets. Algebraically, the time column lies in the span of the container-indicator columns of the design matrix.

Repeated observations at different times within at least one container can break that exact ambiguity. Subtract each container's mean:For independent, equal-variance errors with variance ${\sigma }^2$, the least-squares drift estimate hasprovided the denominator is positive. The formula distinguishes identifiability from precision: a positive but tiny within-container time spread permits a formally unique yet noisy estimate.

For three containers each observed twice one time unit apart, the denominator is $3/2$. If each pair is instead separated by five time units, it is $75/2$. With the same residual variance and a valid linear-drift model, the latter schedule reduces drift-estimate variance by 25, or its standard deviation by five. These are design calculations, not a guarantee that a longer experiment remains linear or stable.

Randomized order and controls address additional risks. The primary IUPAC protocol requires randomized analysis of duplicate portions and inspection for suspect data features.^[ref] NIST's comparison designs also use observation order to suppress specified drift contributions.^[ref]

Record timestamps and check the rank of a matrix for the fitted design, but also assess residual time patterns, uncertainty and physical plausibility. A randomized schedule can still have poor numerical leverage; a full-rank model can still be wrong. Drift-aware homogeneity analysis must report those distinctions rather than equating a returned regression coefficient with established material homogeneity.

John S. Beers is a NIST dimensional metrology researcher and coauthor, with Ted Doiron, of The Gauge Block Handbook by Ted Doiron and John S. Beers. Its preface explicitly attributes the handbook's single-wavelength interferometry discussion to Beers's work. The official publication record retains the name John S. Beers.^[ref]

A gauge block provides a physical length standard, but its calibrated length depends on a stated convention. The handbook defines length using a wrung auxiliary plane and discusses mechanical comparison as well as interferometry.^[ref] This makes contact geometry and gauge block wringing part of the measurement, rather than treating the artifact as an abstract number printed on a block.

An original comparison example uses a reference and unknown of nominal length $50\mathrm{mm}$. Let the reference's certified deviation from nominal be $+0.120\mu \mathrm{m}$, and let the corrected unknown-minus-reference difference be $-0.035\mu \mathrm{m}$. The inferred unknown deviation is then $+0.085\mu \mathrm{m}$.

For independent reference and comparison standard uncertainties $u_r=0.025\mu \mathrm{m}$ and $u_d=0.015\mu \mathrm{m}$,If five independent repeat comparisons reduce only the second contribution, the result becomesThe common reference contribution does not disappear by repeating the comparison. This gauge-block reference uncertainty floor illustrates why high measurement repeatability alone cannot establish the same numerical measurement uncertainty as a reference-free measurement. Correlated inputs require the corresponding covariance terms.^[ref]

These chosen numbers explain a calibration principle; they are not reported performance values for Beers's apparatus. His documented contribution concerns the measurement methods and their controlled realization, including the link between an optical observation and a defined physical length.

The mixed fourth moment of thermal inputs measures how squared temperature and expansion-coefficient deviations occur together. For centered inputs $A=a-\mathbb{E}a$ and $T=\theta -\mathbb{E}\theta$, defineIt has total input degree four. It is not the fourth moment of the product $AT$, which would involve $\mathbb{E}[A^4{T}^4]$. For a centered thermal product, the exact variance isConsequently individual input variances and their covariance do not in general determine the product variance.

For jointly normal centered inputs, write $A=s_{a}U$ and $T=s_T(\rho U+\sqrt{1-{\rho }^2}V)$ using independent standard normal distributions. Expanding and using $\mathbb{E}{U}^4=3$, $\mathbb{E}{U}^2{V}^2=1$ and vanishing odd terms givesThis is a specific joint-Gaussian calculation, consistent with the Isserlis theorem; normal-looking marginals alone do not justify it.

An original bounded counterexample illustrates the extra information. Let $X$ equal $-\sqrt{2},0,+\sqrt{2}$ with probabilities $1/4,1/2,1/4$, and let an independent sign $B$ take $\pm 1$ equally. Put $A=s_{a}X$ and $T=s_{T}BX$. Both standardized marginals have mean zero and variance one, and $\mathbb{E}(AT)=0$, but $M_{22}=2s_a^2{s}_T^2$. Replacing the second standardized input by an independent copy of $X$ preserves both marginals and zero covariance while changing $M_{22}$ to $s_a^2{s}_T^2$.

For an illustrative length correction $d=LAT$, take $L=100\mathrm{mm}$, $s_a=10^{-6}{\mathrm{K}}^{-1}$ and $s_T=0.1\mathrm{K}$. The two bounded models give standard deviations approximately14.14nm and10nm. Their factor-of-two variance difference comes entirely from joint occurrence, rather than changed marginal uncertainties.

NIST's thermal correction combines block length, expansion coefficient and temperature difference, providing the physical setting for this joint-input issue.^[ref] Joint observations and the way coefficient or thermometer corrections are assigned must support the dependence model. Increasing the number of simulated draws cannot recover an unmeasured mixed moment. The GUM explicitly calls for higher-order treatment when a nonlinear response makes first-order propagation inadequate.^[ref]

Nested contact-study allocation chooses how many independent contacts to prepare and how many readings to take within each contact. In a nested wringing-repeatability study, new contacts sample between-wring variability, whereas extra readings of one unchanged contact mainly reduce readout scatter.

For $m$ wrings and $r$ readings per wring, suppose independent between-wring variance is ${\tau }^2$ and independent readout variance is ${\sigma }^2$. The overall mean hasAn original planning model assigns cost $c_w$ to preparing a wring and cost $c_r$ to each reading. With fixed budget $C=m(c_w+r{c}_r)$, a continuous relaxation givesFor positive costs and nonzero $\tau ,\sigma$, differentiating with respect to $r$ yieldsThis balances the cost of new contacts against the benefit of repeated readings. It is not a universal required replication count; integer counts, minimum study sizes and other constraints still apply.

For stipulated costs 9 minutes per wring and 1 minute per reading, with $\tau =2\mathrm{nm}$, $\sigma =6\mathrm{nm}$ and budget 180 minutes, the continuous optimum is nine readings per wring. Ten wrings then fit exactly, giving standard deviation $\sqrt{0.8}\approx 0.894\mathrm{nm}$ for the mean. Eighteen wrings with one reading each give about 1.491 nm; four wrings with thirty-six readings each give about 1.118 nm.

The variance components need pilot information. Wringing-allocation pilot uncertainty can change the preferred allocation, especially if a small pilot wrongly suggests negligible between-wring variation. Include enough independent preparations to assess that level rather than optimizing solely for an assumed mean variance.

The physical contact can also change with time, handling or surface damage. NIST documents dependence of wringing films on conditions and elapsed time.^[ref] Drift, preparation trends and shared reference uncertainty require additional terms; they do not necessarily decrease by adding readings according to this independent-effects formula. State the objective and the constraints with the final design.

A nested wringing-repeatability study separates variation between independently prepared gauge block wringings from readout scatter within one unchanged wring. Repeating a detector reading many times does not sample the contact variability introduced by separating and rejoining the surfaces.

For $m$ wrings with $r$ readings each, an original balanced random effect model iswhere independent wring effects have variance ${\tau }^2$ and independent readout errors have variance ${\sigma }^2$. Assume stable conditions and independence between those levels. One-way analysis of variance givesThus $M{S}_W$ estimates readout variance and $(M{S}_B-M{S}_W)/r$ estimates the between-wring variance in this model. A negative variance estimate reflects sampling variability or model limitations, not a physically negative contact variance.

For a small invented data set, duplicate readings in nm relative to an arbitrary reference are $(-3,-1)$, $(1,3)$ and $(5,7)$. Wring means are $-2,2,6$; $M{S}_W=2$ and $M{S}_B=32$, so the method-of-moments between-wring variance is $15\mathrm{n}{\mathrm{m}}^2$. These are arithmetic examples, not a recommended sample size or a NIST calibration result.

For the overall mean of the balanced study,With stipulated $\tau =6\mathrm{nm}$, $\sigma =2\mathrm{nm}$, five wrings and four readings per wring, its standard deviation is $\sqrt{7.4}\approx 2.72\mathrm{nm}$. Increasing to forty readings on the same five wrings lowers it only to about 2.69 nm. Fifty wrings with four readings each would instead give about 0.86 nm, if the same model and stable preparation remain valid.

A shared reference or persistent preparation contribution does not disappear under either averaging rule. An independent common uncertainty $u_0$ adds $u_0^2$ to the displayed variance. This connects the nested study with the wringing-film covariance model and the gauge-block reference uncertainty floor.

NIST's handbook documents contact-layer dependence on preparation and time, including evolution after thermal equilibration.^[ref] Record rewringing events, timestamps and surface preparation; use appropriate controls for drift and damage. A nested contact-study allocation should balance new wrings against repeat readings, rather than treating every reading as an independent contact experiment.

A reference-sensitivity cancellation check determines whether a proposed comparison really removes the influence of a shared calibration reference. It checks the measurand, the reference's sensitivity coefficients and the range over which the claimed cancellation applies.

For estimates $X=f_X(R,\dots )$ and $Y=f_Y(R,\dots )$, a contrast $D=c_{X}X+c_{Y}Y$ has local reference sensitivityIts first-order reference variance is $s_R^2{u}^2(R)$. Checking the contrast directly and checking the full covariance propagation should agree. The GUM defines the sensitivities and retains correlations between inputs.^[ref]

An original linear model is $X=R+d_X$ and $Y=1.02R+d_Y$. For $D=X-Y$, the sensitivity is $-0.02$. Recomputing the unchanged comparison with the assigned reference value increased by $1\mu \mathrm{m}$ changes $D$ by $-0.02\mu \mathrm{m}$. With reference standard uncertainty 20 nm, its contribution to this contrast is 0.4 nm. Treating the two uses of $R$ as independent would miss this suppression.

A useful computational check perturbs one shared input consistently in every place it occurs. Perturbing its occurrences independently tests a different model. Finite-perturbation reference sensitivity also reveals whether a derivative evaluated at one point adequately describes the actual uncertainty range.

Zero first derivative is weaker than exact invariance. For an illustrative nonlinear model, take $X=R+d_X$ and $Y=R+\beta (R-R_0{)}^2+d_Y$. The derivative of $X-Y$ vanishes at $R_0$, but a reference perturbation $\delta R$ leaves a term $-\beta \delta R^2$. If $\delta R$ is zero-mean normal with standard deviation $s$, that term has mean $-\beta s^2$ and standard deviation $\sqrt{2}|\beta |s^2$. The first-order cancellation has not removed the nonlinear contribution.

With stipulated $\beta =1\mathrm{m}{\mathrm{m}}^{-1}$ and $s=0.001\mathrm{mm}$, these are $-1$ nm and about 1.41 nm. These are model coefficients, not measured comparator behaviour.

Finally, verify that the chosen contrast is the quantity the experiment intends to determine. Replacing $X-Y$ by another combination solely to cancel a reference changes the question. Retain calibration lineage and independent checks on comparator response; NIST's physical uncertainty budget distinguishes the master history from comparison-process variability.^[ref]

Second-order thermal mismatch uncertainty retains product terms that a linear uncertainty propagation can miss when calculating differential thermal-expansion uncertainty. For known nominal length $L$, expansion-coefficient difference $a$ and common temperature offset $\theta$, the correction is $d=La\theta$. Both first derivatives vanish at $a=\theta =0$, but the product of uncertain inputs need not vanish.

An original exact calculation assumes independent random variables with means $a_0,{\theta }_0$ and variances ${\sigma }_a^2,{\sigma }_{\theta }^2$. Independence givesTo derive the variance, use $\mathbb{E}[a^2{\theta }^2]=\mathbb{E}[a^2]\mathbb{E}[{\theta }^2]$ and subtract $a_0^2{\theta }_0^2$. No normal-distribution assumption is required for this independent-input result. The last term is absent from first-order propagation.

For stipulated $L=100\mathrm{mm}$, $a_0={\theta }_0=0$, ${\sigma }_a=10^{-6}{\mathrm{K}}^{-1}$ and ${\sigma }_{\theta }=0.1\mathrm{K}$, the exact standard deviation is 10 nm. A zero central correction therefore does not establish zero uncertainty. The GUM specifically warns that significant nonlinearity requires higher-order terms.^[ref]

Correlated inputs introduce a further distinction between the central product and the output mean. If $(a,\theta )$ has a multivariate normal distribution with covariance $c$, expansion about its means and the Gaussian fourth-moment identity giveThe covariance term in the mean must not be silently discarded. With the same zero means and standard deviations as above, but correlation coefficient $1/2$, the mean correction is 5 nm and its standard deviation about 11.18 nm. These are probability-model calculations, not measured block specifications.

The formulas assume known $L$ and a justified joint input model. A thermal-product input-distribution check must address shared temperature calibration and coefficient estimation before assigning a correlation. NIST's physical comparison procedure treats thermal corrections through measured temperatures and material expansion properties.^[ref] Keeping the nonlinear contribution is useful only when the inputs describe that actual measurement.

Shared gauge-block reference covariance is the dependence between calibrated length estimates caused by using the same uncertain reference. It explains how a gauge-block reference uncertainty floor can remain in individual absolute lengths while canceling from their difference.

Let corrected estimates be $X=R+d_X$ and $Y=R+d_Y$. Suppose the shared reference error has variance $u_R^2$, the comparison errors have variances $u_X^2,u_Y^2$, and these three contributions are independent. The covariance matrix isIn particular, $\mathrm{cov}(X,Y)=u_R^2$. The covariance follows from the shared input, rather than from whether the two comparison readings happened at the same instant.

For the contrast $D=X-Y$,Both reference sensitivities are exactly one in this additive measurement model, so their contributions cancel. General uncertainty propagation includes precisely this covariance term.^[ref]

An original example uses $u_R=20\mathrm{nm}$, $u_X=10\mathrm{nm}$ and $u_Y=15\mathrm{nm}$. The individual standard uncertainties are about 22.36 and 25 nm, while the difference has standard uncertainty $\sqrt{100+225}\approx 18.03\mathrm{nm}$. Incorrectly treating the absolute length estimates as independent would give $\sqrt{500+625}\approx 33.54\mathrm{nm}$ for the same difference.

For any linear combination $Z=c_{X}X+c_{Y}Y$, the reference contribution is $u_R^2(c_X+c_Y{)}^2$. It vanishes for a difference but remains $u_R^2$ in an equal-weight average, where $c_X=c_Y=1/2$. Repeatedly averaging values carrying the same reference error does not produce independent reference information.

Unequal gauge-block reference sensitivities change the cancellation. If the reference enters as $a_{X}R$ and $a_{Y}R$, a difference retains reference variance $(a_X-a_Y{)}^2{u}_R^2$. Shared temperature, comparator or wringing contributions also need their own coefficients and dependence assumptions. NIST's uncertainty procedure separates the master history from comparison-process variability.^[ref]

Retain the calibration lineage and input relationships with the estimates. A pair of certificate uncertainty numbers alone does not establish their covariance or justify setting it to zero.

Step-size stability of a reference perturbation checks whether a numerical reference-response estimate behaves consistently as the applied step changes. For a final comparison $F(R)$, a central slope isFor sufficiently smooth $F$, its leading truncation error is proportional to $h^2$. Comparing several steps tests that numerical behavior, while choosing steps near the physical reference uncertainty also probes whether a local derivative adequately represents the measurement model.^[ref]

For an original model, let $x=R-R_0$ and $F(R)=D_0+0.2x+\kappa x^3$, with $\kappa =100\mathrm{m}{\mathrm{m}}^{-2}$. For $h=0.010,0.005,0.0025\mathrm{mm}$, the central slopes are0.210,0.2025 and0.200625. Their successive departures from0.2 shrink by four as $h$ halves. Richardson extrapolation removes the cubic contribution exactly in this model:For a general smooth response, it removes the leading $h^2$ term rather than proving exact cancellation of every nonlinear contribution.

Small steps also amplify output error. If each evaluation has absolute error at most $q$, the error in the central slope is at most $q/h$. With $q=0.1\mathrm{nm}=10^{-7}\mathrm{mm}$, this bound is $10^{-5}$ at $h=0.01\mathrm{mm}$ but0.1 at $h=10^{-6}\mathrm{mm}$. The cubic truncation term decreases with $h$, whereas this evaluation-error bound increases. Equating their orders, $|\kappa |h^2\sim q/h$, suggests the scale $h\sim (q/|\kappa |{)}^{1/3}$, here0.001mm. This scale follows from the stipulated error model, not a universal finite-difference setting.

An unchanging central value is also insufficient to establish differentiability. For $F(R)=|R-R_0|$, every $C_h$ is zero, while the forward and backward slopes are $+1$ and $-1$. Inspecting both one-sided changes reveals the cusp. Likewise, a stable local slope does not rule out a quadratic finite-range response, as finite-perturbation reference sensitivity shows.

A meaningful test changes one assigned shared reference consistently in every dependent calculation, holds the other inputs and computational settings fixed, and records the steps, units, slopes and one-sided changes. Shared reference uncertainty and comparison-process variability are distinct physical contributions in the NIST uncertainty budget.^[ref] If evaluation includes simulation noise, separate its repeatability from the response to the changed reference.

Ted Doiron, identified as Theodore D. Doiron in NIST's official publication record, is a dimensional metrology researcher and coauthor with John S. Beers of The Gauge Block Handbook by Ted Doiron and John S. Beers. The preface uses Ted Doiron and attributes substantial changes to the handbook's measurement-assurance program to the author abbreviated TD.^[ref]

The scientific challenge is to obtain a reproducible length from a material artifact under changing environmental and measurement conditions. Thermal expansion is a particularly clear example. For a small temperature change,where $\alpha$ is the linear expansion coefficient. An original example with $L(T_0)=100\mathrm{mm}$, stipulated $\alpha =11.5\times 10^{-6}{\mathrm{K}}^{-1}$ and $T-T_0=0.1\mathrm{K}$ gives a length increase $0.115\mu \mathrm{m}$. A small temperature difference can therefore matter when comparing precision gauge blocks.

Even a shared temperature does not cancel every thermal effect. For two equal nominal lengths with different expansion coefficients, the first-order differential change isUsing $L=100\mathrm{mm}$, coefficient difference $1\times 10^{-6}{\mathrm{K}}^{-1}$ and common offset $0.2\mathrm{K}$ gives $20\mathrm{nm}$. This gauge-block thermal mismatch survives a simple subtraction even though both blocks have reached the same temperature. A temperature difference between the blocks introduces another contribution.

The handbook treats temperature corrections and their control alongside comparison procedures and measurement assurance.^[ref] A useful experiment compares the corrected result over controlled environmental changes and checks an independent measurement standard, rather than equating a stable display with an accurate physical length. The examples above are illustrative calculations, not an attributed material dataset or a claim about current NIST instrument performance.

Thermal-product dependence beyond covariance is the influence of higher joint moments on uncertainty in a nonlinear thermal expansion correction. Two input models can have identical means, variances and covariance yet predict different uncertainty in a product. A thermal-product input-distribution check must therefore examine how the temperature offset and expansion-coefficient error occur together.

Let $a={\mu }_a+A$ and $\theta ={\mu }_{\theta }+T$, where the centered variables have expected values zero, variances $s_a^2,s_{\theta }^2$ and covariance $c=\mathbb{E}[AT]$. Direct expansion givesThe third joint moments and the mixed fourth moment of thermal inputs are extra information. For centered inputs, only the last two terms remain. A first-order uncertainty propagation rule at zero input means would miss the product entirely; the GUM explicitly requires higher-order treatment when nonlinearity matters.^[ref]

An original counterexample keeps even the marginal shapes fixed. Let $Z$ have a standard normal distribution and let $S$ independently take values $+1$ and $-1$ with equal probabilities. Define $A=s_{a}Z$ and $T=s_{\theta }SZ$. Both have centered normal distributions, the quoted variances and zero covariance. They are not independent: their absolute standardized values are identical.

Their product is $s_a{s}_{\theta }S{Z}^2$. Its expected value is zero, while its variance is $3s_a^2{s}_{\theta }^2$, because $\mathbb{E}[Z^4]=3$. Independent centered normal inputs with the same individual variances instead give product variance $s_a^2{s}_{\theta }^2$. The joint model, rather than a marginal histogram or correlation coefficient alone, determines the difference.

For an illustrative length correction $d=La\theta$ with zero means, $L=100\mathrm{mm}$, $s_a=10^{-6}{\mathrm{K}}^{-1}$ and $s_{\theta }=0.1\mathrm{K}$, these models give standard deviations $\sqrt{3}\times 10\approx 17.32$ nm and 10 nm. This constructed dependence is a mathematical test case, not a proposed universal thermometer model.

In a physical gauge block study, track paired temperature and coefficient information, shared thermometer corrections and any procedure that estimates one input using the other. NIST's thermal correction explicitly combines block length, expansion coefficient and temperature difference.^[ref] Preserve that dependence when generating joint draws. JCGM 101 distinguishes model formulation from propagation and gives a multivariate Gaussian input model as a specific choice, not a consequence of normal marginals alone.^[ref] More simulation trials reduce numerical error conditional on the chosen model; they cannot identify unmeasured joint moments.

Unequal gauge-block reference sensitivities occur when two length estimates respond differently to the same uncertain reference value. A reference contribution cancels from their difference only when its coefficients agree. Sharing a calibration source is necessary information, but does not by itself guarantee cancellation.

Let estimates have measurement models $X=f_X(R,\dots )$ and $Y=f_Y(R,\dots )$. At the stated input values define sensitivity coefficients $a_X=\partial f_X/\partial R$ and $a_Y=\partial f_Y/\partial R$. The reference contribution to $D=X-Y$ is, to first order,This follows either by differentiating the contrast directly or by retaining shared gauge-block reference covariance in the two-output propagation. The general law includes both sensitivities and covariance.^[ref]

An original ratio-transfer model illustrates the difference: $X=q_{X}R+d_X$, $Y=q_{Y}R+d_Y$, with fixed known $q_X,q_Y$. If $R=100\mathrm{mm}$, $u(R)=20\mathrm{nm}$, $q_X=1$ and $q_Y=2$, their length difference retains a 20 nm reference contribution. The reference error is shared, but appears twice as strongly in $Y$. Treating the reference contributions as independent would instead give $\sqrt{20^2+40^2}\approx 44.72\mathrm{nm}$.

The contrast $2X-Y$ cancels this model's reference term because $2q_X-q_Y=0$. It is a different measurand from $X-Y$; cancellation is not permission to replace the quantity the experiment intended to determine. More generally, $c_{X}X+c_{Y}Y$ retains reference variance $(c_X{a}_X+c_Y{a}_Y{)}^2{u}^2(R)$.

Nearly matched sensitivities can provide strong suppression. With $a_X=1$, $a_Y=0.98$ and the same 20 nm reference uncertainty, the residual reference contribution is 0.4 nm. That does not imply the total uncertainty is 0.4 nm: comparison scatter, coefficient uncertainty, temperature and wringing contributions remain.

If the transfer ratios are themselves uncertain, their derivatives introduce terms proportional to the reference length, along with any shared covariance. Determine them from the actual calibration model rather than assigning ideal coefficients solely to obtain cancellation. NIST separates reference-master history from comparator-process variability in its uncertainty calculation.^[ref] A reference-sensitivity cancellation check should accompany a claimed low-uncertainty contrast.