REVIEW 4 major objections 6 minor 17 references
Surface-Encoded Partial Coherence Transformation: Modeling Source Coherence Effects in Wave Optics
T0 review · 4 major / 6 minor · reviewed 2026-08-07 · deepseek-v4-flash
Pith's one-line read This paper argues that a single surface operator can imprint partial coherence on a wavefront, after which standard Fresnel propagation reproduces Van Cittert–Zernike theory.
desk verdict A clean restatement of textbook coherence results wrapped around a deterministic-coherence claim that fails on its own equations; not a new result and not ready for referees. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The central object is the surface coherence operator $C_S$, a linear integral operator with kernel $K_S(r,r',\lambda)$, composed with the Fresnel propagation operator $P_z$ as $U_d=P_z(C_S(R(U_i)))$. The identity doing the work is the claimed surface correlation $\Gamma_S(r_1,r_2)=|B_0|^2K_S(r_1-r_2,\lambda)$, which converts the kernel into a coherence function, and the kernel choice $K_S=\mathcal{F}\{I_s(\rho/\lambda z)\}$ that links the construction to the Van Cittert–Zernike theorem. The proof machinery is the coherence-function formalism, the Fresnel propagator, and the convolution theorem applied to the double propagation integral.
What would settle it
For a constant input $U_i=A_0$ and any spatially stationary kernel $K_S(r-r')$, the surface operator output is $U_s(r)=B_0\int K_S(r-r')\,d^2r'$, which is constant in $r$; the resulting correlation $\Gamma_S(r_1,r_2)=|B_0|^2\,|\int K_S|^2$ does not depend on $r_1-r_2$, whereas Theorem 2.4 claims it equals $|B_0|^2K_S(r_1-r_2,\lambda)$. Computing this two-point correlation for a constant input settles the central claim.
Extended reading notes
Core claim
The paper's discovery claim is that the mutual coherence function of a partially coherent field can be built in two deterministic steps rather than tracked as a four-dimensional object through the entire system. At the surface, the coherence operator $C_S$ with kernel $K_S(r,r',\lambda)$ acts on the reflected field; for a uniform plane-wave input the claimed result is $\Gamma_S(r_1,r_2)=|B_0|^2 K_S(r_1-r_2,\lambda)$, so the kernel itself becomes the mutual coherence function. Propagation to the detector is carried by the Fresnel operator $P_z$, and the full composition $U_d=P_z(C_S(R(U_i)))$ is claimed to be equivalent to the Van Cittert–Zernike theorem when $K_S$ is the Fourier transform of the source intensity. The paper derives explicit Airy and Gaussian kernels for circular and Gaussian sources, extends the framework to curved surfaces and polychromatic light, and claims the multi-surface computation scales as $O(M N^2\log N)$ instead of $O(N^{2M})$.
Load-bearing premise
The load-bearing premise is that a single deterministic linear operation at the surface can produce the same second-order correlations as a genuinely random partially coherent field; with a constant coherent input, the output field has no fluctuation, so the required correlation statistics are effectively assumed rather than derived.
Editorial extensions
If this is right
- For systems with multiple surfaces, the framework claims the cost of a coherence-aware simulation drops from $O(N^{2M})$ to $O(M N^2\log N)$, because each surface–propagation pair is applied sequentially and stationary kernels can be implemented with FFTs.
- The circular and Gaussian source corollaries give direct recipes for modeling stellar and Gaussian-Schell sources: choose the Airy or Gaussian kernel at the surface, then use ordinary Fresnel propagation.
- The separation of surface and propagation effects yields a dimensionless parameter $\eta$ that identifies whether surface roughness or propagation distance dominates the coherence change, guiding where to spend computational effort.
- The combined operator $K_{\text{tot}}$ acts as a spatial filter, so partial coherence effects can be computed and interpreted in the spatial-frequency domain.
- Polychromatic and temporal coherence enter by spectral weighting and a tensor-product spatial/temporal operator, so the framework claims to cover broadband sources without changing its structure.
Reading between the lines
- Read constructively, the equivalence proof suggests that $C_S$ is best viewed as a definition of the target coherence function at the surface; a follow-up development would be to formalize $C_S$ as a map on the mutual coherence function itself rather than on field amplitudes.
- The convergence theorem is stated as a limit with no rate, so a practical extension is to bound the finite-$N_c$ error and use it to set sampling requirements in simulations.
- For segmented-mirror telescopes, the framework suggests a pipeline in which source coherence is imprinted once at the pupil and segment phase errors are applied deterministically; this could be tested against existing stellar interferometry measurements.
- Replacing the Fresnel kernel with the exact Euclidean-distance kernel would extend the framework to non-paraxial and curved-surface systems, at the cost of losing the FFT speedup; the paper identifies this trade-off but does not quantify it.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript proposes a 'Surface-Encoded Coherence Transformation' (SECT) framework for modeling partial coherence in wave optics. The framework consists of a surface operator CS, defined by a kernel KS, and a free-space propagation operator Pz; the central claim is that applying CS to a fully coherent incident field deterministically produces a field with the second-order statistics of a partially coherent field, and that the composition Pz∘CS is mathematically equivalent to the Van Cittert–Zernike (VCZ) theorem when KS is chosen as the Fourier transform of the source intensity. The authors further claim an O(M N^2 log N) complexity for M-surface systems versus O(N^{2M}) for conventional mutual-coherence propagation, and provide convergence, energy-conservation, and operator-property results. The paper contains no numerical experiments, no code, and no comparison with existing partial-coherence simulations.
Significance. If the central claim were correct, the paper would offer a genuinely new viewpoint: encoding partial coherence through a deterministic surface operator rather than through ensembles, together with an exponential complexity reduction for multi-surface systems. The manuscript also correctly reviews standard coherence formalism and includes a correct positive-semidefiniteness argument for Fourier-transform kernels (Lemma A.1). However, the load-bearing assertion—that a deterministic linear operator can turn a coherent field into a field with genuine partially coherent statistics—is not established and is in fact contradicted by the manuscript's own equations. The claimed equivalence to VCZ is circular because the appendix assumes the desired correlation function rather than deriving it. Consequently, the advertised conceptual and computational advantages are not realized by the presented framework. The paper would need a fundamentally different stochastic or ensemble-based construction to support the stated claims, which is outside the current manuscript's scope.
major comments (4)
- [Theorem 2.4, Eqs. (29)–(32)] The reduction from the double integral in Eq. (31) to Eq. (32) is unjustified. For a deterministic linear operator with stationary kernel KS applied to the constant field B0, one obtains ΓS(r1,r2)=|B0|^2 [∫ KS(r1−r') d^2r']* [∫ KS(r2−r'') d^2r''], which is a constant in (r1,r2) for any integrable stationary kernel, not |B0|^2 KS(r1−r2). No 'appropriate normalization' can convert a constant into a function of r1−r2 unless KS itself is a delta function, in which case the output field is still B0 and has no fluctuations. Thus Theorem 2.4 does not establish that the SECTS component produces a partially coherent field.
- [Appendix A, Eq. (94)–(95)] The proof of SECT-VCZ equivalence assumes precisely the statistics it is supposed to derive. The field at the surface is written as Us(r)=B0+ΔUs(r) with ⟨ΔUs⟩=0 and ⟨ΔUs*(r1)ΔUs(r2)⟩=|B0|^2[γS(r1,r2)−1], where γS is already taken to be the Fourier-transformed source intensity. This is the desired second-order correlation, inserted by hand. The subsequent Fresnel propagation only propagates this assumed correlation; it does not show that the deterministic operator CS produces such a correlation. Theorem A.2 is therefore circular.
- [Section 2.7, Eqs. (46)–(48)] The same unjustified factorization appears in the combined formulation. Equation (47) gives Γd(r1,r2)=|A0|^2∬ Ktot*(r1,r1')Ktot(r2,r2') d^2r1' d^2r2', which is a rank-one kernel in (r1,r2) for any deterministic Ktot. Equation (48) then claims this reduces to |A0|^2 Ktot(r1−r2). That reduction is valid only if the double integral collapses to a single kernel value, which in general it does not; the step silently assumes that the product of the two integrals factorizes into Ktot(r1−r2), which is false for any non-delta kernel. This invalidates the claimed unified operator equivalence.
- [Theorem 4.1, proof surrounding Eqs. (63)–(70)] The convergence theorem asserts that the SECT framework converges to the ensemble average, but the proof invokes randomness that the framework does not contain. The statement that 'each grid point effectively represents an independent random sample from the ensemble' is unsupported because CS is a deterministic convolution; no probability space or ensemble of realizations is defined. In the ρc→0 limit claimed to produce independent samples, the kernel KS approaches a delta function and the output becomes B0, with no fluctuations at all. The central limit theorem and law of large numbers are therefore inapplicable, and the claimed convergence is not established.
minor comments (6)
- [Eq. (6) and Appendix A] Several equations have inconsistent dimensional scalings: Eq. (6) writes KS(r−r')∝F{Is(ρ/λz)}(r−r'), but the Fourier transform of Is(ρ/λz) contains factors of (λz)^2 that are not tracked, and the same issue recurs in Eqs. (9)–(11) and (14)–(16). The notation should be made dimensionally consistent.
- [Section 3.2, Table 2] The η values listed for each scenario are asserted without any derivation, simulation, or measurement. The claims such as 'η≈0.1 (SECTP dominant)' are unsupported and should be either derived from a concrete model or removed.
- [Figure 4 and caption] The figure claims coherence length increases from approximately 1–2 units at the surface to 2.5–5.5 units at the detection plane, but no simulation parameters, numerical method, or error bars are given. If this is an illustrative schematic, the caption should say so explicitly.
- [Proposition 4.2] The energy conservation claim ∫|Ctot(U)|^2 = ∫|U|^2 for arbitrary convolution kernels is not generally true; a Gaussian low-pass kernel, for example, reduces total energy. No proof is provided, and the proposition should either be restricted to unitary (phase-only) kernels or corrected.
- [Section 5.3] The complexity claim O(M N^2 log N) assumes stationary kernels and FFT-based convolution for both CS and Pz. For curved surfaces and non-planar detectors, the manuscript itself notes that FFT shortcuts are unavailable (Section 5.1.3), so the stated complexity reduction does not apply to the generalized framework.
- [Throughout] There are numerous typos and unfinished expressions, e.g., 'e ffects', 'di fferent', and the dangling 'R' in 'and Appendix A begins by assuming' do not appear in the manuscript, but the text contains repeated spacing errors and inconsistent notation such as σθ used both as an angular width and a radial coordinate in Corollary 2.3. A thorough copyedit is needed.
Circularity Check
The SECT-VCZ 'equivalence' is circular by construction: Eq. (6) defines the surface kernel as the VCZ coherence function, and Appendix A inserts the desired surface correlation by hand; the deterministic operator never produces the claimed ensemble statistics.
-
self definitional
[Section 2.2.1, Theorem 2.1, Eq. (6)]
"Theorem 2.1 (SECT-VCZ Equivalence). Given a fully coherent incident field Ui(r) = A0 (uniform plane wave illumination) impinging on a surface with coherence operator CS characterized by kernel KS(r, r′, λ), followed by propagation operator Pz, the mutual coherence function at the detection plane is mathematically equivalent to that obtained from the VCZ theorem for an incoherent source with intensity distribution Is(ρ), provided that: KS(r, r′, λ) ∝ F{Is(ρ/λz)}(r−r′). (6)"
The VCZ theorem states exactly that the spatial coherence at distance z from an incoherent source is the Fourier transform of the source intensity Is at scaled separation Δr/λz. Equation (6) therefore defines the surface kernel to be the VCZ mutual coherence function. The theorem's hypothesis contains its conclusion; the corollaries merely Fourier-transform chosen source profiles and label the result an equivalence. Any propagation algebra that follows re-derives what was already placed into KS.
-
self definitional
[Appendix A, Eqs. (93)-(95)]
"We represent this by considering the field at the surface as: Us(r) = B0 + ΔUs(r) where ⟨ΔUs(r)⟩ = 0 and ⟨ΔU∗s(r1)ΔUs(r2)⟩ = |B0|2[γS(r1, r2)−1]. Here, γS(r1, r2) is the normalized mutual coherence function at the surface, which for a spatially stationary kernel depends only on the separation r1−r2. For the specific kernel form stated in the theorem: γS(r1, r2) = F{Is(ρ/λz)}(r1−r2). (95)"
This is the result to be proved, not a derived property. With Ui = A0 and deterministic R, the output of the linear operator CS is a deterministic field, so the ensemble correlation of ΔUs cannot acquire VCZ statistics unless those statistics are postulated. The rest of Appendix A propagates this pre-assigned correlation through the Fresnel kernel and rediscovers VCZ. The equivalence is therefore assumed by construction.
2 more flagged steps
-
self definitional
[Section 2.5, Theorem 2.4, Eqs. (30)-(32)]
"For a fully coherent incident field, ⟨U∗i(r′1)Ui(r′2)⟩ = U∗i(r′1)Ui(r′2). Assuming a uniform plane wave illumination Ui(r) = A0 and a uniform reflection coefficient R(r) = R0 for simplicity: ΓS(r1,r2)=|R0|2|A0|2∫∫K∗S(r1,r′1,λ)KS(r2,r′2,λ) ei[ϕs(r′2)−ϕs(r′1)]d2r′1d2r′2. For a flat surface, where ϕs(r)=ϕ0 is constant, and assuming a spatially stationary kernel KS(r,r′;λ)=KS(r−r′;λ) with appropriate normalization, this reduces to: ΓS(r1,r2)=|R0|2|A0|2KS(r1−r2;λ)."
For deterministic Ui = A0, the field after CS is Us(r) = B0∫KS(r−r′)d2r′, so Eq. (31)'s double integral factorizes into a product of constants for a stationary kernel. No 'appropriate normalization' can turn that rank-one constant product into the non-constant coherence kernel KS(r1−r2) appearing in Eq. (32). The equality is imposed by fiat, making the claimed transformation circular rather than derived.
-
other
[Section 4.1, proof of Theorem 4.1]
"In the SECT approach, as Nc→∞ (and ρc→0), each grid point effectively represents an independent random sample from the ensemble of possible field values. The propagation operator Pz then combines these independent samples with appropriate weights."
The SECT implementation uses deterministic linear operators CS and Pz applied to a deterministic plane wave; there are no random variables and no ensemble from which to sample. The assertion that grid points are 'independent random samples' assumes the very statistical ensemble that the deterministic framework claims to create, so the convergence to an ensemble average is assumed rather than derived.
full rationale
The central derivation chain is circular in two linked places. First, Theorem 2.1 defines the surface kernel KS as the Fourier transform of the source intensity, which is precisely the VCZ spatial coherence function; the theorem's 'provided that' condition already contains the conclusion. Second, Appendix A decomposes the deterministic surface field as Us = B0 + ΔUs and assigns ⟨ΔUs*(r1)ΔUs(r2)⟩ = |B0|²[γS(r1,r2)−1] with γS = F{Is}, i.e., it inserts by hand the very correlation the paper claims to derive. Theorem 2.4 repeats the same move under the phrase 'with appropriate normalization': for a constant deterministic input the double integral factorizes to a constant and cannot produce a non-constant coherence kernel. Theorem 4.1 similarly declares deterministic grid points to be independent random samples, again assuming the ensemble statistics. The paper does contain standard, independent mathematical material—Fresnel propagation of coherence functions, positive-semidefiniteness of the Fourier transform of an intensity, and FFT complexity estimates—but these do not rescue the load-bearing claim that a deterministic linear operator encodes partial coherence or that the framework is equivalent to VCZ. The self-citation [16] on figure/roughness separation is not load-bearing for this circularity. Because the main 'prediction' is assumed by construction, the circularity score is 9.
Assumptions & free parameters
free parameters (2)
- Scaling exponents for rho_s/rho_i and rho_d/rho_s
- eta values per scenario in Table 2 =
0.1, 0.9, 0.4-0.6, 0.8, 0.2, 0.7
assumptions (6)
- domain assumption Paraxial approximation: max(L_surface, L_detector) << z
- domain assumption Statistical stationarity and separability of spatial and temporal coherence: Gamma(r1,r2,tau)=gamma(r1,r2) phi(tau)
- ad hoc to paper The surface kernel KS is defined as the Fourier transform of the source intensity (Eq. 6)
- ad hoc to paper A deterministic operator CS can encode partial coherence statistics
- ad hoc to paper The surface field can be decomposed as Us = B0 + Delta Us with <Delta Us>=0 and <Delta Us*(r1) Delta Us(r2)> = |B0|^2 [gamma_S(r1,r2)-1]
- ad hoc to paper Spatial averaging over coherence cells converges to ensemble averaging
Cite this review
Pith. "Pith review of Surface-Encoded Partial Coherence Transformation: Modeling Source Coherence Effects in Wave Optics." pith.science (2026). https://pith.science/paper/VNPMR2SJ
@misc{pith2026250517754,
author = {Pith},
title = {Pith review of: Surface-Encoded Partial Coherence Transformation: Modeling Source Coherence Effects in Wave Optics},
year = {2026},
howpublished = {\url{https://pith.science/paper/VNPMR2SJ}},
note = {Machine review of arXiv:2505.17754}
}
read the original abstract
We present a new mathematical framework for incorporating partial coherence effects into wave optics simulations through a comprehensive surface-to-detector approach. Unlike traditional ensemble averaging methods, our dual-component framework models partial coherence through: (1) a surface-encoded transformation implemented via a linear integral operator with a spatially-dependent kernel that modifies coherence properties at the reflection interface, followed by (2) a propagation component that evolves these coherence properties to the detection plane. This approach differs fundamentally from conventional models by explicitly separating surface interactions from propagation effects, while maintaining a unified mathematical structure. We derive the mathematical foundation based on the coherence function formalism, establish the connection to the Van Cittert-Zernike theorem, and prove the equivalence of our framework to conventional partial coherence theory. The method reduces the dimensional complexity of coherence calculations and offers potential computational advantages, particularly for systems involving multiple surfaces and propagation steps. Applications include optical testing and astronomical instrumentation. We provide rigorous mathematical proofs, demonstrate the convergence properties, and analyze the relative importance of surface and propagation effects across different optical scenarios.
Figures
Figures from the paper (1 more)
Reference graph
Works this paper leans on
-
[1]
P. H. van Cittert, “Die wahrscheinliche Schwingungsverteilung in einer von einer Lichtquelle direkt oder mittels einer Linse beleuchteten Ebene,” Physica, vol. 1, pp. 201–210, 1934. https://doi.org/10.1016/S0031-8914(34)90026-1
-
[2]
The concept of degree of coherence and its application to optical problems,
F. Zernike, “The concept of degree of coherence and its application to optical problems,” Physica, vol. 5, pp. 785–795, 1938. https://doi.org/10.1016/ S0031-8914(38)80075-4
work page 1938
-
[3]
L. Mandel and E. Wolf, Optical Coherence and Quantum Optics , Cambridge Uni- versity Press, 1995. doi:10.1017/CBO9781139644105
-
[4]
Measurement of the diameter of Alpha Orionis with the interferometer,
A. A. Michelson and F. G. Pease, “Measurement of the diameter of Alpha Orionis with the interferometer,”Astrophysical Journal, vol. 53, pp. 249–259, 1921.https: //doi.org/10.1086/142620
-
[5]
A test of a new type of stellar interferometer on Sirius,
R. Hanbury Brown and R. Q. Twiss, “A test of a new type of stellar interferometer on Sirius,” Nature, vol. 178, pp. 1046–1048, 1956. https://doi.org/10.1038/ 1781046a0
work page 1956
-
[6]
New Theory of Partial Coherence in the Space-Frequency Domain,
E. Wolf, “New Theory of Partial Coherence in the Space-Frequency Domain,”JOSA, vol. 72, no. 3, pp. 343–351, 1982. doi:10.1364/JOSA.72.000343
-
[7]
J. W. Goodman, Statistical Optics , 2nd ed., Wiley, 2015. doi:10.1002/9781119009719
-
[8]
Wigner functions in optics: describing beams as ray bundles and pulses as particle ensembles,
M. A. Alonso, “Wigner functions in optics: describing beams as ray bundles and pulses as particle ensembles,” Advances in Optics and Photonics , vol. 3, no. 4, pp. 272–365, 2011
work page 2011
Show all 17 references
-
[9]
Spatial correlation properties of focused partially coherent light,
D. G. Fischer and T. D. Visser, “Spatial correlation properties of focused partially coherent light,” Journal of the Optical Society of America A , vol. 21, no. 11, pp. 2097–2102, 2004
2004
-
[10]
Coherence and the spatial distribution of intensity,
F. Gori, M. Santarsiero, R. Borghi, and G. Guattari, “Coherence and the spatial distribution of intensity,”Journal of the Optical Society of America A, vol. 20, no. 10, pp. 2107–2119, 2003
2003
-
[11]
Propagation parameters of Gaussian Schell- model beams,
A. T. Friberg and R. J. Sudol, “Propagation parameters of Gaussian Schell- model beams,” Optics Communications , vol. 41, no. 6, pp. 383–387, 1982. doi:10.1016/0030-4018(82)90161-4
1982 doi
-
[12]
Interferometry and aperture synthesis in astronomy,
D. F. Buscher, “Interferometry and aperture synthesis in astronomy,” Reports on Progress in Physics , vol. 65, pp. 267–311, 2002. https://doi.org/10.1088/ 0034-4885/65/2/202 49
2002
-
[13]
Born, M., & Wolf, E. (1999). Principles of optics: electromagnetic theory of propa- gation, interference and diffraction of light. Cambridge University Press
1999
-
[14]
Goodman, J. W. (2005). Introduction to Fourier optics. Roberts and Company Pub- lishers
2005
-
[15]
Electromagnetic di ffraction in optical systems. II. Struc- ture of the image field in an aplanatic system,
B. Richards and E. Wolf, “Electromagnetic di ffraction in optical systems. II. Struc- ture of the image field in an aplanatic system,” Proceedings of the Royal Society of London. Series A. Mathematical and Physical Sciences 253, 358–379 (1959)
1959
-
[16]
Bridging Statistical Scattering and Aberration Theory: Ray-Deflection Function I: Theoretical Framework,
N. Moriya, “Bridging Statistical Scattering and Aberration Theory: Ray-Deflection Function I: Theoretical Framework,” arXiv:2505.01019 [physics.optics] (2025)
2025 arXiv
-
[17]
Optical interferometry in astronomy,
J. D. Monnier, “Optical interferometry in astronomy,” Reports on Progress in Physics, vol. 66, no. 5, pp. 789–857, 2003, doi: 10.1088/0034-4885/66/5/203. 50
2003 doi
Reviewed August 7, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.