Physical applications of higher terms of Taylor series

Question

Depressingly many of the physical "applications" of Taylor series that I can find in textbooks and online are actually just applications of linear approximation, since they only take the constant and linear term of a Taylor series. For instance, this how you get the Newtonian limit of the kinetic energy of special relativity, the $1/r^3$ behavior of an electric dipole far away, and the standard approximation to the period of a pendulum. Since the whole point of a Taylor series is that we can go on from the linear approximation to get quadratic and higher approximations, I find this unsatisfying.

So what are some good physical applications of the higher order terms in a Taylor expansion? (I'm looking for something more than just "we can get a better approximation than the linear one by including more terms" — ideally, a situation where the higher terms give some new qualitative insight.)

Answer 1

Here's one that I just thought of, by modifying a problem in the ODEs section of my textbook.

Question: In the presence of air resistance, does a thrown ball take longer to go up or to come down?

We assume that the force of air resistance is proportional to velocity, with constant of proportionality $p$. Thus, we have

$$ F = m a = - p v - m g. $$

Since $a=v'$, this is a differential equation for $v$. For simplicity, let's divide out by $m$ and write $q=p/m$, so it becomes

$$ v' = - q v - g. $$

The differential equation is separable, so we can solve it to get

$$ v = \left(v_0 + \frac{g}{q}\right) e^{-q t} - \frac{g}{q}. $$

Integrating again, we get the height as a function of time:

$$ y = \left(\frac{v_0}{q} + \frac{g}{q^2}\right)(1-e^{-q t}) - \frac{g t}{q}. $$

We can solve $v = 0$ to find the time $t_{\mathrm{max}}$ at which the ball reaches its maximum height, but it's not so easy to solve $y=0$ to find the time at which it reaches the ground. The book answers the question by using a clever argument to conclude that $y(2t_{\mathrm{max}})>0$. But we can also approximate $y$ by a power series in the small parameter $q$.

Let's substitute the second-degree Taylor polynomial for $e^{-q t}$, namely

$$1 - q t + \frac{1}{2} q^2 t^2,$$

into $v$. We get

$$ v \approx \left(v_0 + \frac{g}{q}\right) \left(1 - q t + \frac{1}{2} q^2 t^2\right) - \frac{g}{q} $$

Multiplying this out, we get

$$ v\approx v_0 + \frac{g}{q} - (v_0 q + g) t + \frac{1}{2} (v_0 q^2 + g q) t^2 - \frac{g}{q}. $$

Canceling the $\frac{g}{q}$s and discarding the term involving $q^2$ (since $q$ is small), we obtain

$$ v\approx v_0 + (v_0 q + g) t + \frac{1}{2} g q t^2 . $$

to first order in $q$. Note that all the $q$s in the denominators canceled out, and this is a linear approximation to $v$ as a function of $q$. In principle, we could have derived it by simply taking the derivative of $v$ with respect to $q$ at $q=0$. However, it's not at all obvious from the formula for $v$ that it's even differentiable at $q=0$! (It's not even necessarily obvious that it has a limit as $q\to 0$, but of course the physical considerations imply that it must.) The power series expansion is much nicer.

We can also integrate $v$ to get the height

$$y \approx v_0 t + \frac{1}{2} (v_0 q + g) t^2 + \frac{1}{6} g q t^3. $$

Now we can answer the original question. Putting $v=0$ and solving for $t$ with the quadratic formula, we get

$$t_{\mathrm{max}} \approx \frac{1}{q} + \frac{v_0}{g} \pm \frac{1}{q}\left(1+\frac{v_0^2 q^2}{g^2}\right)^{1/2}.$$

Using the linear Taylor polynomial for the binomial series, we get

$$t_{\mathrm{max}} \approx \frac{1}{q} + \frac{v_0}{g} \pm \frac{1}{q}\left(1+\frac{v_0^2 q^2}{2g^2}\right).$$

We must use the minus sign to cancel out the impossible $1/q$s, yielding

$$t_{\mathrm{max}} \approx \frac{v_0}{g} - \frac{v_0^2 q}{2 g}.$$

Of course, the first term, $\frac{v_0}{g}$, is the obvious value in the no-air-resistance case $v = v_0 - g t$.

We can also set $y=0$ and solve for $t$ to find the time when the ball lands. The solution $t=0$ (when it was thrown) factors out and we can use the quadratic formula again:

$$ t_{\mathrm{lands}} \approx \frac{3}{2q} + \frac{3v_0}{2g} \pm \frac{3}{2q}\left(1 - \frac{2qv_0}{3g} + \frac{v_0^2 q^2}{g^2}\right)^{1/2}.$$

Using the second-degree Taylor polynomial for the binomial series this time (since there is a power of $q$ rather than just $q^2$ inside the square root) but discarding the resulting $q^4$ term, we get

$$ t_{\mathrm{lands}} \approx \frac{3}{2q} + \frac{3v_0}{2g} \pm \frac{3}{2q}\left(1 - \frac{qv_0}{3g} + \frac{v_0^2 q^2}{2g^2} - \frac{q^2 v_0^2}{18 g^2}\right).$$

Again we must take the minus sign to cancel the impossible $\frac{3}{2q}$, and we get

$$ t_{\mathrm{lands}} \approx \frac{2 v_0}{g} - \frac{2 v_0^2 q}{3g^2}. $$

Observe that $t_{\mathrm{lands}} > 2 t_{\mathrm{max}}$.

Answer 2

First f all, check out this great paper. It has some interesting examples, especially in Appendix A.

A typical example is relativistic mechanics: the usual first order approximation yields the unsatisfactory Newtonian mechanics (see for example here), which is in many applications not enough to explain what is happening.

ADDED: Since OP asked for an example, calculations in special relativity usually boil down to the use of the binomial series, hence they are mathematically suitable for a calculus class. Whether people understand the underlying physics is a different question.

For example, Section 15.5 in these great notes explains how in the relativistic harmonic oscillator the period depends on the amplitude (Section 15.5.1, page 311), Formula (15.132).

Answer 3

Much of quantum mechanics is only known in a perturbative framework. Essentially, the basic objects of interest are series. Terms are calculated by Feynman diagrams of increasingly complex diagrams for higher order corrections. So, the series description is pragmatically fundamental. I suppose in principle, non-pertubative solutions exist, however, if no one can hack the math then for a physical purposes (literally in my answer) the series generated by Feynman diagrams represents the best view we have of quantum field theory.

Also, low order applications are interesting. We use $F=mg$ to approximate $F = \frac{GmM}{(R_E+h)^2}$ at $h=0$. In fact, $F=mg$ is just the constant term in the Taylor series for gravity. I submit to you this is still interesting, and, in view of the final I just graded, still sufficiently challenging for many undergraduates.

Applications where you use all the terms in a series are perhaps nicely found in math itself. For example, the alternating harmonic series convergence to $\ln(2)$ as seen from the expansion of $\ln(1+x)$. This is a nice example which uses all the terms in the series. Many examples of this nature are found in every textbook, they make formidable homework problems when taken out of context.

Answer 4

I don't know if this is the sort of thing you're thinking about -- it seems to me like a too-obvious example, but it seems to fit your question pretty well. Suppose you have a particle moving along some path $r(t)$, and suppose we expand $r(t)$ in a Taylor series as $r(t) = a_0 + a_1 t + (1/2) a_2 t^2 + (1/6) a_3 t^3 + ....$ Then $a_0$ tells us the initial location of the particle, $a_1$ its initial velocity, $a_2$ its initial acceleration, $a_3$ its initial jerk, $a_4$ the initial "jounce" or "snap" (beyond that the derivatives don't have standard names).

A second example, one that maybe isn't "physical" enough but also might give insight into the higher-order terms: Compound interest. If a bank account balance is growing according to $B(t)=A e^{rt}$, and we expand the RHS as a Taylor series, the constant term tells you the initial balance, and the linear term tells you how the balance would grow if we were only considering simple interest; the higher-order terms take into account the effects of compounding. (You can adapt this example to deal with any kind of exponential growth or decay.)

Answer 5

While this isn't a 'mainstream' topic, I feel like homotopy analysis methods are worth a glance here. The idea is to approximate some non-linear equation by a linear one, take a continuous connection between them, solve the linear problem, and 'slide' the solution over to the non-linear case. If $L\{u(t)\}=0$ is the linear operator and $N\{u(t)\}$ is the non-linear operator, we introduce a new parameter $q$ and consider the equation

$$(1-q)L\{u(t;q)\}+qN\{u(t;q)\}=0$$

where the solution now depends on $q$. $q=0$ is the linear case and $q=1$ is the full non-linear equation you're trying to solve. What we do is expand the dependence on $q$ as a Taylor series $u(x;q)=\sum u_n(x)q^n$. This leads to a recursion relation for $u_n$ in terms of $u_{n-1}$ that can be solved iteratively. Fiddling around with parameters to ensure that the radius of convergence includes $q=1$ then gives you the non-linear solution as an infinite series. While not physics per se, it has major applications when dealing with nonlinear waves and similar things. The wikipedia page says it better than I can. It might not be appropriate to go through it in detail (it could be done with only an ODE course, but it would probably take as much time to understand as a major topic of the curriculum), but it's a good example of where arbitrarily high terms in the Taylor series become important. It's not just a formality; the first or second order approximations are rarely of any meaningful accuracy.

Another major example (though I again apologize for it not being purely physics) would be in numerical approximations, where Taylor series are bread and butter. Being able to expand Taylor series is fundamental in deriving many methods for numerically solving ODEs, and is crucial for being able to estimate bounds and convergence orders. If you only take linear approximations you'll have a horrible time of everything.

Answer 6

I'm not sure if this is a Taylor series, but I know in physical chemistry, the higher virial coefficients of pressure have an intuitive explanation, for the deviations from ideal gas law. Like the first term is just the law itself (PV=nRT) and then the second and third are deviations based on (a) volume of the molecules themselves (not point masses bouncing around) and (b) attraction between the moleculues when close.

I forget which is which, but see:

https://chem.libretexts.org/Bookshelves/Physical_and_Theoretical_Chemistry_Textbook_Maps/Physical_Chemistry_(LibreTexts)/16%3A_The_Properties_of_Gases/16.05%3A_The_Second_Virial_Coefficient

Or any physical chemistry, or engineering or physics thermodynamics text.

Anyhow, hope it helps since I always felt it was cool that these extra terms had an intuitive physical meaning!