Pith. sign in

REVIEW 4 major objections 4 minor 32 references

Statistics of Non-Rayleigh Speckles Generated from Nonlinear Media

T0 review · 4 major / 4 minor · reviewed 2026-08-07 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

Pith's one-line read This paper derives an analytical expression for the probability density of speckle intensity in a Kerr nonlinear medium, predicting that focusing nonlinearity stretches the high-intensity tail while defocusing compresses it.

desk verdict The numerical correlation observations are worth a look, but the analytic PDF rides on an unjustified intensity mapping that a plane-wave solution of the validating NLSE falsifies. read the letter →

arxiv 2506.04816 v3 pith:PJGMNSBH submitted 2025-06-05 physics.optics

Deependra Singh Gaur , Akanksha Gautam , Rakesh Kumar Singh , Akhilesh Kumar Mishra This is my paper
classification physics.optics
keywords speckleKerrnonlinearityintensityprobabilitydensityfunctionnon-RayleighstatisticscorrelationfieldnonlinearSchrödingerequationSiegertrelation
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

The paper aims to give a closed-form formula for the intensity probability density function (PDF) of a speckle field propagating through a cubic Kerr medium. It claims that focusing nonlinearity enhances the probability of bright speckles, producing a longer tail in the PDF, while defocusing nonlinearity suppresses bright speckles and shortens the tail. The authors validate the formula against numerical solutions of the nonlinear Schrödinger equation. If correct, the result supplies a simple analytic tool for predicting non-Rayleigh intensity statistics and their dependence on the sign of the nonlinearity in complex scattering media.

What carries the argument

The central machinery is the variable-transformation method for probability densities, applied to the pointwise intensity mapping $I_{NL}=I+n_2 I^2$. The inverse $I=(\sqrt{1+4n_2 I_{NL}}-1)/(2n_2)$ combined with the derivative $dI/dI_{NL}=1/\sqrt{1+4n_2 I_{NL}}$ converts the linear negative-exponential PDF into the closed form of Eq. (6). The paper additionally uses the nonlinear Schrödinger equation (NLSE) in normalized units to numerically generate and propagate speckles, supplying the comparison data that the analytic PDF is fitted against.

What would settle it

Measure the one-point intensity PDF of a fully developed speckle after it propagates through a Kerr medium with known $n_2$ and compare quantitatively with Eq. (6). The formula is falsified if the observed tail shape deviates beyond numerical error, and especially in the defocusing case, where the model predicts an imaginary exponential for $4|n_2|I_{NL}>1$ and thus cannot predict any real PDF at all.

Watch

Extended reading notes

Core claim

The paper's central claim is that, when a speckle field's intensity in a linear medium is exponentially distributed with PDF $\rho(I)=\frac{1}{2\sigma^2}e^{-I/(2\sigma^2)}$, the intensity in a Kerr medium follows an exact transformed distribution obtained from the algebraic mapping $I_{NL}=I+n_2 I^2$. Using the standard variable-transformation rule for monotonic functions, the authors arrive at Eq. (6): $\rho(I_{NL})=\frac{1}{2\sigma^2\sqrt{1+4n_2 I_{NL}}}\exp\left[\frac{1-\sqrt{1+4n_2 I_{NL}}}{4n_2\sigma^2}\right]$, choosing the branch with a physical inverse. This expression shows that focusing nonlinearity ($n_2>0$) extends the tail of the intensity PDF, whereas defocusing nonlinearity ($n_2<0$) shortens it, matching the numerical NLSE simulations. The paper also reports that the first-order field correlation remains unchanged while the second-order intensity correlation is modified by the nonlinearity, breaking the linear-media Siegert relation in a direction that depends on the sign of $n_2$.

Load-bearing premise

The entire derivation rests on the claim that the intensity in the nonlinear medium is given by the simple algebraic formula $I_{NL}=I+n_2 I^2$, which is stated without being derived from the wave equation; if that pointwise relation is not an accurate description of Kerr propagation, the predicted PDF does not describe the actual intensity statistics.

Editorial extensions

