Mathieu equation

Mathieu equation
Name	Mathieu equation
Type	Linear second-order ordinary differential equation with periodic coefficients
Introduced by	Émile Léonard Mathieu
Year	1868

Mathieu equation The Mathieu equation is a linear second-order ordinary differential equation with periodic coefficients introduced in the 19th century by Émile Léonard Mathieu for problems in elliptic coordinates, vibrating membranes, and stability of parametric oscillators. It appears in analyses associated with separable problems in Laplace's equation, Helmholtz equation, and in stability studies related to Hill's equation, Floquet theory, and the Mathieu instability phenomena observed in engineering and physics. The equation and its solutions—Mathieu functions—play roles in contexts ranging from quantum mechanics models, electromagnetism in elliptical cavities, to periodic potentials in solid state physics.

Introduction

The Mathieu equation arose historically when Émile Léonard Mathieu studied small vibrations of an elliptic membrane and when later researchers such as Lord Rayleigh, George William Hill, and Augustin-Louis Cauchy examined periodic differential operators. It is a canonical example in the theory developed by Gaston Floquet and later synthesized in texts by E. T. Whittaker, George B. Arfken, and N. N. Lebedev. The problem connects to classical topics handled by Joseph Fourier and later exploited in applied work by figures like Erwin Schrödinger in quantum models and Paul Dirac in wave mechanics.

Definition and canonical forms

The standard, canonical form of the equation is y'' + (a − 2q cos 2x) y = 0, where a and q are parameters. Equivalent parameterizations and scaled variants are used in the literature by authors such as Morse and Feshbach and by contributors to the NIST Digital Library of Mathematical Functions compendia. Alternate forms replace cos 2x with a periodic coefficient of period π or 2π or introduce complex shifts used by researchers in complex analysis and spectral theory contexts. Transformations linking the canonical form to Hill's equation and to the Mathieu–Hill family are standard in monographs by Harold Jeffreys and Austin Hobart.

Mathieu functions and characteristic values

Solutions that are bounded and periodic for given parameters a and q are the even and odd Mathieu functions, often denoted ce_r(x,q) and se_r(x,q), with r indexing characteristic exponents. Their associated eigenvalues, called characteristic values, are typically labeled a_r(q) and b_r(q). These functions appear in treatises by E. T. Whittaker, G. N. Watson, and tables compiled by NIST and by H. Bateman in the Handbook of Mathematical Functions. The properties of orthogonality, completeness, and normalization are established analogously to results for Bessel functions, Legendre polynomials, and Hermite functions, and are exploited in expansions used by Lord Rayleigh and James Clerk Maxwell-related studies of resonances.

Stability, Floquet theory, and Mathieu chart

Floquet theory, developed by Gaston Floquet, provides the framework for analyzing stability of solutions via characteristic exponents (Floquet multipliers). Stability regions in the (a,q) parameter plane—commonly called the Mathieu chart or Ince–Strutt diagram—were delineated by investigators such as Edward L. Ince and William J. Strutt, 3rd Baron Rayleigh. The boundaries between stable and unstable regions connect to concepts examined by George William Hill in celestial mechanics and by A. M. Lyapunov in stability theory. Resonance tongues, parametric resonance, and band-gap structures echo phenomena studied in Anderson localization and in the band theory of solid state physics.

Special cases and limiting behaviors

Several important limits connect Mathieu functions to other special functions. For q → 0 the equation reduces to simple harmonic motion with solutions related to sines and cosines studied by Jean le Rond d'Alembert and Leonhard Euler. For large q and scaled variables, asymptotics approach forms involving Airy functions treated by George B. Airy and WKB approximations associated with Hermann Weyl and Harold Jeffreys. Degenerate or integer-related parameter choices yield periodicities and orthogonality analogous to classical expansions in Fourier series and to eigenproblems in the theory advanced by David Hilbert.

Applications in physics and engineering

Mathieu-type equations model vibrating elliptical membranes analyzed historically by Émile Léonard Mathieu and later used in studies by Lord Rayleigh and George Green. They arise in wave propagation in elliptical waveguides and cavities investigated by James Clerk Maxwell and in stability of parametrically driven oscillators used in engineering contexts traced to Galileo Galilei and developed through André-Marie Ampère-era dynamics. In quantum mechanics Mathieu equations describe particles in cosine potentials akin to the Kronig–Penney model and to approximations for cold-atom optical lattices considered by contemporary researchers including teams at institutions such as CERN and MIT. Applications extend to ion trap stability in Paul trap designs by Wolfgang Paul and to Mathieu-type parametrically excited systems in aerospace engineering and electrical engineering.

Numerical methods and computation of solutions

Computation of Mathieu functions and characteristic values is implemented in numerical libraries and software packages produced by institutions like National Institute of Standards and Technology, and in mathematical software developed by groups at Wolfram Research, MathWorks, and authors of libraries used at Los Alamos National Laboratory. Methods include matrix truncation of Hill determinants, Fourier series expansion, continued fractions techniques popularized by J. Meixner and F. W. J. Olver, and spectral collocation strategies in the spirit of John P. Boyd. Stability charts and band-structure computations employ eigenvalue solvers developed in numerical linear algebra traditions of John von Neumann and Alan Turing.

Category:Ordinary differential equations