博客

  • 从支撑接触到 Hausdorff 分层:凸 k-Hessian 解的四步证明机制

    这篇 note 介绍论文第3–6章的主证明。 目标是说明:若 $u$ 是开凸域 $\Omega$ 上满足

    $$\sigma_k(D^2u)=1,\qquad 2\le k\le n,$$

    的凸黏性解,则其非局部 $C^2$ 集满足

    $$\mathcal H^{n-1}_{\mathrm{loc}}\bigl(\operatorname{Sing}(u)\bigr)=0.$$

    证明把支撑接触集的维数、John 椭球的短轴、严格 section 的体积衰减和 Mooney 覆盖定理连接成一条闭合链。临界接触维数是 $q=n-k+1$。

    1. 接触分层与主证明链

    对 $x\in\Omega$ 和 $p\in\partial u(x)$,令 $\ell_{x,p}$ 为相应的支撑仿射函数,并用 $d_u(x,p)$ 表示局部接触集 $\{u=\ell_{x,p}\}$ 在小尺度稳定后的仿射维数。定义

    $$C_m(u)=\left\{x:d_u(x,p)\ge m\ \text{for every }p\in\partial u(x)\right\}.$$

    “对每一个支撑斜率”是必要量词,因为第6章的 covering theorem 要求同一点的每一个 section 都满足体积衰减。

    1. 第3章:若某个支撑斜率的接触维数不超过 $n-k$,则 $u$ 在该点附近光滑。
    2. 第4章:正的 $k$-Hessian 密度控制小次水平集 John 椭球最短 $k$ 根轴的乘积。
    3. 第5章:$m$ 维接触与二次局部化共同产生 $|S_h|\lesssim h^{(m+k)/2}$。
    4. 第6章:Mooney 覆盖把 section 体积衰减转化为 Hausdorff 分层。

    $$\operatorname{Sing}(u)\subset C_q(u),\qquad q=n-k+1,\qquad \mathcal H^{2n-m-k}_{\mathrm{loc}}\bigl(C_m(u)\bigr)=0.$$

    最后取 $m=q$,便有 $2n-q-k=n-1$。

    2. 第一个引擎:小接触集推出局部光滑

    这一节中,标准逼近、内点二阶估计和 Evans–Krylov 都是熟悉的步骤;真正决定论证能否启动的是一个马鞍面。减去一个支撑仿射函数后,归一化为 $u(0)=0$、$u\ge0$。若 $\{u=0\}$ 在小球内落入至多 $n-k$ 维的子空间 $Z$,取正交分解 $\mathbb R^n=Y\oplus Z$。此时 $\dim Y\ge k$;只需在 $Y$ 中选出 $k$ 个方向,并构造

    $$Q(y,z)=A|y|^2-|z|^2+\alpha r^2.$$

    这个马鞍面的两组方向承担相反而互补的任务:$Y$ 中至少 $k$ 个强正曲率方向保证 $D^2Q\in\Gamma_k$;沿接触空间 $Z$ 的负曲率则把 $Q$ 在球面边界压到 $u$ 的下方。选择 $A$、$\alpha$ 并适当缩放后,可以同时做到“边界上 $Q<u$”与“中心附近 $Q>u$”。这正是单纯凸抛物面无法提供的几何。

    对 CNS 给出的光滑 $k$-admissible 逼近 $u_j$,选择正则值 $t_j$,并取包含原点的连通分支

    $$D_j=\{Q-u_j>t_j\}.$$

    边界与中心的符号差保证 $D_j$ 是困在小球内部的移动区域,而正则值保证其边界光滑。于是可以在 $D_j$ 上应用 Chou–Wang 局部二阶估计,得到与 $j$ 无关的 Hessian 界。随后对凹算子 $F=\sigma_k^{1/k}$ 使用 Evans–Krylov 与 Schauder,才得到极限 $u$ 的局部光滑性。

    一致收敛 $u_j\to u$ 本身不产生 $C^2$ 极限;真正的正则化步骤是马鞍面切出的变化域和 Chou–Wang Hessian 界。

    图1。 马鞍面的两个方向:接触方向中的负曲率负责封住边界,至少 $k$ 个横向正曲率负责保持 $k$-admissibility。图面可用鼠标拖动旋转,并用滚轮缩放。

    3. 第二个引擎:John 椭球与内部黏性测试

    严格 section $K_h=\{u\le h\}\cap\overline B_R$ 虽然是凸集,却可能有复杂边界、随尺度旋转,并且各方向尺度相差很大;证明并没有假定它本身是椭球。John 定理的作用,是用一个内接椭球把这种高复杂度压缩成 $n$ 个半轴和一组主方向:

    $$E_h=z_h+A_hB_1,\qquad E_h\subset K_h\subset z_h+nA_hB_1,$$

    把 $A_h$ 的半轴按 $a_1\ge\cdots\ge a_n$ 排列。于是 section 的全部凸几何复杂度,在固定维数常数 $n$ 的损失内,变成了这 $n$ 根有序轴的账本。由 $E_h$ 构造测试函数

    $$P_h(x)=h\left|A_h^{-1}(x-z_h)\right|^2.$$

    令 $\mu_h=\min(P_h-u)$。平移后的 $P_h-\mu_h$ 在某个辅助内点 $y_h$ 从上方接触 $u$。因为 $D^2P_h=2hA_h^{-2}>0$,黏性不等式给出

    $$a_q(h)\cdots a_n(h)\lesssim \lambda^{-1/2}h^{k/2},\qquad q=n-k+1.$$

    这是最短 $k$ 根 John 轴乘积的上界。方程控制的是乘积,而不是每一根轴。

    图2。 第4章内部黏性测试的几何结构。抬高平面显示 graph cap;$K_h$ 与 $E_h$ 位于自变量空间,测试点 $y_h$ 由最小化 $P_h-u$ 自动产生。

    4. 第三个引擎:接触维数与 section 体积

    这一节的核心不是把 $K_h$ 想象成某个规则形状,而是把上一节得到的有序 John 轴分成三段。记 $q=n-k+1$。测试函数与黏性方程给出最短 $k$ 根轴的乘积上界

    $$a_q\cdots a_n\lesssim \lambda^{-1/2}h^{k/2}.$$

    另一方面,若支撑接触集含有半径 $\rho$ 的 $m$ 维相对球,John 上包含和宽度的极小极大表征给出 $a_m\gtrsim\rho$。由于 $a_1\ge\cdots\ge a_n$,这实际意味着

    $$a_1,\ldots,a_m\gtrsim\rho.$$

    当 $m\ge q$ 时,两段索引在中间确实重合:$a_q,\ldots,a_m$ 同时属于“接触给下界”的区间和“PDE 控制乘积”的区间。把这些已有下界的重合因子从短轴乘积中除去,便得到

    $$a_{m+1}\cdots a_n\lesssim \lambda^{-1/2}\rho^{-(m-q+1)}h^{k/2}.$$

    这一步迫使重合区间之后的最短 $n-m$ 根轴总体快速衰减,但还留下一个漏洞:前 $m$ 根轴虽然有下界,却完全可能非常长;所以仅靠上式还不能控制整个 section 的体积。修补这个漏洞的技巧是令

    $$v(x)=u(x)+\frac{|x|^2}{2}.$$

    加不加这个光滑二次函数,奇异集完全相同:$u$ 在一点附近属于 $C^2$ 当且仅当 $v$ 属于 $C^2$。这里也不对 $v$ 使用 $k$-Hessian 方程;所有 PDE 轴乘积估计仍来自 $u$。二次项只承担几何局部化。对相应严格 section,$v(x)\le h$ 立即蕴含 $|x|^2/2\le h$,因此

    $$S_h^v\subset K_h\cap B_{\sqrt{2h}}.$$

    这只是一项光滑抬升,却把可能任意长的前 $m$ 个接触方向全部截到 $O(\sqrt h)$,贡献 $h^{m/2}$;重合消去后剩余短轴的乘积贡献 $h^{k/2}$。两部分合在一起得到

    $$|S_h^v|\lesssim \lambda^{-1/2}\rho^{-(m-q+1)}h^{(m+k)/2}.$$

    图3。 半轴按 $a_1\ge\cdots\ge a_n$ 排列。上括号覆盖 $q\le i\le n$,下括号覆盖 $1\le i\le m$,所以 $q\le i\le m$ 是两种控制的重合区间;着色的 $m+1\le i\le n$ 正是消去重合后要处理的短轴。
    图4。 Proposition 5.1 的几何顺序:接触片、$K_h$、John 轴下界以及 $S_h^v\subset K_h\cap B_{\sqrt{2h}}$。

    5. 第四个引擎:从体积衰减到 Hausdorff 分层

    对 $x\in C_m(u)$ 以及 $v$ 的每一个支撑斜率,前一节给出

    $$|S_h^v(x)|\le C_{x,p}h^{(m+k)/2}.$$

    令 $s=m+k-n$,则 $(m+k)/2=(n+s)/2$。Mooney 的 section-covering theorem 因而推出

    $$\mathcal H^{n-s}\bigl(C_m(u)\bigr)=\mathcal H^{2n-m-k}\bigl(C_m(u)\bigr)=0.$$

    覆盖定理允许常数与高度阈值依赖 $(x,p)$,所以接触球半径 $\rho$ 没有统一下界并不妨碍论证。另一方面,满维接触会产生曲率任意小的上方二次测试,与 $\sigma_k\ge\lambda$ 矛盾,因此 $C_n(u)=\varnothing$。

    第3章给出 $\operatorname{Sing}(u)\subset C_q(u)$。在 $m=q=n-k+1$ 处应用分层结论,便得到 $\mathcal H^{n-1}_{\mathrm{loc}}(\operatorname{Sing}(u))=0$。

    6. 模型核算:$n=4$、$k=3$

    此时 $q=2$。对临界层 $m=2$,接触几何给出 $a_2\gtrsim\rho$,而 PDE 给出 $a_2a_3a_4\lesssim h^{3/2}$。加入二次项后,

    $$|S_h|\lesssim h^{5/2},\qquad \mathcal H^3(C_2)=0.$$

    对更深的 $m=3$ 层,

    $$|S_h|\lesssim h^3,\qquad \mathcal H^2(C_3)=0.$$

    图5。 $n=4$、$k=3$ 时,$m=2$ 与 $m=3$ 两层的 John 轴乘积、section 体积指数和 Hausdorff 指数。

    7. 端点与限制

    • $k=n$:$q=1$,结论恢复 Monge–Ampère 情形的 $\mathcal H^{n-m}(C_m)=0$。
    • $k=2$:本方法给出临界接触层的余维一零测;全光滑还需要 strict $2$-convexity 的额外输入。
    • 轴的控制:方程只控制短轴乘积,主轴可以随尺度旋转并重新分配,所以当前结论是 Hausdorff nullity,而不是统一 Minkowski packing。

    关于本文。 这是一篇介绍 Xiyu Hu, Sharp Hausdorff Bounds for the Interior Singular Set of Convex k-Hessian Solutions 第3–6章的证明导读 note。论文第8章的 mean-value 公式提供概念动机,但不作为上述主证明链的输入。


  • From Sum-Free Sets to Strongly 2-Primitive Sets: Localization, Patching, and the Search for Stability

    Two extremal problems can share a proof architecture without sharing a proof. Benjamin Bedert’s 2025 breakthrough on large sum-free subsets is additive and Fourier-analytic. A July 2026 manuscript on strongly $2$-primitive sets is multiplicative and hypergraph-theoretic. Their common language is not a transferable sieve, but a four-step design: localize, build a strong object in each local block, prevent leakage between blocks, and add the gains. That comparison points toward the right stability question – and also reveals why the most naive version of stability is false.

    Status and scope. Bedert’s theorem is the February 2025 arXiv preprint cited below. The $27/2$ theorem is contained in a five-page manuscript by Przemek Chojecki supplied for this post in July 2026; no public preprint or peer-reviewed version was found, and the Erdős Problems page for #793 still labels the problem open at the time of writing. The written proof was checked here for internal coherence, but that is not independent peer review. The factor-graph and normalized-stability results later in this article are the additional analysis developed in the accompanying long note. Conditional and open statements are labelled explicitly.

    Here is the headline before the details. The supplied manuscript claims

    $$
    F(n)=\pi(n)+\left(\frac{27}{2}+o(1)\right)
    \frac{n^{2/3}}{(\log n)^2},
    $$

    where $F(n)$ is the largest size of a strongly $2$-primitive subset of $[1,n]$. One might expect every near-extremizer to be close to the construction behind this formula: almost all primes, plus products of three primes near $n^{1/3}$ arranged as a nearly saturated linear $3$-uniform hypergraph. That literal statement is false. There is, however, an exact and useful replacement: after contracting harmless private-prime lifts, every near-extremizer has a canonical prime layer and a small exceptional hub carrying the entire second-order excess. If the residual hub is mostly made of squarefree triples, then pair saturation and $n^{1/3}$-scale stability follow. Proving – or refuting – that triple-core hypothesis is the remaining inverse problem.

    1. Two problems that look similar only from far away

    1.1 The sum-free extraction problem

    A set $B\subset\mathbb Z$ is sum-free if it contains no $x,y,z$, not necessarily distinct, with

    $$x+y=z.$$

    For a finite set $A$ of integers, write

    $$
    S(A)=\max\{|B|:B\subseteq A,\ B\text{ is sum-free}\},
    $$

    and, for $N$-element sets of positive integers,

    $$
    S(N)=\min_{|A|=N}S(A).
    $$

    This is an extraction problem. An adversary gives us an arbitrary host set $A$, and we must find a large structured subset inside it. Erdős’s middle-third argument gives $S(A)\ge |A|/3$. Alon and Kleitman improved this to $(N+1)/3$, and Jean Bourgain proved $S(N)\ge (N+2)/3$ in 1997. The long-standing qualitative question was whether the additive improvement can tend to infinity:

    $$S(N)\ge \frac N3+\omega(N),\qquad \omega(N)\longrightarrow\infty.$$

    Bedert answered yes, proving that an absolute $c\gt0$ exists such that

    $$
    S(A)\ge \frac{|A|}{3}+c\log\log |A|.
    $$

    1.2 The strongly $2$-primitive packing problem

    A set $A\subseteq[1,n]$ is strongly $2$-primitive when

    $$
    a\nmid bc
    \qquad
    (a,b,c\in A,\ a\notin\{b,c\}),
    $$

    where $b=c$ is allowed. The word strongly matters. Under a more recent convention, a $2$-primitive set only forbids witnesses $b,c$ that are distinct. For example, $\{4,5,6\}$ passes that weaker test but fails the strong one because $4\mid6^2$.

    Now define

    $$
    F(n)=\max\{|A|:A\subseteq[1,n],\ A\text{ is strongly }2\text{-primitive}\}.
    $$

    This is a packing problem. The host interval is fixed, and we directly construct the largest possible forbidden-divisibility family. The all-primes set gives the leading term $\pi(n)$. In 1938 Erdős proved upper and lower bounds of the form

    $$
    \pi(n)+c_1\frac{n^{2/3}}{(\log n)^2}
    \le F(n)\le
    \pi(n)+c_2\frac{n^{2/3}}{(\log n)^2},
    $$

    and later asked whether the second-order term has an asymptotic constant. This is the modern Erdős Problem #793. The supplied 2026 manuscript proposes that the constant is $27/2$; it is important not to say that Erdős himself conjectured this numerical value.

    Comparison of the sum-free extraction problem and the strongly 2-primitive packing problem
    Figure 1. The quantifiers already separate the two questions. Sum-free theory extracts a subset from an arbitrary host; the multiplicative problem packs a forbidden-divisibility family into a fixed interval.

    1.3 Equality versus order

    The relation $a\nmid bc$ is not ordinary product-freeness. The set $\{6,10,15\}$, for instance, has no internal equality $xy=z$, yet $6\mid10\cdot15$. Prime valuations expose the real geometry:

    $$
    a\mid bc
    \quad\Longleftrightarrow\quad
    v_p(a)\le v_p(b)+v_p(c)
    \quad\text{for every prime }p.
    $$

    Thus the forbidden relation is a coordinatewise domination inequality in the divisor lattice. For squarefree integers, if $E_a=\{p:p\mid a\}$, it becomes

    $$E_a\subseteq E_b\cup E_c.$$

    So the natural combinatorial object is a $2$-cover-free family, not the solution set of a linear equation. Fourier characters are superb at detecting equations such as $x+y=z$. They do not come with an evident contractive projection that detects the partial order $\nu(a)\le\nu(b)+\nu(c)$. This is the first reason Bedert’s proof cannot simply be copied into the multiplicative setting.

    2. The $27/2$ upper bound: every element needs a private factor

    Set

    $$
    y=n^{1/3},\qquad
    M=\frac{y}{\log n},\qquad
    \Sigma_n=M^2=\frac{n^{2/3}}{(\log n)^2}.
    $$

    The upper bound begins with a small lemma that contains more information than the inequality it proves.

    Private-factor lemma. Let $\mathcal B$ be a set of positive integers, and choose for every $a\in A$ a factorization $a=u_av_a$ with $u_a,v_a\in\mathcal B$. If $A$ is strongly $2$-primitive, then $|A|\le|\mathcal B|$.

    Think of $E_a=\{u_a,v_a\}$ as a two-element multiset. If $u_a\ne v_a$ and neither coordinate is private to $a$, another chosen pair contains $u_a$ and another contains $v_a$; the corresponding two elements have product divisible by $a$. If $u_a=v_a=x$, failure of privacy would give another pair containing $x$ twice, hence another element equal to $x^2=a$. Therefore every $a$ has a private coordinate, and these coordinates are distinct. That is the injection $A\hookrightarrow\mathcal B$.

    The manuscript chooses the multiplicative $2$-basis

    $$
    \begin{aligned}
    \mathcal B_0&=[1,n^{3/5}],\\
    \mathcal B_1&=\{p\text{ prime}:n^{3/5}\lt p\le n\},\\
    \mathcal B_2&=\{pq:p,q\le y\text{ prime}\},\\
    \mathcal B_3&=\{qr:y\lt q\le n^{2/5},\ r\le n/q^2,\ q,r\text{ prime}\}.
    \end{aligned}
    $$

    Every $m\le n$ can be factored into two elements of $\mathcal B=\mathcal B_0\cup\mathcal B_1\cup\mathcal B_2\cup\mathcal B_3$. The case split is elementary but carefully tuned. Small $m$ can be balanced into two factors below $n^{3/5}$, a prime factor above $n^{2/5}$ can be split off, and the remaining difficult case groups two of the three largest prime factors into an element of $\mathcal B_2$ or $\mathcal B_3$.

    The main prime layer comes from $\mathcal B_1$. The genuinely second-order counts are

    $$
    |\mathcal B_2|
    =\binom{\pi(y)+1}{2}
    =\left(\frac92+o(1)\right)\Sigma_n,
    $$

    because $\pi(n^{1/3})\sim3n^{1/3}/\log n=3M$, and

    $$
    |\mathcal B_3|
    =\sum_{y\lt q\le n^{2/5}}\pi\!\left(\frac n{q^2}\right)
    =(9+o(1))\Sigma_n.
    $$

    The contribution of $\mathcal B_0$ and all relevant overlaps is $o(\Sigma_n)$. Hence

    $$
    |\mathcal B|
    =\pi(n)+\left(\frac92+9+o(1)\right)\Sigma_n
    =\pi(n)+\left(\frac{27}{2}+o(1)\right)\Sigma_n.
    $$

    The private-factor injection then gives the upper bound. A caution that becomes crucial for stability: $\mathcal B_2$ and $\mathcal B_3$ are labels in an upper-bound certificate. Near-saturation of these labels does not immediately say that the original elements of $A$ are themselves products of two or three primes.

    3. The lower bound: turn properly coloured edges into prime triples

    The lower bound lives in a linear $3$-uniform hypergraph. Let $\mathcal H$ be a family of triples of distinct primes such that any two triples share at most one prime and every edge product is at most $n$. Define

    $$
    A_{\mathcal H}
    =\{p\le n:p\text{ prime and }p\notin V(\mathcal H)\}
    \cup
    \left\{\prod_{p\in E}p:E\in\mathcal H\right\}.
    $$

    Then

    $$|A_{\mathcal H}|=\pi(n)-|V(\mathcal H)|+|\mathcal H|.$$

    Why is this strongly $2$-primitive? A target triple product has three distinct prime coordinates. Each other hyperedge supplies at most one of them, so two other elements supply at most two. The same argument still works when the two witnesses coincide. Primes outside the vertex set remain singleton elements and cannot divide any other chosen element.

    3.1 Logarithmic prime bins

    Fix a small mesh $h\gt0$ and divide primes near $y=n^{1/3}$ into bins

    $$
    P_i=\{p\text{ prime}:ye^{ih}\lt p\le ye^{(i+1)h}\},
    \qquad
    \Delta_i=e^{(i+1)h}-e^{ih}.
    $$

    For a cell satisfying

    $$i\le j,\qquad i+2j\le-4,$$

    put $k=-i-j-3$. Then $i\le j\lt k$, and any $p\in P_i$, $q\in P_j$, $r\in P_k$ obeys $pqr\le n$.

    If $i\lt j$, properly edge-colour the complete bipartite graph between $P_i$ and $P_j$, using primes of $P_k$ as colours. When $i=j$, do the same with the complete graph on $P_i$. The lower pair $\{p,q\}$, coloured by $r$, becomes the hyperedge $\{p,q,r\}$.

    A complete bipartite graph between two prime bins properly edge-coloured by a third prime bin
    Figure 2. Within a cell, every colour class is a matching, so two produced triples never share a lower pair. Across cells, the sorted signature $(i,j,k)$ has fixed sum $-3$; two shared bin indices force the third and hence force the same cell.

    The proper colouring is the local no-interference mechanism. The constant-sum signature is the global one. Together they make the union over all cells a linear hypergraph.

    3.2 The cell weight

    For any fixed finite collection of bins, the prime number theorem gives

    $$|P_i|=(3+o(1))M\Delta_i.$$

    Off the diagonal, a cell contributes approximately $9M^2\Delta_i\Delta_j$ edges. A diagonal cell contributes half as many unordered pairs. The exact geometric-series identity is

    $$
    \sum_{\substack{i\lt j\\i+2j\le-4}}\Delta_i\Delta_j
    +\frac12\sum_{i\le-2}\Delta_i^2
    =e^{-h}+\frac12e^{-2h}.
    $$

    Letting $h\to0$, the weight tends to $3/2$, so

    $$
    |\mathcal H|
    =\left(9\cdot\frac32-o(1)\right)\Sigma_n
    =\left(\frac{27}{2}-o(1)\right)\Sigma_n.
    $$

    The number of vertices used is only $o(\Sigma_n)$. Replacing those primes by the triple products therefore yields the matching lower bound.

    The logarithmic feasible pair region split into the 9 over 2 and 9 counting regimes
    Figure 3. The same constant is visible in the canonical feasible pair space. The labels $9/2$ and $9$ are prime-counting contributions, not Euclidean areas of the drawing.

    4. From Erdős’s middle third to Bedert’s $c\log\log N$

    We now return to the additive problem. Let $\mathbb T=\mathbb R/\mathbb Z$, and let $\phi=\mathbf1_{(1/3,2/3)}$. The middle third of the circle is sum-free: two points in it cannot add, modulo $1$, to another point in it. Therefore, for every $x\in\mathbb T$,

    $$A_x=\{a\in A:ax\pmod1\in(1/3,2/3)\}$$

    is sum-free. Averaging $|A_x|$ over $x$ gives $|A|/3$. The problem is to force a positive fluctuation above that mean.

    After a harmless normalization, the relevant Fourier series has the form

    $$
    F_A(x)=\sum_{a\in A}\sum_{m\ge1}
    \frac{\chi(m)}m\cos(2\pi max),
    $$

    where $\chi$ is the non-principal real character modulo $3$. A lower bound for $\|F_A\|_1$ gives a one-sided large value and hence a large sum-free subset.

    4.1 Where Bourgain’s sieve loses a logarithm

    The historical correction is worth making explicit. The relevant predecessor is Jean Bourgain’s 1997 paper, not “John Bogan, 1993.” Also, the Littlewood $L^1$ conjecture had already been proved in 1981, independently by McGehee-Pigno-Smith and Konyagin. Bourgain’s argument combines Fourier analysis with a Möbius operation that filters harmonic indices.

    Schematically, if $R_Q$ denotes the positive $Q$-rough integers, then

    $$
    \sum_{d\mid\prod_{p\le Q}p}
    \frac{\mu(d)\chi(d)}d F_A(dx)
    =\sum_{a\in A}\cos(2\pi ax)+\mathcal R_Q(x).
    $$

    The left side costs

    $$
    \prod_{p\le Q}\left(1+\frac1p\right)\asymp\log Q
    $$

    under the triangle inequality. If $Q$ is large enough to make the rough remainder negligible directly, this factor consumes the logarithm delivered by a Littlewood-type lower bound. This explains why merely “using the Littlewood theorem harder” does not produce the desired unbounded gain.

    4.2 Bedert’s split: medium primes are sieved, small primes are projected

    Bedert chooses

    $$Q_1=(\log N)^{1/2},\qquad Q=(\log N)^{20}.$$

    He Möbius-sieves only the medium primes $Q_1\le p\le Q$. Their reciprocal sum is $O(1)$, so the norm loss is a constant rather than a logarithm. The small primes $p\le Q_1$ are handled by a different operation.

    For every small prime, record the exact valuation $\nu_p(a)$ and the unit residue after removing that prime power. Chinese remaindering packages all this data into a joint fibre $A(r,\nu)$. Fourier projection to a residue class,

    $$
    \operatorname{Proj}(H;\rho\bmod q)(x)
    =\sum_{m\equiv\rho\ (q)}\widehat H(m)e(mx),
    $$

    is an $L^1$-contraction:

    $$
    \|\operatorname{Proj}(H;\rho\bmod q)\|_1\le\|H\|_1.
    $$

    At the main valuation level, the projection removes small-prime harmonic contamination without paying the full Möbius product. After truncation, one sees a large exponential sum plus an $L^2$-small error, to which a robust McGehee-Pigno-Smith test function applies.

    At a general valuation level, lower fibres can still leak into the projected channel. Bedert proves an isolation-or-descent alternative: either the desired $L^1$ lower bound is already present, or a large fibre has a strictly lower fibre losing at most a controlled polylogarithmic factor. Iterating and reversing this descent produces a chain

    $$
    \nu^{(1)}\prec\nu^{(2)}\prec\cdots\prec\nu^{(J)},
    \qquad
    J\gg\frac{\log N}{\log\log N},
    $$

    whose fibre sizes grow geometrically.

    4.3 The non-Archimedean no-leakage lemma

    Associate to these levels the nested moduli

    $$
    q_i=\prod_{p\le Q_1}p^{\nu_p^{(i)}+1},
    \qquad q_1\mid q_2\mid\cdots\mid q_J.
    $$

    Let $g_i$ be the normalized indicator of the $i$-th residue block and put $Q_i(x)=\exp(-|\widehat g_i(x)|)$. The non-Archimedean MPS test function is assembled from

    $$
    \Phi_J=
    \widehat g_J+\widehat g_{J-1}Q_J+\cdots+
    \widehat g_1Q_2\cdots Q_J.
    $$

    The decisive observation is that $|\widehat g_i|$ is $1/q_i$-periodic, hence $\widehat Q_i$ is supported on $q_i\mathbb Z$. Multiplying an earlier block by a later $Q_k$ shifts its frequencies only by multiples of $q_k$. Since $q_i\mid q_k$, the earlier block cannot escape its residue class modulo $q_i$. Main inner products contribute one controlled unit per block, while geometric fibre growth makes the cross terms summable.

    That is Bedert’s genuine no-leakage mechanism. It is more precise than the slogan “use many congruence classes and patch them.”

    One further qualification prevents a common misstatement. Bedert obtains a suitable $F_4$-isomorphic model $B$ of the original set and proves the required $L^1$ bound there. The norm $\|F_A\|_1$ is not claimed to be invariant under an $F_4$-isomorphism. What is preserved is the four-term additive information needed for sum-freeness, so $S(A)=S(B)$.

    Bedert’s inverse output is also stronger than the numerical lower bound. If $S(A)\le N/3+C$, his results force low additive dimension, a dense $F_4$-model, large additive energy in every substantial subset, and a “99% Structure Theorem” decomposing all but a small exceptional set into large small-doubling pieces. This is a genuine inverse theorem for host sets resisting sum-free extraction, but not an edit-distance classification by one canonical extremizer.

    5. The real bridge: localize, build, prevent leakage, sum

    Parallel flowcharts for Bedert's p-adic Fourier proof and the logarithmic hypergraph construction
    Figure 4. The analogy is architectural. Bedert preserves Fourier residue lanes through nested moduli; the multiplicative construction preserves low codegrees through matching colours and unique scale signatures.

    The dictionary is now clean:

    Role Bedert’s sum-free proof Strongly $2$-primitive proof
    Underlying relation Linear equation in an abelian group Coordinatewise domination of prime valuations
    Local coordinates Exact small-prime valuations and unit residues Logarithmic sizes of prime factors
    Local block Joint $p$-adic fibre and MPS block Scale cell and properly coloured graph
    No interference Nested moduli preserve Fourier support Colour matchings and unique cell signatures preserve linearity
    Accumulation One controlled inner product per fibre One triple per admissible lower pair
    Inverse output Additive dimension, energy, small-doubling pieces Private factors and a cover-free exponent core

    This is a substantial connection, but it is a connection of proof design. Bedert’s Fourier projection has no direct analogue for the divisor partial order. The most plausible transfer is therefore not a line-by-line proof but a research program: find a multiplicative localization, an isolation-or-descent alternative, and a no-leakage invariant that survives repeated prime powers and larger supports.

    6. The first stability guess is false

    The lower construction suggests a tempting statement: every set with

    $$
    |A|\ge
    \pi(n)+\left(\frac{27}{2}-o(1)\right)\Sigma_n
    $$

    should differ in only $o(\Sigma_n)$ places from almost all primes plus a near-optimal linear family of triple products. Two examples show why this is too rigid.

    6.1 A private-prime lift

    Start with a near-optimal linear triple family $\mathcal H_n$ using odd primes below $n/2$. Keep primes above $n/2$, replace every unused prime $p\le n/2$ by $2p$, and retain the triple products. Each $2p$ has the private prime divisor $p$; no other chosen element contains that prime. The triple products are still protected by linearity. The new set remains strongly $2$-primitive and has the same second-order asymptotic size.

    But essentially every prime below $n/2$ has been replaced by a composite. The symmetric-difference distance from the all-primes model is

    $$
    \Theta(\pi(n/2))=\Theta\!\left(\frac n{\log n}\right),
    $$

    and

    $$
    \frac{\pi(n/2)}{\Sigma_n}
    \asymp n^{1/3}\log n\longrightarrow\infty.
    $$

    So literal edit-distance stability fails by much more than the scale of the second-order term.

    6.2 A nonlinear sunflower at almost no cost

    Even if we insist on squarefree triple products, linearity need not hold in the original representation. Add a sunflower

    $$\{\{2,3,r\}:r\in R\}$$

    for many private petals $r$. Remove the singleton primes $2,3,r$ and add the products $6r$. The net cardinality loss is only $2$, while the support family contains $\asymp n/\log n$ edges sharing the pair $\{2,3\}$. It is $2$-cover-free because every target has its own private petal, but any linear subfamily contains at most one sunflower edge.

    The two examples have the same cause: degree-one prime coordinates carry a huge amount of neutral decoration. Stability can only become true after quotienting out this freedom.

    7. Near equality in the upper bound is an exact star forest

    Fix the chosen factorizations $a=u_av_a$ from the multiplicative basis proof. Make a graph $G_A$ with vertex set $\mathcal B$, one edge $\{u_a,v_a\}$ for each $a\in A$, and allow loops.

    Exact factor-graph theorem. Every nonloop edge has a degree-one endpoint, and every loop is an isolated component. Thus $G_A$ is a disjoint union of stars and isolated loops. If $U$ is the number of unused vertices and $c(G_A)$ is the number of nonloop star components, then

    $$|\mathcal B|-|A|=U+c(G_A).$$

    The proof is the private-factor lemma read without discarding information. If a nonloop edge $\{u,v\}$ had both endpoints in other edges, the corresponding two elements would have product divisible by $uv=a$. If a loop at $u$ met another edge, then $u^2$ would divide the square of the other element.

    A factor graph decomposed into stars, an isolated loop and unused vertices
    Figure 5. The upper-bound defect is counted exactly: an unused basis coordinate costs one, and each nonloop star component costs one.

    More quantitatively,

    $$
    \bigl|\{b\in\mathcal B\setminus\mathcal B_0:
    \deg_{G_A}(b)\ne1\}\bigr|
    \le |\mathcal B|-|A|.
    $$

    If

    $$
    |A|\ge
    \pi(n)+\left(\frac{27}{2}-\delta_n\right)\Sigma_n,
    $$

    then the defect $D_n=|\mathcal B|-|A|$ is at most $(\delta_n+o(1))\Sigma_n$. Hence almost every secondary coordinate in both $\mathcal B_2$ and $\mathcal B_3$ occurs in exactly one chosen factor pair. This is a rigorous stability theorem for the upper-bound certificate. It is not yet a classification of the original integers.

    8. Private-prime compression gives the correct normal form

    Suppose a prime $p$ divides $a\in A$ and divides no other element of $A$. Replacing $a$ by $p$ preserves the cardinality and strong $2$-primitivity. The new target $p$ cannot divide a product of two other elements because neither contains $p$; for any unchanged target, replacing $a$ by a divisor only decreases the product available to cover it.

    For primes $p\gt n^{3/5}$, the basis factorization can be chosen to expose $p$ in every multiple. Thus a degree-one large-prime coordinate is genuinely a globally private prime and can be contracted. Contract all such coordinates simultaneously.

    Normalized stability modulo private-prime lifts. From every strongly $2$-primitive $A\subseteq[1,n]$ satisfying the near-extremal bound above, private-prime contractions produce a strongly $2$-primitive set $\widetilde A$ of the same size with

    $$\widetilde A=(\mathbb P\cap[1,n]\setminus V)\cup C,\qquad \operatorname{supp}(c)\subseteq V\quad(c\in C).$$

    where

    $$|V|\le(\delta_n+o(1))\Sigma_n,\qquad |C|-|V|=|A|-\pi(n).$$

    In particular, if $\delta_n=o(1)$, then

    $$
    |V|=o(\Sigma_n),
    \qquad
    |C|=\left(\frac{27}{2}+o(1)\right)\Sigma_n.
    $$

    The geometry is now transparent. Before compression, the leading $\pi(n)$ coordinates may be represented by arbitrary private multiples. After compression, almost every prime appears literally. Every residual composite is supported entirely on the small exceptional prime hub $V$, and its excess over the missing primes is exactly the second-order gain.

    Private-prime lifts contracted to a normal form with singleton primes and a small exceptional hub
    Figure 6. Literal stability fails on the left. The quotient by private-prime contractions produces the exact normal form on the right.

    There is one more unconditional consequence. At most $|V|$ members of $C$ have support of size at most two. Indeed, any one- or two-prime residual element must have a valuation coordinate in which it is a strict global maximum; assign the element to such a prime. Two elements cannot receive the same prime. Therefore, in a normalized near-extremizer, all but $o(\Sigma_n)$ residual composites have at least three distinct prime divisors.

    This is close to the hoped-for triple picture, but it does not prove that the typical element has exactly three prime factors, that it is squarefree, or that its prime factors lie near $n^{1/3}$.

    9. Conditional stability when the residual core is made of triples

    Now impose an additional hypothesis:

    Squarefree-triple hypothesis. All but $o(\Sigma_n)$ elements of $C$ are products of three distinct primes.

    Let $\mathcal H$ be the support family of those triples. Strong $2$-primitivity says precisely that it is $2$-cover-free:

    $$
    E\nsubseteq F\cup G
    \qquad(E,F,G\in\mathcal H,\ E\notin\{F,G\}),
    $$

    with $F=G$ allowed.

    9.1 Pruning to a linear core

    If two triples $E,F$ share a pair and $x$ is the third vertex of $E$, then $x$ has degree one. Otherwise another edge $G$ containing $x$, together with $F$, would cover $E$. Delete such a petal edge whenever a repeated pair occurs. Every deletion removes at least one vertex as well as one edge, so the excess $|\mathcal H|-|V(\mathcal H)|$ does not decrease.

    The result is a linear subfamily $\mathcal L$ with

    $$
    |\mathcal L|-|V(\mathcal L)|
    \ge |\mathcal H|-|V(\mathcal H)|,
    $$

    and only $o(\Sigma_n)$ edges are lost in the normalized near-extremal setting. Notice the nuance: linearity is recovered after pruning private petals; it is not asserted for every cover-free representation.

    9.2 Saturating the feasible pair space

    Define

    $$
    \mathcal D_n=
    \{(p,q):p\lt q\text{ prime and }pq^2\le n\}.
    $$

    Sort an edge of $\mathcal L$ as $p\lt q\lt r$ and map it to $(p,q)$. Linearity makes this map injective. Since $r\gt q$ and $pqr\le n$, the image lies in $\mathcal D_n$. Direct counting gives

    $$
    |\mathcal D_n|
    =\left(\frac{27}{2}+o(1)\right)\Sigma_n.
    $$

    The range $q\le n^{1/3}$ contributes $(9/2+o(1))\Sigma_n$; the range $n^{1/3}\lt q\le n^{2/5}$ contributes $(9+o(1))\Sigma_n$; the tail is negligible. Since the linear core already has $(27/2-o(1))\Sigma_n$ edges, its lower-pair map misses only $o(\Sigma_n)$ feasible pairs.

    Moreover, for every fixed $\varepsilon\gt0$, all but $o(\Sigma_n)$ edges satisfy

    $$
    n^{1/3-\varepsilon}
    \le p,q,r\le
    n^{1/3+\varepsilon}.
    $$

    This is the desired scale stability. It does not imply uniqueness. Different one-factorizations or different proper edge-colourings can change $\Theta(\Sigma_n)$ triples while preserving the same pair occupancy and extremal count. The stable object is the saturated pair space, not an individual colouring.

    10. The remaining inverse problem is weighted and cover-free

    After normalization, write every $c\in C$ as an exponent vector

    $$
    \nu(c)=(v_p(c))_{p\in V}\in\mathbb Z_{\ge0}^{V}.
    $$

    Strong $2$-primitivity becomes the weighted cover-free condition

    $$
    \nu(c)\nleq\nu(c_1)+\nu(c_2)
    \quad\text{coordinatewise}
    $$

    whenever $c\notin\{c_1,c_2\}$. We know that

    $$
    |C|=\left(\frac{27}{2}+o(1)\right)\Sigma_n,
    \qquad |V|=o(\Sigma_n),
    $$

    and that almost every vector has support at least three. What remains unknown is the sharper assertion

    $$
    \#\{c\in C:c\text{ is not a squarefree product of three primes}\}
    =o(\Sigma_n).
    $$

    This gap may conceal genuine alternative near-extremizers. If a triple $pqr$ has product slack, one can contemplate replacing it by $2pqr$ while deleting the singleton prime $2$. For a linear support family, the underlying three private coordinates still prevent coverage. It is not known whether a positive density of such bounded-hub multiplier layers can coexist globally near the $27/2$ threshold, or whether cross-layer collisions force them to be negligible.

    Open normalized stability problem. Classify weighted $2$-cover-free exponent families under the product constraint $c\le n$ whose excess is $(27/2-o(1))\Sigma_n$. Decide whether the residual core is predominantly squarefree of support three, or whether bounded-hub multiplier layers produce genuinely different normalized near-extremizers.

    10.1 What a Bedert-style descent would need

    The additive proof suggests three design requirements, not three ready-made lemmas.

    1. A local order projection. One needs to isolate a valuation or support layer while controlling contamination from lower exponent vectors. A divisor-lattice zeta or Möbius transform is a possible language, but no analogue of Bedert’s $L^1$-contractive residue projection is presently available.
    2. An isolation-or-descent alternative. If the secondary basis slots do not have predominantly prime cofactors, the argument should descend to a structured hub of repeated factors with a quantitative gain or a summable loss.
    3. A no-leakage invariant. Low codegree and scale signatures work for squarefree triples. A general invariant must survive repeated exponents and supports of size at least four.

    The star-forest theorem is already a first inverse statement of this kind: near equality forces almost every basis coordinate to be a leaf attached to a small collection of hubs. The hard step is to turn that certificate-level hub structure into a classification of the original weighted exponent vectors.

    11. What is proved, conditional, and open

    • Published/preprint additive result: Bedert proves $S(A)\ge |A|/3+c\log\log|A|$, together with inverse information involving additive dimension, dense $F_4$-models, energy and a $99\%$ decomposition into large small-doubling pieces.
    • Supplied 2026 manuscript: the claimed $27/2$ asymptotic follows from the multiplicative-basis upper bound and the logarithmic-cell hypergraph construction described above. The status caveat at the beginning remains in force.
    • Unconditional stability developed in the accompanying note: literal edit stability is false; the upper-bound factor graph is an exact star forest; private-prime compression gives the normal form $(\mathbb P\setminus V)\cup C$; and almost every residual composite has at least three distinct prime divisors.
    • Conditional stability: if almost all residual composites are squarefree triples, pruning gives a near-maximal linear core whose lower pairs saturate $\mathcal D_n$, and almost every prime factor lies at scale $n^{1/3+o(1)}$.
    • Open: prove that the normalized residual core is predominantly squarefree of support three, or construct a different near-extremal weighted cover-free core.

    The conceptual moral is simple but useful. Bedert’s $p$-adic fibres and the multiplicative logarithmic cells are not the same mathematical object. Yet both proofs win by choosing coordinates in which local constructions can be made strong and then finding an exact invariant that prevents those constructions from interfering. In the stability problem, the same philosophy says to quotient the neutral directions first. Once private-prime lifts are removed, the true obstruction becomes visible: a weighted cover-free family on a tiny prime hub. That is the right object for the next theorem.

    References and source trail

    1. P. Erdős, On sequences of integers no one of which divides the product of two others and on some related problems, 1938.
    2. P. Erdős, On some applications of graph theory to number-theoretic problems, 1969; see also Erdős Problem #793.
    3. O. Carruth McGehee, L. Pigno and B. Smith, Hardy’s inequality and the $L^1$-norm of exponential sums, Annals of Mathematics 113 (1981), 613-618.
    4. J. Bourgain, Estimates related to sumfree subsets of sets of integers, Israel Journal of Mathematics 97 (1997), 71-92.
    5. B. Bedert, Large sum-free subsets of sets of integers via $L^1$-estimates for trigonometric series, arXiv:2502.08624v1, 12 February 2025.
    6. P. Chojecki, The Second Term for Strongly 2-Primitive Sets, user-supplied five-page manuscript, July 2026; no public identifier located at the time of writing.
  • Ramanujan-Nagell theorem:平方数与 2^n 相差 7 的有限性

    旧博客原文

    原题:The Ramanujan-Nagell Theorem: Understanding the Proof

    The Ramanujan-Nagell Theorem: Understanding the Proof


    补充说明

    以下是新整理的中文说明;上方旧博客原文保持不变。

    Ramanujan-Nagell theorem 研究的是一个看起来非常小的指数丢番图方程:

    $$x^2+7=2^n.$$

    它的结论是整数解只有有限个,而且正整数解恰好为

    $$(x,n)=(1,3),(3,4),(5,5),(11,7),(181,15).$$

    Ramanujan-Nagell theorem:平方数与 2^n 相差 7 的有限性
    Ramanujan-Nagell 方程把初等同余、二次域分解和 Lucas 序列约束压缩在同一个指数丢番图问题中。

    1. 初等筛选

    先看奇偶性。若 $x$ 为偶数,则左边 $x^2+7$ 为奇数,不可能等于 $2^n$;所以 $x$ 必为奇数。再看模 $8$,奇数平方恒为 $1$,所以 $x^2+7\equiv0\pmod 8$,这只说明 $n\ge3$,但已经排除了很多无意义情形。

    2. 代数数论中的分解

    真正有力的观察是把方程写成

    $$(x+\sqrt{-7})(x-\sqrt{-7})=2^n.$$

    在 $\mathbb Q(\sqrt{-7})$ 的整数环中,$2$ 可以分解成两个共轭因子。由于这个二次域的类数很小,理想层面的分解可以被提升为元素层面的约束,于是 $x+\sqrt{-7}$ 必须接近某个基本元素的 $n$ 次幂。

    3. Lucas 序列的出现

    把共轭相减,得到的不是任意等式,而是一个 Lucas 型序列项必须等于很小的数:

    $$\frac{\alpha^n-\bar\alpha^n}{\alpha-\bar\alpha}=\pm 1\quad\text{or}\quad \pm 7.$$

    这种递推序列增长很快,同时在模意义下有强限制。少数小 $n$ 需要直接检查,大 $n$ 则被递推结构和同余条件排除。

    4. 为什么这条定理有代表性

    这类问题的典型形状是:初等同余给出粗过滤,二次域分解给出结构,最后用 Lucas 序列或线性形式估计把无限可能压成有限检查。Ramanujan-Nagell 方程的漂亮之处在于,所有这些工具都集中在一个非常短的公式里。

  • BSD 猜想入门:椭圆曲线、L 函数与 Mordell-Weil rank

    旧博客原文

    原题:A crash introduction to BSD conjecture

    The pdf version is A crash introduction to BSD conjecture .

    We begin with the Weierstrass form of elliptic equation, i.e. look it as an embedding cubic curve in {\mathop{\mathbb P}^2}.

    Definition 1 (Weierstrass form) {E \hookrightarrow \mathop{\mathbb P}^2 }, In general the form is given by,

    \displaystyle E: y^2+a_1xy+a_3y=x^3+a_2x^2+a_4x+a_6 \ \ \ \ \ (1)

    If {char F \neq 2,3}, then, we have a much more simper form,

    \displaystyle y^2=x^3+ax+b, \Delta:=4a^3+27b^2\neq 0. \ \ \ \ \ (2)

    Remark 1

    \displaystyle \Delta(E)=\prod_{1\leq i,\neq j\leq 3}(z_i-z_j)

    Where {z_i^3+az_i+b=0, \forall 1\leq i\leq 3}.

    We have two way to classify the elliptic curve {E} living in a fix field {F}. \paragraph{j-invariant} The first one is by the isomorphism in {\bar F}. i.e. we say two elliptic curves {E_1,E_2} is equivalent iff

    \displaystyle \exists \rho:\bar F\rightarrow \bar F

    is a isomorphism such that {\rho(E_1)=E_2}.

    Definition 2 (j-invariant) For a elliptic curve {E}, we have a j-invariant of {E}, given by,

    \displaystyle j(E)=1728\frac{4a^3}{4a^3+27b^2} \ \ \ \ \ (3)

    Why j-invariant is important, because j-invariant is the invariant depend the equivalent class of {E} under the classify of isomorphism induce by {\bar F}. But in one equivalent class, there also exist a structure, called twist.

    Definition 3 (Twist) For a elliptic curve {E:y^2=x^3+ax+b}, all elliptic curve twist with {E} is given by,

    \displaystyle E^{(d)}:y^2=x^3+ad^2x+bd^3 \ \ \ \ \ (4)

    So the twist of a given elliptic curve {E} is given by:

    \displaystyle H^1(Gal(\bar F/ F), Aut(E_{\bar F})) \ \ \ \ \ (5)

    Remark 2 Of course a elliptic curve {E:y^2=x^3+ax+b} is the same as {E:y^2=x^3+ad^2x+bd^4}, induce by the map {\mathop{\mathbb P}^1\rightarrow \mathop{\mathbb P}^1, (x,y,1)\rightarrow (x,dy,1)}.

    But this moduli space induce by the isomorphism of {F} is not good, morally speaking is because of the abandon of universal property. see \cite{zhang}. \paragraph{Level {n} structure} We need a extension of the elliptic curve {E}, this is given by the integral model.

    Definition 4 (Integral model) {s:=Spec(\mathcal{O}_F)}, {E\rightarrow E_s}. {E_s} is regular and minimal, the construction of {E_s} is by the following way, we first construct {\widetilde{E_s} } and then blow up. {\widetilde E_s} is given by the Weierstrass equation with coefficent in {\mathcal{O}_F}.

    Remark 3 The existence of integral model need Zorn’s lemma.

    Definition 5 (Semistable) the singularity of the minimal model of {E} are ordinary double point.

    Remark 4 Semistable is a crucial property, related to Szpiro’s conjecture.

    Definition 6 (Level {n} structure)

    \displaystyle \phi: ({\mathbb Z}/n{\mathbb Z})_s^2\longrightarrow E[N] \ \ \ \ \ (6)

    {P=\phi(1,0), Q=\phi(o,1)} The weil pairing of {P,Q} is given by a unit in cycomotic fields, i.e. {<P,Q>=\zeta_N\in \mu_{N}(s)}

    What happen if {k={\mathbb C}}? In this case we have a analytic isomorphism:

    \displaystyle E({\mathbb C})\simeq {\mathbb C}/\Lambda \ \ \ \ \ (7)

    Given by,

    \displaystyle {\mathbb C}/\Lambda \longrightarrow \mathop{\mathbb P}^2 \ \ \ \ \ (8)

    \displaystyle z\longrightarrow (\mathfrak{P}(z), \mathfrak{P}'(z), 1 ) \ \ \ \ \ (9)

    Where {\mathfrak{P(z)}=\frac{1}{z^2}+\sum_{\lambda\in \Lambda,\lambda\neq 0}(\frac{1}{(z-\lambda)^2}-\frac{1}{\lambda^2})}, and the Weierstrass equation {E} is given by {y^2=4x^3-60G_4(\Lambda)x-140G_6(\Lambda)}. The full n tructure of it is given by {{\mathbb Z}+{\mathbb Z}\lambda} and the value of {P,Q}, i.e.

    \displaystyle P=\frac{1}{N}, Q=\frac{\tau}{N} \ \ \ \ \ (10)

    Where {\tau} is induce by

    \displaystyle \Gamma(N):=ker(SL_2({\mathbb Z})\rightarrow SL_2({\mathbb Z}/n{\mathbb Z})) \ \ \ \ \ (11)

    The key point is following:

    Theorem 7 {k={\mathbb C}}, the moduli of elliptic curves with full level n-structure is identified with

    \displaystyle \mu_N^*\times H/\Gamma(N) \ \ \ \ \ (12)

    Now we discuss the Mordell-Weil theorem.

    Theorem 8 (Mordell-Weil theorem)

    \displaystyle E(F)\simeq {\mathbb Z}^r\oplus E(F)_{tor}

    The proof of the theorem divide into two part:

    1. Weak Mordell-Weil theorem, i.e. {\forall m\in {\mathbb N}}, {E(F)/mE(F)} is finite.
    2. There is a quadratic function,

      \displaystyle \|\cdot\|: E(F)\longrightarrow {\mathbb R} \ \ \ \ \ (13)

      {\forall c\in {\mathbb R}}, {E(F)_c=\{P\in E(F), \|P\|<c\}} is finite.

    Remark 5 The proof is following the ideal of infinity descent first found by Fermat. The height is called Faltings height, introduce by Falting. On the other hand, I point out, for elliptic curve {E}, there is a naive height come from the coefficient of Weierstrass representation, i.e. {\max\{|4a^3|,|27b^2|\}}.

    While the torsion part have a very clear understanding, thanks to the work of Mazur. The rank part of {E({\mathbb Q})} is still very unclear, we have the BSD conjecture, which is far from a fully understanding until now.

    But to understanding the meaning of the conjecture, we need first constructing the zeta function of elliptic curve, {L(s,E)}.

    \paragraph{Local points} We consider a local field {F_v}, and a locally value map {F\rightarrow F_{\nu}}, then we have the short exact sequences,

    \displaystyle 0\longrightarrow E^0(F_{\nu})\longrightarrow E(F_{\nu})=E_s(\mathcal{O}_F)\longrightarrow E_s(K_0)\longrightarrow 0 \ \ \ \ \ (14)

    Topologically, we know {E(F_{\nu})} are union of disc indexed by {E_s(k_{\nu})},

    \displaystyle |E_s(k_{\nu})| \sim q_{\nu}+1=\# \mathop{\mathbb P}^1(k_{\nu})

    . Define {a_{\nu}=\# \mathop{\mathbb P}^1(k_{\nu})-|E_s(k_{\nu})|}, then we have Hasse principle:

    Theorem 9 (Hasse principle)

    \displaystyle |a_{\nu}|\leq 2\sqrt{q_{\nu}} \ \ \ \ \ (15)

    Remark 6 I need to point out, the Hasse principle, in my opinion, is just a uncertain principle type of result, there should be a partial differential equation underlying mystery.

    So count the points in {E(F)} reduce to count points in {H^1(F_{\nu},E(m))}, reduce to count the Selmer group {S(E)[m]}. We have a short exact sequences to explain the issue.

    \displaystyle 0\longrightarrow E(F)/mE(F) \longrightarrow Sha(E)[m] \longrightarrow E(F)/mE(F)\longrightarrow 0 \ \ \ \ \ (16)

    I mention the Goldfold-Szipiro conjecture here. {\forall \epsilon>0}, there {\exists C_{\epsilon}(E)} such that:

    \displaystyle \# (E)\leq c_{\epsilon}(E)N_{E/{\mathbb Q}}(N)^{\frac{1}{2}+\epsilon} \ \ \ \ \ (17)

    \paragraph{L-series} Now I focus on the construction of {L(s,E)}, there are two different way to construct the L-series, one approach is the Euler product.

    \displaystyle L(s,E)=\prod_{\nu: bad}(1-a_{\nu}q_{\nu}^{-s})^{-1}\cdot \prod_{\nu:good}(1-a_{\nu}q_{\nu}^{-s}+q_{\nu}^{1-2s})^{-1} \ \ \ \ \ (18)

     

    Where {a_{\nu}=0,1} or {-1} when {E_s} has bad reduction on {\nu}.

    The second approach is the Galois presentation, one of the advantage is avoid the integral model. Given {l} is a fixed prime, we can consider the Tate module:

    \displaystyle T_l(E):=\varprojlim_{l^n} E[l^n] \ \ \ \ \ (19)

    Then by the transform of different embedding of {F\hookrightarrow \bar F}, we know { T_{l}(E)/Gal(\bar F/F)}, decompose it into a lots of orbits, so we can define {D_{\nu}}, the decomposition group of {w}(extension of {\nu} to {\bar F}). We define {I_{\nu}} is the inertia group of {D_{\nu}}.

    Then {D_{\nu}/I_{\nu}} is generated by some Frobenius elements

    \displaystyle Frob{\nu}x\equiv x^{q_{\nu}} (mod w),\forall x\in \mathcal{O}_{\bar Q} \ \ \ \ \ (20)

    So we can define

    \displaystyle L_{\nu}(s,E)=(1-q_{\nu}^{-s}Frob_{\nu}|T_{l}(E)^{I_{\nu}})^{-1} \ \ \ \ \ (21)

    And then {L(s,E)=\prod_{\nu}L_{\nu}(s,E)}.

    Faltings have proved {L_{\nu}(s,E)} is the invariant depending the isogenous class in the follwing meaning:

    Theorem 10 (Faltings) {L_{\nu}(s,E)} is an isogenous ivariant, i.e. {E_1} isogenous to {E_2} iff {\forall a.e. \nu}, {L_{\nu}(s,E_1)=L_{\nu}(s,E_2)}.

    \displaystyle L(s,E)=L(s-\frac{1}{2},\pi ) \ \ \ \ \ (22)

    Where {\pi} come from an automorphic representation for {GL_2(A_F)}. Now we give the statement of BSD onjecture. {R} is the regulator of {E}, i.e. the volume of fine part of {E(F)} with respect to the Neron-Tate height pairing. {\Omega} be the volume of {\prod_{v|\infty}F(F_v)} Then we have,

    1. {ord_{s=1}L(s,E)=rank E(F)}.
    2. {|Sha(E)|<\infty}.
    3. {\lim_{s\rightarrow 0}L(s,E)(s-1)^{-rank(E)}=c\cdot \Omega(E)\cdot R(E)\cdot |Sha(E)|\cdot |E(F)_{tor}|^{-2}}

    Here {c} is an explictly positive integer depending only on {E_{\nu}} for {\nu} dividing {N}.

     


    补充说明

    以下是新整理的中文说明;上方旧博客原文保持不变。

    Birch and Swinnerton-Dyer 猜想把椭圆曲线的有理点群和它的 $L$ 函数在 $s=1$ 处的零点阶联系起来。它是数论中最核心的桥之一:一边是 Diophantine 方程的解,另一边是解析函数的特殊值。

    BSD 猜想入门:椭圆曲线、L 函数与 Mordell-Weil rank
    BSD 猜想连接椭圆曲线有理点群的 rank 与 $L$ 函数在 $s=1$ 的零点阶。

    1. Weierstrass 形式

    在特征不是 $2,3$ 的域上,椭圆曲线常写成

    $$E:y^2=x^3+ax+b,$$

    并要求判别式

    $$\Delta=-16(4a^3+27b^2)\ne0.$$

    非零判别式保证曲线光滑。椭圆曲线不仅是代数曲线,它的点还带有 Abelian group 结构。

    2. $j$-invariant 与 twist

    $j$-invariant 分类椭圆曲线在代数闭包上的同构类:

    $$j(E)=1728\frac{4a^3}{4a^3+27b^2}.$$

    但在固定域 $K$ 上,同一个 $j$-invariant 可能对应不同的 twists。twist 说明“几何上同构”和“在基域上同构”之间有差别。

    3. Mordell-Weil theorem

    若 $E$ 定义在 $\mathbb Q$ 上,则 Mordell-Weil theorem 说

    $$E(\mathbb Q)\cong E(\mathbb Q)_{\operatorname{tors}}\oplus\mathbb Z^r.$$

    整数 $r$ 称为 Mordell-Weil rank。求 $r$ 是理解有理点结构的核心问题。

    4. Hasse-Weil L 函数

    对每个好素数 $p$,令

    $$a_p=p+1-\#E(\mathbb F_p).$$

    椭圆曲线的 $L$ 函数由 Euler product 组成:

    $$L(E,s)=\prod_p(1-a_pp^{-s}+p^{1-2s})^{-1}$$

    再在坏素数处加入修正因子。模性定理保证这个 $L$ 函数有解析延拓和函数方程。

    5. BSD 猜想

    BSD 猜想最核心的断言是

    $$\operatorname{rank}E(\mathbb Q)=\operatorname{ord}_{s=1}L(E,s).$$

    更精细的版本还给出 $L(E,s)$ 在 $s=1$ 处首项系数,涉及 regulator、Tate-Shafarevich group、Tamagawa numbers 和 torsion subgroup。

    这个猜想的意义在于,它把有理点这个离散、代数的问题,转化成了 $L$ 函数特殊值这个解析问题。许多数论现代方法都在这座桥上来回移动。

  • $SL_2(\mathbb Z)$ 与 congruence subgroups:分式线性作用和 modular curves

    旧博客原文

    原题:SL_2(Z) and its congruence subgroups

    The pdf version is SL_2(Z) and its congruence subgroups.

    We know we can always do the following thing:

    \displaystyle R\ commutative\ ring \longrightarrow \ "general\ linear\ group" \ GL_2(R) \ \ \ \ \ (1)

    Where

    \displaystyle GL_2(R):=\{\begin{pmatrix} a & b \\ c & d \end{pmatrix}: det \begin{pmatrix} a & b \\ c & d \end{pmatrix}=R^*, a,b,c,d\in R\} \ \ \ \ \ (2)

     

    Remark 1 Why it is {R^*} but not 1? if it is 1, then the action {R/GL_2(R)} distribute is not trasitive on {R}, i.e. every element in unite group present a connected component

    Now we consider the subgroup {SL_2(R)\subset GL_2(R)}.

    \displaystyle SL_2(R)=\{\begin{pmatrix} a & b\\ c & d \end{pmatrix}: det \begin{pmatrix} a & b\\ c & d \end{pmatrix}=1, a,b,c,d\in R \} \ \ \ \ \ (3)

    We are most interested in the case {R={\mathbb Z}, {\mathbb Z}/n{\mathbb Z}}. So how to investigate {SL_2({\mathbb R})}? We can look at the action of it on something, for particular, we look at the action of it on Riemann sphere, i.e. { \hat {\mathbb C}/({\mathbb R})} given by fraction linear map:

    \displaystyle g(z):=\frac{az+b}{cz+d}, g(\infty)=\frac{a}{c} \ \ \ \ \ (4)

    Remark 2 What is fraction linear map? This action carry much more information than the action on vector, thanks for the exist of multiplication in {{\mathbb C}} and the algebraic primitive theorem. Due to I always looks the fraction linear map as something induce by the permutation of the roots of polynomial of degree 2, this is true at least for fix points, and could natural extension. So how about the higher dimension generate? consider the transform of {k-1} tuples induce by polynomial with degree {k}?

    Remark 3

    1. {SL_2({\mathbb R})/\pm I:=PSL_2({\mathbb R})}, then {PSL_2({\mathbb R}) } action faithful on {\hat C}, i.e. except identity, every action is nontrivial. This is easy to be proved, observed,

      \displaystyle \frac{az+b}{cz+d}=z,\forall z\in \hat{\mathbb C}\Longrightarrow \begin{pmatrix} a & b\\ c & d \end{pmatrix}=\begin{pmatrix} 1 & 0\\ 0&1 \end{pmatrix} or \begin{pmatrix} -1&0\\ 0&-1 \end{pmatrix} \ \ \ \ \ (5)

    2. Up half plane {H} is invariant under the action of {PSL_2({\mathbb R})}, i.e. {\forall g\in PSL_2({\mathbb R})}, {gH=H}. The proof is following,
    3. \displaystyle \begin{array}{rcl} Im(\frac{az+b}{cz+d}) & = & Im(\frac{(az+b)(c\bar z+d)}{|cz+d|^2})\\ & = & Im(\frac{ac|z|^2+bc\bar z+adz+bd}{|cz+d|^2})\\ & > & 0,\ due\ to\ ad=bc+1. \end{array}

     

    Now we focus on {SL_2({\mathbb Z})} or the same,{PSL({\mathbb Z})}. All the argument for {SL_2(R)} make sense for

    \displaystyle \Gamma:= SL_2({\mathbb Z}), \bar \Gamma:=SL_2({\mathbb Z})/\pm I \ \ \ \ \ (6)

    Fix {N\in {\mathbb N}}, define,

    \displaystyle \Gamma(N):=\{\begin{pmatrix} a &b\\ c&d \end{pmatrix}, a,d \equiv 1(mod N), b,c\equiv 0(mod N)\} \ \ \ \ \ (7)

    Then {\Gamma(N)} is the kernel of map {SL_2({\mathbb Z})\rightarrow SL_2({\mathbb Z}/n{\mathbb Z})}, i.e. we have short exact sequences,

    \displaystyle 0\longrightarrow \Gamma(N)\longrightarrow \Gamma\longrightarrow SL_2({\mathbb Z}/n{\mathbb Z})\longrightarrow 0 \ \ \ \ \ (8)

    Remark 4 The relationship of {\Gamma(N)\subset \Gamma} is just like {N{\mathbb Z}+1\subset {\mathbb Z}}.

    Definition 1 (Congruence group) A subgroup of {\Gamma} is called a congruence group iff {\exists n\in {\mathbb N}}, {\Gamma(N)\subset G}.

    Example 1 We give two examples of congruence subgroups here.

    1. \displaystyle \Gamma_1(N)=\{\begin{pmatrix} 1& *\\ 0 &1 \end{pmatrix} mod N\} \ \ \ \ \ (9)

    2. \displaystyle \Gamma_0(N)=\{\begin{pmatrix} * & *\\ 0 & * \end{pmatrix}mod N\} \ \ \ \ \ (10)

     

    Definition 2 (Fundamental domain)

    \displaystyle F=\{z\in H:-\frac{1}{2}\leq Re(z)\leq \frac{1}{2}\ and |z|\geq 1\} \ \ \ \ \ (11)

    Now here is a theorem charistization the fundamental domain.

    Theorem 3 This domain {F} is a fundamental domain of {\hat {\mathbb H}/({\mathbb R})}

    Proof: Ths key point is {SL_2({\mathbb Z})} have two generators,

    1. {\tau_a: z\rightarrow z+a, \forall a\in {\mathbb Z}}.
    2. {s: z\rightarrow \frac{1}{z}}.

    Thanks to this two generator exactly divide the action of {\Gamma} on {H} into a lots of scales, then {\Omega} is a fundamental domain is a easy corollary. \Box

    Remark 5 This is not rigorous, {H} need be replace by {\hat H}, but this is very natural to get a modification to a right one.

    Remark 6 {z_1,z_2\in \partial F} are {\Gamma} equivalent iff {Re(z)=\pm \frac{1}{2}} and {z_2=z_1\pm 1} or if {z_1} on the unit circle and {z_2=-\frac{1}{z_1}}

    Remark 7 If {z\in F}, then {\Gamma_z=\pm I} expect in the following three case:

    1. {\Gamma=\pm \{\tau,s\}} if {z=i}.
    2. {\Gamma=\pm\{ I,s\tau, (s\tau)^2\}} if {z=w=-\frac{1}{2}+\frac{\sqrt{-3}}{2}}.
    3. {\Gamma=\pm\{I,\tau s, (\tau s)^2\}} if {z=-\bar w=\frac{1}{2}+\frac{\sqrt{-3}}{2}}.

    Where {\tau=\tau_1}.

    Remark 8 The group {\bar \Gamma=SL_2({\mathbb Z})/\pm I} is generated by the two elements {s}, {\tau}. In other word, any fraction linear transform is a “word” induce by {s,\tau,s^{-1}.\tau^{-1}}. But not free group, we have relationship {s^2=-I,(s\tau)^3=-I}.

    The natural function space on {F} is the memorphic function, under the map: {H\rightarrow D-\{0\}}, it has a {q}-expension,

    \displaystyle f(q)=\sum_{k\in {\mathbb Z}}a_kq^k \ \ \ \ \ (12)

    And there are only finite many negative {k} such that {a_k\neq 0}.


    补充说明

    以下是新整理的中文说明;上方旧博客原文保持不变。

    $SL_2(\mathbb Z)$ 是 modular forms 和 modular curves 的基本群。它通过分式线性变换作用在上半平面上,而 congruence subgroups 则把模 $N$ 的算术信息放进几何商空间。

    $SL_2(\mathbb Z)$ 与 congruence subgroups:分式线性作用和 modular curves
    $SL_2(\mathbb Z)$ 通过分式线性变换作用在上半平面,congruence subgroups 对应 modular curves 的 level structure。

    1. 分式线性作用

    对

    $$\gamma=\begin{pmatrix}a&b\\ c&d\end{pmatrix}\in SL_2(\mathbb Z),$$

    定义

    $$\gamma z=\frac{az+b}{cz+d}.$$

    若 $\operatorname{Im}z>0$,则 $\operatorname{Im}(\gamma z)>0$,所以上半平面 $\mathbb H$ 在作用下不变。

    2. 为什么是 $PSL_2$

    $I$ 和 $-I$ 给出同一个分式线性变换。因此真正忠实作用的是

    $$PSL_2(\mathbb Z)=SL_2(\mathbb Z)/\{\pm I\}.$$

    3. Congruence subgroups

    主同余子群定义为

    $$\Gamma(N)=\ker\bigl(SL_2(\mathbb Z)\to SL_2(\mathbb Z/N\mathbb Z)\bigr).$$

    更常见的还有 $\Gamma_0(N)$ 和 $\Gamma_1(N)$,它们分别对矩阵的某些项加模 $N$ 条件。

    4. Modular curves

    商空间

    $$Y(\Gamma)=\Gamma\backslash\mathbb H$$

    是 modular curve 的开部分。补上 cusps 后得到紧化 $X(\Gamma)$。不同 congruence subgroup 对应不同 level structure。

    5. 算术与几何

    $SL_2(\mathbb Z)$ 的作用把矩阵、分式线性变换、椭圆曲线的 level structure 和 modular forms 连接起来。研究 congruence subgroups,就是研究模 $N$ 算术信息如何改变上半平面商的几何。

  • Dirichlet hyperbola method:用双曲线拆分除数和

    旧博客原文

    原题:Dirichlet hyperbola method

    A pdf version is Dirichlet hyperbola method.

    1. Introduction

    Theorem 1

    \displaystyle \sum_{1\leq n\leq x}d(n)=\sum_{1\leq n\leq x}[\frac{x}{n}]=xlogx+(2\gamma-1) x+O(\sqrt{x}) \ \ \ \ \ (1)

     

    Remark 1 I thought this problem initial 5 years ago, cost me several days to find a answer, I definitely get something without the argument of Dirchlet hyperbola method and which is weaker but morally the same camparable with the result get by Dirichlet hyperbola method.

    Remark 2 How to get the formula:

    \displaystyle \sum_{1\leq n\leq x}d(x)=\sum_{1\leq n\leq x}[\frac{x}{n}]? \ \ \ \ \ (2)

    In fact,

    \displaystyle \sum_{1\leq n\leq x}d(x)=\sum_{1\leq ab\leq x}1=\sum_{1\leq n\leq x}[\frac{x}{n}]\ \ \ \ \ (3)

    Which is the integer lattices under or lying on the hyperbola {\{(a,b)|ab=x\}}.

    Remark 3 By trivial argument, we can bound the quantity as following way,

    \displaystyle \begin{array}{rcl} \sum_{1\leq n\leq x}[\frac{x}{n}] & = & \sum_{1\leq ab\leq x}1\\ & = & x\sum_{i=1}^x\frac{1}{i}-\sum_{i=1}^x\{\frac{x}{i}\}\\ & = &xlnx+\gamma x+O(x) \end{array}

    The error term is {O(x)}, which is too big. But fortunately we can use the symmetry of hyperbola to improve the error term.

    Proof:

    \displaystyle \begin{array}{rcl} \sum_{1\leq n\leq x}d(n) & = & \sum_{ab\leq x}1\\ & = & \sum_{a\geq \sqrt{x}}[\frac{x}{b}]+\sum_{b\geq \sqrt{x}}[\frac{x}{a}]-\sum_{1\leq a,b\leq \sqrt{x}}1\\ & = & xlogx+(2\gamma-1)x+O(\sqrt{x}) \end{array}

    \Box

    Theorem 2

    Given a natural number k, use the hyperbola method together
    with induction and partial summation to show that

    \displaystyle \sum_{n\leq x}d_k(n) = xP_k(log x) + O(x^{1-\frac{1}{k}+\epsilon}), n\leq x \ \ \ \ \ (4)

    where {P_k(t)} denotes a polynomial of degree {k-1} with leading term {\frac{t^{k-1}}{(k-1)!}}.

    Remark 4 {P_k(x)} is the residue of {\zeta(s)^kx^ss^{-1}} at {s=1}.

    Proof:

    We can establish the dimension 3 case directly, which is the following asymptotic formula,

    \displaystyle \sum_{1\leq xy\leq n}[\frac{n}{xy}]=xP_2(logx)+O(x^{1-\frac{1}{3}+\epsilon}) \ \ \ \ \ (5)

    The approach is following, we first observe that

    \displaystyle \sum_{1\leq xy\leq n}[\frac{n}{xy}]=\sum_{xyz\leq n}1 \ \ \ \ \ (6)

    The problem transform to get a asymptotic formula for the lattices under 3 dimension hyperbola. The first key point is, morally {([n^{\frac{1}{3}}],[n^{\frac{1}{3}}],[n^{\frac{1}{3}}])} is the central point under the hyperbola.

    Then we can divide the range into 3 parts, and try to get a asymptotic formula for each part then add them together. Assume we have:

    1. {A_x=\sum_{1\leq r\leq [n^{\frac{1}{3}}]}\sum_{1\leq yz\leq [n^{\frac{2}{3}}]}[\frac{r}{yz}]}.
    2. {A_y=\sum_{1\leq r\leq [n^{\frac{1}{3}}]}\sum_{1\leq xz\leq [n^{\frac{2}{3}}]}[\frac{r}{yz}]}.
    3. {A_z=\sum_{1\leq r\leq [n^{\frac{1}{3}}]}\sum_{1\leq xy\leq [n^{\frac{2}{3}}]}[\frac{r}{yz}]}.

    Then the task transform to get a asymptotic formula,

    \displaystyle A_x=A_y=A_z=xQ_2(logx)+O(x^{1-\frac{1}{3}+\epsilon}) \ \ \ \ \ (7)

    But we can do the same thing for {\sum_{1\leq yz\leq [n^{\frac{2}{3}}]}[\frac{r}{yz}]} and then integral it. This end the proof. For general {k\in {\mathbb N}}, the story is the same, by induction.

    Induction on {k} and use the Fubini theorem to calculate {\sum_{x_1...x_r\leq n}\frac{n}{x_1...x_r},\forall 1\leq r\leq k}. \Box

    There is a major unsolved problem called Dirichlet divisor problem.

    \displaystyle \sum_{n\leq x}d(n) \ \ \ \ \ (8)

    What is the error term? The conjecture is the error term is {O(x^{\theta}), \forall \theta>\frac{1}{4}}, it is known that {\theta=\frac{1}{4}} is not right.

    Remark 5

    To beats this problem, need some tools in algebraic geometry.

    2. Several problems

    {\forall k\in {\mathbb N}}, is there a asymptotic formula for {\sum_{t=1}^n\{\frac{kn}{t}\}} ?

    {\forall k\in {\mathbb N}}, {f(n)} is a polynomial with degree {k}, is there a asymptotic formula for {\sum_{t=1}^n\{\frac{f(n)}{t}\}} ?

    {\forall k\in {\mathbb N}}, {g(n)} is a polynomial with degree {k}, is there a asymptotic formula for {\sum_{t=1}^n\{\frac{n}{g(t)}\}} ?

    Theorem 3 {k\in {\mathbb N}}, then we have

    \displaystyle \lim_{n\rightarrow \infty}\frac{\{\frac{kn}{1}\}+\{\frac{kn}{2}\}+...+\{\frac{kn}{n}\}}{n}=k(\sum_{i=1}^k\frac{1}{i}-lnk-\gamma) \ \ \ \ \ (9)

    Proof:

    \displaystyle \begin{array}{rcl} \frac{\{\frac{kn}{1}\}+\{\frac{kn}{2}\}+...+\{\frac{kn}{n}\}}{n} & = & \frac{\sum_{i=1}^k\frac{kn}{i}-\sum_{i=1}^n[\frac{kn}{i}]}{n}\\ & = & k(lnn+\gamma +\epsilon_n)-\frac{\sum_{i=1}^{kn}[\frac{kn}{i}]-\sum_{i=n+1}^{kn}[\frac{kn}{i}]}{n} \end{array}

    \Box

    Now we try to estimate

    \displaystyle S_k(n)=\sum_{i=1}^{kn}[\frac{kn}{i}]-\sum_{i=n+1}^{kn}[\frac{kn}{i}] \ \ \ \ \ (10)

    In fact, we have,

    \displaystyle \begin{array}{rcl} S_k(n) & = & (2\sum_{i=1}^{[\sqrt{kn}]}[\frac{kn}{i}]-[\sqrt{kn}]^2)-(\sum_{i=1}^k[\frac{kn}{i}]-kn)\\ & = & 2\sum_{i=1}^{[\sqrt{kn}]}\frac{kn}{i}-\sum_{i=1}^k\frac{kn}{i}+2\{\sqrt{kn}\}[\sqrt{kn}]+\{\sqrt{kn}\}^2-2\sum_{i=1}^{[\sqrt{kn}]}\{\frac{kn}{i}\}+\sum_{i=1}^k\{\frac{kn}{i}\}\\ & = & 2kn(ln[\sqrt{kn}]+\gamma+\epsilon_{[\sqrt{kn}]})-kn\sum_{i=1}^k\frac{1}{i}+r(n)\\ & = & knln(kn)+kn(2\gamma-\sum_{i=1}^k\frac{1}{i})+r'(n)\\ & = & knln n+kn(2\gamma+lnk-\sum_{i=1}^k\frac{1}{i})+r'(n) \end{array}

    Where {-3\sqrt{n}<r(n)<3\sqrt{n}}, {-3\sqrt{n}<r'(n)<3\sqrt{n}}.

    So by 1 we know,

    \displaystyle \begin{array}{rcl} \frac{\{\frac{kn}{1}\}+...+\{\frac{kn}{n}\}}{n} & = & k(lnn+\gamma+\epsilon_n)-klnn-k(2\gamma+lnk-\sum_{i=1}^k\frac{1}{i})+\frac{r'(n)}{n}\\ & = & k(\sum_{i=1}^k\frac{1}{i}-lnk-\gamma)+\frac{r'(n)}{n}+k\epsilon_n \end{array}

    So we have,

    \displaystyle \lim_{n\rightarrow \infty}\frac{\{\frac{kn}{1}\}+...+\{\frac{kn}{n}\}}{n} =k(\sum_{i=1}^k\frac{1}{i}-lnk-\gamma)=k\epsilon_k \ \ \ \ \ (11)

    Remark 6 In fact we can get {0<k\epsilon_k<\frac{1}{2}, \forall k\in {\mathbb N}}, by combining the theorem 3 and 1.

    3. Lattice points in ball

    Gauss use the cube packing circle get a rough estimate,

    \displaystyle \sum_{n\leq x}r_2(n)=\pi x+O(\sqrt{x}) \ \ \ \ \ (12)

     

    In the same way one can obtain,

    \displaystyle \sum_{n\leq x}r_k(n)=\rho_kx^{\frac{k}{2}}+O(x^{\frac{k-1}{2}}) \ \ \ \ \ (13)

    Remark 7 Where {\rho_k=\frac{\pi^{\frac{k}{2}}}{\Gamma(\frac{k}{2}+1)}} is the volume of the unit ball in {k} dimension.

    Dirchlet’s hyperbola method works nicely for the lattic points in a ball of dimension {k\geq 4}. Langrange proved that every natural number can be represented as the sum of four squares, i.e. {r_4(n)>0}, and Jacobi established the exact formula for the number of representations

    \displaystyle r_4(n)=8(2+(-1)^n)\sum_{d|n,d\ odd}d. \ \ \ \ \ (14)

    Hence we derive,

    \displaystyle \begin{array}{rcl} \sum_{n\leq x}r_4(n) & = & 8\sum_{m\leq x}(2+(-1)^m)\sum_{dm\leq x, d\ odd}d\\ & = & 8\sum_{m\leq x}(2+(-1)^m)(\frac{x^2}{4m^2}+O(\frac{x}{m}))\\ & = & 2x^2\sum_1^{\infty}(2+(-1)^m)m^{-2}+O(xlogx)\\ & = & 3\zeta(2)x^2+O(xlogx) = \frac{1}{2}(\pi x)^2+O(xlogx) \end{array}

    This result extend easily for any {k\geq 4}, write {r_k} as the additive convolution of {r_4} and {r_{k-4}}, i.e.

    \displaystyle r_k(n)=\sum_{0\leq t\leq n}r_4(t)r_{k-4}(n-t) \ \ \ \ \ (15)

    Apply the above result for {r_4} and execute the summation over the remaining {k-4} squares by integration.

    \displaystyle \sum_{n\leq x}r_k(n)=\frac{(\pi x)^{\frac{k}{2}}}{\Gamma(\frac{k}{2}+1)}+O(x^{\frac{k}{2}-1}logx) \ \ \ \ \ (16)

     

    Remark 8 Notice that this improve the formula 12 which was obtained by the method of packing with a unit square. The exponent {\frac{k}{2}-1} in 16 is the best possible because the individual terms of summation can be as large as the error term (apart from {logx}), indeed for {k=4} we have {r_4(n)\geq 16n} if {n} is odd by the Jacobi formula. The only case of the lattice point problem for a ball which is not yet solved (i.e. the best possible error terms are not yet established) are for the circle({k=2}) and the sphere ({k=3}).

    Theorem 4

    \displaystyle \sum_{n\leq x}\tau(n^2+1)=\frac{3}{\pi}xlogx+O(x) \ \ \ \ \ (17)

    4. Application in finite fields

    Suppose {f(x)\in {\mathbb Z}[x]} is a irreducible polynomial. And for each prime {p}, let

    \displaystyle \rho_f(p)=\# \ of \ solutions\ of f(x)\equiv 0(mod\ p) \ \ \ \ \ (18)

    By Langrange theorem we know {\rho_f(p)\leq deg(f)}. Is there a asymptotic formula for

    \displaystyle \sum_{p\leq x}\rho_f(p)? \ \ \ \ \ (19)

    A general version, we can naturally generated it to algebraic variety.

    \displaystyle \rho_{f_1,...,f_k}(p)=\#\ of \ solutions\ of f_i(x)\equiv 0(mod\ p) ,\ \forall 1\leq i\leq k \ \ \ \ \ (20)

    Is there a asymptotic formula for

    \displaystyle \sum_{p\leq x}\rho_{f_1,...,f_k}(p)? \ \ \ \ \ (21)

    Example 1 We give an example to observe what is involved. {f(x)=x^2+1}. We know {x^2+1\equiv 0 (mod \ p)} is solvable iff {p\equiv 1 (mod\ 4)} or {p=2}. One side is easy, just by Fermat little theorem, the other hand need Fermat descent procedure, which of course could be done by Willson theorem. In this case,

    \displaystyle \sum_{p\leq n}\rho_f(p)=\# \ of \{primes \ of \ type\ 4k+1 \ in \ 1,2,...,n\} \ \ \ \ \ (22)

    Which is a special case of Dirichlet prime theorem.

    Let {K} be an algebraic number field, i.e. the finite field extension of rational numbers, let

    \displaystyle \mathcal{O}_K=\{\alpha\in K, \alpha \ satisfied \ a\ monic \ polynomial\ in\ {\mathbb Z}[x]\} \ \ \ \ \ (23)

     

    Dedekind proved that,

    Theorem 5

    1. {\mathcal{O}_K} is a ring, we call it the ring of integer of {K}.
    2. He showed further every non-zero ideal of {\mathcal{O}_K} could write as the product of prime ideal in {\mathcal{O}_k} uniquely.
    3. the index of every non-zero ideal {I} in {\mathcal{O}_K} is finite, i.e. {[\mathcal{O}_K:I]<\infty}, and we can define the norm induce by index.

      \displaystyle N(I):=[\mathcal{O}_K:I] \ \ \ \ \ (24)

      Then the norm is a multiplication function in the space of ideal, i.e. {N(IJ)=N(I)N(J), \forall I,J \in \ ideal\ class\ group\ of\ \mathcal{O}_K}.

    4. Now he construct the Dedekind Riemann zeta function,

      \displaystyle \zeta_K(s)=\sum_{N(I)\neq 0}\frac{1}{N(I)^s}=\prod_{J\ prime \ ideal\ }\frac{1}{1-\frac{1}{N(J)^s}},\ \forall Re(s)>1 \ \ \ \ \ (25)

     

    Now we consider the analog of the prime number theorem. Let {\pi_K(x)=\{I,N(I)<x\}}, does the exist a asymptotic formula,

    \displaystyle \pi_K(x)\sim \frac{x}{ln x}\ as\ x\rightarrow \infty? \ \ \ \ \ (26)

    Given a prime {p}, we may consider the prime ideal

    \displaystyle p\mathcal{O}_K=\mathfrak{P}_1^{e_1}\mathfrak{P}_2^{e_2}...\mathfrak{P}_k^{e_k} \ \ \ \ \ (27)

    Where {\mathfrak{P}_i } is different prime ideal in {\mathcal{O}_K}. But the question is how to find these {\mathfrak{P}_i}? For the question, there is a satisfied answer.

    Lemma 6 (existence of primitive element) There always exist a primetive elements in {K}, such that,

    \displaystyle K={\mathbb Q}(\theta) \ \ \ \ \ (28)

    Where {\theta} is some algebraic number, which’s minor polynomial {f(x)\in {\mathbb Z}[x]}.

    Theorem 7 (Dedekind recipe) Take the polynomial {f(x)}, factorize it in the polynomial ring {{\mathbb Z}_p[x]},

    \displaystyle f(x)\equiv f_1(x)^{e_1}...f_{r}(x)^{e_r}(mod \ p) \ \ \ \ \ (29)

    Consider {\mathfrak{P}_i=(p, f_i(\theta)) \subset \mathcal{O}_K}. Then apart from finite many primes, we have,

    \displaystyle p\mathcal{O}_K=\mathfrak{P}_1^{e_1}\mathfrak{P}_2^{e_2}...\mathfrak{P}_k^{e_k} \ \ \ \ \ (30)

    Where {N(\mathfrak{P}_i)=p^{deg{f_i}}}.

    Remark 9 The apart primes are those divide the discriminant.

    Now we can argue that 4 is morally the same as counting the ideals whose norm is divide by {p} in a certain algebraic number theory.

    And we have following, which is just the version in algebraic number fields of 2.

    Theorem 8 (Weber) {\#} of ideals of {\mathcal{O}_K} with norm {\leq x} equal to,

    \displaystyle \rho_k(X)+O(x^{1-\frac{1}{d}}), where \ d=[K:Q] \ \ \ \ \ (31)

     


    补充说明

    以下是新整理的中文说明;上方旧博客原文保持不变。

    Dirichlet hyperbola method 是解析数论中处理卷积和的基本工具。它的几何图像非常直接:把求和区域看成双曲线 $ab\le x$ 下方的格点,再利用双曲线关于 $\sqrt x$ 的对称性减少误差。

    Dirichlet hyperbola method:用双曲线拆分除数和
    Dirichlet hyperbola method 用双曲线 $ab=x$ 的对称性改进卷积和的误差。

    1. 除数函数的例子

    除数函数满足

    $$d(n)=\sum_{ab=n}1.$$

    因此

    $$\sum_{n\le x}d(n)=\sum_{ab\le x}1.$$

    这就是双曲线 $ab=x$ 下方的整数格点数。

    2. 直接估计的问题

    若对每个 $a$ 求 $\lfloor x/a\rfloor$,得到

    $$\sum_{a\le x}\left\lfloor\frac xa\right\rfloor.$$

    平凡地把 floor 换成 $x/a$ 会产生太大的误差,因为项数有 $x$ 个。hyperbola method 的关键是只在 $a\le\sqrt x$ 和 $b\le\sqrt x$ 的短范围内精确处理。

    3. 基本公式

    由对称性可得

    $$\sum_{n\le x}d(n)=2\sum_{a\le\sqrt x}\left\lfloor\frac xa\right\rfloor-\lfloor\sqrt x\rfloor^2.$$

    于是

    $$\sum_{n\le x}d(n)=x\log x+(2\gamma-1)x+O(\sqrt x).$$

    这个误差已经比直接方法好很多。

    4. 一般卷积

    若 $h=f*g$,则

    $$\sum_{n\le x}h(n)=\sum_{ab\le x}f(a)g(b).$$

    可以选择参数 $Y$,把区域分成 $a\le Y$、$b\le x/Y$ 和重叠部分。合适的 $Y$ 取决于 $f,g$ 的平均阶和可用误差估计。

    5. 高维推广

    对 $k$ 重除数函数 $d_k(n)$,问题变成

    $$a_1a_2\cdots a_k\le x$$

    下方的格点计数。通过归纳、partial summation 和多维双曲面拆分,可以得到

    $$\sum_{n\le x}d_k(n)=xP_{k-1}(\log x)+\text{error},$$

    其中 $P_{k-1}$ 是次数 $k-1$ 的多项式。这个方法的力量在于:它把乘法卷积的求和问题变成可视化的几何区域拆分。

  • Diophantine approximation:Dirichlet 定理、抽屉原理与最佳逼近

    旧博客原文

    原题:Diophantine approximation

    I explain some general ideal in the theory of diophantine approximation, some of them is original by myself, begin with a toy model, then consider the application on folklore Swirsing-Schmidt conjecture.

    \tableofcontents

    1. Dirichlet theorem, the toy model

    The very basic theorem in the theory of Diophantine approximation is the well known Dirichlet approximation theorem, the statement is following.

    Theorem 1 (Dirichlet theorem) for all {\alpha} is a irrational number, we have infinity rational number {\frac{q}{p}} such that:

    \displaystyle |\alpha-\frac{q}{p}|<\frac{1}{p^2} \ \ \ \ \ (1)

    Remark 1 It is easy to see the condition of irrational is crucial. There is a best constant version of it, said, instead of {1}, the best constant in the suitable sense for the theorem 1 should be {\frac{1}{\sqrt{5}}} and arrive by {\frac{\sqrt{5}+1}{2}} at least. The strategy of the proof of the best constant version involve the Frey sequences.

    Now we begin to explain the strategies to attack the problem.

    \paragraph{Argument 1, boxes principle} We begin with a easiest one, i.e. by the argument of box principle, the box principle is following,

    Theorem 2 (Boxes principle) Given {n\in {\mathbb N}} and two finite sets {A={a_1,a_2,...,a_n,a_{n+1}}}, set {B={b_1,...,b_{n}}}, if we have a map:

    \displaystyle f:A\longrightarrow B \ \ \ \ \ (2)

    Then there exists a element {b_k\in B} such that there exist at least two element {a_i,a_j\in A}, {f(a_i)=f(a_j)=b_k}.

    Proof: The proof is trivial. \Box

    Now consider, {\forall N\in {\mathbb N}}, the sequences {x,2x,...,Nx}, then {\{ix\}\in [0,1], \forall i\in \{1,2,...,n\}}. Divide {[0,1]} in an average way to {N} part: {[\frac{k-1}{N},\frac{k}{N}]}. Then the linear structure involve (which, in fact play a crucial role in the approach). And the key point is to look at {\{nx\}} and integers.

    \paragraph{Argument 2, continue fractional} We know, for irrational number {x}, {x} have a infinite long continue fractional:

    \displaystyle x=q_0+\frac{1}{q_1+\frac{1}{q_2+\frac{1}{q_3+....+\frac{1}{q_k+...}}}} \ \ \ \ \ (3)

    Then

    \displaystyle |x-q_0+\frac{1}{q_1+\frac{1}{q_2+\frac{1}{q_3+....+\frac{1}{q_k}}}}|\sim \frac{1}{(q_1q_2...q_{k-1})^2q_k} \ \ \ \ \ (4)

    And we have,

    \displaystyle \frac{1}{q_1+\frac{1}{q_2+\frac{1}{q_3+....+\frac{1}{q_k}}}}=\frac{a_n}{b_n}, (a_n,b_n)=1 \ \ \ \ \ (5)

    Then {b_n=O(q_1...q_k)}.

    \paragraph{Argument 3, Bohr set argument} We begin with some kind of Bohr set:

    \displaystyle B_p=I-\cup_{q\in \{0,1,...,p-1\}}(\frac{q}{p}-\frac{1}{p^2},\frac{q}{p}+\frac{1}{p^2}) \ \ \ \ \ (6)

    The key point is the shift of Bohr set, on the vertical line i.e. {|B_p\cap B_{p+1}|} is very slow, and can be explained by

    \displaystyle \frac{k}{p+1}+\frac{1}{(p+1)^2}>\frac{k}{p}-\frac{1}{p^2} \ \ \ \ \ (7)

    So:

    \displaystyle \frac{1}{p^2}+\frac{1}{(p+1)^2}>\frac{k}{p(p+1)} \ \ \ \ \ (8)

    in {|B_p \cap B_{p+1}|\sim \frac{1}{p(p+1)}} But in fact they are not really independent, as the number of Bohr sets increase, then you can calculate the correlation, thanks to the harmonic sires increasing very slowly, wwe can get something non trivial by this argument, but it seems not enough to cover the whole theorem 1.

    \paragraph{Argument 4, mountain bootstrap argument} This argument is more clever than 3, although both two arguments try to gain the property we want in 1 from investigate the whole space {[0,1]} but not {x}, this argument is more clever.

    Now I explain the main argument, it is nothing but sphere packing, with the set of balls

    \displaystyle \Omega=\{B_{p,q}:=(\frac{q}{p}-\frac{1}{q^2},\frac{q}{p}+\frac{1}{p^2})| \forall p\in {\mathbb N}, 1\leq q\leq p-1 \} \ \ \ \ \ (9)

    and define its subset

    \displaystyle \Omega_l=\{B_{p,q}:=(\frac{q}{p}-\frac{1}{q^2},\frac{q}{p}+\frac{1}{p^2})| \forall 1\leq p\leq l, 1\leq q\leq p-1 \} \ \ \ \ \ (10)

    Then {\Omega_l\subset \Omega}, and {\Omega =\cup_{l\in {\mathbb N}}\Omega_l}. If we can proof,

    Lemma 3 For all {l\in {\mathbb N}}, there is a subset {A_l} of {\Omega-\Omega_l} such that {\cup_{i\in A_l}B_i=[0,1]}.

    Remark 2 If we can proof 3, it is easy to see the theorem 1 follows.

    Proof: The proof follows very standard in analysis, may be complex analysis? Key point is we start with a ball {B_{p,q}}, whatever it is, this is not important, the important thing is we can take some ball {B_{p',q'}} with the center of {B_{p',q'}} in {B_{p,q}}, then try to consider {B_{p',q'}\cup B_{p,q}} to extension {B_{p,q}} and then we find the boudary is also larger then we can extension again, step by step just like mountain bootstrap argument. So we involve in two possible ending,

    1. The extension process could extension {B{p,q}} to whole space.
    2. we can not use the extension argument to extension to the whole space.

    If we are in the first situation, then we are safe, there is nothing need proof. If we are in second case, anyway we take a ball {B_{p,q}=(\frac{q}{p}-\frac{1}{p^2},\frac{q}{p}+\frac{1}{p^2})}. Then try to find good ball {B_{p',q'}} to approximate {B_{p,q}}, but this is difficult… \Box

    Remark 3 Argument 1 is too clever to be true in generalization, argument 2 is standard, by the power of renormalization. argument 3 and argument 4 have gap… I remember I have got a proof similar to argument 4 here many years ago, but I forgot how to get it…

    2. Schimidt conjecture

    The Schimidt conjecture could be look as the generalization of Dirchlet approximation theorem 1 to algebraic number version, to do this, we need define the height of a algebraic number.

    Definition 4 We say a number {\alpha\in {\mathbb C}} is a {k-}order algebraic number if and only is the minimal polynomial of {\alpha}, {f(x)=a_nx^n+...+a_1x+a_0, a_n\neq 0} have degree {deg(f)=n, f\in {\mathbb Z}[x]}.

    Definition 5 (Height) Now we define the height of a {k-}th order algebraic number as {H(\alpha):=\max\{\|a_n\|_h,\|a_{n-1}\|_h,...,\|a_0\|_h\}}, Where

    \displaystyle h(a_i)=\|a_i\|_{\infty} \ \ \ \ \ (11)

    Now we state the conjecture:

    Theorem 6 (Swiring-Schimidt conjecture) For all transendental number {x\in {\mathbb C}}, there is infinitely {\alpha} are {k-}th algebraic number such that:

    \displaystyle |x-\alpha|<\frac{c_k}{H(\alpha)^{k+1}} \ \ \ \ \ (12)

    Where {c_k} is a constant only related to {k} but not {x}.

    I point out the conjecture is very related to the map:

    \displaystyle F:(x_1,...,x_n) \longrightarrow (\sigma_1(x_1,...,x_n),\sigma_2(x_1,...,x_n),...,\sigma_n(x_1,...,x_n)) \ \ \ \ \ (13)

    Where {\sigma_k(x_1,...,x_n)=\sum_{1\leq i_1<...<i_k\leq n}\Pi_{j=1}^kx_{i_1}x_{i_2}...x_{i_k}} is the {k-}th symmetric sum.

    Remark 4 {F} is a map {{\mathbb C}^n\rightarrow {\mathbb C}^n}, what we consider is its inverse, {G=F^{-1}}, but {G} is not smooth, it occur singularity when {x_i=x_j} for some {i\neq j}. And the map, as we know, the singularity depend on the quantity {\Pi_{1\leq i< j\leq n}(x_i-x_j)}.

    Remark 5 I then say something about the geometric behaviour of the map {G}, as we know, what we have in mind is consider the map {G} as a distortion {{\mathbb C}^n\rightarrow {\mathbb C}^n}, Then {H(\alpha)} is just the pullback of the canonical metric on {{\mathbb C}}(morally) to {{\mathbb C}}.


    补充说明

    以下是新整理的中文说明;上方旧博客原文保持不变。

    Diophantine approximation 的基本问题是:一个无理数 $\alpha$ 能被有理数 $p/q$ 逼近到什么程度?Dirichlet 定理给出最基本的答案,而连分数告诉我们哪些分母是真正的最佳逼近。

    Diophantine approximation:Dirichlet 定理、抽屉原理与最佳逼近
    Dirichlet 定理的抽屉原理证明依赖小数部分的差分线性结构。

    1. Dirichlet 定理

    对任意无理数 $\alpha$,存在无穷多个有理数 $p/q$,使得

    $$\left|\alpha-\frac pq\right|<\frac1{q^2}.$$

    一个有限版本是:给定 $Q$,存在 $1\le q\le Q$ 和整数 $p$,使

    $$|q\alpha-p|<\frac1Q.$$

    2. 抽屉原理证明

    看 $Q+1$ 个数的小数部分

    $$0,\{\alpha\},\{2\alpha\},\ldots,\{Q\alpha\}.$$

    把 $[0,1]$ 分成 $Q$ 个长度 $1/Q$ 的区间。两个小数部分落在同一区间,于是它们的差给出某个 $q\alpha$ 距离整数小于 $1/Q$。这就是 Dirichlet 定理最干净的证明。

    3. 线性结构在哪里

    抽屉原理本身只是计数,但这里真正起作用的是线性结构:两个点 $\{a\alpha\}$ 和 $\{b\alpha\}$ 接近,差就变成 $\{(a-b)\alpha\}$ 接近整数。没有这个差分结构,抽屉原理不会自动给出有理逼近。

    4. 连分数与最佳逼近

    连分数展开

    $$\alpha=[a_0;a_1,a_2,\ldots]$$

    给出 convergents $p_k/q_k$。它们满足

    $$\left|\alpha-\frac{p_k}{q_k}\right|<\frac1{q_kq_{k+1}}.$$

    这些分母 $q_k$ 是最佳逼近的自然尺度。若 $a_k$ 有界,则 $\alpha$ 是 badly approximable;若某些 $a_k$ 很大,就会出现异常好的逼近。

    5. 更高维和 Schmidt 猜想的方向

    高维 Diophantine approximation 会把一个数的逼近问题变成向量、线性形式或流形上的逼近问题。此时抽屉原理仍然给出基准结果,但最佳常数、例外集维数和代数数逼近会变得更深。许多问题最后会进入 geometry of numbers、dynamical systems on homogeneous spaces 或 Schmidt game 的语言。

  • 离散调和函数:格点 Laplacian、坏点集合与 Liouville 型问题

    旧博客原文

    原题:Discrete harmonic function in Z^n

    There is some gap, in fact I can improve half of the argument of Discrete harmonic function , the pdf version is Discrete harmonic function in Z^n, but I still have some gap to deal with the residue half…

     

    1. The statement of result

    First of all, we give the definition of discrete harmonic function.

    Definition 1 (Discrete harmonic function) We say a function {f: {\mathbb Z}^n \rightarrow {\mathbb R}} is a discrete harmonic function on {{\mathbb Z}^n} if and only if for any {(x_1,...,x_n)\in {\mathbb Z}^n}, we have:

    \displaystyle f(x_1,...,x_n)=\frac{1}{2^n}\sum_{(\delta_1,...,\delta_n )\in \{-1,1\}^n}f(x_1+\delta_1,...,x_n+\delta_n ) \ \ \ \ \ (1)

     

    In dimension 2, the definition reduce to:

    Definition 2 (Discrete harmonic function in {{\mathbb R}^2}) We say a function {f: {\mathbb Z}^2 \rightarrow {\mathbb R}} is a discrete harmonic function on {{\mathbb Z}^2} if and only if for any {(x_1,x_2)\in {\mathbb Z}^2}, we have:

    \displaystyle f(x_1,x_2)=\frac{1}{4}\sum_{(\delta_1,\delta_2)\in \{-1,1\}^2}f(x_1+\delta_1,x_2+\delta_2 ) \ \ \ \ \ (2)

     

    The result establish in \cite{paper} is following:

    Theorem 3 (Liouville theorem for discrete harmonic functions in {{\mathbb R}^2}) Given {c>0}. There exists a constant {\epsilon>0} related to {c} such that, given a discrete harmonic function {f} in {{\mathbb Z}^2} satisfied for any ball {B_R(x_0)} with radius {R>R_0}, there is {1-\epsilon} portion of points {x\in B_R(x_0)} satisfied {|f(x)|<c}. then {f} is a constant function in {{\mathbb Z}^2}.

    Remark 1 This type of result contradict to the intuition, at least there is no such result in {{\mathbb C}}. For example. the existence of poisson kernel and the example given in \cite{paper} explain the issue.

    Remark 2 There are reasons to explain why there could not have a result in {{\mathbb C}} but in {{\mathbb Z}^2},

    1. The first reason is due to every radius {R} there is only {O(R^2)} lattices in {B_R(x)} in {{\mathbb Z}^2} so the mass could not concentrate very much in this setting.
    2. The second one is due to there do not have infinite scale in {{\mathbb Z}^2} but in {{\mathbb C}}.
    3. The third one is the function in {{\mathbb Z}^2} is automatically locally integrable.

     

    The generation is following:

    Theorem 4 (Liouville theorem for discrete harmonic functions in {{\mathbb R}^n}) Given {c>0,n\in {\mathbb N}}. There exists a constant {\epsilon>0} related to {n,c} such that, given a discrete harmonic function {f} in {{\mathbb Z}^n} satisfied for any ball {B_R(x_0)} with radius {R>R_0}, there is {1-\epsilon} portion of points {x\in B_R(x_0)} satisfied {|f(x)|<c}. then {f} is a constant function in {{\mathbb Z}^n}.

    In this note, I give a proof of 4, and explicit calculate a constant {\epsilon_n>0} satisfied the condition in 3, this way could also calculate a constant {\epsilon_n} satisfied 4. and point the constant calculate in this way is not optimal both in high dimension and 2 dimension.

    2. some element properties with discrete harmonic function

    We warm up with some naive property with discrete harmonic function. The behaviour of bad points could be controlled, just by isoperimetric inequality and maximum principle we have following result.

    Definition 5 (Bad points) We divide points of {{\mathbb Z}^n} into good part and bad part, good part {I} is combine by all point {x} such that {|f(x)|<c}, and {J} is the residue one. So {A\amalg B={\mathbb Z}^n}.

    For all {B_R(0)}, we define {J_R:=J\cap B_R(0), I_R=I\cap B_R(0)} for convenient.

    Theorem 6 (The distribution of bad points) For all bad points {J_R} in {B_R(0)}, they will divide into several connected part, i.e.

    \displaystyle J_R=\amalg_{i\in S_R}A_i \ \ \ \ \ (3)

    and every part {A_i} satisfied {A_i\cap \partial B_R(0)\neq \emptyset}.

    Remark 3 We say {A} is connected in {{\mathbb Z}^n} iff there is a path in {A} connected {x\rightarrow y, \forall x,y\in A}.

    Remark 4 the meaning that every point So the behaviour of bad points are just like a tree structure given in the gragh.

    Proof: A very naive observation is that for all {\Omega\subset {\mathbb Z}^n} is a connected compact domain, then there is a function

    \displaystyle \lambda_{\Omega}: \partial \longrightarrow {\mathbb R} \ \ \ \ \ (4)

    such that {\lambda_{\Omega}(x,y)\geq 0, \forall (x,y)\in \Omega \times \mathring{\Omega}}. And we have:

    \displaystyle f(x)=\sum_{y\in\partial \Omega}\lambda_{\Omega}(x,y)(y) \ \ \ \ \ (5)

     

    This could be proved by induction on the diameter if {\Omega}. Then, if there is a connected component of {\Omega} such that contradict to theorem 6 for simplify assume the connected component is just {\Omega}, then use the formula 5we know

    \displaystyle \begin{array}{rcl} \sup_{x\in \Omega}|f(x)| & = & \sup_{x\in \Omega}\sum_{y\in\partial \Omega}\lambda_{\Omega}(x,y)(y) \\ & \leq & \sup_{\partial \Omega}|f(x)| \\ & \leq & c \end{array}

    The last line is due to consider around {\partial \Omega}. But this lead to: {\forall x\in \Omega, |f(x)|<c} which is contradict to the definition of {\Omega}. So we get the proof. \Box

    Now we begin another observation, that is the freedom of extension of discrete harmonic function in {{\mathbb Z}^n} is limited.

    Theorem 7 we can say something about the structure of harmonic function space of {Z^n}, the cube, you will see, if add one value, then you get every value, i.e. we know the generation space of {Z^n}

    Proof: For two dimension case, the proof is directly induce by the graph. The case of {n} dimensional is similar. \Box

    Remark 5 The generation space is well controlled. In fact is just like n orthogonal direction line in n dimensional case.

    3. sktech of the proof for \ref

    }

    The proof is following, by looking at the following two different lemmas establish by two different ways, and get a contradiction.

    \paragraph{First lemma}

    Lemma 8 (Discrete poisson kernel) the poisson kernel in {{\mathbb Z}^n}. We point out there is a discrete poisson kernel in {{\mathbb Z}^n}, this is given by:

    \displaystyle f(x)=\sum_{y\in \partial B_R(z)}\lambda_{B_R(z)}(x,y)(y) \ \ \ \ \ (6)

    And the following properties is true:

    1. {\lambda_{B_R(z+h)}(x+h,y+h)=\lambda_{B_R(z)}(x,y)} , {\forall x\in \Omega, h\in {\mathbb Z}^n}.
    2. \displaystyle \lambda_{B_R(z)}(x,y)\rightarrow \rho_R(x,y) \ \ \ \ \ (7)

    Remark 6 The proof could establish by central limit theorem, brown motion, see the material in the book of Stein \cite{stein}. The key point why this lemma 8 will be useful for the proof is due to this identity always true {\forall x\in B_R(0)}, So we will gain a lots of identity, These identity carry information which is contract by another argument.

    \paragraph{Second lemma} The exponent decrease of mass.

    Lemma 9 The mass decrease at least for exponent rate.

    Remark 7 the proof reduce to a random walk result and a careful look at level set, reduce to the worst case by brunn-minkowski inequality or isoperimetry inequality.

    \paragraph{Final argument} By looking at lemma 1 and lemma 2, we will get a contradiction by following way, first the value of {f} on {\partial B_R(0)} increasing too fast, exponent increasing by lemma2, but on the other hand, it lie in the integral expresion involve with poisson kernel, but the pertubation of poisson kernel is slow, polynomial rate in fact…

    \newpage

    {99} \bibitem{paper} A DISCRETE HARMONIC FUNCTION BOUNDED ON A LARGE PORTION OF Z2 IS CONSTANT

    \bibitem{stein} Functional analysis

     


    补充说明

    以下是新整理的中文说明;上方旧博客原文保持不变。

    离散调和函数是连续调和函数在格点 $\mathbb Z^n$ 上的版本。它保留了平均值性质、最大值原理和随机游走解释,但也带来新的组合几何问题:如果一个离散调和函数在很多点上都不大,是否能推出它必须是常数?

    离散调和函数:格点 Laplacian、坏点集合与 Liouville 型问题
    离散调和函数的平均值性质可以用随机游走解释,坏点集合的几何受到最大值原理控制。

    1. 定义

    函数 $u:\mathbb Z^n\to\mathbb R$ 称为离散调和,如果对每个 $x\in\mathbb Z^n$,

    $$u(x)=\frac1{2n}\sum_{y\sim x}u(y),$$

    其中 $y\sim x$ 表示 $y$ 与 $x$ 相邻。等价地,离散 Laplacian

    $$\Delta_d u(x)=\sum_{y\sim x}(u(y)-u(x))$$

    满足 $\Delta_d u=0$。

    在二维中,这就是

    $$u(i,j)=\frac14\bigl(u(i+1,j)+u(i-1,j)+u(i,j+1)+u(i,j-1)\bigr).$$

    2. 最大值原理

    离散调和函数满足最大值原理:如果 $u$ 在有限连通区域内部调和,那么最大值和最小值出现在边界上。原因很简单:一个点的值是邻点平均,若内部点达到严格最大值,则所有邻点也必须取同样的值,连通性迫使整个区域常数。

    这使得坏点集合的形状受到限制。设

    $$G=\{x:|u(x)|\le 1\},\qquad B=\mathbb Z^n\setminus G.$$

    如果某个坏点连通块完全被好点包围,那么最大值原理会迫使它不能真正坏。因此坏点必须以某种方式连接到边界或无穷远。

    3. Poisson kernel 与随机游走

    在有限区域 $\Omega\subset\mathbb Z^n$ 上,离散调和函数由边界值决定:

    $$u(x)=\sum_{z\in\partial\Omega}P_\Omega(x,z)u(z).$$

    这里 $P_\Omega(x,z)$ 是从 $x$ 出发的简单随机游走第一次离开 $\Omega$ 时落在 $z$ 的概率。这就是离散 Poisson kernel。

    这个表示把分析问题转成概率问题:若边界上大值点所占比例很小,那么内部点看到大值的概率也会受到控制。

    4. Liouville 型命题

    经典 Liouville theorem 说,有界的整调和函数必须是常数。离散版本也有类似结论。更细的问题是:如果 $u$ 不假设处处有界,但在每个大球里都有固定比例的点满足 $|u|\le 1$,是否仍能推出 $u$ 是常数?

    这类命题的证明通常要比较两个方向。第一,离散调和性和 Poisson 表示迫使内部值由边界平均控制。第二,若存在越来越大的坏点连通块,那么等周不等式或随机游走逃逸概率会给出大值传播。二者冲突时,就只能得到常数解。

    5. 为什么维数和尺度重要

    在 $\mathbb Z^n$ 中,每个半径球只有有限个格点,边界体积和内部体积之间有明确关系。坏点集合若想在所有尺度上保持稀疏,就很难同时支撑一个非平凡调和函数的增长。

    这和连续情形的差别在于,离散空间把局部传播路径变成了组合对象。一个值要从边界影响内部,必须沿随机游走路径进入;而路径数量、逃逸概率和等周结构都会参与估计。

    6. 证明图像

    可以把证明想成两条 lemma 的冲突。一个 lemma 来自 Poisson kernel:内部值是边界值的随机平均,因此不能随意增长。另一个 lemma 来自坏点集合的几何:如果坏点在每个尺度都存在足够结构,它的质量会以某种速率传播。若假设“好点比例”在所有尺度上都足够大,这两种趋势最终矛盾。

    这个问题的有趣之处在于,它把调和分析、随机游走和离散等周不等式绑在一起。离散调和函数不是连续理论的机械翻译,而是一种真正带有格点几何味道的分析对象。

  • 对数平均 Sarnak 猜想:从 BSZ 准则到熵下降

    旧博客原文

    原题:Log average sarnak conjecture

     

    This is a note concentrate on the log average Sarnak conjecture, after the work of Matomaki and Raziwill on the estimate of multiplication function of short interval. Given a overview of the presented tools and method dealing with this conjectue.

     

    1. Introduction

    Sarnak conjecture \cite{Sarnak} assert that for any obersevable {\{f(T^n(x_0))\}_{n=1}^{\infty}} come from a determination systems {(T,X),T:X\rightarrow X}, where {h(T)=0}, {x_0\in X, f\in C(X)}. The correlation of it and the Liuvillou function is 0, i.e. they are orthongonal to each other, more preseicesly it is just to say,

    \displaystyle \sum_{n<x}\mu(x)f(T^n(x_0))=o(x) \ \ \ \ \ (1)

     

    This is a very natural raised conjecture, Liuville function is the presentation of primes, due to we always believe the distribution of primes in {\mathbb N} should be randomness.

    It has been known as observed by Landau \cite{Laudau} that the simplest case,

    \displaystyle \sum_{n<x}\mu(n)=o(x)

    already equivalent to the prime number theorem. It is not difficult to deduce the spetial case of Sarnak conjecture when with the obersevation in $latex {(1)}&fg=000000$ come from finite dynamic system is equivalent to the prime number theorem in athremetic progress by the similar argument. Besides this two classical result, may be the first new result was established by Davenport,

    Theorem 1 Let {T:S_1\rightarrow S_1, T(x)=x+\alpha}, {\alpha} is a inrational, then the obersevation come from {(T,S_1)} is orthogonal to Mobius function. due to {\{e^{2\pi ikx}\}_{k\in \mathbb Z}} is a basis of {C(S_1)}, suffice to proof,

    \displaystyle \sum_{n<x}e^{2\pi ikn\alpha}\mu(n)=o(x), \forall k\in \mathbb N

    There is a lots of spetial situations of Sarnak’s conjecture have been established, The parts I mainly cared is the following:

    1. Interval exchange map.
    2. Skew product flow.
    3. Obersevable come from One dimensional zero entropy flow.
    4. Nilsequences.

    But in this note, I do not want to explain the tecnical and tools to establish this result, but considering an equivalent conjecture of Sarnak conjecture, named Chowla conjecture, and explain the underlying insight of the suitable weak statement, i.e. the log average Chowla conjecture and the underlying insight of it.

    The note is organized as following way, in the next section $latex {(2)}&fg=000000$, we give a self-contained introduction on the tools called Bourgain-Sarnak-Ziegler critation, explain the relationship of this critation and the sum-product phenomenon, also given some more general critation along the philosephy use in establish the Bourgain-Sarnak-Ziegler critation, which maybe useful in following development combine with some other tools. The key point is transform the sum from linear sum to bilinear sum and decomposition the bilinear sum into diagonal part and off-diagonal part, use the assume in the critation to argue the off-diagonal part is small and on the orther hand the diagonal part is also small by the trivial estimate and the volume of diogonal is small, this is very similar to a suitable Caderon-Zugmund decomposition.

    In section $latex {(4)}&fg=000000$, I try to give a proof sketch of the result of Matomaki and Raziwill, which is also a key tools to understanding the Sarnak conjecture, or equivalent the Chowla conjecture. The key points of the proof contains following:

    1. Find a suitable fourier indentity
    2. Construct a multiplication-addition dense subset {S}, and proof that the theorem MR hold we need only to proof it hold for {S\cap [1,2,...,n]} instead of {[1,2,...,n]}
    3. Involve the power of euler product formula. divide the whole interval into a lot of small interval with smaller and smaller scale and a residue part. We look the part come from every small scale as a major term and look the residue part as minor term.
    4. Deal with the major term at every scale, by a combitorios identity and second moments method.
    5. find a enough decay estimate from a scale to the next smaller scale.
    6. Deal with the minor term by the H… lemma.

    Due to the theorem of MR do not exausted the method they developed, we trying to make some more result with their method, Tao and Matomaki attain the average version of Chowla conjecture is true by this way, and combine this argument and the entropy decresment argument they established the 2 partten of the log average Chowla conjecture is true. Very recently Tao and his coperator proved the odd partten case of log average chowla conjecture is true, combine an argument of frustenberg crresponding principle and entopy decresment argument. But it seems the even and large than 2 case is much difficult and seems need something new to combine with the method of MR and entropy decresment and frunstenberg corresponfing principle to make some progress.

    So, in section $latex {(5)}&fg=000000$, we give a self-contain introduction to the entropy decresment argument of Tao, and combine with the frustenberg corresponding principle.

    In the last section $latex {(6)}&fg=000000$, I state some result and method and phylosphy of them I get on nilsequences and wish to combine them with the previous method to make some progress on log average Chowla conjecture on the even partten case.

    \newpage

    2. Bourgain-Sarnak-Zieglar creation

    We begin with the easiest one, this is the main result established in \cite{BSZ}, I try to give the main ideal under the proof, but with a no quantitative version is the following,

    Theorem 2 (Bourgain-Sarnak-Zieglar creation, not quantitative version) if for all primes {p,q>>1} we have:

    \displaystyle \sum_{n=1}^Nf(T^{pn}(x))\overline{f(T^{qn}(x))}=o(N) \ \ \ \ \ (2)

     

    Then for multiplication function {g(n)} we have

    \displaystyle \sum_{n=1}^Ng(n)\overline{ f(T^n(x))}=o(N) \ \ \ \ \ (3)

     

    Remark 1 For simplify we identify {f(T^n(x)):=F(n)}.

    Remark 2

    The idea is following, break the sum into a bilinear one, so, of course, we multiplication it with itself. i.e. we consider to control,

    \displaystyle |\sum_{i=1}^Ng(n)\overline{ F(n)}|^2=\sum_{n=1}^N\sum_{m=1}^Ng(n)g(m)\overline{F(n)F(m)} \ \ \ \ \ (4)

     

    To control 4, we need exhausted the mutiplication property of {g(n)}, we have {g(mn)=g(n)g(m),\forall\ m,n\in {\mathbb N}}. We can not get good estimate for all term,

    \displaystyle g(n)g(m)\overline{F(n)F(m)} \ \ \ \ \ (5)

    The condition in our hand if following,

    \displaystyle \sum_{n=1}^NF(pn)\overline{F(qn)}=o(N), \forall \ p,q\in \mathop{\mathbb P} \ \ \ \ \ (6)

    So, just like the situation of Cotlar-Stein lemma \cite{Cotlar-Stein lemma}, we wish to estimate like following:

    \displaystyle \begin{array}{rcl} |\sum_{p\in W}\sum_{n\in V}F(pn)g(pn)| & \leq & \sum_{n\in V}|g(n)|\cdot |\sum_{p\in W}F(pn)g(p)| \\ & \leq &\sum_{n\in V}|\sum_{p \in W}F(pn)g(p)|\\ & \overset{Cauchy-Schwarz}\leq & |V|^{\frac{1}{2}}[\sum_{n\in V}|\sum_{p\in W}F(pn)g(p)|^2]^{\frac{1}{2}}\\ & = & |V|^{\frac{1}{2}}[\sum_{p_1,p_2\in W}\sum_{n\in V}F(p_1n)\overline{F(p_2n)}g(p_1)\overline{g(p_2)}]^{\frac{1}{2}}\\ \end{array}

    Then we consider divide the sum into diagonal part and non-diagonal part, as following,

    \displaystyle |V|^{\frac{1}{2}}[\sum_{p_1\neq p_2\in W}\sum_{n\in V}F(p_1n)\overline{F(p_2n)}g(p_1)\overline{g(p_2)}]^{\frac{1}{2}}+|V|^{\frac{1}{2}}[\sum_{p\in W}\sum_{n\in V}|F(pn)|^2]^{\frac{1}{2}} \ \ \ \ \ (7)

    But the first part is small, i.e.

    \displaystyle |V|^{\frac{1}{2}}[\sum_{p_1\neq p_2\in W}\sum_{n\in V}F(p_1n)\overline{F(p_2n)}g(p_1)\overline{g(p_2)}]^{\frac{1}{2}} =o(|W||V|) \ \ \ \ \ (8)

    Because of

    \displaystyle \sum_{n\in V}F(p_1n)\overline{F(p_2n)}=o(V), \forall p_1\neq p_2\in W \ \ \ \ \ (9)

    and the second part is small, i.e.

    \displaystyle |V|^{\frac{1}{2}}[\sum_{p\in W}\sum_{n\in V}|F(pn)|^2]^{\frac{1}{2}}=o(|W||V|) \ \ \ \ \ (10)

    Because diagonal part is small in {W\times W} and trivial inequality

    \displaystyle \sqrt{\sum_{n\in V}|F(pn)|^2}\leq |V|^{\frac{1}{2}} \ \ \ \ \ (11)

    But the method in remark 2 is not always make sense in any situation, we need to construct two suitable sets {W,V} and then break up {\{1,2,...,n{\mathbb N}\}} into {W\times V}, this mean,

    \displaystyle \{1,2,...,N\}\sim W\times V+o(N) \ \ \ \ \ (12)

    But this {W,V} could be construct in this situation, thanks to the prime number theorem,

    Theorem 3 (Prime number theorem)

    \displaystyle \pi(n)\sim \frac{n}{ln(n)} \ \ \ \ \ (13)

    Morally speaking, this is the statement that the primes, which is the generator of multiplication function, is not very sparse.

    3. Van der curpurt trick

    There is the statement of Van der carport theorem:

    Theorem 4 (Van der curpurt trick) Given a sequences { \{x_n\}_{n=1}^{\infty}} in { S_1}, if { \forall k\in N^*}, { \{x_{n+k}-x_n\}} is uniformly distributed, then { \{x_n\}_{n=1}^{\infty}} is uniformly distributed.

    I do not know how to establish this theorem with no extra condition, but this result is true at least for polynomial flow. \newpage Proof:

    \displaystyle \begin{array}{rcl} |\sum_{n=1}^Ne^{2\pi imQ(n)}|& = &\sqrt{(\sum_{n=1}^Ne^{2\pi imQ(n)})(\overline{\sum_{n=1}^Ne^{2\pi imQ(n)}})}\\ & = &\sqrt{\sum_{h_1=1}^N\sum_{n=1}^{N-h_1}e^{2\pi imQ(n+h_1)-Q(n)}}\\ & = &\sqrt{\sum_{h_1=1}^N\sum_{n=1}^{N-h_1}e^{2\pi im \partial^1_{h_1}Q(n)}}\\ & \leq & \sqrt{\sum_{h_1=1}^N|\sum_{n=1}^{N-h_1}e^{2\pi \partial^1_{h_1}Q(n)}|}\\ & = &\sqrt{\sum_{h_1=1}^N\sqrt{ (\sum_{n=1}^{N-h_1}e^{2\pi \partial^1_{h_1}Q(n)} )(\overline{\sum_{n=1}^{N-h}e^{2\pi \partial^1_{h_1}Q(n)})}}}\leq\sqrt{\sum_{h_1=1}^N\sqrt{ \sum_{h_2=1}^N|\sum_{n=1}^{N-h_1}e^{2\pi\partial^1_{h_2} \partial^1_hQ(n)} |}}\\ & \leq ....\leq & \\ & = & \sqrt{\sum_{h_1=1}^N\sqrt{ \sum_{h_2=1}^N \sqrt{....\sqrt{\sum_{h_{k-1}=1}^{N-h_{k-2}}|\sum_{n=1}^{N-h_{k-1}}e^{2\pi\partial_{h_1h_2...h_{k-1}Q(n)}}|}}}} =o(1) \end{array}

    \Box

    This type of trick could also establish the following result, which could be understand as a discretization of the Vinegradov lemma.

    Remark 3

    Uniformly distribution result of { F_p}: Given {Q(n)=a_kn^k+...+a_1n+a_0}, {\{Q(0),Q(1),...,Q(p-1)\}} coverages to a uniformly distribution in {\{0,1,...,p-1\}} as {p \rightarrow \infty}.

    Remark 4 But I definitely do not know how to establish the similar result when {Q(n)=n^{-1}}.

    Remark 5

    This trick could also help to establish estimate of correlation of low complexity sequences and multiplicative function, such as result:

    \displaystyle S(x)=\sum_{n\le x}\left(\frac{n}{p}\right)\mu(n)=o(n)

    Maybe with the help of B-Z-S theorem.

    \newpage

    4. Matomaki and Raziwill’s work

    In this section we explain the main idea underlying the paper \cite{KAISA MATOMA 虉KI AND MAKSYM RADZIWILL}. But play with a toy model, i.e. the corresponding corollary of the original result on Liouville鈥檚 function.

    Definition 5 (Lioville’s function)

    \displaystyle \lambda(n)=(-1)^{\alpha_1+\alpha_2+...+\alpha_k}, \forall \ n=p_1^{\alpha_1}...p_k^{\alpha_k}. \ \ \ \ \ (14)

    Remark 6

    \displaystyle |\int_{X}^{2X}\lambda(n)dx|=o(x) \ \ \ \ \ (15)

    is equivalent to the prime number theorem 3.

    The most important beakgrouth of analytic number theory is the new understanding of multiplication function on share interval, this result is established by Kaisa Matom盲ki and Maksym Radziwill. Two very young and intelligent superstars.

    The main theorem in them article is :

    Theorem 6 (Matomaki,Radziwill) As soon as {H\rightarrow \infty} when {x\rightarrow \infty}, one has:

    \displaystyle \sum_{x\leq n\leq x+H}\lambda(n)= o(H) \ \ \ \ \ (16)

    for almost all {1\leq x\leq X} .

    In my understanding of the result, the main strategy is:

    1. Parseval indetity, transform to Dirchelet polynomial.
    2. Involved by multiplication property, spectral decomposition.
    3. From linear to multilinear , Cauchy schwarz inequality.
    4. major term estimate.
    5. Estimate the contribution of area which is not filled.

    4.1. Parseval indetity, transform to Dirchelet polynomial

    We wish to establish the equality,

    \displaystyle \frac{1}{X}\int_{X}^{2X}|\sum_{x\leq n\leq x+H}\lambda(n)|dx=o(H) \ \ \ \ \ (17)

    This is the {L^1} norm, by Chebyschev inequality, this could be control by {L^2} norm, so we only need to establish the following,

    \displaystyle \frac{1}{X}\int_X^{2 X}|\sum_{x\leq n\leq x+H}\lambda(n)|^2dx=o(H^2) \ \ \ \ \ (18)

     

    We wish to transform from the discretization sum to a continue sum, that is,

    \displaystyle \int_{{\mathbb R}}|\sum_{xe^{-\frac{1}{T}}\leq n\leq xe^{\frac{1}{T}}}\lambda(n)1_{X\leq n\leq 2X}|^2\frac{dx}{x} \ \ \ \ \ (19)

     

    Remark 7 There are two points to understand why 19 and 18 are the same.

    1. {[xe^{-\frac{1}{T}},xe^{\frac{1}{T}}]\sim [x-H,x+H]}.
    2. {1_{x\leq n\leq 2x}} and {\frac{1}{x}} is to make that {x=O(X)}.

    So the Magnitude of 18 and 19 are the same. i.e.

    \displaystyle \int_{{\mathbb R}}|\sum_{xe^{-\frac{1}{T}}\leq n\leq xe^{\frac{1}{T}}}\lambda(n)1_{X\leq n\leq 2X}|^2\frac{dx}{x}\sim \frac{1}{X}\int_X^{2 X}|\sum_{x\leq n\leq x+H}\lambda(n)|^2dx \ \ \ \ \ (20)

    Now we try to transform 19 by Parseval indetity, this is something about the {L^2} norms of the quality we wish to charge. It is just trying to understanding 19 as a quantity in physical space by a more chargeable quality in frequency space. Image,

    \displaystyle \int_{{\mathbb R}}|\sum_{xe^{-\frac{1}{T}}\leq n\leq xe^{\frac{1}{T}}}\lambda(n)1_{X\leq n\leq 2X}|^2\frac{dx}{x}:=\int_{{\mathbb R}}|f_X(x)|^2dx \ \ \ \ \ (21)

    Then {f_X(x)=\int_{xe^{-\frac{1}{T}}\leq n\leq xe^{\frac{1}{T}}}\lambda(x)1_{X\leq n\leq 2X}}. Note that,

    \displaystyle \begin{array}{rcl} \widehat{f_X(\xi)} & = & \int_{{\mathbb R}}f_X(x)e^{2\pi ix\xi}dx\\ & = & \sum_{x\leq n\leq 2x}\lambda(x)\int_{logn-\frac{1}{T}}^{logn+\frac{1}{T}}e^{2\pi ix\xi}dx, \ T=\frac{X}{H}\\ & = & \sum_{X\leq n\leq 2X}\lambda(x)e^{2\pi ilog(n)\cdot \xi}\cdot\frac{e^{2\pi i\frac{\xi}{T}}-e^{2\pi i-\frac{\xi}{T}}}{2\pi i\xi}\\ \end{array}

    So by Parseval identity, we have,

    \displaystyle \begin{array}{rcl} \int_{{\mathbb R}}|f_X(x)|^2dx & = & \int_{{\mathbb R}}|\widehat{f_X(\xi)}|^2d\xi \\ & = & \int_{{\mathbb R}}|\sum_{X\leq n\leq 2X}\lambda(n)\cdot n^{2\pi i\xi}|^2(\frac{e^{2\pi i\frac{\xi}{T}}-e^{2\pi i\frac{-\xi}{T}}}{2\pi i\xi})^2d\xi\\ & \sim & \int_{{\mathbb R}}|\sum_{X\leq n\leq 2X}\lambda(n)\cdot n^{2\pi i\xi}|^2\frac{1}{T^2}1_{|\xi|^2\leq T}\\ \end{array}

    Remark 8 We know the Fejer kernel satisfied,

    \displaystyle (\frac{e^{2\pi i\frac{\xi}{T}}-e^{2\pi i\frac{-\xi}{T}}}{2\pi i\xi})^2\sim \frac{1}{T^2}1_{|\xi|\leq T} \ \ \ \ \ (22)

    So morally speaking, we get the following identity.

    \displaystyle \frac{1}{X}\int_{X}^{2X}|\sum_{x\leq n\leq x+H}\lambda(n)|^2dx\sim \frac{1}{(x/H)^2}\int_{0}^{\frac{X}{H}}|\sum_{x\leq n\leq 2x}\lambda(x)x^{2\pi i\xi}|^2d\xi \ \ \ \ \ (23)

    In fact we do a cutoff, the quality we really consider is just:

    \displaystyle \frac{1}{X^2}\int_{|log(X)|^{100}}^{\frac{X}{H}}|\sum_{n\leq X}\lambda(n)n^{it}|^2dt \ \ \ \ \ (24)

    established the monotonically inequality:

    Theorem 7 (Paserval type identity)

    \displaystyle \frac{1}{X}\int_{X}^{2X}|\frac{1}{H}\sum_{x\leq n\leq x+H}\lambda(n)|^2dx \sim聽\frac{1}{X^2}\int_{|log(X)|^{100}}^{\frac{X}{H}}|\sum_{n\leq X}\lambda(n)n^{it}|^2dt \ \ \ \ \ (25)

     

    Remark 9

    In my understanding, This is a perspective of the quality, due to the quality is a multiplicative function integral on a domain { \mathbb N^*} with additive structure, it could be looked as a lots of wave with the periodic given by primes, so we could do a orthogonal decomposition in the fractional space, try to prove the cutoff is a error term and we get such a monotonically inequality.

    But at once we get the monotonically inequality, we could look it as a聽compactification process and this process still carry most of the information so lead to the inequality.

    It seems something similar occur in the attack of the moments estimate of zeta function by the second author. And it is also could be looked as something similar to the 聽spectral decomposition with some basis come from multiplication generators, i.e. primes.

    4.2. Involved by multiplication property, spectral decomposition

    I called it is “spectral decomposition”, but this is not very exact. Anyway, the thing I want to say is that for multiplication function {\lambda(n)}, we have Euler-product formula:

    \displaystyle \Pi_{p,prime}(\frac{1}{1-\frac{\lambda(p)}{p^s}})=\sum_{n=1}^{\infty} \frac{\lambda(n)}{n^s} \ \ \ \ \ (26)

     

    But anyway, we do not use the whole power of multiplication just use it on primes, i.e. {\lambda(pn)=\lambda(p)\lambda(n)} leads to following result:

    \displaystyle \lambda(n)=\sum_{n=pm,p\in I}\frac{\lambda(p)\lambda(m)}{\# \{p|m, p\in I\}+1}+\lambda(n)1_{p|n;p\notin I} \ \ \ \ \ (27)

    This is a identity about the function {\lambda(n)}, the point is it is not just use the multiplication at a point,i.e. {\lambda(mn)=\lambda(m)\lambda(n)}, but take average at a area which is natural generated and compatible with multiplication, this identity carry a lot of information of the multiplicative property. Which is crucial to get a good estimate for the quality we consider about.

    4.3. From linear to multilinear , Cauchy schwarz

    Now, we do not use one sets {I}, but use several sets {I_1,...,I_n } which is carefully chosen. And we do not consider [X,2X] with linear structure anymore , instead reconsider the decomposition:

    {[X,2X]=\amalg_{i=1}^n (I_i\times J_i) \amalg U}

    On every {I_i\times J_i} it equipped with a bilinear structure. And {U} is a very small set, {|U|=o(X)} which is in fact have much better estimate.

    {\int_{|log(X)|^{100}}^{\frac{X}{H}}|\sum_{n\leq X}\lambda(n)n^{it}|^2dt =\sum_{i=1}^n\int_{I_i\times J_i}聽聽\frac{1}{X^2}\int_{|log(X)|^{100}}^{\frac{X}{H}}|\sum_{n\leq X}\lambda(n)n^{it}|^2dt +\int_N |\sum_{n\leq X}\lambda(n)n^{it}|^2dt}

    Now we just use a Cauchy-Schwarz:

    {\sum_{i=1}^n\int_{I_i\times J_i}聽聽\frac{1}{X^2}\int_{|log(X)|^{100}}^{\frac{X}{H}}|\sum_{n\leq X}\lambda(n)n^{it}|^2dt +\int_N |\sum_{n\leq X}\lambda(n)n^{it}|^2dt}

    4.4. major term estimate

    {=\sum_{i=1}^n\int_{I_i\times J_i}聽聽\frac{1}{X^2}\int_{|log(X)|^{100}}^{\frac{X}{H}}|\sum_{n\leq X}\lambda(n)n^{it}|^2dt}

    {\int_N |\sum_{n\leq X}\lambda(n)n^{it}|^2dt}

    4.5. estimate the contribution of area which is not filled

    \newpage

    5. Entropy dcrement argument

    \newpage

    6. Correlation with nilsequences

    I wish to establish the following estimate: {\lambda(n)} is the liouville function we wish the following estimate is true.

    \displaystyle \int_{0\leq x\leq X}|\sup_{f\in \Omega^m}\sum_{x\leq n\leq x+H}\lambda(n)e^{2\pi if(x)}|dx =o(XH). \ \ \ \ \ (28)

    Where we have { H\rightarrow \infty} as { x\rightarrow \infty},

    \displaystyle \Omega^m=\{a_mx^m+a_{m-1}x^{m-1}+...+a_1x+a_0 | a_m,...,a_1,a_0\in [0,1]\}

    is a compact space.

    I do not know how to prove this but this is result is valuable to consider, because by a Fourier identity we could transform the difficulty of (log average) Chowla conjecture to this type of result.

    There is some clue to show this type of result could be true, the first one is the result established by Matomaki and Raziwill in 2015:

    Theorem 8 (multiplication function in short interval)

    {f(n): \mathbb N\rightarrow \mathbb C} is a multiplicative function, i.e. { f(mn)=f(n)f(m), \forall m,n\in \mathbb N}. {H\rightarrow \infty} as {x\rightarrow \infty}, then we have the following result,

    \displaystyle \int_{1\leq x\leq X}|\sum_{x\leq n\leq x+H}f(n)|=o(XH). \ \ \ \ \ (29)

    And there also exists the result which could be established by Vinagrodov estimate and B-S-Z critation :

    Theorem 9 (correlation of multiplication function and nil-sequences in long interval)

    {f(n): \mathbb N\rightarrow \mathbb C} is a multiplicative function, i.e. { f(mn)=f(n)f(m), \forall m,n\in \mathbb N}. {g(n)=a_n^m+...+a_1n+a_0} is a polynomial function then we have the following result,

    \displaystyle \int_{1\leq n \leq X}|f(n)e^{2\pi i g(n)}|=o(X) \ \ \ \ \ (30)

    \newpage {9} \bibitem{Sarnak} Peter Sarnak, Mobius Randomness and Dynamics.

    \texttt{https://publications.ias.edu/sites/default/files/Mahler }. \bibitem{Laudau} JA 虂NOS PINTZ (BUDAPEST). LANDAU鈥橲 PROBLEMS ON PRIMES.

    \texttt{https://users.renyi.hu/~pintz/pjapr.pdf} \bibitem{BSZ} Knuth: Computers and Typesetting,

    \texttt{http://www-cs-faculty.stanford.edu/\~{}uno/abcde.html}

    \bibitem{Cotlar-Stein lemma} Almost orthogonality

    \texttt{https://hxypqr.wordpress.com/2017/12/18/almost-orthogonality/}

    \bibitem{KAISA MATOMA 虉KI AND MAKSYM RADZIWILL} KAISA MATOMA 虉KI AND MAKSYM RADZIWIL, MULTIPLICATIVE FUNCTIONS IN SHORT INTERVALS.

    \texttt{https://arxiv.org/abs/1501.04585v4/}.

     


    补充说明

    以下是新整理的中文说明;上方旧博客原文保持不变。

    这篇笔记想整理的是 Sarnak 猜想的一条现代路线:从莫比乌斯函数和零熵动力系统的正交性出发,经过 Bourgain-Sarnak-Ziegler 准则,把线性相关和转成双线性相关;再借助 Matomaki-Radziwill 的短区间乘法函数估计,以及 Tao 的 entropy decrement,把问题推向对数平均 Chowla 猜想。

    对数平均 Sarnak 猜想:从 BSZ 准则到熵下降
    Sarnak/Chowla 的对数平均路线:线性相关先转成双线性相关,再用短区间估计和熵下降选择合适尺度。

    1. Sarnak 猜想的基本形状

    Sarnak 猜想说,如果 $(X,T)$ 是零拓扑熵动力系统,$f\in C(X)$,那么对任意 $x\in X$,莫比乌斯函数与观测序列 $f(T^n x)$ 应该正交:

    $$\frac1N\sum_{n\le N}\mu(n)f(T^n x)\to 0.$$

    这里的哲学是:$\mu(n)$ 携带素数分布中的振荡,而零熵系统产生的是低复杂度的确定序列;二者不应该长期对齐。

    最简单的有限动力系统情形已经包含素数定理和算术级数中的素数定理。圆周无理旋转情形则接近 Davenport 型估计:对任意无理数 $\alpha$,指数序列 $e(n\alpha)$ 与 $\mu(n)$ 的相关和具有消失。

    2. 从 Sarnak 到 Chowla

    Chowla 猜想更像是对 $\mu$ 或 Liouville 函数自身随机性的表述。一个典型的 $k$ 点相关形式是

    $$\frac1N\sum_{n\le N}\lambda(n+h_1)\cdots \lambda(n+h_k)\to 0,$$

    其中 $h_i$ 两两不同。Sarnak 关心的是乘法函数与低复杂度确定序列的相关,Chowla 关心的是乘法函数自身不同平移之间的相关。二者之间可以用 Furstenberg correspondence principle 和动力系统模型联系起来。

    对数平均版本把普通平均换成

    $$\frac1{\log N}\sum_{n\le N}\frac{a(n)}n.$$

    这个权重让尺度选择更加稳定,也更适合短区间分析。很多情况下,对数平均结论比普通平均结论先被证明,因为它容许把不同尺度上的误差以更柔和的方式叠加。

    3. Bourgain-Sarnak-Ziegler 准则

    BSZ 准则的核心是一个线性到双线性的转换。设 $a_n$ 是有界序列,如果对很多不同素数 $p\ne q$,都有

    $$\frac1N\sum_{n\le N}a_{pn}\overline{a_{qn}}\to 0,$$

    那么可以推出

    $$\frac1N\sum_{n\le N}\mu(n)a_n\to 0.$$

    直观上,$\mu$ 的乘法结构允许我们把原来的线性相关拆成不同素数伸缩后的相关。对角部分 $p=q$ 用平凡估计和体积小来处理;非对角部分则由上面的假设控制。这一点和 Calderon-Zygmund 分解的精神很像:把难对象拆成一个小的对角坏集和一个可估计的非对角主体。

    4. Matomaki-Radziwill 的短区间输入

    Matomaki-Radziwill 的工作说明,乘法函数在大多数短区间里的平均行为可以被控制。粗略地说,对许多短区间 $[x,x+H]$,有

    $$\frac1H\sum_{x

    接近它在长区间中的平均。这使得我们可以把一个全局相关和切成许多短尺度,再用欧拉乘积、二阶矩和组合恒等式逐层比较。

    短区间估计真正有用的地方在于:它让乘法函数的局部随机性可以被拿来攻击动力系统中的相关问题,而不是只停留在平均阶的数论命题。

    5. Entropy decrement 的作用

    entropy decrement 的思想是寻找一个尺度,使得随机变量 $n$ 与它的素数倍 $pn$ 之间的条件信息变少。换句话说,在合适尺度上,系统看到的结构不会因为乘一个小素数而增加太多复杂度。

    这一步的意义是把“乘法平移”转成“动力系统里可比较的两个观测”。当这个信息损失足够小的时候,短区间估计、BSZ 双线性结构和 correspondence principle 就可以接上。

    6. Nilsystem 方向

    nilsequence 是零熵系统中非常重要的一类模型。它既有足够丰富的几何结构,又保留了可计算的 Fourier 分析。把 Sarnak、Chowla、短区间乘法函数和 nilsequence 放在一起看,真正的问题是:乘法函数的随机性如何穿过 nilmanifold 上的低复杂度轨道。

    这条路线目前最有价值的地方,不是把所有情形一次性解决,而是提供了一张方法地图:线性相关转双线性相关;短区间估计提供局部随机性;entropy decrement 选择尺度;Furstenberg 原理把数论相关放回动力系统。

  • Laplace 谱、Rayleigh quotient 与无限维 Morse 图像

    旧博客原文

    原题:Two stupid question

    The story of the infinite dimensional space of $\Delta$ is following, we eliminate ourself with compact smooth non-boundary manifold $M$ with metric $g$, then we have Betrami-Laplace operator $\Delta_g$. We could instead $\Delta_g$ by hodge laplace $dd^*+d^*d$, but let we consider $\Delta_g$ the eigenvalue problem:
    $$\Delta_g u=\lambda u$$
    A classical way to investigate the eigenvalue problem is according to consider variational principle and max-min principle. We equip the path integral on the function space $C^{\infty}(M)$:
    $$E(f)=\frac{\int_M |\nabla u|^2}{\int_M |u|^2 }$$
    Then it have a sequences of eigenvalue, negative of course: $$0<-\lambda_1<-\lambda_2<…<\lambda_k<…$$

    Then things become interesting, the morse theory of infinite space involve, called the infinite space as $X$, so at least, shrink the far place of $X$ as a point, in physics, this mean, cut off at fix scale. And we can take the scale to infinite small, we use the cutoff one to approximation the real one. What I can do is the following, I can proof the eigenvalue function is uniformly distributed in $L^2(M_g)$ (after rescaling of course) and the classical weyl law(although can not give a good error term estimate), but thing become more complicated when I try to consider the infinite space $X_{M_g}$’s topology, at finite scale at least, i.e. $X_{M_g}^{h}$ which is the cut off at scale $h$. Among the other thing, I believe the following issue is true, but without ability to proof it:

    >**Problem**
    for every manifold $M$ and metric $g$ on $M$, the topology of infinite space $X_{M_g}$ is the same, beside this, the inverse could be true, i.e. If $X_{M_1},X_{M_2}$ is not homomorphism for some scale $h$ then $M_1,M_2$ is not homomorphism.

     

     

     

     

     

    By intuition, I think it is depend by the underling manifold’s topology. But I do not have a rigorous proof, I definitely have a non-rigorous one, if ignore the coverage…

    As I find this problem when I try to give a proof of weyl law, I do not check the reference, may be this problem is a classical one? As always, I will appreciate to any interesting comments and answers, thanks a lots!A

     

     

     

     

     

     

     

    2.

    We begin with our favorite situation, the Dirchlet problem on bounded simple-connected domain $\Omega$ in $\mathbb R^n$. Let $\lambda_1$ be the first eigenvalue of $$\Delta u=\lambda u \ in\ \Omega$$
    $$u=0\ \ on\ \partial\Omega$$
    Rescaling $u$ such that $\sup_{\Omega} u=1$, I think the following property of the first eigenvalue is true.
    >**Problem**
    We have, the Minkowski functional of $\Omega$, called $M_{\Omega}$ and the Minkowski functional with the ball $B$ such that $vol(B)=vol(\Omega)$, then along the level set of $u$, i.e. the fiber: $$\Omega=\cup_{t\in [0,1]}l_t, l_t:=\{t|x\in \Omega, u(x)=t\}$$
    We pretend for the isolate point $l_1$ to be a ball with radius 0, so equipped it with the uniformly density at every direction in $S^1$, i.e. the mass distribution given by $M_B$ and the total mass coincide with the total mass induce by $M_{\Omega}$ in $l_0$, i.e.
    $$\int_{e\in S_1}M_{\Omega}(e)d\mu=\int_{e\in S_1}M_{B}(e)d\mu$$
    The measure $d\mu$ equipped on $S^1$ is the natural Haar measure. And the cost function is given by $c(x,y)=\|x-y\|^2$. Then, among this setting,
    I wish the following property to be true:
    Along the direction $1\to 0$, the transport of density $\partial_{t_0} M_{\cup_{t=t_0}^1l_t}$ given the unique optimal transport of the natural measure induce by $M(\Omega)$ and $M(B)$.

    **Remark 1** As point out by SebastianGoette, the multiplicity of the first eigenvalue must be one, thanks to the eigenfunction never change the symbol, so we are in the best case.

    **Remark 2**:I am not very sure this property could always true, there may be a center example when $\Omega$ is not convex, but I tend to believe it is true at least when $\Omega$ is convex.

    **Remark 3**: As point out by Dirk, when you try to consider the optimal transport problem, you always need to point out the cost function $c(x,y)$ defined on $\Omega \times \Omega$, for there, I think the naive choice is $c(x,y)=\|x-y\|^2$

     

    The thing I can proof is the following, the level set of $u$ should be convex by brunn-minkowski inequality, and some type of monotonically property, i.e. more and more like a ball when the level set is more and more shirking smaller form $\partial \Omega$ to the point $f$ arrive maximum.

    I will appreciate for any relevant comments and answer, thanks!

     

     


    补充说明

    以下是新整理的中文说明;上方旧博客原文保持不变。

    Laplace 特征值问题可以看成无限维函数空间上的变分问题。Rayleigh quotient 的临界点给出 eigenfunctions,min-max principle 给出 eigenvalues。这种图像很像 Morse theory,只是空间变成了无限维。

    Laplace 谱、Rayleigh quotient 与无限维 Morse 图像
    Laplace 特征值可以由 Rayleigh quotient 的 min-max 原理得到,形成无限维 Morse 图像。

    1. Laplace 特征值问题

    在紧无边界 Riemannian manifold 上,考虑

    $$-\Delta_g u=\lambda u.$$

    谱离散,特征值可以排成

    $$0=\lambda_0<\lambda_1\le\lambda_2\le\cdots\to\infty.$$

    2. Rayleigh quotient

    定义

    $$E(u)=\frac{\int_M|\nabla u|^2\,dV_g}{\int_M|u|^2\,dV_g}.$$

    在 $L^2$ 单位球上,$E$ 的临界点正是 Laplace eigenfunctions。第一非零特征值是与常数正交的函数中 $E$ 的最小值。

    3. Min-max principle

    第 $k$ 个特征值可由

    $$\lambda_k=\inf_{\dim V=k+1}\sup_{u\in V\setminus\{0\}}E(u)$$

    给出。这个公式把谱问题转成函数空间中的拓扑/变分问题。

    4. 无限维 Morse 图像

    在有限维 Morse theory 中,临界点和拓扑变化有关。Laplace 特征函数也可以看成 Rayleigh quotient 的临界点;其 Morse index 与低于该特征值的谱空间维数相关。

    5. 与 Weyl law 的关系

    Weyl law 告诉我们高能临界点的数量增长率。变分图像解释 eigenfunctions 从哪里来,Weyl law 则解释这些临界点在高频区域有多密。