If this is right

  • Focusing Kerr nonlinearity yields super-Rayleigh speckle statistics with enhanced intensity contrast and a higher probability of bright spots, while defocusing yields sub-Rayleigh statistics with reduced contrast.
  • The Siegert relation for speckle intensity fluctuations is violated in nonlinear media: focusing makes the intensity correlation exceed the squared field correlation, and defocusing makes it fall below.
  • The complex field amplitude retains Gaussian statistics and its spatial correlation is unaffected by the nonlinearity, so all deviations from Rayleigh behavior appear only in intensity moments.
  • In the limit $n_2=0$, Eq. (6) reduces exactly to the familiar negative-exponential intensity PDF of fully developed Rayleigh speckle.
  • The derived PDF provides a predictive analytic tool for tailored speckle statistics in nonlinear scattering experiments, potentially useful for imaging and metrology that rely on controlled intensity fluctuations.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The assumed mapping $I_{NL}=I+n_2 I^2$ is a pointwise algebraic relation, not directly derived from the NLSE; the analytic PDF is therefore best interpreted as describing the local self-action of a thin nonlinear slice, while full propagation effects captured by the NLSE may introduce additional corrections beyond this simple transformation.
  • The same variable-transformation route could be applied to other monotonic intensity-dependent refractive-index models, such as saturable or higher-order Kerr nonlinearities, yielding closed-form PDFs that could be tested numerically with the same NLSE machinery.
  • The paper's admitted breakdown for defocusing media when $|4n_2 I_{NL}|>1$ suggests that for strong defocusing the quadratic mapping becomes non-invertible over the physical domain; a two-branch treatment or a different intensity definition would be needed to cover that regime.
  • The observed enhancement of the intensity tail for focusing nonlinearity is qualitatively similar to the statistics of rogue wave or filamentation events, hinting that the formula might capture early-stage intensity spike formation in disordered nonlinear media.
Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Referee Report

4 major / 4 minor

Summary. The paper claims to derive an analytical expression for the probability density function (PDF) of speckle intensity in a Kerr nonlinear medium. The derivation is based on an assumed pointwise algebraic relation between the linear intensity I and a so-called nonlinear intensity I_NL, Eq. (1): I_NL = I + n2 I^2. Using standard variable transformation, the authors obtain Eq. (6), a modified exponential PDF whose tail is enhanced for focusing nonlinearity (n2 > 0) and reduced for defocusing nonlinearity (n2 < 0). The paper also presents numerical simulations of the nonlinear Schrödinger equation (Eq. 7) and reports that the field correlation function is unchanged while the intensity correlation function is modified. The central claim is that Eq. (6) provides an analytic description of non-Rayleigh speckle statistics in nonlinear media.

Significance. If Eq. (6) were correct and properly derived, it would offer a compact analytic description of first-order statistics for speckle in Kerr media, which could be useful for imaging and light-matter interaction studies. However, the central assumption Eq. (1) is not derived from the underlying wave equation and is actually inconsistent with the very NLSE used for validation: for a plane wave, the NLSE predicts constant |E|^2, whereas Eq. (1) predicts an intensity increase. The paper also contains an internal contradiction between the claim of Gaussian field statistics and a non-exponential intensity PDF. Because the main theoretical result rests on an invalid premise, the paper's contribution as it stands is not a reliable prediction; the claimed numerical agreement is presented only as a 'best fit' without quantitative metrics. No machine-checked proofs, reproducible code, or parameter-free derivations are provided.

major comments (4)
  1. [Theory, Eqs. (1)-(4) and Eq. (7)] Eq. (1), I_NL = I + n2 I^2, is the load-bearing assumption of the entire derivation, but it is not derived from Maxwell's equations or from the NLSE used for validation. In the NLSE (Eq. 7), the physical intensity is always |E|^2; the nonlinearity modifies the field E through the term n2|E|^2 E, not the definition of intensity. A direct counterexample is the plane-wave solution E(z) = A exp(i n2 |A|^2 z) of Eq. (7) with zero transverse Laplacian, for which |E(z)|^2 = |A|^2 is unchanged, whereas Eq. (1) would predict |A|^2 + n2 |A|^4. Thus Eq. (6), which follows from Eq. (1), does not describe the statistics of the validating model.
  2. [Introduction and Eq. (6)] The paper states that 'the field distribution of the generated speckle pattern follows Gaussian statistics and remains invariant in all nonlinear conditions', while simultaneously claiming that the intensity PDF is the non-exponential Eq. (6). For a circular complex Gaussian field, the intensity PDF is necessarily the negative exponential distribution; a non-exponential intensity PDF implies non-Gaussian field statistics. This is an internal mathematical inconsistency, not merely a deviation from consensus.
  3. [Results and Discussion, Fig. 4] The numerical validation of Eq. (6) is described only as 'best fitted', with no statement of how n2 and sigma were chosen, no error bars, no goodness-of-fit measure, and no comparison on the same plot between simulation and theory. Since Eq. (1) is itself used to construct the fit, the agreement cannot be considered an independent test of the analytical prediction; the procedure is circular.
  4. [Results and Discussion, paragraph after Fig. 4] The paper admits that for defocusing nonlinearity with |4 n2 I_NL| > 1 the exponential in Eq. (6) becomes imaginary and the expression fails. Because the support of a physically meaningful intensity PDF extends to arbitrarily large intensities, this failure affects a significant part of the distribution. Consequently, Eq. (6) cannot be claimed to describe defocusing media, contradicting the abstract's assertion that the tail is reduced in that case.
