pith. sign in

arxiv: 2603.07855 · v2 · submitted 2026-03-09 · ⚛️ physics.optics · cond-mat.mes-hall· math-ph· math.MP

Explicit Construction of Floquet-Bloch States from Arbitrary Solution Bases of the Hill Equation

Gregory V Morozov This is my paper

Pith reviewed 2026-05-15 14:30 UTC · model grok-4.3

classification ⚛️ physics.optics cond-mat.mes-hallmath-phmath.MP
keywords Hill equationFloquet-Bloch statesmonodromy matrixphotonic crystalsperiodic systemstransfer matrixband edges
0
0 comments X

The pith

Explicit closed-form formulas map any pair of independent Hill-equation solutions to the Floquet-Bloch basis via the monodromy matrix.

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

The paper develops a direct method to obtain Floquet-Bloch states, which describe wave behavior in one-dimensional periodic systems, from any two linearly independent solutions of the Hill equation rather than from specially normalized ones. It supplies explicit expressions that use the monodromy matrix to perform the mapping and handles the generic case at band edges where the matrix takes a Jordan form. The construction is illustrated with one-dimensional photonic crystals and can be rewritten in terms of the transfer matrix to expose remaining freedom in the basis choice. This supplies an implementation-ready route for both analytical and numerical work on periodic media.

Core claim

For the Hill equation with periodic coefficients, an arbitrary fundamental system of two linearly independent solutions can be converted to the corresponding Floquet-Bloch basis by closed-form formulas that involve only the monodromy matrix; the formulas remain valid in the generic Jordan-block case at band edges and do not require canonical normalization of the input solutions.

What carries the argument

The monodromy matrix, which records the linear transformation experienced by any solution pair after one spatial period, acts as the explicit bridge that converts an arbitrary fundamental system into the Floquet-Bloch states.

If this is right

  • The same formulas apply directly to one-dimensional photonic crystals for band-structure calculations.
  • The transfer-matrix form of the construction makes residual basis freedom explicit for numerical codes.
  • The method covers both ordinary and degenerate band-edge cases without separate handling.
  • The resulting framework supports direct implementation in symbolic or numerical software for periodic systems.

Where Pith is reading between the lines

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

  • Numerical solvers for periodic problems could adopt the construction to bypass the need to locate specially normalized solutions at each parameter value.
  • The approach may generalize to other linear periodic differential equations whose fundamental matrices can be tracked over one period.
  • Dispersion relations could be extracted from the monodromy eigenvalues without first constructing canonical Bloch waves.

Load-bearing premise

The input pair must be linearly independent solutions of the Hill equation and the monodromy matrix must be well-defined from the periodic coefficient.

What would settle it

Compute the constructed states from two arbitrary independent solutions, propagate them over one period, and check whether they fail to reproduce the original differential equation or the expected Floquet multiplier.

Figures

Figures reproduced from arXiv: 2603.07855 by Gregory V Morozov.

