pith. sign in

arxiv: 2604.08129 · v1 · submitted 2026-04-09 · 🧮 math.PR

Polarity of points for Gaussian random fields in critical dimension

Youssef Hakiki , Cheuk Yin Lee , Yimin Xiao This is my paper

Pith reviewed 2026-05-10 17:57 UTC · model grok-4.3

classification 🧮 math.PR
keywords gaussian random fieldspolaritycritical dimensionhitting pointssojourn timeshausdorff measurestationary incrementsvariance function
0
0 comments X

The pith

In the critical dimension, points are polar for Gaussian random fields exactly when the log-variance exponent gamma is at most 1/d.

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

The paper determines when fixed points are almost surely not hit by continuous Gaussian random fields that have stationary increments and a regularly varying variance of index H. It first gives a sufficient condition for polarity in terms of the variance growth. Then, focusing on the critical dimension d equals N over H and the specific variance sigma of r to the H times log to the gamma, it proves polarity holds if and only if gamma is at most one over d. This is the same as when the integral of r to the N minus one over sigma to the d diverges. The necessity of this integral condition holds under general assumptions, and the sufficiency is shown by extending a covering argument from Talagrand using sojourn time estimates to get Hausdorff measure bounds.

Core claim

In the critical dimension d = N/H, for sigma(r) = r^H (log(1/r))^gamma, points are polar if and only if gamma ≤ 1/d, or equivalently when the integral from 0+ of r^{N-1} / sigma^d(r) dr diverges. This integral condition is necessary for polarity under general assumptions on the variance. The sufficiency proof in the specific case extends Talagrand's covering argument based on sojourn time estimates to obtain the required Hausdorff measure bounds.

What carries the argument

The divergence of the integral ∫ r^{N-1} / σ^d(r) dr near zero, which is the test for polarity of points in the critical dimension for this class of fields.

If this is right

  • When the integral diverges, no fixed point is hit almost surely.
  • When gamma > 1/d the integral converges and points are not polar.
  • The integral condition is necessary for polarity for general regularly varying sigma.
  • The sojourn time covering method yields Hausdorff measure estimates sufficient to prove non-hitting.

Where Pith is reading between the lines

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

  • The integral test could be checked for polarity of higher-dimensional sets such as curves or surfaces hit by the field.
  • Results like this may help classify the sample path properties of Gaussian fields in terms of their modulus of continuity and dimension.
  • Discretized versions of the field could be simulated to observe the hitting behavior near the critical value of gamma.

Load-bearing premise

Extending Talagrand's covering argument with sojourn time estimates is sufficient to obtain the Hausdorff measure bounds needed for proving that points are polar when the integral diverges.

What would settle it

Numerical simulation of many sample paths of the field with gamma slightly larger than 1/d to check if the probability of hitting a fixed point is positive, which would be inconsistent with the polarity claim if it occurs when the integral diverges.

read the original abstract

We study the property of hitting points for a class of $\mathbb{R}^d$-valued continuous Gaussian random fields on $\mathbb{R}^N$ with stationary increments, i.i.d. coordinates, and a regularly varying variance function $\sigma$ of index $0<H<1$. We first prove that if \[ \lim_{r\to 0^+} \frac{r^N}{\sigma^d\left(r\left( \log\log\frac{1}{r}\right)^{-1/N}\right)} = \infty, \] then every fixed point is polar (i.e., not hit almost surely). In general, this criterion may not be optimal in the critical dimension $d=N/H$. To aim for an optimal condition, we consider the specific case $\sigma(r) = r^H (\log(1/r))^\gamma$ and prove that, in the critical dimension $d=N/H$, points are polar if and only if $\gamma \le 1/d$, or equivalently in this specific case, \[ \int_{0^+} \frac{r^{N-1}}{\sigma^d(r)} dr = \infty. \] This integral condition is also necessary for points to be polar under general assumptions. Our main contribution lies in the proof of sufficiency of this condition in the specific case, where we extend a covering argument of Talagrand (1998) based on sojourn time estimates to obtain Hausdorff measure bounds and solve polarity of points in the critical dimension.

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.

