[motivation] ## Motivation The method of [matched asymptotic expansions](/pages/1162) handles problems where the solution has different structure in different *spatial regions* — a [boundary](/page/Boundary) layer here, a smooth outer solution there. But there is a second, fundamentally different way that a regular perturbation expansion can fail: the small perturbation may have a **cumulative effect over long times**, even though the solution looks well-behaved at any fixed instant. The prototypical setting is a **weakly perturbed oscillator**: \begin{align*} \ddot{y} + y + \epsilon\,f(y, \dot{y}) = 0, \qquad \epsilon \ll 1. \end{align*} For $\epsilon = 0$, the solution oscillates forever with fixed amplitude and frequency. For small $\epsilon$, the perturbation is weak at any given instant — but over $O(1/\epsilon)$ oscillation periods, its cumulative effect becomes $O(1)$. A regular expansion $y \sim y_0 + \epsilon y_1 + \cdots$ produces terms that grow without bound in time (**secular terms**), causing the expansion to break down at $t = O(1/\epsilon)$. The resolution is that the solution depends on $t$ and $\epsilon$ through the combination $\epsilon t$ (as well as through $t$ itself) — the secular term is a Taylor expansion of a function of $\epsilon t$ that has been artificially broken apart. The **method of multiple scales** addresses this by introducing $\tau = t$ and $T = \epsilon t$ as independent variables, constructing a Poincaré expansion of a two-variable extension $Y(\tau, T;\epsilon)$ in which $\epsilon t$ is never broken apart, and then evaluating the result on the diagonal $(\tau, T) = (t, \epsilon t)$. This page develops the method from first principles: the failure of regular perturbation, the two-time-scale extension, the removal of secular terms, and two fully worked examples. [/motivation]

## The Problem: Weakly Perturbed Oscillators

We consider equations of the form

\begin{align*} \ddot{y} + y + \epsilon\,f(y, \dot{y}) = 0, \qquad \epsilon \searrow 0, \end{align*}

where dots denote [derivatives](/page/Derivative) with respect to $t$. We have normalised so that the unperturbed oscillator ($\epsilon = 0$) has angular frequency $\omega_0 = 1$ (period $2\pi$).

A standard Poincaré expansion $y(t;\epsilon) \sim y_0(t) + \epsilon y_1(t) + \cdots$ (as $\epsilon \to 0$ with $t$ fixed) is valid for any fixed $t$, but produces secular terms — terms that grow without bound in $t$ — making the expansion useless for $t = O(1/\epsilon)$. As we will show, these secular terms arise because the true solution depends on $\epsilon t$ as a unit, and expanding in $\epsilon$ with $t$ fixed breaks this combination apart.

The method of multiple scales resolves this as follows:

1. Define a two-variable extension $Y(\tau, T;\epsilon)$ (where $\tau = t$, $T = \epsilon t$) such that $y(t;\epsilon) = Y(t, \epsilon t;\epsilon)$. 2. Find a Poincaré expansion $Y(\tau, T;\epsilon) \sim Y_0(\tau, T) + \epsilon Y_1(\tau, T) + \cdots$ as $\epsilon \to 0$ with $\tau$ and $T$ **both fixed**, subject to the condition that each $Y_k$ is bounded in $\tau$. 3. Evaluate the result on the **physical diagonal** $(\tau, T) = (t, \epsilon t)$ to obtain an approximation to $y(t;\epsilon)$.

The output of step 3 is *not* a Poincaré expansion of $y(t;\epsilon)$ — it is the restriction to a diagonal of a Poincaré expansion of a different function $Y$. The boundedness condition in step 2 is what prevents secular terms from appearing; it determines the previously free functions in the leading-order solution.

## Failure of Regular Perturbation: Secular Terms