minor comments (4)
  1. [Theory, Eqs. (1)-(2)] Eqs. (1) and (2) are identical; the duplication appears to be an editing artifact and should be removed.
  2. [Theory, Eq. (6)] The notation I_Nl in the square root is inconsistent with I_NL used elsewhere.
  3. [Results and Discussion, Fig. 2] The definition of C_I(delta r) as a sum over M realizations with normalization by M is not standard; the normalized correlation (DOC) should be explicitly defined, including how the mean is subtracted and how the zero-lag value is normalized.
  4. [Introduction] The sentence 'The available literature reports on the experimental and numerical studies of the of non-Rayleigh speckles' contains a duplicated 'of'; also 'analytically driven results' in the abstract is a typo for 'analytically derived results'.

Circularity Check

2 steps flagged · score 6.0 of 10

Eq. (6) is the pushforward of the Rayleigh PDF through the assumed map I_NL = I + n2 I^2 (Eq. 1), and the numerical validation is described as a best fit; the non-Rayleigh prediction reduces to its own input.

  1. self definitional [Theory, Eqs. (1)-(6)]
    "the intensity in such a medium can be represented as [30,31] I_NL = |E|^2 + n2|E|^4 = I + n2 I^2. (1) ... Here, we use the variable transformation method to derive the PDF of the intensity in nonlinear media. ... rho(I_NL) = 1/(2 sigma^2 sqrt(1+4 n2 I_NL)) exp[(1 - sqrt(1+4 n2 I_NL))/(4 n2 sigma^2)] (6)"

    Equation (6) is obtained by applying the change-of-variables formula (Eq. 5) to the assumed pointwise monotonic map (Eq. 1). The non-Rayleigh character of the result is entirely contained in the term n2 I^2 chosen in Eq. (1); no physical derivation of that map from Maxwell's equations or from the NLSE used in Eq. (7) is given. The 'analytically derived' PDF is therefore the assumed nonlinear-intensity relation restated in probability form, not an independent prediction.

  2. fitted input called prediction [Results and Discussion, Fig. 4]
    "These numerical outcomes were best fitted by our analytical expression derived in Eqn. (6). The plots of the analytically derived intensity probability distributions are shown in Fig. 4 (b)"

    The numerical simulation is presented as validation, but the paper states the analytical expression is 'best fitted' to the numerical outcomes. Since Eq. (6) contains the scale sigma and the same n2 used in the NLSE, fitting it to the simulated histogram does not test the model independently; it demonstrates that the assumed map can be adjusted to match the data. The paper also concedes that for defocusing nonlinearity with |4 n2 I_NL| > 1 the expression becomes imaginary and 'fails to predict the PDF', further limiting the claimed agreement.

full rationale

The central claim, Eq. (6), follows mathematically from the linear exponential intensity PDF and the assumed relation I_NL = I + n2 I^2 (Eq. 1). That relation is the sole source of the non-Rayleigh tail; it is asserted with citations to textbooks rather than derived from the NLSE (Eq. 7), whose physical intensity remains |E|^2 while nonlinearity enters the propagation operator. Thus the PDF is a change of variables on an extra input, not a consequence of the nonlinear wave equation used for validation. The numerical agreement is described as a 'best fit' of Eq. (6) to the simulated outcomes, so it is not an independent prediction. The paper's own caveat that the expression fails for defocusing nonlinearity when |4 n2 I_NL| > 1 further undercuts the claimed unified description. There is no load-bearing self-citation here; the circularity is internal to the assumed map and the fitted comparison. The variable-transformation algebra itself is internally consistent, so the paper is not a complete tautology, but the main 'analytical derivation' reduces by construction to its assumed nonlinear-intensity map.