Referee Report

2 major / 2 minor

Summary. The paper studies polarity of fixed points for R^d-valued continuous Gaussian random fields on R^N with stationary increments, i.i.d. coordinates, and regularly varying variance σ of index 0 < H < 1. It proves a general sufficient condition for polarity via the limit lim r→0+ r^N / σ^d(r (log log 1/r)^{-1/N}) = ∞. In the critical dimension d = N/H, for the specific case σ(r) = r^H (log(1/r))^γ, it establishes that points are polar if and only if γ ≤ 1/d (equivalently, the integral ∫_{0+} r^{N-1}/σ^d(r) dr diverges). Necessity of the integral condition is claimed under general assumptions, while sufficiency in the specific case is obtained by extending Talagrand (1998)'s covering argument based on sojourn-time estimates to derive the required Hausdorff measure bounds.

Significance. If the extension of the covering argument holds, the work supplies an optimal if-and-only-if characterization of point polarity precisely in the critical dimension, a longstanding boundary case in the theory of Gaussian random fields and their sample-path properties. The adaptation of sojourn-time estimates to absorb logarithmic corrections would constitute a technical contribution with potential applicability to other critical regimes in hitting probabilities and Hausdorff dimensions for Gaussian processes.

major comments (2)
  1. [Abstract / sufficiency proof section] Abstract and main sufficiency argument: the claim that the integral condition implies polarity (i.e., the Hausdorff measure of the level set is zero) in the critical dimension for σ(r) = r^H (log(1/r))^γ rests on extending Talagrand (1998)'s covering argument via sojourn-time estimates. The manuscript must explicitly verify that the extra (log(1/r))^γ factors are absorbed uniformly in the covering numbers and probability bounds; at the boundary γ = 1/d any unabsorbed logarithmic growth would prevent the measure from vanishing and invalidate the sufficiency direction of the iff statement.
  2. [Abstract / necessity argument] General necessity claim: the integral condition is asserted to be necessary for polarity under general assumptions on σ, yet the sufficiency proof is restricted to the specific log-power case. The manuscript should delineate the precise hypotheses under which necessity holds and confirm that the general criterion does not leave a gap at the critical dimension where the limit condition is known to be non-optimal.
minor comments (2)
  1. Notation: ensure consistent use of the regularly varying index H and the critical dimension d = N/H throughout; the relation between the general limit criterion and the integral test should be stated more explicitly for the specific σ.
  2. References: complete the citation to Talagrand (1998) and clarify which results from that paper are invoked verbatim versus adapted.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the detailed and constructive report. The comments highlight important points regarding the clarity of our proofs and statements, which we will address through targeted revisions and clarifications. Below we respond point by point to the major comments.

