previous index next Some Useful Integrals of Exponential Functions Michael Fowler We’ve shown that differentiating the exponential function just multiplies it by the constant in the exponent, that is to say, d ax e = ae ax . dx Integrating the exponential function, of course, has the opposite effect: it divides by the constant in the exponent: ∫e ax dx = 1 ax e , a as you can easily check by differentiating both sides of the equation. An important definite integral (one with limits) is ∞ ∫e 0 − ax dx = 1 . a Notice the minus sign in the exponent: we need an integrand that decreases as x goes towards infinity, otherwise the integral will itself be infinite. To visualize this result, we plot below e-x and e-3x. Note that the green line forms the hypotenuse of a right-angled triangle of area 1, and it is very plausible from the graph that the total area under the e-x curve is the same, that is, 1, as it must be. The e-3x curve has area 1/3 under it, (a = 3). 2 Now for something a bit more challenging: how do we evaluate the integral I (a) = ∞ ∫e − ax 2 dx ? −∞ (a has to be positive, of course.) The integral will definitely not be infinite: it falls off equally fast in both positive and negative directions, and in the positive direction for x greater than 1, it’s smaller than e-ax, which we know converges. To see better what this function looks like, we plot it below for a = 1 (red) and a = 4 (blue). 3 Notice first how much faster than the ordinary exponential e-x this function falls away. Then note that the blue curve, a = 4, has about half the total area of the a = 1 curve. In fact, the area goes as 1/ a . The green lines help see that the area under the red curve (positive plus negative) is somewhat less than 2, in fact it’s π = 1.77 approximately. But—it’s not so easy to evaluate! There is a trick: square it. That is to say, write ( I ( a )) 2 ∞ = ∫e − ax 2 −∞ ∞ dx ∫ e − ay dy 2 −∞ Now, this product of two integrals along lines, the x-integral and the y-integral, is exactly the same as an integral over a plane, the (x, y) plane, stretching to infinity in all directions. We can rewrite it ∞ ∫e −∞ − ax 2 ∞ dx ∫ e −∞ − ay 2 ∞ ∞ dy = ∫ ∫e −∞ −∞ − ax 2 e − ay 2 ∞ ∞ dxdy = ∫ ∫e − ar 2 dxdy −∞ −∞ where r 2 = x 2 + y 2 , r is just the distance from the origin (0, 0) in the (x, y) plane. The plane is divided up into tiny squares of area dxdy, and doing the integral amounts to finding the value of 2 e − ar for each tiny square, multiplying by the area of that square, and adding the contributions from every square in the plane. In fact, though, this approach is no easier than the original problem—the trick is to notice that 2 the integrand e − ar has a circular symmetry: for any circle centered at the origin (0, 0), it has the same value anywhere on the circle. To exploit this, we shouldn’t be dividing the (x, y) plane up 4 into little squares at all, we should be dividing it into regions having all points the same distance from the origin. These are called “annular” regions: the area between two circles, both centered at the origin, the inner one of radius r, the slightly bigger outer one having radius r + dr. We take dr to be very small, so this is a thin circular strip, of length 2πr (the circumference of the circle) and breadth dr, and therefore its total area is 2πrdr (neglecting terms like dr2, which become negligible for dr small enough). So, the contribution from one of these annular regions is e − ar 2π rdr , and the complete integral over the whole plane is: 2 ( I ( a )) 2 ∞ = π ∫ e − ar 2rdr. 2 0 This integral is easy to evaluate: make the change of variable to u = r2, du = 2rdr giving ∞ π ( I ( a ) ) = π ∫ e− au du = 2 a 0 so taking square roots I (a) = ∞ ∫e −∞ − ax 2 dx = π a . 5 Some Integrals Useful in the Kinetic Theory of Gases We can easily generate more results by differentiating I(a) above with respect to the constant a! Differentiating once: ∞ ∞ 2 2 d d π 1 π =− e − ax dx = − ∫ x 2 e − ax dx = ∫ da −∞ da a 2a a −∞ so we have ∞ ∫xe 2 − ax 2 dx = −∞ 1 π 2a a and differentiating this result with respect to a gives ∞ ∫xe 4 − ax 2 dx = −∞ 3 4a 2 π a . The ratio of these two integrals comes up in the kinetic theory of gases in finding the average kinetic energy of a molecule with Maxwell’s velocity distribution. ∞ ∫xe dx ∫xe dx 4 − ax 2 −∞ ∞ = 2 − ax 2 −∞ 3 π 4a 2 a 1 π 2a a = 3 . 2a Finding this ratio without doing the integrals: It is interesting to note that this ratio could have been found with much less work, in fact without evaluating the integrals fully, as follows: Make the change of variable y 2 = ax 2 , so dy = adx and ∞ ∫xe −∞ 2 − ax 2 dx = ( a ) −3/ 2 ∞ ∫ye 2 − y2 dy = Ca −3/ 2 −∞ where C is a constant independent of a, because a has completely disappeared in the integral over y. (Of course, we know C = π / 2 , but that took a lot of work.) Now, the integral with x4 in place of x2 is given by differentiating the x2 integral with respect to a, and multiplying by -1, as discussed above, so, differentiating the right hand side of the above equation, the x4 integral is just ( 3 / 2 ) Ca −5/ 2 , and the C cancels out in the ratio of the integrals. 6 However, we do need to do the integrals at one point in the kinetic theory: the overall normalization of the velocity distribution function is given by requiring that ∞ ∞ 1 = ∫ f (v)dv = ∫ 4π A3v 2 e − mv 0 2 / 2 kT 0 and this in fact determines A, using the results we found above, giving ⎛ m ⎞ f (v) = 4π ⎜ ⎟ ⎝ 2π kT ⎠ 3/ 2 v 2 e − mv 2 / 2 kT . One last trick… We didn’t need this in the kinetic theory lecture, but is seems a pity to review exponential integrals without mentioning it. It’s easy to do the integral I ( a, b ) = ∞ − ax bx ∫ e e dx 2 −∞ It can be written I ( a, b ) = ∞ − a ( x −b / 2 a ) b e ∫e 2 ∞ 2 / 4a dx = eb −∞ 2 / 4a − a ( x −b / 2 a ) dx = eb ∫e 2 −∞ where to do the last step just change variables from x to y = x - b/2a. This can even be used to evaluate for example ∞ ∫e − ax 2 −∞ by writing the cosine as a sum of exponentials. previous index next cos bx dx 2 / 4a π /a