Incoprea research · leads K · soft matter and fluid measurement · 3 November 2026

Two hundred leads on soft matter and fluid measurement

This page exists because of one sentence that appears, in some form, in most soft-matter papers and in almost none of their error budgets: the sample was measured in a single sweep. A frequency sweep, a stress ramp, a force–distance curve, a mobility, a sessile-drop angle. Each of them is a one-dimensional cut through a space that has more dimensions than the cut, and the model fitted to the cut is usually chosen from a shelf of three or four standard forms that the cut cannot separate.

The four sections below take four such cuts. In rheology the ambiguity is between a spectrum and a single stretched mode. At an interface it is the coordinate itself — the separation is known only up to an offset, and the offset is not measured. In a colloid it is a conversion factor between what the instrument reports (a mobility) and what the paper quotes (a potential). In wetting it is a convention: an angle that depends on the direction in which the contact line was last moved. Each section is fifty leads, and each lead is a sentence that could be shown false.

Part 1 · the leads

The four sections, and what the tags mean

test means the sentence could be settled by computation or by an examination of published data and standard models, without a new experiment. note means it could be settled only by a study of practice — a survey, a round robin, or an audit of published tables — which this programme has not run and, in most cases, could not run from a desk.

M. Rheology: what a single sweep cannot distinguish

Fifty leads on the gap between what a rheometer reports and what the constitutive statement behind it would require, and on the model families that a single frequency or rate sweep leaves tied.

N. Interfacial forces and the unknown zero of separation

Fifty leads on the quantity that every force measurement presupposes and none reports: the separation at which the two surfaces are actually zero apart.

O. Colloids: charge, aggregation and the meaning of a zeta potential

Fifty leads on the conversion that turns a measured mobility into a quoted potential, and on what that potential does and does not predict about aggregation.

P. Wetting and capillarity: conventions, hysteresis and contact-line pinning

Fifty leads on an angle that is quoted to a tenth of a degree and defined to no better than ten, and on the pinning that makes the definition unstable in the first place.

Part 2 · tests run

Tests run

Four computations, all synthetic. Each has a known generating truth, so the failure modes can be measured rather than argued about. Code: incoprea_nov3/code/softmatter_leads_compute.py. Output: incoprea_nov3/results/softmatter_leads.json. Seed 20261103, numpy 2.5.4. Re-run with python3 softmatter_leads_compute.py.

T1. A relaxation curve fitted by two standard models

A stress-relaxation curve generated as a two-mode Maxwell spectrum, G(t)/G₀ = ½e^−t/τ₁ + ½e^−t/τ₂, with τ₁ = 1 s and τ₂ = Rτ₁, sampled at 61 points over four decades with 1% proportional noise, 200 replicates. Two standard fits: a Kohlrausch–Williams–Watts stretched exponential (three parameters) and the correct biexponential (four). The question is not whether the models differ but what each one reports.

The models are separable much earlier than the folklore says. At R = 3 the biexponential wins by ΔAIC = 35.6, at R = 10 by 123.9, and at R ≥ 3 no replicate in 200 prefers the stretched form. Degeneracy survives only at R = 1.5, where 94% of replicates put the two within 2 AIC units of each other. Four decades and 1% noise is a great deal of information about a relaxation, and a single sweep of that span is not blind.

What the sweep cannot do is report a time. At R = 10 the fitted stretched exponential gives β = 0.593 and τ = 3.43τ₁, and the true modes are at 1 and 10. The reported time is the geometric mean of the two modes, not either of them, and it is the only time the model has. An author who quotes ‘a relaxation time of 3.4 s’ has quoted a parameter of the fit, not a property of the material, and the fit that produced it was rejected by its own AIC.

The exponent depends on the window, which is the part that matters. Holding R = 10 and changing only how many decades were recorded:

β moves from 0.493 at two decades to 0.629 at six, a drift of 0.136 in a parameter that is routinely reported to two decimal places as a material constant. Neither fit is wrong; the material is the same in both columns. The number describes the measurement window.

The same failure in creep, with no model competition at all

A Burgers creep compliance J(t) = J₀ + J₁(1 − e^−t/τ) + t/η with J₀ = 1.0×10⁻⁴ Pa⁻¹, J₁ = 4.0×10⁻⁴, τ = 2 s and η = 5×10⁵ Pa·s, fitted by a straight line in log–log coordinates over successive one-decade windows. No competing model, no noise — one curve, eight windows:

The exponent runs from 0.0656 to 0.4450, a factor of 6.8, and every window returns R² ≥ 0.93 — the fit looks good everywhere and means something different every time. The control is the same procedure applied to a genuine power law (J = 3×10⁻⁴ t^0.25): it returns 0.2500 in every window to four decimals, so the estimator is not the problem and neither is the arithmetic. The failure is that a power law is being asked to describe a curve with a retardation time in it.

What a single sweep cannot distinguish, in one line. It separates the models (AIC does that by R = 3) but not the parameters: β drifts by 0.136 with the window, the KWW time is the geometric mean of two modes, and a power law fitted to a creep curve returns any exponent between 0.07 and 0.44 depending on which decade was fitted.