Assumptions & free parameters 1 free parameters · 5 assumptions · 0 invented entities

The central claim rests on a nonstandard and unjustified relation between nonlinear and linear intensities. The numerical validation is a curve fit, so the formula's predictive power is untested. The internal inconsistency about Gaussian field statistics further undermines the physical basis.

free parameters (1)
  • sigma (sigma) = not specified
    The linear intensity PDF is assumed exponential with parameter 1/(2 sigma^2). If sigma is adjusted to fit the numerical PDF, which the 'best fitted' phrasing suggests is possible, it acts as a free parameter.
assumptions (5)
  • standard math The intensity PDF in a linear medium is exponential: rho(I) = 1/(2 sigma^2) exp(-I/(2 sigma^2)).
    Well-known result for fully developed speckle with Gaussian complex amplitude, used as the starting point of the transformation.
  • ad hoc to paper The intensity in a nonlinear medium is I_NL = I + n2 I^2.
    Stated in Eq. (1) without derivation; standard nonlinear optics does not add a term to the intensity itself, rather the nonlinearity modifies the refractive index. This is the foundational assumption of the paper.
  • standard math The variable transformation formula rho_I_NL(I_NL) = rho(f^{-1}(I_NL)) |dI/dI_NL| is valid.
    Standard probability theory for monotonic transformations, used in Eq. (5).
  • domain assumption The NLSE (Eq. 7) accurately models speckle propagation in the nonlinear medium.
    The numerical validation relies on this standard model for paraxial nonlinear propagation.
  • ad hoc to paper The field distribution follows Gaussian statistics and remains invariant in all nonlinear conditions.
    Stated in the text and used to argue field correlation invariance, but it contradicts the reported non-exponential intensity PDF.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Statistics of Non-Rayleigh Speckles Generated from Nonlinear Media." pith.science (2026). https://pith.science/paper/PJGMNSBH

@misc{pith2026250604816,
  author       = {Pith},
  title        = {Pith review of: Statistics of Non-Rayleigh Speckles Generated from Nonlinear Media},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/PJGMNSBH}},
  note         = {Machine review of arXiv:2506.04816}
}
read the original abstract

We analytically derive an expression for a speckle field's intensity probability density function (PDF) in a nonlinear medium. The analytically driven results are in good agreement with the numerical outcomes. In a focusing nonlinear medium, the local intensity of the speckle is enhanced as manifested through the longer tail of the PDF. In contrast, the local intensity of speckle is reduced in the presence of a defocusing nonlinearity, and the tail of the probability density function also reduces. This change in local intensity of the speckles arises due to the cubic Kerr nonlinearity, which eventually modifies the second-order statistics. Hence, the intensity correlation is altered as per the nature of the associated nonlinearity while the field correlation remains invariant of both types of the nonlinear conditions.

Discussion (0). Sign in to comment.

Reference graph

Works this paper leans on

32 extracted references · 32 canonical work pages

  1. [1]

    J. W. Goodman, Second Edi., (Roberts and Company Publishers, Greenwood Village (Colorado), 2007)

    work page 2007

  2. [2]

    J. C. Dainty, V ol. 9 (Springer Science & Business Media, New York, 2013)

    work page 2013

  3. [3]

    Sutton, S

    M. Sutton, S. G.J. Mochrie, T. Greytak, S. E. Nagler, L. E. Berman, G.A. Held, and G.B. Stephenson, Nature 352, 608 (1991)

    work page 1991

  4. [4]

    Wagner, S

    R. Wagner, S. Smith, J. Sandrik, and H. Lopez, IEEE Trans. SonicsUltrason.30,156(1983)

    work page 1983

  5. [5]

    J. W. Goodman, Second Edi., (John Wiley & Sons, 2015)

    work page 2015

  6. [6]

    Apostol and A

    A. Apostol and A. Dogariu, Phys. Rev. Lett. 91 (2003)

    work page 2003

  7. [7]

    Apostol and A

    A. Apostol and A. Dogariu, Phys. Rev. E 72, 025602R (2005)

    work page 2005

  8. [8]

    M. P. Cagigal and V . F. Canales, Opt. Eng. 40, 2690 (2001)

    work page 2001

