Valley-Peak Modulation in Phase Space: an Exposure-Invariant VPM and its Theta-Function Structure

Aaron J. Hendrickson; David P. Haefner

arxiv: 2603.01199 · v2 · submitted 2026-03-01 · ⚛️ physics.ins-det

Valley-Peak Modulation in Phase Space: an Exposure-Invariant VPM and its Theta-Function Structure

Aaron J. Hendrickson , David P. Haefner This is my paper

Pith reviewed 2026-05-15 18:01 UTC · model grok-4.3

classification ⚛️ physics.ins-det

keywords valley-peak modulationread noisewrapped GaussianJacobi theta functionsphase spaceelliptic integralsDSERNCMOS sensors

0 comments

The pith

A phase mapping to wrapped Gaussians makes valley-peak modulation an exposure-invariant measure of read noise, expressed as a normalized Jacobi theta ratio.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper shows how to remove the dependence on exposure from the valley-peak modulation (VPM) metric used to quantify read noise in deep sub-electron read noise CMOS sensors. By applying a phase mapping that wraps the output to quotient out the integer electron count, the sensor statistics become a wrapped Gaussian whose properties depend only on the read noise. This wrapped density admits a representation in terms of Jacobi theta functions, and the exposure-invariant core is the ratio R(σ) = ϑ₄(q)/ϑ₃(q). The measured VPM is then a contrast normalization of this ratio, and the read noise can be recovered in closed form using elliptic integrals. This matters because it allows accurate noise characterization without controlling the exact number of photons or electrons incident on the sensor.

Core claim

Starting from the standard Poisson-Gaussian model, a phase mapping is applied that quotients out the integer electron count, yielding a wrapped-Gaussian density parameterized only by read noise. This density has both lattice-sum and Jacobi theta-function representations. The fundamental exposure-invariant quantity is the theta ratio R(σ)=ϑ₄(q)/ϑ₃(q), of which any VPM is a contrast normalization. A closed-form inverse expressing read noise in terms of VPM is obtained using elliptic integrals.

What carries the argument

The phase-mapped wrapped-Gaussian density and its representation as a ratio of Jacobi theta functions ϑ₄/ϑ₃, which carries the exposure-invariant contrast information.

If this is right

Read noise can be estimated directly from observed VPM without knowledge of exposure level.
Existing exposure-independent approximations are recovered as low-order truncations of the lattice-sum expansion of the theta ratio.
The elliptic-integral inversion provides an exact method to convert measured modulation back to read noise sigma.
Practical estimation is illustrated via simulation in phase space for DSERN sensors.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

Precomputing the elliptic-integral inverse could enable fast, real-time read-noise monitoring in camera firmware without additional calibration exposures.
This phase-space approach may extend to other periodic or wrapped noise metrics in sensor arrays or imaging systems.
Connections to theta-function techniques in other domains, such as lattice models in physics, could yield further analytic tools for noise analysis.

Load-bearing premise

The Poisson-Gaussian model for sensor output remains accurate after the phase mapping, and the wrapping fully removes the integer electron count without biasing the noise metric.

What would settle it

If independent measurements of read noise from controlled low-exposure images deviate systematically from the values obtained by inverting the observed VPM using the elliptic-integral formula, the invariance or the model would be falsified.

read the original abstract

Valley-peak modulation (VPM) was introduced as a metric for quantifying read noise in deep sub-electron read noise (DSERN) CMOS sensors. In the original amplitude-domain definition, VPM depends on both read noise and quanta exposure, yet Starkey & Fossum demonstrated exposure-independent approximations that hold in the DSERN regime. In this note we identify the exposure-invariant object those approximations probe. Starting from the standard Poisson-Gaussian model, we apply a phase mapping that quotients out the integer electron count, yielding a wrapped-Gaussian density parameterized only by read noise and admitting both lattice-sum and Jacobi theta-function representations. The fundamental exposure-invariant quantity is shown to be the theta ratio $R(\sigma)=\vartheta_4(q)/\vartheta_3(q)$, of which any VPM is a contrast normalization; the existing exposure-independent approximations are then recovered as low-order truncations of the lattice-sum representation of $R$. A closed-form inverse expressing read noise in terms of VPM is obtained using elliptic integrals, and a short simulation example illustrates practical estimation of read noise from the VPM in phase space.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

This note pins the exposure-invariant core of VPM to the Jacobi theta ratio R(σ) and supplies an elliptic-integral inversion for read noise.

read the letter

The main takeaway is that the exposure-invariant part of valley-peak modulation is the theta-function ratio R(σ) = ϑ₄(q)/ϑ₃(q), with a closed-form elliptic-integral expression that recovers the read-noise width σ from a measured VPM value. The paper starts from the standard Poisson-Gaussian model, maps to the fractional part to produce a wrapped Gaussian whose density depends only on σ, and uses Poisson summation to reach the theta representation. The ratio R(σ) is independent of exposure μ, any contrast-normalized VPM is a monotonic function of it, and the earlier approximations by Starkey and Fossum appear as low-order truncations of the lattice sum. The elliptic-integral inversion is the practical addition for estimating σ directly from data. This identification and the explicit inverse are not in the prior VPM literature, so the reformulation is new. The derivation stays within classical results, the circularity is low because σ enters as the physical parameter rather than a fitted quantity, and the steps are traceable. The note is short and the simulation is only an illustration, so it lacks extensive validation across sample sizes or departures from the Poisson-Gaussian assumption. Those are minor gaps rather than load-bearing problems; the central argument holds without contradiction. The work is aimed at the narrow community that characterizes deep sub-electron read noise CMOS sensors. A reader who needs a theoretical handle on VPM or a direct formula for read noise will find it useful and reproducible. It deserves a serious referee because the math is grounded, the result is falsifiable with sensor measurements, and the contribution clarifies an existing practical metric without overreach.

