Fourier-Poisson duality is a framework that exchanges spatial functions and their frequency counterparts through lattice summation and dual measures.
It extends classical Poisson summation to discrete-support and nonperiodic settings, offering new insights in spectral analysis, tiling, and diffraction.
The theory generalizes to weighted, geometric, and quantum formulations, providing computational tools for solving PDEs and exploring quantum group dynamics.
Fourier-Poisson duality is a family of correspondences in which Fourier transform exchanges a spatial object and a frequency-side object linked by lattice summation, periodicity, or a dual measure. Its classical prototype is the Poisson Summation Formula (PSF), which in one dimension may be written as
n∈Z∑f(n)=n∈Z∑f(n),f(ξ)=∫Re−2πiξxf(x)dx,
or distributionally as
(n∈Z∑δn)∧=n∈Z∑δn.
In the literature considered here, that template reappears in several distinct but related forms: discrete-support Fourier pairs on R, formal duality for periodic configurations, diffraction of model sets, spectral-pair measures for unbounded tilings, weighted and geometric Fourier transforms, symbol-level formulations of the Poisson equation, and local Fourier dualities in quantum and Poisson-Lie settings (Kolountzakis, 2015, Cohn et al., 2013, Chakraborty et al., 27 Oct 2025, Tudoran, 2024, Luu et al., 2015, Massar, 1 Dec 2025).
1. Classical lattice duality and the Poisson-summation template
and Λ∗={y:⟨x,y⟩∈Z for all x∈Λ}. In one dimension, the identity
(n∈Z∑δn)∧=n∈Z∑δn
gives a self-dual discrete measure, and by translating, modulating, and dilating this identity one obtains many other discrete-support Fourier pairs (Kolountzakis, 2015).
For radial potentials, the Gaussian is the standard analytic model. If
Gc(r)=e−πcr2,
then
Gc(r)=c−n/2e−πr2/c,
so the parameter transformation (n∈Z∑δn)∧=n∈Z∑δn.0 is built directly into Fourier transform. In this form, Poisson summation expresses a relation between a lattice and its dual lattice, and also between the energies of dual periodic configurations (Cohn et al., 2013).
The significance of this prototype is its rigidity. The PSF generates Fourier pairs whose supports inherit lattice or periodic structure, and much of the later literature can be read either as an extension of this pattern or as a demonstration that the PSF-generated class is not exhaustive.
A central nonclassical result is the existence of a translation-bounded measure
(n∈Z∑δn)∧=n∈Z∑δn.1
with (n∈Z∑δn)∧=n∈Z∑δn.2 discrete, such that
(n∈Z∑δn)∧=n∈Z∑δn.3
with (n∈Z∑δn)∧=n∈Z∑δn.4 discrete, while neither (n∈Z∑δn)∧=n∈Z∑δn.5 nor (n∈Z∑δn)∧=n∈Z∑δn.6 is contained in a finite union of arithmetic progressions. Here “translation bounded” means that the total variation in any unit interval is uniformly bounded, and “discrete” means that every bounded interval contains only finitely many support points. Since measures obtained from finitely many applications of the PSF, together with translation, modulation, and dilation, have support contained in a finite union of arithmetic progressions, this pair cannot be derived from finitely many applications of the PSF (Kolountzakis, 2015).
The construction proceeds in two stages. First, for any prescribed (n∈Z∑δn)∧=n∈Z∑δn.7, one builds a nonzero Fourier pair (n∈Z∑δn)∧=n∈Z∑δn.8 of discrete-support measures such that both (n∈Z∑δn)∧=n∈Z∑δn.9 and R0 vanish on R1. The starting point is a finite cyclic group R2, where one chooses a nonzero function R3 vanishing on
Second, one combines infinitely many such periodic pairs with increasing gaps. Choosing Λ⊂Rn5, a sequence Λ⊂Rn6 that is linearly independent over Λ⊂Rn7, and normalizing factors Λ⊂Rn8, one defines
Λ⊂Rn9
The weights ensure convergence and translation boundedness, the expanding zero neighborhoods ensure discreteness of the supports, and the $\sum_{x \in \Lambda} f(x) = \frac{1}{\covol(\Lambda)} \sum_{y \in \Lambda^*} \widehat f(y),$0-linearly independent shifts obstruct containment in any finite union of arithmetic progressions (Kolountzakis, 2015).
This result broadens Fourier-Poisson duality decisively. It shows that a measure and its Fourier transform can both be purely atomic without arising from any finite PSF construction and without possessing global periodic structure.
for every Schwartz function $\sum_{x \in \Lambda} f(x) = \frac{1}{\covol(\Lambda)} \sum_{y \in \Lambda^*} \widehat f(y),$6. In the finite-group reformulation, with $\sum_{x \in \Lambda} f(x) = \frac{1}{\covol(\Lambda)} \sum_{y \in \Lambda^*} \widehat f(y),$7 for a finite abelian group $\sum_{x \in \Lambda} f(x) = \frac{1}{\covol(\Lambda)} \sum_{y \in \Lambda^*} \widehat f(y),$8, formal duality is equivalent to
This makes formal duality a combinatorial statement about difference sets and character sums. The paper gives concrete examples such as TITO, with f(y)=∫Rnf(x)e−2πi⟨x,y⟩dx1 and f(y)=∫Rnf(x)e−2πi⟨x,y⟩dx2, and a family arising from quadratic Gauss sums in f(y)=∫Rnf(x)e−2πi⟨x,y⟩dx3 (Cohn et al., 2013).
A different but related realization occurs in cut-and-project theory. For a lattice f(y)=∫Rnf(x)e−2πi⟨x,y⟩dx4 in a cut-and-project scheme and a weight f(y)=∫Rnf(x)e−2πi⟨x,y⟩dx5, the weighted Dirac comb is
f(y)=∫Rnf(x)e−2πi⟨x,y⟩dx6
The generalized PSF states that, under the transformability hypotheses of the paper,
f(y)=∫Rnf(x)e−2πi⟨x,y⟩dx7
For Riemann integrable f(y)=∫Rnf(x)e−2πi⟨x,y⟩dx8, the autocorrelation is
f(y)=∫Rnf(x)e−2πi⟨x,y⟩dx9
and the diffraction measure is
Λ∗={y:⟨x,y⟩∈Z for all x∈Λ}0
For regular model sets, this yields pure point diffraction. The paper proves that the diffraction formula for regular model sets is equivalent to the PSF for the underlying lattice (Richard et al., 2015).
Taken together, these results shift Fourier-Poisson duality away from the single question of “lattice versus dual lattice.” In formal duality it becomes a relation between pair correlations and Fourier magnitudes; in model-set diffraction it becomes a statement that Bragg structure is precisely the PSF of a hidden higher-dimensional lattice.
4. Spectral pairs and unbounded tiling sets
In the spectral-pair framework, a measurable set Λ∗={y:⟨x,y⟩∈Z for all x∈Λ}1 and a positive Radon measure Λ∗={y:⟨x,y⟩∈Z for all x∈Λ}2 form a spectral pair if the Fourier transform gives an isometric isomorphism
Λ∗={y:⟨x,y⟩∈Z for all x∈Λ}3
When Λ∗={y:⟨x,y⟩∈Z for all x∈Λ}4 has finite measure, this reduces to the classical notion of a spectrumΛ∗={y:⟨x,y⟩∈Z for all x∈Λ}5, where Λ∗={y:⟨x,y⟩∈Z for all x∈Λ}6 is counting measure on Λ∗={y:⟨x,y⟩∈Z for all x∈Λ}7 and Λ∗={y:⟨x,y⟩∈Z for all x∈Λ}8 is an orthogonal basis for Λ∗={y:⟨x,y⟩∈Z for all x∈Λ}9, with
A concrete one-dimensional theorem identifies a tiling condition with an exact dual measure. If (n∈Z∑δn)∧=n∈Z∑δn2 is open and (n∈Z∑δn)∧=n∈Z∑δn3, then
(n∈Z∑δn)∧=n∈Z∑δn4
if and only if (n∈Z∑δn)∧=n∈Z∑δn5 is spectral with pair measure
(n∈Z∑δn)∧=n∈Z∑δn6
The corresponding Plancherel identity is
(n∈Z∑δn)∧=n∈Z∑δn7
for (n∈Z∑δn)∧=n∈Z∑δn8, and the Fourier transform is onto: (n∈Z∑δn)∧=n∈Z∑δn9
The proof constructs finite approximants
The reverse implication derives structural constraints from the spectral side. For Gc(r)=e−πcr2,1, the periodization
Gc(r)=e−πcr2,2
has period Gc(r)=e−πcr2,3, which implies
Gc(r)=e−πcr2,4
and therefore
Gc(r)=e−πcr2,5
In this formulation, Fourier-Poisson duality is not between two discrete sets. It is a duality between an unbounded tiling set and a periodic Lebesgue measure, with Fourier transform acting as an onto isometry between the two Hilbert spaces (Chakraborty et al., 27 Oct 2025).
5. Weighted, geometric, and algebraic generalizations
One line of generalization replaces the classical Fourier transform by a weighted operator built from dilation sums. For a sequence Gc(r)=e−πcr2,6,
Gc(r)=e−πcr2,7
and with convolution inverse Gc(r)=e−πcr2,8, the Fourier-Poisson operator is defined by
Gc(r)=e−πcr2,9
Under the decay hypothesis
Gc(r)=c−n/2e−πr2/c,0
for some Gc(r)=c−n/2e−πr2/c,1, and the same condition for Gc(r)=c−n/2e−πr2/c,2, the paper proves that Gc(r)=c−n/2e−πr2/c,3 is unitary and is uniquely determined by
Gc(r)=c−n/2e−πr2/c,4
After logarithmic-Fourier conjugation it becomes multiplication by
Gc(r)=c−n/2e−πr2/c,5
so the operator is diagonalized by a ratio of Dirichlet series on the critical line (Faifman, 2011).
A second line starts from a nondegenerate bilinear form Gc(r)=c−n/2e−πr2/c,6 on Gc(r)=c−n/2e−πr2/c,7. The associated left and right geometric Fourier transforms are
Gc(r)=c−n/2e−πr2/c,8
For the geometric pair Gc(r)=c−n/2e−πr2/c,9, they satisfy
(n∈Z∑δn)∧=n∈Z∑δn.00
and they are exactly the transforms that place Poisson summation for the full lattices (n∈Z∑δn)∧=n∈Z∑δn.01 and (n∈Z∑δn)∧=n∈Z∑δn.02 into the canonical forms
(n∈Z∑δn)∧=n∈Z∑δn.03
(n∈Z∑δn)∧=n∈Z∑δn.04
When (n∈Z∑δn)∧=n∈Z∑δn.05 is positive definite, the same framework yields a geometric fractional Laplacian defined by the multiplier (n∈Z∑δn)∧=n∈Z∑δn.06 (Tudoran, 2024).
Other papers isolate algebraic aspects of the same theme. One interprets the continuous Fourier transform as a limiting CRT decomposition: sampled values of (n∈Z∑δn)∧=n∈Z∑δn.07 are packaged into a polynomial (n∈Z∑δn)∧=n∈Z∑δn.08, sampled values of (n∈Z∑δn)∧=n∈Z∑δn.09 appear as residues modulo the factors of (n∈Z∑δn)∧=n∈Z∑δn.10, and the Poisson summation formula
(n∈Z∑δn)∧=n∈Z∑δn.11
emerges from the finite cyclic structure in the limit (n∈Z∑δn)∧=n∈Z∑δn.12 (Xu, 2018). Another shows that for the Gaussian (n∈Z∑δn)∧=n∈Z∑δn.13, Poisson summation and cosine Fourier expansion on (n∈Z∑δn)∧=n∈Z∑δn.14 are equivalent once the image terms are negligible for sufficiently large (n∈Z∑δn)∧=n∈Z∑δn.15, specifically (n∈Z∑δn)∧=n∈Z∑δn.16 (Abrarov et al., 2012).
These generalizations preserve the defining feature of Fourier-Poisson duality: an operator or summation identity exchanges two structures that are not identical on the nose but are tied by a transform law, often with an explicit reconstruction mechanism.
6. Poisson equations, symbols, and computational formulations
In numerical PDE analysis, the phrase denotes the diagonalization of the Poisson operator in frequency space. For
(n∈Z∑δn)∧=n∈Z∑δn.17
the Fourier differentiation rule gives
(n∈Z∑δn)∧=n∈Z∑δn.18
hence
(n∈Z∑δn)∧=n∈Z∑δn.19
and therefore
(n∈Z∑δn)∧=n∈Z∑δn.20
The singularity at (n∈Z∑δn)∧=n∈Z∑δn.21 requires the solvability correction
(n∈Z∑δn)∧=n∈Z∑δn.22
so that (n∈Z∑δn)∧=n∈Z∑δn.23. In the discrete setting, the DFT
(n∈Z∑δn)∧=n∈Z∑δn.24
and its FFT implementation reduce the cost from (n∈Z∑δn)∧=n∈Z∑δn.25 to (n∈Z∑δn)∧=n∈Z∑δn.26. The method is naturally suited to periodic grids and is constrained by the Nyquist frequency
A more recent symbol-level formulation removes the need for an explicit Green’s function. For the scalar Poisson equation, the inverse symbol is
(n∈Z∑δn)∧=n∈Z∑δn.28
and the key identity is
(n∈Z∑δn)∧=n∈Z∑δn.29
Splitting the time integral at (n∈Z∑δn)∧=n∈Z∑δn.30 gives
(n∈Z∑δn)∧=n∈Z∑δn.31
so the solution decomposes into a localized part (n∈Z∑δn)∧=n∈Z∑δn.32 and a smooth nonlocal history part (n∈Z∑δn)∧=n∈Z∑δn.33. The localized component admits asymptotic expansions in powers of (n∈Z∑δn)∧=n∈Z∑δn.34 for volume, single-layer, and double-layer potentials, with coefficients expressed in terms of local geometry and derivatives of the source data. The same Fourier-domain construction extends to coupled strongly elliptic systems by replacing (n∈Z∑δn)∧=n∈Z∑δn.35 with a matrix symbol (n∈Z∑δn)∧=n∈Z∑δn.36 and using (n∈Z∑δn)∧=n∈Z∑δn.37 (Fryklund, 13 Apr 2026).
In this computational setting, Fourier-Poisson duality is no longer primarily a statement about dual supports or lattices. It is the statement that the inverse Poisson operator becomes an algebraic multiplier in Fourier space, and that this multiplier can be regularized, discretized, or asymptotically decomposed while remaining entirely at the level of symbols.
7. Local Fourier duality in quantum curves and Poisson-Lie quantum groups
For quantum curves, the basic object is a pair of ordinary differential operators
(n∈Z∑δn)∧=n∈Z∑δn.38
with bi-degree (n∈Z∑δn)∧=n∈Z∑δn.39. There are two natural KP deformations: one obtained by normalizing (n∈Z∑δn)∧=n∈Z∑δn.40, the other by normalizing (n∈Z∑δn)∧=n∈Z∑δn.41. The paper associates to the quantum curve a companion-matrix connection (n∈Z∑δn)∧=n∈Z∑δn.42 on the formal punctured disc and proves that the two KP orbits are related by local Fourier transform. If (n∈Z∑δn)∧=n∈Z∑δn.43 are coprime and both (n∈Z∑δn)∧=n∈Z∑δn.44 are normalized, then
(n∈Z∑δn)∧=n∈Z∑δn.45
where the new KP times are determined by compositional inversion relations. In 2D gravity, this yields a conceptual proof of the (n∈Z∑δn)∧=n∈Z∑δn.46 duality: the local Fourier transform exchanges the connections attached to the two models, and therefore exchanges the corresponding KP data (Luu et al., 2015). A closely related treatment realizes the (n∈Z∑δn)∧=n∈Z∑δn.47 duality of 2D gravity as a local Fourier duality of irregular (n∈Z∑δn)∧=n∈Z∑δn.48-modules attached to Kac-Schwarz operators (Luu, 2014).
In Poisson-Lie and quantum-group theory, the duality takes yet another form. For semidirect-product Poisson-Lie groups
(n∈Z∑δn)∧=n∈Z∑δn.49
with (n∈Z∑δn)∧=n∈Z∑δn.50 abelian, suitable abelian approximations
(n∈Z∑δn)∧=n∈Z∑δn.51
allow the construction of an explicit unitary operator
(n∈Z∑δn)∧=n∈Z∑δn.52
the Poisson-Fourier transform. It implements the quantum duality principle in the form
(n∈Z∑δn)∧=n∈Z∑δn.53
by inducing an isomorphism of locally compact quantum groups (Massar, 1 Dec 2025). In a related QP-manifold formulation of Poisson-Lie T-duality, canonical transformations on a doubled correspondence space lead to a proposed Fourier-Mukai-type integral transform with kernel (n∈Z∑δn)∧=n∈Z∑δn.54, where
(n∈Z∑δn)∧=n∈Z∑δn.55
so that the kernel intertwines the twisted differential complexes on the dual sides (Arvanitakis et al., 2021).
These quantum-geometric formulations retain the recognizable Fourier-Poisson pattern—an explicit transform exchanges dual data—but the objects are now connections, KP orbits, bicrossed products, or twisted flux complexes rather than functions on (n∈Z∑δn)∧=n∈Z∑δn.56. The common structure is that Fourier transform, or a Fourier-type operator, is the mechanism that converts one side of the duality into the other.