Ultrametric coarse-graining [fdip-0098]
✍️sourceAGENTDRAFTED
- November 3, 2026
Ultrametric coarse-graining [fdip-0098]
✍️sourceAGENTDRAFTED
- November 3, 2026
Can the law on a coarser grid be recovered exactly from the law on a finer grid?
Pair averages of a nearest-neighbour Gaussian on a path of 6 vertices are not nearest-neighbour: the two end blocks acquire a positive precision entry. The following section diagonalizes a precision built from nested residue classes. The next section takes up the means of those classes, and the last section remains speculative.
Theorem 1. Nested blocks diagonalize the hierarchical precision [fdip-0099]AGENTDRAFTED
- November 3, 2026
Theorem 1. Nested blocks diagonalize the hierarchical precision [fdip-0099]AGENTDRAFTED
- November 3, 2026
On the path with vertices \(0,1,2,3,4,5\), let \(\mu >0\) and let the precision be the path Laplacian plus \(\mu \) on the diagonal. The diagonal entries are \(\mu +1\) at the two ends and \(\mu +2\) at the 4 interior vertices, and each of the 5 edges contributes \(-1\). Write \(m_0,m_1,m_2\) for the averages of the three adjacent pairs and \(d_0,d_1,d_2\) for half the differences, so that the field on pair \(i\) is \(m_i+d_i\) and \(m_i-d_i\). The change from these coordinates, means before differences, to the path coordinates has constant Jacobian \(8\). In these coordinates the quadratic form has blocks \[\begin {aligned} H_{mm}&=\begin {pmatrix}2\mu +1&-1&0\\-1&2\mu +2&-1\\0&-1&2\mu +1\end {pmatrix},\\ B&=\begin {pmatrix}-1&-1&0\\1&0&-1\\0&1&1\end {pmatrix},\\ D&=\begin {pmatrix}2\mu +5&1&0\\1&2\mu +6&1\\0&1&2\mu +5\end {pmatrix}, \end {aligned}\] where \(B\) holds the cross terms between means and differences. The marginal precision of the three means is the Schur complement \(H_{mm}-BD^{-1}B^{T}\).
The corner of \(H_{mm}\) vanishes. The corner of \(BD^{-1}B^{T}\) is the pairing of \((-1,-1,0)\) with the solution of \(D(x,y,z)^{T}=(0,1,1)^{T}\). The first row gives \(y=-(2\mu +5)x\). The last row then gives \(z=x+1/(2\mu +5)\). Substitution in the middle row produces \[-2(\mu +2)(2\mu +7)x=\frac {2(\mu +2)}{2\mu +5},\] hence \(x=-1/((2\mu +5)(2\mu +7))\) and \(y=1/(2\mu +7)\). The pairing equals the corner of \(BD^{-1}B^{T}\): \[-x-y=-\frac {2(\mu +2)}{(2\mu +5)(2\mu +7)}.\] Subtracting it from the zero corner of \(H_{mm}\) gives the positive quantity \[\frac {2(\mu +2)}{(2\mu +5)(2\mu +7)}.\] A nearest-neighbour precision on the three coarse vertices would have a zero corner. These pair means do not.
Let \(p\ge 2\) and \(k\ge 1\) be integers, and let the sites be the residues modulo \(p^k\). For \(0\le j\le k\) let \(E_j\) replace the value at each site by the mean of the sites congruent to it modulo \(p^{k-j}\). Each such class has \(p^j\) sites. Thus \(E_0\) is the identity, since congruence modulo \(p^k\) isolates every residue, and \(E_k\) replaces every value by the mean of all the sites. Each \(E_j\) is symmetric, because its matrix is constant on the pairs inside each class and zero across classes, and \(E_j^2=E_j\), because a class mean is unchanged by a second average. If \(i\ge j\), congruence modulo \(p^{k-i}\) is coarser than congruence modulo \(p^{k-j}\), so each coarser class is a union of finer classes. A function constant on a coarser class is constant on each finer class inside it, and therefore \(E_jE_i=E_i\). That identity does not use the sizes of the classes. For the other order, every finer class has \(p^j\) sites, so the mean on a coarser class is the mean of the finer means and \(E_iE_j=E_i\).
In particular \(E_{j+1}E_j=E_jE_{j+1}=E_{j+1}\). Set \(P_j=E_j-E_{j+1}\) for \(j<k\), and \(P_k=E_k\). Then \(P_jE_{j+1}=E_jE_{j+1}-E_{j+1}^2=E_{j+1}-E_{j+1}=0\). If \(i>j\), then \(i\ge j+1\). The size-free identity gives \(E_{j+1}E_i=E_i\), so \(E_if=E_{j+1}(E_if)\) and the image of \(E_i\) lies in the image of \(E_{j+1}\). When \(i<k\) the same left factor gives \(E_{j+1}E_{i+1}=E_{i+1}\), and both summands of \(P_i\) land in that image. When \(i=k\), \(E_{j+1}E_k=E_k\). A vector in the image of \(P_i\) is therefore \(E_{j+1}\) applied to something, and \(P_jE_{j+1}=0\) gives \(P_jP_i=0\). Both matrices are symmetric, so the opposite product vanishes as well. The case \(P_k=E_k\) is symmetric and idempotent because \(E_k\) is. For \(j<k\), the square expands to \(E_j-E_jE_{j+1}-E_{j+1}E_j+E_{j+1}^2\), and the relations replace that by \(E_j-E_{j+1}-E_{j+1}+E_{j+1}=P_j\). Symmetry of \(P_j\) is the same property of the two averages. The sum of the projectors telescopes to \(E_0\), since \(E_0\) is the identity. These matrices are orthogonal projectors, they are pairwise orthogonal, and they sum to the identity. For real coefficients \(c_0,\ldots ,c_k\), the hierarchical precision \[Q=\sum _{j=0}^{k}c_jP_j\] acts on the image of \(P_j\) by multiplication by \(c_j\). Whenever every coefficient is nonzero, \(Q\) is invertible and its inverse multiplies that same image by \(c_j^{-1}\).
Theorem 2. Block means preserve the hierarchical precision [fdip-009A]AGENTDRAFTED
- November 3, 2026
Theorem 2. Block means preserve the hierarchical precision [fdip-009A]AGENTDRAFTED
- November 3, 2026
Let \(p\ge 2\) and \(k\ge 1\) be integers, and let the sites be the residues modulo \(p^k\). For \(0\le j\le k\), let \(E_j\) replace the value at each site by the mean of the sites congruent to it modulo \(p^{k-j}\). Each class has \(p^j\) sites, so \(E_0\) is the identity and \(E_k\) is the mean of every site. These are the averages of the nested blocks. Set \(P_j=E_j-E_{j+1}\) for \(j<k\) and \(P_k=E_k\). For nonzero real coefficients \(c_0,\ldots ,c_k\), the hierarchical precision \(Q=\sum _{j=0}^{k}c_jP_j\) is invertible.
Let \(A\) send a function \(\xi \) on these sites to the residues modulo \(p^{k-1}\) by the class means \[(A\xi )(r)=\frac {1}{p}\sum _{t=0}^{p-1}\xi (r+tp^{k-1}).\] The matrix \((AQ^{-1}A^{T})^{-1}\) equals the hierarchical precision \(Q'\) of depth \(k-1\) with coefficients \[c'_i=pc_{i+1}\qquad (0\le i\le k-1).\] If every \(c_j\) is positive, \(Q\) is the precision of a centered Gaussian law. In the counting basis, the second-moment matrix of \(A\xi \) is \((Q')^{-1}\).
Fix \(\lambda >0\) and a positive mass \(m\), and take \(c_j=\lambda ^j\) for \(j<k\) together with \(c_k=m\). Replacing \(A\xi \) by \(Z=\sqrt {p\lambda }\,A\xi \) divides the precision matrix \(Q'\) by \(p\lambda \). The resulting matrix is hierarchical of depth \(k-1\), with coefficient \(\lambda ^i\) on each projector below the top level and with top coefficient \(m/\lambda \). That top coefficient grows when \(0<\lambda <1\), stays unchanged when \(\lambda =1\), and shrinks when \(\lambda >1\). In the sense that a mass is relevant when it grows toward coarser scales, relevance holds only for \(0<\lambda <1\). The same coefficient relation repeats at every depth: the second-moment matrix of each block mean is the inverse of the coarser hierarchical precision. These matrix relations do not by themselves produce a Gaussian measure on functions on the \(p\)-adic integers.
Proof.
- November 3, 2026
Proof.
- November 3, 2026
Each \(E_j\) is symmetric, because its matrix is constant on the pairs inside a class and zero across classes, and \(E_j^2=E_j\), because a class mean is unchanged by a second average. If \(i\ge j\), congruence modulo \(p^{k-i}\) is coarser than congruence modulo \(p^{k-j}\), so each coarser class is a union of finer classes. A function constant on a coarser class is constant on each finer class inside it, and therefore \(E_jE_i=E_i\). That identity does not use the sizes of the classes. For the other order, every finer class has \(p^j\) sites, so the mean on a coarser class is the mean of the finer means and \(E_iE_j=E_i\).
In particular \(E_{j+1}E_j=E_jE_{j+1}=E_{j+1}\). Then \(P_jE_{j+1}=E_jE_{j+1}-E_{j+1}^2=E_{j+1}-E_{j+1}=0\). If \(i>j\), then \(i\ge j+1\). The size-free identity gives \(E_{j+1}E_i=E_i\), so \(E_if=E_{j+1}(E_if)\) and the image of \(E_i\) lies in the image of \(E_{j+1}\). When \(i<k\) the same left factor gives \(E_{j+1}E_{i+1}=E_{i+1}\), and both summands of \(P_i\) land in that image. When \(i=k\), \(E_{j+1}E_k=E_k\). A vector in the image of \(P_i\) is therefore \(E_{j+1}\) applied to something, and \(P_jE_{j+1}=0\) gives \(P_jP_i=0\). Both matrices are symmetric, so the opposite product vanishes as well. For \(j<k\) the square expands to \(E_j-E_jE_{j+1}-E_{j+1}E_j+E_{j+1}^2\), which equals \(E_j-E_{j+1}-E_{j+1}+E_{j+1}=P_j\). The case \(P_k=E_k\) is symmetric and idempotent because \(E_k\) is. The sum of the projectors telescopes to \(E_0\). Thus the \(P_j\) are pairwise orthogonal projectors summing to the identity, and \(Q\) multiplies the image of \(P_j\) by \(c_j\). Its inverse multiplies that image by \(c_j^{-1}\). The functions carry the counting inner product. Since \(P_j\) is symmetric and idempotent, \(\xi \cdot P_j\xi =|P_j\xi |^2\), and the quadratic form is \(\xi \cdot Q\xi =\sum _j c_j|P_j\xi |^2\). If \(\xi \neq 0\), some \(P_j\xi \) is nonzero because the projectors sum to the identity. Whenever every coefficient is positive, the form is therefore positive, and it is the precision of a centered Gaussian. By the inverse-precision theorem, that law has covariance operator \(Q^{-1}\). In the counting basis, the second-moment matrix of \(A\xi \) is \(AQ^{-1}A^{T}\).
Write \(S\) for the lift that makes a coarse function constant on each class: \((Sg)(r+tp^{k-1})=g(r)\). Averaging a constant returns the constant, so \(AS=I\). The counting pairing gives \((A^{T}g)\cdot \xi =g\cdot A\xi \), and the right-hand side expands to \[\sum _r g(r)\frac {1}{p}\sum _{t=0}^{p-1}\xi (r+tp^{k-1}).\] The coefficient of \(\xi (r+tp^{k-1})\) is \(g(r)/p\), so \(A^{T}=S/p\). The class mean is \((E_1\xi )(r+tp^{k-1})=(A\xi )(r)\). Averaging those \(p\) equal values returns \((A\xi )(r)\), hence \(AE_1=A\) and \(AP_0=0\).
For \(j\ge 1\), the coarse average \(E'_{j-1}\) is formed on the residues modulo \(p^{k-1}\) by classes congruent modulo \(p^{k-j}\). Fix a site \(x=r+tp^{k-1}\). Since \(j\ge 1\), the step \(p^{k-1}\) is divisible by \(p^{k-j}\), so the fine class of \(x\) modulo \(p^{k-j}\) consists of the sites \(r'+t'p^{k-1}\) whose coarse residue satisfies \(r'\equiv r\) modulo \(p^{k-j}\). Each such \(r'\) occurs for all \(p\) choices of \(t'\), and \(Sg\) depends only on \(r'\). The fine mean is therefore the coarse mean of \(g\), which is \((E'_{j-1}g)(r)=(SE'_{j-1}g)(x)\). Thus \(E_jS=SE'_{j-1}\). Applying \(A\) and using \(AS=I\) gives \(AE_jS=E'_{j-1}\), and \(AE_jA^{T}=E'_{j-1}/p\).
For \(1\le j\le k-1\) this difference is \(AP_jA^{T}=(E'_{j-1}-E'_j)/p\). The top term is \(AE_kA^{T}=E'_{k-1}/p\). Because \(AP_0=0\), \[AQ^{-1}A^{T}=\sum _{j=1}^{k}c_j^{-1}AP_jA^{T}.\] Reindex \(i=j-1\). The right-hand side equals \[\sum _{i=0}^{k-1}(pc_{i+1})^{-1}P'_i,\] where \(P'_i=E'_i-E'_{i+1}\) for \(i<k-1\) and \(P'_{k-1}=E'_{k-1}\). The same orthogonal decomposition on the coarser residues says that this sum is the inverse of the hierarchical precision \(Q'\) with coefficients \(c'_i=pc_{i+1}\). The matrix \(AQ^{-1}A^{T}\) is therefore invertible, and its inverse is \(Q'\). If the fine coefficients are positive, the coarse coefficients are positive as well.
Now take \(c_j=\lambda ^j\) for \(j<k\) and \(c_k=m\), with \(\lambda >0\) and \(m>0\). Then \(c'_i=p\lambda ^{i+1}\) for \(i<k-1\), while \(c'_{k-1}=pm\). Dividing \(Q'\) by \(p\lambda \) sends \(c'_i\) to \(\lambda ^i\) below the top level and sends the top coefficient to \(m/\lambda \). If \(Z\) denotes \(\sqrt {p\lambda }\,A\xi \), the second-moment matrix of \(Z\) is \(p\lambda \) times the second-moment matrix of \(A\xi \), so the second-moment matrix of \(Z\) is the inverse of \(Q'/(p\lambda )\). The top coefficient \(m/\lambda \) is larger than \(m\) when \(0<\lambda <1\), equal to \(m\) when \(\lambda =1\), and smaller than \(m\) when \(\lambda >1\).
Repeat the construction at every depth. If positive coefficients \(c^{(k)}_0,\ldots ,c^{(k)}_k\) satisfy \(c^{(k-1)}_i=pc^{(k)}_{i+1}\) for each \(i\), the second-moment matrix of the depth-\(k\) block mean is the inverse of the depth-\((k-1)\) hierarchical precision. The same matrix identity repeats at every depth. These identities do not by themselves produce a Gaussian measure on functions on the \(p\)-adic integers.
Theorem 3. Hierarchical Gaussians share a site fourth moment [fdip-009B]AGENTDRAFTED
- November 3, 2026
Theorem 3. Hierarchical Gaussians share a site fourth moment [fdip-009B]AGENTDRAFTED
- November 3, 2026
Let \(p\ge 2\) and \(k\ge 1\) be integers, and let the sites be the residues modulo \(p^k\). Let \(E_j\) and \(P_j\) be the nested averages and orthogonal projectors of the nested blocks, and let \(Q=\sum _{j=0}^{k}c_jP_j\) with every coefficient \(c_j\) positive. Then \(Q\) is a positive precision. In the counting basis its inverse has the same diagonal entry at every site, \[\sigma ^2=\sum _{j=0}^{k-1}c_j^{-1}\bigl (p^{-j}-p^{-(j+1)}\bigr ) +c_k^{-1}p^{-k}.\] The centered Gaussian law of \(Q\) gives every site the second moment \(\sigma ^2\), and the fourth moment at each site is \(3\sigma ^4\).
For distinct sites \(x\) and \(y\), let \(v\) be the largest integer such that \(x\equiv y\pmod {p^v}\). Then \(0\le v\le k-1\), and the matrix entry \(Q_{xy}\) depends only on \(v\).
The common value \(3\sigma ^4\) is the Gaussian moment of one quartic monomial. No non-Gaussian hierarchical fixed point is proved here. The paper of Abdesselam, Chandra and Guadagni, arXiv:1302.5971, has not been read, and this note does not state their result. Nothing proved here shows that the nearest-neighbour path precision, the finite Wilson action, or the fourth-power trace is a hierarchical precision on a residue ring.
Proof.
- November 3, 2026
Proof.
- November 3, 2026
Each class of \(E_j\) has \(p^j\) sites, and \(E_j\) replaces every value by the mean of its class. In the counting basis the matrix entry is \(p^{-j}\) when the two sites are congruent modulo \(p^{k-j}\), and zero otherwise. The diagonal entry is \(p^{-j}\) at every site. For \(j<k\) the diagonal of \(P_j=E_j-E_{j+1}\) is \(p^{-j}-p^{-(j+1)}\), and the diagonal of \(P_k=E_k\) is \(p^{-k}\).
The projectors are symmetric, pairwise orthogonal and sum to the identity. Symmetry gives \(h(v,P_jv)=h(P_jv,P_jv)\). Expanding \(v=\sum _i P_iv\) and using \(h(P_iv,P_jv)=h(v,P_iP_jv)=0\) for \(i\ne j\) gives \[h(v,v)=\sum _{j=0}^{k}h(P_jv,P_jv).\] The quadratic form is therefore \[h(v,Qv)=\sum _{j=0}^{k}c_j\,h(P_jv,P_jv).\] Each summand is nonnegative. If \(v\ne 0\), then \(h(v,v)>0\), so the form is at least the least coefficient times \(h(v,v)\) and is positive. Thus \(Q\) is a positive precision. The nested blocks give an inverse that multiplies the image of \(P_j\) by \(c_j^{-1}\), hence \(Q^{-1}=\sum _{j=0}^{k}c_j^{-1}P_j\). Every diagonal entry equals \(\sigma ^2\). The inverse-precision theorem converts that entry into the second moment of the coordinate.
Fix a site \(x\) and let \(X\) be the line through the counting basis vector \(\delta _x\). Retaining that line, the retained Gaussian marginal is the Gaussian law on \(X\) whose precision \(S\) satisfies \(h(\delta _x,S^{-1}\delta _x)=\sigma ^2\). The coordinate on \(X\) is the scalar \(t\) in \(t\delta _x\). Its law is the scalar Gaussian of precision \(a=1/\sigma ^2\).
Write \(g_a(t)=\exp (-at^2/2)\). The product rule gives \[(t^3g_a(t))'=3t^2g_a(t)-at^4g_a(t).\] Integration on \([-R,R]\) gives \[a\int _{-R}^{R}t^4g_a(t)\,dt =3\int _{-R}^{R}t^2g_a(t)\,dt -2R^3\exp (-aR^2/2).\] For \(u\ge 0\) the exponential series gives \(e^u\ge u^2/2\). At \(u=aR^2/2\) and \(R\ge 1\) the boundary term is at most \(16/(a^2R)\), so it tends to zero. The integrands are nonnegative, and polynomial moments make their full integrals finite, so the compact integrals increase to those limits. Therefore \[a\int _{\mathbb R}t^4g_a(t)\,dt =3\int _{\mathbb R}t^2g_a(t)\,dt.\] The scalar Gaussian integral is \(I(a)\). Dividing by it and using the scalar identity \(\int t^2g_a=I(a)/a\) from the inverse-precision theorem gives the normalized moment \(3/a^2=3\sigma ^4\). The polynomial pullback in Theorem [fdip-000E] carries this moment to the coordinate \(\phi (x)\).
If \(x\ne y\), the integer \(v\) above is the largest for which the two sites are congruent modulo \(p^v\), and \(v\le k-1\) because the sites are distinct residues modulo \(p^k\). Congruence modulo \(p^{k-j}\) holds exactly when \(k-j\le v\). The matrix entry of \(E_j\) at this pair is therefore determined by \(v\). The same is true of every \(P_j\), and of \(Q\).