[example: Secular Terms in the Duffing Equation] Consider the Duffing equation \begin{align*} \ddot{y} + y + \epsilon y^3 = 0, \qquad \epsilon \searrow 0. \end{align*} We attempt a regular expansion $y(t;\epsilon) \sim y_0(t) + \epsilon\,y_1(t)$. Substituting into the ODE: \begin{align*} (\ddot{y}_0 + \epsilon\,\ddot{y}_1 + \cdots) + (y_0 + \epsilon\,y_1 + \cdots) + \epsilon(y_0 + \cdots)^3 = 0. \end{align*} Expanding $(y_0 + \epsilon y_1 + \cdots)^3 = y_0^3 + O(\epsilon)$ and grouping by powers of $\epsilon$: \begin{align*} \underbrace{(\ddot{y}_0 + y_0)}_{O(1)} + \epsilon\underbrace{(\ddot{y}_1 + y_1 + y_0^3)}_{O(\epsilon)} + O(\epsilon^2) = 0. \end{align*} *At $O(1)$:* $\ddot{y}_0 + y_0 = 0$, with general solution $y_0 = a\cos(t + \phi_0)$ for constants $a > 0$ and $\phi_0$. *At $O(\epsilon)$:* \begin{align*} \ddot{y}_1 + y_1 = -y_0^3. \end{align*} Computing the right-hand side using the identity $\cos^3\theta = \frac{3}{4}\cos\theta + \frac{1}{4}\cos 3\theta$: \begin{align*} -y_0^3 = -a^3\cos^3(t + \phi_0) = -\frac{3a^3}{4}\cos(t + \phi_0) - \frac{a^3}{4}\cos 3(t + \phi_0). \end{align*} The equation for $y_1$ is now: \begin{align*} \ddot{y}_1 + y_1 = -\frac{3a^3}{4}\cos(t + \phi_0) - \frac{a^3}{4}\cos 3(t + \phi_0). \end{align*} The $\cos 3(t+\phi_0)$ term on the right is unproblematic — it drives a bounded particular solution proportional to $\cos 3(t+\phi_0)$. But the $\cos(t+\phi_0)$ term is a **resonant forcing**: it solves the homogeneous equation $\ddot{y}_1 + y_1 = 0$. By the theory of forced oscillators, the particular solution for this term is: \begin{align*} y_1^{(\text{resonant})} = -\frac{3a^3}{8}\,t\sin(t + \phi_0). \end{align*} This grows linearly in $t$. When $t = O(1/\epsilon)$, the "correction" $\epsilon y_1 = O(1)$ — it is no longer small compared to $y_0$. The regular expansion has broken down. These unbounded terms are called **secular terms** (from the Latin *saeculum*, meaning "age" — they become important over long time scales). [/example]

[remark: When Do Secular Terms Arise?] Secular terms appear whenever the right-hand side of the $O(\epsilon)$ equation contains terms that solve the corresponding homogeneous equation. For oscillator problems with natural frequency $\omega_0 = 1$, this means terms proportional to $e^{\pm it}$ (in complex notation) or $\cos t$ and $\sin t$ (in real notation). In Fourier-[series](/page/Series) language: secular terms arise from the **first-harmonic** (fundamental frequency) components of the forcing. Higher harmonics produce bounded responses. [/remark]

## The Two-Time-Scale Extension

### Motivation

The secular term tells us something physical: the weak perturbation slowly modulates the oscillation over $O(1/\epsilon)$ time periods. The solution involves two processes operating on different scales:

- Fast oscillation on the scale $\tau = t$ (period $2\pi$) - Slow evolution of amplitude and phase on the scale $T = \epsilon t$

The idea of the method of multiple scales is to treat these as **independent variables** and construct a two-variable function $Y(\tau, T;\epsilon)$ such that

\begin{align*} y(t;\epsilon) = Y(t,\, \epsilon t;\, \epsilon) \end{align*}

on the **physical diagonal** $(\tau, T) = (t, \epsilon t)$. This is called a **two-variable extension** of $y$. Crucially, this extension is **not unique**: given $y(t;\epsilon)$, there are many two-variable [functions](/page/Function) $Y(\tau, T;\epsilon)$ that agree with $y$ when restricted to the diagonal $T = \epsilon\tau$. For instance, adding any function that vanishes on the diagonal (such as $(\tau - T/\epsilon) \cdot g(\tau, T)$) produces a different $Y$ with the same $y$. This non-uniqueness is not a defect — it is the **freedom** that the method exploits. The no-secular-terms condition (Step 5 below) will select a specific $Y$ from all the possibilities by requiring that the expansion coefficients $Y_0, Y_1, \ldots$ remain bounded in $\tau$.

### The chain rule