Figure 1
Figure 1. Figure 1: FIG. 1. Schematic representation of a binary photonic cryst [PITH_FULL_IMAGE:figures/full_fig_p009_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. Transmittance and cos( [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Absolute values of the Floquet–Bloch waves over [PITH_FULL_IMAGE:figures/full_fig_p010_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Absolute values of the Floquet–Bloch waves over [PITH_FULL_IMAGE:figures/full_fig_p011_4.png] view at source ↗
read the original abstract

For the Hill equation describing one-dimensional periodic systems, a constructive formulation is developed for generating Floquet-Bloch states directly from arbitrary pairs of linearly independent solutions. One-dimensional photonic crystals are used as a concrete illustration. Explicit closed-form formulas map an arbitrary fundamental system to the corresponding Floquet-Bloch basis via the monodromy matrix, including the generic Jordan band-edge case, without reliance on canonically normalized solutions. The construction can be expressed directly in terms of the transfer matrix, making the residual representation freedom transparent and providing an implementation-ready framework for analytical and numerical studies of periodic systems.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 3 minor

Summary. The paper develops an explicit constructive mapping from an arbitrary pair of linearly independent solutions of the Hill equation to the corresponding Floquet-Bloch basis. It forms the monodromy matrix M = Y(0)^{-1} Y(T) from the fundamental matrix Y(x), extracts its eigenvectors (or generalized eigenvector in the Jordan case when trace(M) = ±2), and supplies closed-form algebraic expressions for the change-of-basis coefficients. The construction is recast in terms of the transfer matrix, handles the generic band-edge Jordan block without additional normalization, and is illustrated for one-dimensional photonic crystals.

Significance. If the derivations hold, the result supplies a practical, implementation-ready procedure that removes the need for canonically normalized input solutions and makes representation freedom transparent. This is valuable for analytical work and numerical studies of periodic systems governed by Hill's equation, including band-structure calculations and wave propagation in photonic crystals. The explicit closed-form treatment of the Jordan case at band edges is a concrete strength.

minor comments (3)
  1. [§2.2] §2.2, after Eq. (7): the definition of the fundamental matrix Y(x) would benefit from an explicit statement that its columns are the two input solutions y1(x), y2(x) together with their Wronskian normalization.
  2. [§4.1] §4.1, paragraph following Eq. (19): the Jordan-block construction for the generalized eigenvector should include a brief reminder that the linear-in-x term is the standard form required by the nilpotent part of the monodromy matrix.
  3. [Figure 2] Figure 2 caption: the photonic-crystal example would be clearer if the period T and the specific coefficient functions q(x) were stated numerically alongside the plotted dispersion curves.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for the positive evaluation of our manuscript, the recognition of its practical value for periodic systems, and the recommendation of minor revision. We appreciate the emphasis on the closed-form treatment of the Jordan band-edge case and the transparency of representation freedom.

Circularity Check

0 steps flagged

Direct algebraic construction from monodromy matrix; no circularity

full rationale

The paper derives Floquet-Bloch states by forming the monodromy matrix M = Y(0)^{-1} Y(T) from an arbitrary fundamental matrix Y(x) of linearly independent solutions to the Hill equation, then extracting eigenvectors (or generalized eigenvectors for the Jordan block when trace(M)=±2). All steps are explicit algebraic combinations of endpoint values and derivatives; the construction is self-contained within standard Floquet theory, requires no fitted parameters, imposes no self-citation load-bearing premises, and does not rename or smuggle in prior results. The result is a direct change-of-basis that remains independent of the specific input solutions beyond the defining linear independence.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The central claim rests on standard properties of second-order linear ODEs and the definition of the monodromy matrix; no free parameters or invented entities are introduced.

axioms (2)
  • standard math Any second-order linear homogeneous ODE possesses a fundamental system of two linearly independent solutions.
    Invoked when the paper refers to an arbitrary fundamental system of solutions to the Hill equation.
  • standard math The monodromy matrix is obtained by propagating the fundamental solutions over one period of the periodic coefficient.
    Central to the mapping; follows from the existence and uniqueness theorem for linear ODEs.

pith-pipeline@v0.9.0 · 5395 in / 1352 out tokens · 49292 ms · 2026-05-15T14:30:29.891116+00:00 · methodology

discussion (0)

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

Reference graph

Works this paper leans on

41 extracted references · 41 canonical work pages

  1. [1]

    We derive explicit closed-form transformation for- mulas that map an arbitrary fundamental matrix to the corresponding Floquet–Bloch fundamental system via reduction of the monodromy matrix to its canonical (diagonal or Jordan) form

    work page

  2. [2]

    We provide a complete classification of the residual representation freedom of Floquet–Bloch states in- duced by the choice of a fundamental system: inde- pendent rescalings in the nondegenerate case, and upper-triangular mixing of the hybrid mode with the periodic (or antiperiodic) Floquet–Bloch wave in the degenerate (Jordan) case

    work page

  3. [3]

    We reformulate the construction directly in terms of the one-period transfer matrix, an intrinsic prop- agation operator of the differential equation whose eigenvectors correspond to Floquet–Bloch states. This formulation allows Floquet–Bloch states to be obtained directly from the transfer matrix, with- out introducing canonically normalized solutions or ...

    work page

  4. [4]

    Taken together, these results complement the classi- cal spectral theory by rendering the basis transformation underlying Floquet–Bloch states fully explicit and oper- ational

    We organize the results into a systematic step-by- step procedure that translates the abstract Floquet structure into a practical workflow for analytical and numerical implementations. Taken together, these results complement the classi- cal spectral theory by rendering the basis transformation underlying Floquet–Bloch states fully explicit and oper- ation...

    work page

  5. [5]

    We first treat the generic diagonalizable case, characterized by tr[ ˜A] ≡ a11 + a22 ̸= ± 2, for which the multipliers satisfy ρ1 ̸= ρ2

    Bandgaps and allowed bands. We first treat the generic diagonalizable case, characterized by tr[ ˜A] ≡ a11 + a22 ̸= ± 2, for which the multipliers satisfy ρ1 ̸= ρ2. An explicit diagonalizing matrix may be chosen as ˜B =    1 a12 ρ2 − a11 a21 ρ1 − a22 1    . (21) In this case, both Floquet–Bloch states are Floquet– Bloch waves, and in accordance with ...

    work page

  6. [6]

    When tr[ ˜A] ≡ a11 + a22 = ± 2, Eqs

    Band-edge case. When tr[ ˜A] ≡ a11 + a22 = ± 2, Eqs. (18) - (19) immediately yield a double Floquet mul- tiplier ρ1 = ρ2 ≡ ρ = 1 2 (a11 + a22) = ± 1, and consequently cos( µd ) = ± 1, see Eq. (20). The sys- tem is therefore at a band edge of the spectrum. The construction of the Floquet–Bloch basis now depends on whether the off-diagonal entries of ˜A vani...

    work page

  7. [7]

    Compute any pair of linearly independent solutions E1(z) and E2(z) of Eq

    Compute a fundamental system on the first period. Compute any pair of linearly independent solutions E1(z) and E2(z) of Eq. (2) on the interval 0 ≤ z ≤ d

    work page

  8. [8]

    Construct the fundamental matrix ˜E(z) from E1(z) and E2(z) and evaluate ˜A = ˜E− 1(0) ˜E(d), see Eq

    F orm the monodromy matrix. Construct the fundamental matrix ˜E(z) from E1(z) and E2(z) and evaluate ˜A = ˜E− 1(0) ˜E(d), see Eq. (11)

    work page

  9. [9]

    • If tr[ ˜A] ̸= ± 2 (allowed bands or bandgaps), use Eq

    Construct Floquet–Bloch states on the first period. • If tr[ ˜A] ̸= ± 2 (allowed bands or bandgaps), use Eq. (22) to construct the two Floquet– Bloch waves. • If tr[ ˜A] = ± 2 (band edges), distinguish: (a) If at least one off-diagonal element of ˜A is nonzero (generic band edge), use Eqs. (24) and (25) to construct the Floquet–Bloch wave and the hybrid Flo...

    work page

  10. [10]

    For Floquet– Bloch waves, use Eq

    Extend the states to all periods. For Floquet– Bloch waves, use Eq. (29). For the hybrid Flo- quet mode, use Eq. (30). In incipient bands, the Floquet–Bloch states are already periodic or an- tiperiodic, so extension is immediate. Optimized construction (with transfer matrix). Assume that the transfer matrix ˜W (z, 0) is known on the first period 0 ≤ z ≤ d

    work page

  11. [11]

    Choose any invertible 2 × 2 matrix as ˜E(0)

    Specify a fundamental matrix at the origin. Choose any invertible 2 × 2 matrix as ˜E(0). As dis- cussed in Section 5, this uniquely specifies a funda- mental system of solutions of Eq. (2) via its initial data

    work page

  12. [12]

    Compute ˜A = ˜E− 1(0) ˜Wd ˜E(0), see Eq

    F orm the monodromy matrix. Compute ˜A = ˜E− 1(0) ˜Wd ˜E(0), see Eq. (46)

    work page

  13. [13]

    Construct Floquet–Bloch states at the ori- gin. Proceed exactly as in Step 3 of the direct construction above, using the monodromy matrix ˜A just computed, distinguishing between nondegener- ate bands, generic band edges, and incipient bands (case 2b of Section 3), but apply the construction only to determine the Floquet–Bloch initial vectors F1(0) and F2...

    work page

  14. [14]

    For 0 ≤ z ≤ d, evaluate ˜F (z) = ˜W (z, 0) ˜F (0)

    Propagate across the first period. For 0 ≤ z ≤ d, evaluate ˜F (z) = ˜W (z, 0) ˜F (0)

    work page

  15. [15]

    Extend to all subsequent periods. Use Eq. (29) for Floquet–Bloch waves and Eq. (30) for the hybrid Floquet mode. In incipient bands, pe- riodicity or antiperiodicity makes further extension trivial. Comment. The transfer-matrix implementation avoids computing individual solutions E1(z) and E2(z) through- out the first period. Floquet–Bloch states are con- ...

    work page

  16. [16]

    83 µm − 1 in a bandgap (Fig. 4). In full agreement with the analysis of Section IV, the two constructions differ only by constant ( z-independent) complex factors, α 1 = F (pw) 1 (z) F (n) 1 (z) , α 2 = F (pw) 2 (z) F (n) 2 (z) . (61) FIG. 3. Absolute values of the Floquet–Bloch waves over the first three periods of the binary photonic crystal, with paramet...

    work page

  17. [17]

    McLachlan N W 1947 Theory and applications of Math- ieu functions (New York: Oxford Univ. Press)

    work page 1947

  18. [18]

    Stoker J J 1950 Nonlinear Vibrations (New York: Inter- science)

    work page 1950

  19. [19]

    Magnus W and Winkler S 1966 Hill’s Equation (New York: Wiley)

    work page 1966

  20. [20]

    Eastham M S P 1975 The Spectral Theory of Periodic Differential Equations (Edinburgh: Scottish Academic Press)

    work page 1975

  21. [21]

    Yakubovich V A and Starzhinskii V M 1975 Linear Dif- ferential Equations with Periodic Coefficients (New York: Wiley)

    work page 1975

  22. [22]

    Fedoryuk M V 1985 Ordinary Differential Equations [in Russian] (Moscow: Nauka)

    work page 1985

  23. [23]

    Cottey A A 1972 J. Phys. C: Solid State Physics 5 2583– 2590

    work page 1972

  24. [24]

    Nusinsky I and Hardy A A 2006 Phys. Rev. B 73 125104

    work page 2006

  25. [25]

    Lett.) 96 54005

    Morozov G V and Sprung D W L 2011 EPL (Europhys. Lett.) 96 54005

    work page 2011

  26. [26]

    Caffrey S, Morozov G V, Sprung D W L and Martorell J 2017 Opt. Quant. Electron. 49 112

    work page 2017

  27. [27]

    Bukov M, D’Alessio L and Polkovnikov A 2015 Adv. Phys. 64 139–226

    work page 2015

  28. [28]

    Eckardt A 2017 Rev. Mod. Phys. 89 011004

    work page 2017

  29. [29]

    Ikeda T N 2018 Phys. Rev. A 97 063413

    work page 2018

  30. [30]

    Asboth J K, Tarasinski B and Delplace P 2014 Phys. Rev. B 90 125143

    work page 2014

  31. [31]

    Minguzzi J, Zhu Z, Sandholzer K, Walter A S, Viebahn K and Esslinger T 2022 Phys. Rev. Lett. 129 053201

    work page 2022

  32. [32]

    Wang P 2023 Phys. Rev. A 108 063304

    work page 2023

  33. [33]

    Coppini F and Santini P M 2024 J. Phys. A: Math. Theor. 57 075701

    work page 2024

  34. [34]

    Mogilner A I and Loly P D 1992 J. Phys. A: Math. Gen. 25 L855–L860

    work page 1992

  35. [35]

    Ibrahim A, Sprung D W L and Morozov G V 2018 J. Opt. Soc. Am. B 35 1223–1232

    work page 2018

  36. [36]

    Sprung D W L, Wu H and Martorell J 1993 Am. J. Phys. 61 1118–1124

    work page 1993

  37. [37]

    Lekner J 1994 J. Opt. Soc. Am. A 11 2892–2899

    work page 1994

  38. [38]

    Morozov G V, Placido F and Sprung D W L 2011 J. Opt. 13 035102

    work page 2011

  39. [39]

    Morozov G V 2023 Opt. Quant. Electron. 55 1120

    work page 2023

  40. [40]

    Meissner E 1918 Schweizerische Bauzeitung 72 95–98

    work page 1918

  41. [41]

    Wang Z, Cheng Y, Nie Y, Wang X and Gong R 2014 J. Appl. Phys. 116 054905

    work page 2014