Papers
Topics
Authors
Recent
Search
2000 character limit reached

Fibonacci Potential: Structure & Spectral Features

Updated 8 July 2026
  • Fibonacci potential is generated by concatenating elementary blocks following the Fibonacci substitution rule, resulting in a nonperiodic, quasi-crystalline order.
  • Spectral analysis reveals Cantor-like sets, mini-gap hierarchies, and scale-recursive fractal structures through transfer-matrix and trace map techniques.
  • The concept is versatile, featuring in continuum Schrödinger operators, effective polariton models, and golden-periodic complex potentials in planar hydrodynamics.

Searching arXiv for papers on Fibonacci potentials and the provided IDs. arxiv_search(query="Fibonacci potential Schrödinger Fibonacci continuum spectrum polariton", max_results=10, sort_by="relevance") arxiv_search(query="(Fillman et al., 2017)", max_results=5, sort_by="relevance") arxiv_search(query="(Tanese et al., 2013)", max_results=5, sort_by="relevance") arxiv_search(query="(Pashaev, 2021)", max_results=5, sort_by="relevance") A Fibonacci potential is a potential generated by concatenating two elementary constituents according to the Fibonacci substitution rule S(a)=abS(a)=ab, S(b)=aS(b)=a. In the cited literature, this construction appears in several technically distinct settings: continuum Schrödinger operators on R\mathbb{R} whose local potential pieces are arranged by the Fibonacci subshift, effective one-dimensional potentials for polariton gases in Fibonacci quasi-periodic cavities, and golden-periodic complex potentials arising from the method of images in planar hydrodynamics (Fillman et al., 2017, Tanese et al., 2013, Pashaev, 2021). Across these settings, the common structural feature is substitutional order without periodicity, which yields Cantor-like spectra, mini-gap hierarchies, and scale-recursive analytic structures.

1. Substitutional definition and geometric realization

The underlying symbolic construction starts from a two-letter alphabet. In the continuum Schrödinger setting, one takes A={a,b}A=\{a,b\}, fixes two positive lengths la,lb>0l_a,l_b>0, and chooses two building-block potentials

faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).

The Fibonacci substitution is

S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.

Starting from aa and iterating, one obtains the one-sided fixed point

u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.

The associated two-sided Fibonacci subshift ΩFAZ\Omega_F\subset A^{\mathbb Z} consists of all bi-infinite words whose every finite subword appears somewhere in S(b)=aS(b)=a0. The dynamical system S(b)=aS(b)=a1 is strictly ergodic and minimal (Fillman et al., 2017).

For S(b)=aS(b)=a2, the real-line potential is defined by concatenation. Writing

S(b)=aS(b)=a3

one sets

S(b)=aS(b)=a4

An equivalent box notation is

S(b)=aS(b)=a5

where the box marks the location of the origin S(b)=aS(b)=a6 (Fillman et al., 2017).

A related, but not identical, realization appears in the polariton work. There the letters are denoted S(b)=aS(b)=a7 and S(b)=aS(b)=a8, with substitution rules S(b)=aS(b)=a9, R\mathbb{R}0, or equivalently

R\mathbb{R}1

The length of R\mathbb{R}2 is the Fibonacci number R\mathbb{R}3, satisfying R\mathbb{R}4, and R\mathbb{R}5. In the discrete description, the R\mathbb{R}6th segment of length R\mathbb{R}7 carries potential

R\mathbb{R}8

In the continuous reduction used for wire-width modulation, the effective one-dimensional potential is

R\mathbb{R}9

where A={a,b}A=\{a,b\}0 alternates between A={a,b}A=\{a,b\}1 and A={a,b}A=\{a,b\}2 in Fibonacci order (Tanese et al., 2013).

These constructions show that a Fibonacci potential is not restricted to a single formalism. It may be implemented as a continuum concatenation of compactly supported pieces or as an effective quasiperiodic profile derived from geometry, but in both cases the symbolic order is generated by the same inflation rule.

2. Continuum Fibonacci Schrödinger operators

For each A={a,b}A=\{a,b\}3, the continuum Fibonacci Schrödinger operator acts on A={a,b}A=\{a,b\}4 as

A={a,b}A=\{a,b\}5

Because each A={a,b}A=\{a,b\}6, the maximal operator associated with the differential expression A={a,b}A=\{a,b\}7 with domain A={a,b}A=\{a,b\}8 is essentially self-adjoint. Hence

A={a,b}A=\{a,b\}9

and no boundary-condition choices are needed at la,lb>0l_a,l_b>00 (Fillman et al., 2017).