T2. The unknown zero of separation, and what actually breaks the degeneracy

An oscillatory solvation force measured against piezo displacement z, where the true separation is D = z − δ and δ is unknown. The model fitted is F(D) = A e^−D/λ cos(2πD/p + φ) − C/D², with A = 0.80 nN, λ = 1.60 nm, p = 0.85 nm and C = 0.50 nN nm², 120 points from 0.6 to 9.0 nm and 1% noise. The exponential term is exactly invariant under D → D + δ: for any assumed zero d the parameter set (A e^−d/λ, λ, p, φ + 2πd/p) reproduces the curve to machine precision. The 1/D² term is the only thing in the model that is not.

The true offset in this run is 0.20 nm, so the phase at an assumed zero of d = 0 is already 275° and the absolute value of the phase means nothing on its own. The SSR column is the same number to every printed digit for every assumed zero, and the relative change is at the level of 10⁻¹⁴ — floating-point noise. The phase moves by −2πδ/p exactly. A purely exponential force law cannot identify the zero at any noise level, with any number of points, because the model has an exact symmetry and the data cannot break it. This is not a precision limit; it is a rank deficiency, and it is why every one-exponential fit in the literature has a zero that is a convention rather than a measurement.

Adding the van der Waals term makes the zero measurable — conditionally

Keeping the 1/D² term and profiling the fit over the assumed zero, the 1σ uncertainty on the zero is obtained from the curvature of the profile, as a function of how far into the near field the data extend, 30 replicates per row:

With data from 0.6 nm the zero is recovered to 0.0107 nm — a real measurement. With data starting at 3 nm it is recovered to 0.0786 nm, and since the layer spacing is 0.85 nm that is a tenth of a layer. The identifiability comes entirely from the term that decays as a power law, and it degrades by a factor of seven as the fit range is pulled back from the surface.

The cost of getting it wrong, when the van der Waals term has been subtracted first. Subtracting a fitted or literature 1/D² background before fitting the oscillation is standard practice, and it removes exactly the term that carried the information. Generating data with a true zero offset of 0.15 nm and fitting with the zero forced to 0:

Every parameter is recovered to a few per cent, the period is right to 1.3%, and the phase is wrong by 63.5° — which is the whole of the offset, 360° × 0.15/0.85. The fit is not bad; it is the identical fit, translated. If the analyst keeps the 1/D² term, the wrong zero is rejected at Δχ² = 81.8 and cannot be reported. If the analyst subtracts that term — as almost everyone does, because it is not the force under study — the wrong zero is not rejected at all: Δχ² = 0.0, and every parameter except the phase is unchanged.

The result. The unknown zero is not an uncertainty that better instruments will reduce. It is an exact symmetry of any force law that is a function of exponentials in the separation, and it is broken only by the term experimenters routinely subtract — which is why a 0.15 nm error costs Δχ² = 0 with the background removed and 81.8 with it kept.

T3. Contact-angle hysteresis from a random pinning landscape

A one-dimensional contact line on a heterogeneous surface, with the energy E(x)/γ = −cosθ_Y x + u(x), where u is a Gaussian random field of rms amplitude a and correlation length ξ. The line is pinned where u′(x) = cosθ, so it advances when cosθ exceeds the largest u′ along its path and recedes at the smallest. With sd(u′) = a/ξ and the extremes of a correlated field growing as √(2 ln(L/ξ)) over a contact line of length L, the model predicts Δcosθ ≈ 2(a/ξ)√(2 ln(L/ξ)). θ_Y = 60°, 400 realisations per cell, ξ = 1 and L/ξ = 10, 100, 1000.

The last column is a check on the generator: it should be 1.0, and it is 1.00 at L/ξ = 1000 and 1.21–1.24 at L/ξ = 10, where the finite sample has not yet converged. Hysteresis is linear in a/ξ across the whole range and grows only as the square root of the logarithm of the sampled length — a tenfold longer contact line raises a 15° hysteresis to 18.5°, not to 150°.

The predicted scaling against the simulation, over four amplitudes and three lengths:

The largest deviation between the simulation and the closed form is 9.7%, so the scaling law is adequate and the exponent is not in doubt: hysteresis measures a/ξ, the ratio of pinning strength to pinning length, and it measures nothing else.

Why one hysteresis cannot be inverted

Five surfaces with the same ratio a/ξ = 0.04 and the same physical contact-line length L = 100 µm, differing only in the correlation length — that is, in whether the pinning is strong and coarse or weak and fine:

A sixteenfold change in correlation length moves the hysteresis from 12.8° to 17.2° — 34% — while the amplitude changes by the same factor of sixteen to compensate. A measured hysteresis therefore constrains the ratio and not either quantity, and the two limits correspond to physically different surfaces: sparse strong pins against dense weak ones.

The convention, which moves the number more than the physics does

The same Δcos reported as an angular hysteresis, at five base angles:

A surface with Δcos = 0.1 reports a hysteresis of 5.7° at a base angle of 90° and 11.7° at 30°: the same surface, expressed in the units most papers use, differs by a factor of two. The cosine is the quantity the energy balance produces; the angle is the quantity the instrument displays.

The same non-linearity separates the two conventions for a patchy surface. For a mixture of 40° and 80° patches at area fraction f, the Cassie cosine average against the arithmetic angle average:

The two conventions differ by up to 2° at the mixture that maximises the difference, which is small — smaller than the hysteresis the same heterogeneity produces. That is the point: the convention debate about averaging is a 2° effect, while the pinning the same defects produce is the 15° effect in the first table, and the literature spends more words on the first.

T4. What a zeta potential is a measurement of

A fixed measured mobility of 3.0×10⁻⁸ m² V⁻¹ s⁻¹ in water at 298 K, converted to a potential by the Smoluchowski formula and by Henry's function (Ohshima's approximation, f = 1 + ½(1 + 2.5/κa)⁻³, which runs from 1 in the Hückel limit to 1.5 in the Smoluchowski limit):

The Smoluchowski value is the same number in every row, because it does not know how big the particle is. It is wrong by -33.2% at κa = 0.5 and by -16.3% at κa = 10, and the sign of the error is that Smoluchowski understates the potential of a small particle. For a 20 nm particle in 1 mM salt, κa = 2.1 and the error is 30%.

Why the 30% matters: DLVO amplifies it exponentially

Feeding each potential into a DLVO calculation for two 20 nm spheres, A = 5×10⁻²¹ J, 1 mM 1:1 electrolyte, the barrier and the stability ratio W = (2κa)⁻¹ exp(V_max/kT):

The same experiment, analysed two ways, predicts a suspension that is stable for months (10.6 decades of stability ratio) or one that aggregates in minutes (4.7 decades): a factor of 1.43 in the reported potential becomes a factor of 10^5.9 in the predicted stability. Nothing in the measurement changed — only the conversion.

The same amplification applies to the correction nobody makes, surface conductance, parameterised here as an apparent zeta reduced by 1/(1 + Du):

A Dukhin number of 0.5 — ordinary for an oxide at millimolar salt — removes 6.4 decades of predicted stability without moving the measured mobility by a per cent.

And the potential is not the charge

Converting ζ to surface charge through Gouy–Chapman at 1 mM, against the linear Debye–Hückel estimate σ = εε₀κζ:

The ratio is 1.1 at 40 mV and 2.2 at 120 mV: above roughly 60 mV a potential is a compressed report of a charge that is growing faster than it is. Two particles with potentials of 60 and 120 mV differ by a factor of two in potential and by a factor of 3.53 in charge, which is the quantity that enters every adsorption and coagulation model.

The result. A zeta potential is a mobility divided by a model. The model choice moves it by 30% at κa = 2, surface conductance moves it by a factor of two, and DLVO converts either error into 5.9 decades of predicted stability — so the number is worth quoting only with the four items that produced it.

Part 3 · limits

Limits

• Everything in the tests is a simulation. The generating truth is known because I wrote it. That makes the failure modes measurable and says nothing about how often they occur in real instruments. No measured data enters this page.
• The parameter values are chosen, not measured. λ = 1.6 nm, p = 0.85 nm, a/ξ = 0.04 and ζ = 40 mV are ordinary magnitudes, but the specific numbers in every table are consequences of those choices. A different mixture gives different numbers with the same signs.
• The Henry function is an approximation. Ohshima's closed form replaces a numerical solution of the electrokinetic equations, and it ignores surface conduction and polarisation. The 33% and 16% figures are properties of that approximation to within a few per cent, not exact results.
• The pinning model is one-dimensional and quenched. A real contact line can bend, which lowers the effective pinning, so the hysteresis computed here is an upper bound for a given a/ξ. The scaling exponent is the robust part; the prefactor is not.
• Noise is Gaussian and proportional. Rheometer noise near the torque floor is not proportional, and that changes where the tail of a relaxation stops carrying information. The T1 conclusions are for a well-behaved 1% relative noise.
• One seed. The replicate counts (200, 400, 30) are large enough for the medians to be stable to the last printed digit, and the single-run numbers carry sampling error of a few per cent. Seed 20261103 is recorded and re-running reproduces every number.
• The 200 leads are untested by construction. Several are probably known in the rheology, surface-force or wetting literature under other names; no literature search was run. Their value is as a falsifiable checklist, not as findings, and the tags record what kind of study would settle each one.
• No claim that any named instrument, paper or material is wrong. The failures discussed here are failures of reporting and of model choice, which this programme can exhibit without accusing anyone of a measurement error.

If one sentence survives: in each of these four measurements the number reported is a model times a reading, and the model — a stretched exponent, a zero of separation, a conversion factor, an averaging convention — moves the answer by more than the effect being claimed, while the fit quality moves by nothing at all.

Incoprea Research · leads K · 3 November 2026. Two hundred leads, four simulations, no measured data. Code and JSON in incoprea_nov3/code/ and incoprea_nov3/results/. No novelty claimed; no systematic literature search.