Referee Report

1 major / 3 minor

Summary. The manuscript identifies the exposure-invariant core of valley-peak modulation (VPM) for read noise quantification in deep sub-electron read noise (DSERN) CMOS sensors. Starting from the Poisson-Gaussian model, a phase mapping to the fractional part yields a wrapped-Gaussian density depending solely on read noise width σ. This density admits lattice-sum and Jacobi theta-function representations; the ratio R(σ) = ϑ₄(q)/ϑ₃(q) is shown to be the fundamental invariant, with any VPM arising as its contrast normalization. Existing exposure-independent approximations are recovered as low-order truncations of the lattice sum, a closed-form inverse for σ in terms of observed VPM is derived via elliptic integrals, and the approach is illustrated with a simulation example.

Significance. If the mapping and identifications are verified, the work supplies a rigorous theoretical foundation linking empirical VPM metrics to classical special-function identities (Poisson summation and Jacobi theta relations). The parameter-free character of R(σ) (depending only on the physical σ) and the elliptic-integral inverse constitute clear strengths, enabling exposure-independent read-noise estimation without ad-hoc fitting. This could standardize VPM usage in DSERN sensor characterization and calibration. The simulation provides a practical check, though its limited scope leaves room for broader validation.

major comments (1)

[§3] §3 (phase mapping and wrapped density): the assertion that the wrapping operation {n + g} = {g} fully quotients out the integer Poisson count without introducing bias in the subsequent noise metric is central to the exposure-invariance claim, yet the manuscript provides no explicit error analysis or numerical check confirming that the VPM computed from the wrapped density exactly recovers the original amplitude-domain definition of Starkey & Fossum in the DSERN regime.

minor comments (3)

[Notation] The definition of the nome q appearing in ϑ₄(q) and ϑ₃(q) is not stated explicitly in the main text (only implied in the abstract); it should be given as q = exp(−2π²σ²) at first use for clarity.
[Simulation] The simulation example lacks reported exposure values, number of Monte-Carlo trials, and quantitative comparison metrics (e.g., recovered σ versus input σ); adding these would improve reproducibility.
[References] The reference to Starkey & Fossum should be expanded to include full bibliographic details (journal, volume, year, DOI) rather than appearing only by name.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for the careful reading and positive recommendation of minor revision. The single major comment is addressed below; we will incorporate the requested verification in the revised manuscript.

read point-by-point responses

Referee: [§3] §3 (phase mapping and wrapped density): the assertion that the wrapping operation {n + g} = {g} fully quotients out the integer Poisson count without introducing bias in the subsequent noise metric is central to the exposure-invariance claim, yet the manuscript provides no explicit error analysis or numerical check confirming that the VPM computed from the wrapped density exactly recovers the original amplitude-domain definition of Starkey & Fossum in the DSERN regime.

Authors: The equality {n + g} = {g} holds exactly because n is integer-valued; the fractional-part operator is invariant under addition of any integer. Consequently, the distribution of the observed fractional part is precisely the wrapped-Gaussian density parameterized solely by σ, independent of exposure λ. The VPM, being a functional of the valley and peak locations (or their contrast) in this density, therefore coincides exactly with the original amplitude-domain definition, introducing neither bias nor approximation. We will add a short derivation of this exact invariance together with a numerical comparison (wrapped vs. amplitude-domain VPM) in the DSERN regime to §3 of the revised manuscript. revision: yes

Circularity Check

0 steps flagged

No significant circularity; derivation is self-contained

full rationale

The paper starts from the standard Poisson-Gaussian model and applies a phase mapping {n + g} = {g} that isolates the wrapped-Gaussian density depending only on read noise σ. This density is then expressed via standard lattice sums and Jacobi theta functions, yielding the ratio R(σ) = ϑ₄(q)/ϑ₃(q) as the exposure-invariant object by direct mathematical construction. Any VPM appears only as a subsequent contrast normalization, and the elliptic-integral inverse recovers σ from observed VPM without presupposing the result. No steps reduce by definition to fitted VPM data, no self-citations are load-bearing for the central claim, and the approximations are recovered as truncations of the independent lattice representation. The derivation chain is therefore independent of its target metric.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The derivation rests on the standard Poisson-Gaussian noise model for CMOS sensors and the classical theory of wrapped distributions and Jacobi theta functions; no new free parameters or invented entities are introduced beyond the physical read-noise width σ.

free parameters (1)

read noise width σ
Physical parameter of the Gaussian component; the target quantity to be estimated, not fitted ad hoc to produce the invariant.

axioms (2)

domain assumption Sensor output follows the standard Poisson-Gaussian model
Invoked at the start of the phase-mapping step; standard in the DSERN literature but not re-derived here.
standard math Jacobi theta-function identities for the wrapped Gaussian
Used to obtain the closed lattice-sum and theta representations; classical results assumed without proof.

pith-pipeline@v0.9.0 · 5501 in / 1446 out tokens · 48336 ms · 2026-05-15T18:01:50.018260+00:00 · methodology

Valley-Peak Modulation in Phase Space: an Exposure-Invariant VPM and its Theta-Function Structure

Core claim

What carries the argument

If this is right

Where Pith is reading between the lines

Load-bearing premise

What would settle it

discussion (0)