Minimality and unique ergodicity imply the existence of a compact set la,lb>0l_a,l_b>01, called the common spectrum, such that

la,lb>0l_a,l_b>02

Moreover, la,lb>0l_a,l_b>03 is a Cantor set of zero Lebesgue measure whenever la,lb>0l_a,l_b>04 is aperiodic, meaning not both periodic combinations (Fillman et al., 2017).

The polariton model is formulated through an effective one-dimensional Schrödinger equation for the lower-polariton field la,lb>0l_a,l_b>05: la,lb>0l_a,l_b>06 where la,lb>0l_a,l_b>07 is the polariton effective mass la,lb>0l_a,l_b>08, la,lb>0l_a,l_b>09 is the long axis of the wire, and Dirichlet boundary conditions faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).0 mimic the vanishing field at the etched edges (Tanese et al., 2013).

In the infinite Fibonacci limit, the spectral type in that setting is described as neither purely continuous nor pure point but singular continuous, that is, a Cantor-like set of zero total Lebesgue measure (Tanese et al., 2013). A common misconception is that “almost full-dimensional” fractal spectra should become interval-like. The continuum theory shows otherwise: for aperiodic building blocks the spectrum remains a Cantor set of zero Lebesgue measure, even though the local Hausdorff dimension approaches faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).1 in specific asymptotic regimes (Fillman et al., 2017).

3. Trace map, invariant surfaces, and spectral characterization

The fine structure of the continuum spectrum is analyzed through transfer matrices. For each symbol in faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).2, one considers the fundamental solutions of

faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).3

on the appropriate interval, with Dirichlet or Neumann data at faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).4. Denoting them by faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).5 and faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).6, the faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).7 transfer matrix is

faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).8

Since faL2([0,la)),fbL2([0,lb)).f_a\in L^2([0,l_a)),\qquad f_b\in L^2([0,l_b)).9, one has

S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.0

This defines the curve of initial conditions

S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.1

for the trace map

S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.2

(Fillman et al., 2017).

The trace map preserves the Fricke–Vogt invariant

S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.3

Therefore each forward orbit S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.4 lies on the surface

S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.5

A theorem of Damanik–Fillman–Gorodetski gives the exact spectral criterion

S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.6

On the spectrum, S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.7, and the local fractal properties of S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.8 at energy S(a)=ab,S(b)=a.S(a)=ab,\qquad S(b)=a.9 are governed by aa0 (Fillman et al., 2017).

This characterization is central because it converts a spectral problem for a continuum quasiperiodic operator into a dynamical problem on invariant surfaces. The role of aa1 is not merely diagnostic; it parametrizes the local fractal geometry of the spectrum.

4. Local Hausdorff dimension and asymptotic regimes

For each aa2, the local Hausdorff dimension is defined by

aa3

There exists a continuous, real-analytic mapping

aa4

such that

aa5

Its asymptotics are

aa6

and

aa7

(Fillman et al., 2017).

The continuum theory establishes two asymptotic regimes in which aa8, independently of the precise shapes of aa9 and u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.0. In the high-energy regime,

u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.1

and consequently

u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.2

In the small-coupling regime, for

u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.3

one has uniformly in u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.4,

u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.5

so that

u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.6

(Fillman et al., 2017).

The proofs rely on perturbative transfer-matrix approximations. At large u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.7 or small u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.8, the fundamental solutions u=S(u)=abaababaAN.u=S(u)=abaababa\cdots\in A^{\mathbb N}.9 are close to the free solutions ΩFAZ\Omega_F\subset A^{\mathbb Z}0 and ΩFAZ\Omega_F\subset A^{\mathbb Z}1, leading to the estimates

ΩFAZ\Omega_F\subset A^{\mathbb Z}2

ΩFAZ\Omega_F\subset A^{\mathbb Z}3

and hence

ΩFAZ\Omega_F\subset A^{\mathbb Z}4

Substituting these into

ΩFAZ\Omega_F\subset A^{\mathbb Z}5

shows that in the free case ΩFAZ\Omega_F\subset A^{\mathbb Z}6, while perturbative errors yield ΩFAZ\Omega_F\subset A^{\mathbb Z}7 or ΩFAZ\Omega_F\subset A^{\mathbb Z}8 (Fillman et al., 2017).

The function ΩFAZ\Omega_F\subset A^{\mathbb Z}9 arises from a uniformly hyperbolic repeller on S(b)=aS(b)=a00, whose unstable-manifold contraction rates depend smoothly on S(b)=aS(b)=a01. Analytic perturbation theory gives the expansion S(b)=aS(b)=a02 as S(b)=aS(b)=a03 (Fillman et al., 2017). This makes precise the sense in which the continuum Fibonacci spectrum becomes almost full-dimensional in Hausdorff dimension while retaining its Cantor character.