read point-by-point responses

  1. Referee: The claim that the integral condition implies polarity in the critical dimension for σ(r) = r^H (log(1/r))^γ rests on extending Talagrand (1998)'s covering argument. The manuscript must explicitly verify that the extra (log(1/r))^γ factors are absorbed uniformly in the covering numbers and probability bounds; at the boundary γ = 1/d any unabsorbed logarithmic growth would prevent the measure from vanishing.

    Authors: We agree that explicit verification of the absorption of the logarithmic factors is essential for rigor, particularly at the boundary γ = 1/d. In the revised manuscript we will augment the sufficiency proof (in the section extending Talagrand's covering argument) with a dedicated lemma that tracks the (log(1/r))^γ terms through the sojourn-time estimates, covering numbers, and probability bounds. This will include direct computation showing uniform absorption even when γ = 1/d, ensuring the resulting Hausdorff measure vanishes as claimed. revision: yes

  2. Referee: The integral condition is asserted to be necessary for polarity under general assumptions on σ, yet the sufficiency proof is restricted to the specific log-power case. The manuscript should delineate the precise hypotheses under which necessity holds and confirm that the general criterion does not leave a gap at the critical dimension where the limit condition is known to be non-optimal.

    Authors: The necessity result is proved under the standing assumptions of the paper: stationary increments, i.i.d. coordinates, and σ regularly varying of index H. We will revise the introduction and the statement of the necessity theorem to list these hypotheses explicitly. We will also add a clarifying remark noting that the general sufficient limit criterion is not always sharp in the critical dimension d = N/H, which motivates the separate sharp analysis for the log-power case; this removes any apparent gap. revision: yes

Circularity Check

0 steps flagged

No circularity: derivation extends external Talagrand argument independently

full rationale

The paper first proves a general sufficient condition for polarity via the limit lim r^N / σ^d(r (log log 1/r)^{-1/N}) = ∞. For the specific σ(r) = r^H (log(1/r))^γ in critical dimension d = N/H, it establishes the iff statement with γ ≤ 1/d. This equivalence follows directly from substituting the power-log form into the integral ∫ r^{N-1}/σ^d(r) dr and simplifying the exponent (Hd = N yields ∫ dr/(r (log)^{γ d})), which diverges precisely when γ d ≤ 1. Sufficiency of the integral condition is obtained by extending Talagrand (1998)'s covering argument based on sojourn-time estimates to control the extra logarithmic factors and produce the required Hausdorff measure bounds. Necessity of the integral holds under general assumptions on σ. No step reduces by construction to a fitted parameter, self-definition, or load-bearing self-citation; the central proof imports an external technique whose validity is independent of the present results.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The results rest on standard domain assumptions for the class of Gaussian fields studied and the validity of extending Talagrand's covering argument; no free parameters or invented entities are introduced.

axioms (2)
  • domain assumption The random field is continuous, has stationary increments, and i.i.d. coordinates.
    Explicitly stated as the class of fields under study in the abstract.
  • domain assumption The variance function σ is regularly varying of index 0<H<1.
    Assumed for the general polarity criterion.

pith-pipeline@v0.9.0 · 5571 in / 1214 out tokens · 63969 ms · 2026-05-10T17:57:30.833425+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

43 extracted references · 43 canonical work pages

  1. [1]

    Anderson, T. W. (1955). The integral of a symmetric unimodal function over a symmetric convex set and some probability inequalities. Proc. Amer. Math. Soc. 6, 170–176

    work page 1955

  2. [2]

    V., Angulo, J

    Anh, V. V., Angulo, J. M. and Ruiz-Medina, M. D. (1999). Possible long-range dependence in fractional random fields. J. Statist. Plann. Inference, 80, 95–110

    work page 1999

  3. [3]

    Berman, S. M. (1969). Local times and sample function properties of stationary Gaussian processes. Trans. Amer. Math. Soc. 137, 277–299

    work page 1969

  4. [4]

    and Xiao, Y

    Biermé, H., Lacaux, C. and Xiao, Y. (2009). Hitting probabilities and the Hausdorff dimension of the inverse images of anisotropic Gaussian random fields, Bull. Lond. Math. Soc. 41, no. 2, 253–273

    work page 2009

  5. [5]

    H., Goldie, C

    Bingham, N. H., Goldie, C. M. and Teugels, J. L. (1989). Regular variation. Encyclopedia Math. Appl., 27, Cambridge University Press, Cambridge

    work page 1989

  6. [6]

    C., Khoshnevisan, D

    Dalang, R. C., Khoshnevisan, D. and Nualart, E. (2007). Hitting probabilities for systems of non-linear stochastic heat equations with additive noise. ALEA Lat. Am. J. Probab. Math. Stat. 3, 231–271

    work page 2007

  7. [7]

    C., Khoshnevisan, D

    Dalang, R. C., Khoshnevisan, D. and Nualart, E. (2009). Hitting probabilities for systems for non- linear stochastic heat equations with multiplicative noise. Probab. Theory Related Fields 144, no. 3-4, 371–427

    work page 2009

  8. [8]

    C., Khoshnevisan, D

    Dalang, R. C., Khoshnevisan, D. and Nualart, E. (2013). Hitting probabilities for systems of non-linear stochastic heat equations in spatial dimensionk≥1. Stoch. Partial Differ. Equ. Anal. Comput. 1, no. 1, 94–151,

    work page 2013

  9. [9]

    C., Mueller, C

    Dalang, R. C., Mueller, C. and Xiao, Y. (2017). Polarity of points for Gaussian random fields. Ann. Probab. 45, no. 6B, 4700–4751

    work page 2017

  10. [10]

    Dalang, R. C. and Nualart, E. (2004). Potential theory for hyperbolic SPDEs. Ann. Probab. 32, no. 3A, 2099–2148,

    work page 2004

  11. [11]

    Dalang, R. C. and Sanz-Solé, M. (2010). Criteria for hitting probabilities with applications to systems of stochastic wave equations. Bernoulli 16, no. 4, 1343–1368

    work page 2010

  12. [12]

    Dalang, R. C. and Sanz-Solé, M. (2015). Hitting probabilities for nonlinear systems of stochastic waves. Mem. Amer. Math. Soc. 237, no. 1120,

    work page 2015

  13. [13]

    Dudley, R.M.(1967).ThesizesofcompactsubsetsofHilbertspaceandcontinuityofGaussianprocesses. J. Functional Analysis 1, 290–330

    work page 1967

  14. [14]

    and Hakiki, Y

    Erraoui, M. and Hakiki, Y. (2025). Fractional Brownian motion with deterministic drift: how critical is drift regularity in hitting probabilities. Math. Proc. Camb. Philos. Soc. 178, no. 1, 103–132

    work page 2025

  15. [15]

    Falconer, K. (2013). Fractal geometry: mathematical foundations and applications. John Wiley & Sons

    work page 2013

  16. [16]

    (2011), A local-time correspondence for stochastic partial differential equations

    Foondun, M., Khoshnevisan, D., and Nualart, E. (2011), A local-time correspondence for stochastic partial differential equations. Trans. Amer. Math. Soc., 363, no. 5, 2481–2515

    work page 2011

  17. [17]

    and Horowitz, J

    Geman, D. and Horowitz, J. (1980). Occupation densities. Ann. Probab. 8, no. 1, 1–67

    work page 1980

  18. [18]

    Hawkes, J. (1986). Local times as stationary processes. In: From Local Times to Global Geometry. El- worthy, K.D. (Ed.), pp. 111–120, Pitman Research Notes in Mathematics, Vol. 150. Longman, Chicago. 30 Y. HAKIKI, C. Y. LEE, AND Y. XIAO

    work page 1986

  19. [19]

    and Xiao, Y

    Herrell, R., Song, R., Wu, D. and Xiao, Y. (2020). Sharp space-time regularity of the solution to a stochastic heat equation driven by a fractional-colored noise. Stoch. Anal. Appl. 38, 747–768

    work page 2020

  20. [20]

    and Sanz-Solé, M

    Hinojosa-Calleja, A. and Sanz-Solé, M. (2021). Anisotropic Gaussian random fields: criteria for hitting probabilities and applications. Stoch. Partial Differ. Equ. Anal. Comput. 9, no. 4, 984–1030

    work page 2021

  21. [21]

    Kahane, J.-P. (1985). Some Random Series of Functions. 2nd edition, Cambridge University Press

    work page 1985

  22. [22]

    Kesten, H. (1969). Hitting probabilities of single points for processes with stationary independent in- crements. Amer. Math. Soc., Providence RI

    work page 1969

  23. [23]

    and Shi, Z

    Khoshnevisan, D. and Shi, Z. (1999). Brownian sheet and capacity. Ann. Probab. 27, no. 3, 1135–1159

    work page 1999

  24. [24]

    and Zhong, Y

    Khoshnevisan, D., Xiao, Y. and Zhong, Y. (2003a). Measuring the range of an additive Lévy process. Ann. Probab. 31, 1097–1141

    work page

  25. [25]

    and Zhong, Y

    Khoshnevisan, D., Xiao, Y. and Zhong, Y. (2003b). Local times of additive Lévy processes. Stoch. Process. Appl. 104, 193–216

    work page

  26. [26]

    and Matlak, D

    Latała, R. and Matlak, D. (2017). Royen’s proof of the Gaussian correlation inequality, Geometric Aspects of Functional Analysis, Lecture Notes in Math., vol. 2169, Springer, Cham, pp. 265–275

    work page 2017

  27. [27]

    Y., Song, J., Xiao, Y

    Lee, C. Y., Song, J., Xiao, Y. and Yuan, Y. (2023). Hitting probabilities of Gaussian random fields and collision of eigenvalues of random matrices. Trans. Amer. Math. Soc. 376, no. 6, 4273–4299

    work page 2023

  28. [28]

    Lee, C. Y. and Xiao, Y. (2026). Hitting probability, thermal capacity, and Hausdorff dimension results for the Brownian sheet. Electron. J. Probab. 31, no. 26, 1–31

    work page 2026

  29. [29]

    Marcus, M. B. and Rosen, J. (2006). Markov Processes, Gaussian Processes, and Local Times. Cam- bridge University Press, Cambridge

    work page 2006

  30. [30]

    Mattila, P. (1995). Geometry of sets and measures in Euclidean spaces. Fractals and rectifiability. Cambridge Stud. Adv. Math., 44, Cambridge University Press

    work page 1995

  31. [31]

    Pitt, L. D. (1978). Local times for Gaussian random fields. Indiana Univ. Math. J., 27, 309–330

    work page 1978

  32. [32]

    Pitman, E. J. G. (1968). On the behavior of the characteristic function of a probability distribution in the neighbourhood of the origin. J. Australian Math. Soc. Series A, 8, 422–443

    work page 1968

  33. [33]

    Royen, T. (2014). A simple proof of the Gaussian correlation conjecture extended to some multivariate gamma distributions, Far East J. Theor. Stat. 48, no. 2, 139–145

    work page 2014

  34. [34]

    Talagrand, M. (1995). Hausdorff measure of trajectories of multiparameter fractional Brownian motion. Ann. Probab. 23, no. 2, 767–775

    work page 1995

  35. [35]

    Talagrand, M. (1998). Multiple points of trajectories of multiparameter fractional Brownian motion. Probab. Theory Relat. Fields 112, 545–563

    work page 1998

  36. [36]

    Testard, F. (1986). Polarité, points multiples et géométrie de certain processus gaussiens. Publ. du Laboratoire de Statistique et Probabilités de l’U.P.S., Toulouse, 01 - 86

    work page 1986

  37. [37]

    Xiao, Y. (1996). Hausdorff measure of the sample paths of Gaussian random fields. Osaka J. Math. 33, no. 4, 895–913

    work page 1996

  38. [38]

    Xiao, Y. (1997). Hölder conditions for the local times and the Hausdorff measure of the level sets of Gaussian random fields. Probab. Theory Related Fields 109, no. 1, 129–157

    work page 1997

  39. [39]

    Xiao, Y. (1999). Hitting probabilities and polar sets for fractional Brownian motion. Stochastics Stochas- tics Rep. 66, no. 1-2, 121–151

    work page 1999

  40. [40]

    Xiao, Y. (2007). Strong local nondeterminism and the sample path properties of Gaussian random fields. In: Asymptotic Theory in Probability and Statistics with Applications (Tze Leung Lai, Qiman Shao, Lianfen Qian, editors), pp. 136–176, Higher Education Press, Beijing

    work page 2007

  41. [41]

    Xiao, Y. (2009). Sample path properties of anisotropic Gaussian random fields. In: A minicourse on stochastic partial differential equations. pp. 145–212. Lecture Notes in Math., vol. 1962, Springer, Berlin

    work page 2009

  42. [42]

    Yaglom, A. M. (1957), Certain types of random fields in n-dimensional space similar to stationary stochastic processes. Teor. Veroyatnost. i Primenen. 2, 292–338

    work page 1957

  43. [43]

    Yaglom, A. M. (1987), Correlation theory of stationary and related random functions. Vol. I. Basic results. Springer Ser. Statist. Springer-Verlag, New York. POLARITY OF POINTS FOR GAUSSIAN RANDOM FIELDS 31 Department of Mathematics, Purdue University, West Laf ayette, IN 47907, United States Email address:yhakiki@purdue.edu School of Science and Engineer...

    work page 1987