Show all 32 references
  1. [9]

    F. L. Metz, C. Bonatto, and S. D. Prado, Phys. Rev. A 109, 013501(2024)

    2024

  2. [10]

    Strudley, D

    T. Strudley, D. Akbulut, W. L. V os, A. Lagendijk, A. P. Mosk, and O. L. Muskens, Opt. Lett. 39, 6347-6350 (2014)

    2014

  3. [11]

    Th. M. Nieuwenhuizen and M. C.W. van Rossum, Phys. Rev. Lett. 74 (1995)

    1995

  4. [12]

    Wolf, (Cambridge university press, 2007)

    E. Wolf, (Cambridge university press, 2007)

    2007

  5. [13]

    S. Liu, Z. Liu, C. Hu, E. Li, X. Shen, and S. Han, Opt. Commun. 472, 126017 (2020)

    2020

  6. [14]

    Zhang, W

    S. Zhang, W. Wang, R. Yu, and X. Yang, Laser Phys. 26, 055007 (2016)

    2016

  7. [15]

    Z. Li, X. Nie, F. Yang, X. Liu, D. Liu, X. Dong, X. Zhao, T. Peng, M. S. Zubairy, and M. O. Scully, Opt. Express 29, 19621 (2021)

    2021

  8. [16]

    Bromberg, and H

    Y . Bromberg, and H. Cao, Phys. Rev. Lett. 112, 213904 (2014)

    2014

  9. [17]

    S. Y . Han, N. Bender, and H. Cao, Phys. Rev. Lett. 130, 093802 (2023)

    2023

  10. [18]

    Bender, H

    N. Bender, H. Yılmaz, Y . Bromberg, and H. Cao, APL Photon. 4, 110806 (2019)

    2019

  11. [19]

    Bender, H

    N. Bender, H. Yılmaz, Y . Bromberg, and H. Cao, Optica 5, 595-600 (2018)

    2018

  12. [20]

    Bender, Henry Haig, Demetrios N

    N. Bender, Henry Haig, Demetrios N. Christodoulides, and F. W. Wise, Optica 10, 1260-1268 (2023)

    2023

  13. [21]

    J. P. Amaral, E. J. S. Fonseca, and A. J. Jesus-Silva, Phys. Rev. A 92,063851 (2015)

    2015

  14. [22]

    R. Liu, B. Qing, S. Zhao, P . Zhang, H. Gao, S. Chen and F. Li, Phys. Rev. Lett. 127, 180601 (2021)

    2021

  15. [23]

    J. G. M. N. Neto, E. J. S. Fonseca, and A. J. Jesus-Silva, Phys. Rev. A 106, 053519 (2022)

    2022

  16. [24]

    Bromberg, Y

    Y . Bromberg, Y . Lahini, E. Small and Y . Silberberg, Nphoton. 4, 721-726 (2010)

    2010

  17. [25]

    Derevyanko and E

    S. Derevyanko and E. small, Phys. Rev A 85,053816 (2012)

    2012

  18. [26]

    Frostig, E

    H. Frostig, E. Small, A. Daniel, P. Oulevey, S. Derevyanko and Y . Silberberg, Optica 4, 1073-1079 (2017)

    2017

  19. [27]

    S. B. Alves, H. L. D. de S. Cavalcante, G. F. de Oliveira Jr., T. P. de Silans, I. Vidal, M. Chevrollier, and M. Oriá, Phys. Rev. A 99,033838 (2019)

    2019

  20. [28]

    I. L. Rodrigues, J. P. Lopez, and T. P. de Silans, Phys. Rev. A 110,023525 (2024)

    2024

  21. [29]

    M. J. B. Crispim, C. C. S. Pereira, N. T. C. Oliveira, M. Chevrollier, R. A. de Oliveira, W. S. Martins and A. S. Reyna, Sci. Rep. 7239 (2023)

    2023

  22. [30]

    G. P. Agrawal, Fourth Edi., (Academic Press, 2001)

    2001

  23. [31]

    D. L. Mills, Second Edi., Springer (2012)

    2012

  24. [33]

    Mishra, S

    S. Mishra, S. K. Gautam, D. N. Naik, Z. Chen, J. Pu and R. K. Singh, J. Opt. 20, 125605 (2018)

    2018

Pith tools

Reviewed August 7, 2026 · model on record in the stance chip above.