5. Fractal energy spectrum, self-similarity, and gap labeling

In the polariton realization, the fractal character of a Fibonacci potential is described through the integrated density of states

S(b)=aS(b)=a04

About any reference energy S(b)=aS(b)=a05, one has the discrete scaling relation

S(b)=aS(b)=a06

with scaling factors S(b)=aS(b)=a07 depending on S(b)=aS(b)=a08. Equivalently,

S(b)=aS(b)=a09

where the scaling function S(b)=aS(b)=a10 is periodic of period S(b)=aS(b)=a11. This gives power-law envelopes with exponent S(b)=aS(b)=a12, modulated by log-periodic oscillations (Tanese et al., 2013).

A leading-order equivalent form is

S(b)=aS(b)=a13

which explicitly exhibits log-periodic modulation of period S(b)=aS(b)=a14 in S(b)=aS(b)=a15 (Tanese et al., 2013).

The same work states the gap-labeling theorem for one-dimensional quasi-periodic Schrödinger operators: S(b)=aS(b)=a16 with S(b)=aS(b)=a17. Equivalently, the Fourier expansion of S(b)=aS(b)=a18 has Bragg peaks at

S(b)=aS(b)=a19

and each peak opens a mini-gap whose integrated density of states is exactly S(b)=aS(b)=a20 modulo S(b)=aS(b)=a21 (Tanese et al., 2013).

In weak-potential perturbation theory, each Bragg peak S(b)=aS(b)=a22 couples degenerate plane waves S(b)=aS(b)=a23, opening a gap of width

S(b)=aS(b)=a24

The Fourier or structure-factor form of the potential is

S(b)=aS(b)=a25

(Tanese et al., 2013).

These formulas sharpen the spectral meaning of Fibonacci order. The spectrum is not merely fragmented; its gaps are topologically indexed, and its integrated density of states displays discrete self-similarity with log-periodic corrections.

6. Golden-periodic complex potentials in planar hydrodynamics

A distinct usage of Fibonacci structure occurs in planar hydrodynamics through the PQ-calculus of Fibonacci divisors. For any fixed integer S(b)=aS(b)=a26, the S(b)=aS(b)=a27th Fibonacci divisors are

S(b)=aS(b)=a28

where the Fibonacci numbers satisfy S(b)=aS(b)=a29, S(b)=aS(b)=a30, S(b)=aS(b)=a31, and Binet’s formula

S(b)=aS(b)=a32

The S(b)=aS(b)=a33th Golden derivative is defined by

S(b)=aS(b)=a34

and satisfies

S(b)=aS(b)=a35

In the limit S(b)=aS(b)=a36, this reduces to the ordinary S(b)=aS(b)=a37 (Pashaev, 2021).

The corresponding S(b)=aS(b)=a38th golden exponential is

S(b)=aS(b)=a39

with S(b)=aS(b)=a40, and it obeys

S(b)=aS(b)=a41

The golden translation operator is

S(b)=aS(b)=a42

A function is golden-periodic of order S(b)=aS(b)=a43 exactly when

S(b)=aS(b)=a44

and one has

S(b)=aS(b)=a45

(Pashaev, 2021).

In the annulus

S(b)=aS(b)=a46

the primary vortex potential is

S(b)=aS(b)=a47

which fails the no-normal-flow condition on S(b)=aS(b)=a48 and S(b)=aS(b)=a49. By the classical two-circle theorem, one adds two doubly infinite image sequences: images by inversion in the inner circle at points S(b)=aS(b)=a50, and images by inversion in the outer circle at points S(b)=aS(b)=a51, S(b)=aS(b)=a52. After summation, the closed-form complex potential is

S(b)=aS(b)=a53

and the complex velocity is

S(b)=aS(b)=a54

(Pashaev, 2021).

This hydrodynamic usage differs from the Schrödinger and polariton settings: the word “potential” here denotes a complex potential rather than a scalar potential in an eigenvalue problem. The shared element is golden-periodic recursion. A plausible implication is that Fibonacci order functions as a transferable analytic template across spectral theory and image-based boundary constructions, even when the governing equations are different.

Topic to Video (Beta)

No one has generated a video about this topic yet.

Whiteboard

No one has generated a whiteboard explanation for this topic yet.

Follow Topic

Get notified by email when new papers are published related to Fibonacci Potential.