Exponential-Type Integrals

  • POTW
  • Thread starter Euge
  • Start date
  • Tags
    Integrals
In summary, exponential-type integrals are mathematical integrals that involve exponential functions and are used in various fields to model phenomena. They can be solved using techniques like substitution or integration by parts, and have applications in physics, engineering, economics, and other areas. These integrals have special properties and can be evaluated numerically for more accurate results.
  • #1
Euge
Gold Member
MHB
POTW Director
2,054
211
Evaluate, with proof, the definite integral $$\int_{-\infty}^\infty \frac{e^{ax}}{1 + e^x}\, dx$$ where ##0 < a < 1##.
 
  • Like
Likes Greg Bernhardt and topsquark
Physics news on Phys.org
  • #2
Firstly, for [itex]\operatorname{Re}(p) > 0[/itex] and [itex]\operatorname{Re}(q) > 0[/itex] we have the following integral representations of the Beta function:
[tex]
B(p,q) = \int_0^1 t^{p-1}(1 - t)^{q-1}\,dt = \int_0^\infty \frac{u^{p-1}}{(1 + u)^{p+q}}\,du.
[/tex] Secondly, from the reflection formula of the Gamma function we have [tex]
\Gamma(z)\Gamma(1-z) = B(z, 1-z) = \pi \csc(\pi z).[/tex]

I imagine, however, that is not sufficient to just quote these, and that it is intended that one should prove the second result.
 
  • #3
The integral

$$
\int_{-\infty}^\infty \frac{e^{a (x+iy)}}{1 + e^{x+iy}} dx
$$

for ##y = 2 \pi## is

$$
\int_{-\infty}^\infty \frac{e^{a x} e^{i 2 \pi a}}{1 + e^x e^{i 2 \pi}} dx = e^{i 2 \pi a} \int_{-\infty}^\infty \frac{e^{a x}}{1 + e^x} dx
$$

The integral

\begin{align*}
\lim_{R \rightarrow \infty} \left| \int_0^Q \frac{e^{a R} e^{i a y}}{1 + e^R e^{iy}} dy \right| & = \lim_{R \rightarrow \infty} e^{- (1 - a) R} \left| \int_0^Q \frac{e^{i a y}}{e^{-R} + e^{iy}} dy \right|
\nonumber \\
& \leq \lim_{R \rightarrow \infty} e^{- (1 - a) R} \int_0^Q \frac{1}{\sqrt{e^{-2R} + 2 e^{-R} \cos y + 1}} dy
\nonumber \\
& \leq \lim_{R \rightarrow \infty} e^{- (1 - a) R} \int_0^Q \frac{1}{\sqrt{e^{-2R} - 2 e^{-R} + 1}} dy
\nonumber \\
& \leq \lim_{R \rightarrow \infty} e^{- (1 - a) R} \frac{1}{1-e^{-R}} Q
= 0 \quad (*)
\end{align*}

because ##a < 1##. The integral

\begin{align*}
\lim_{R \rightarrow \infty} \left| \int_0^Q \frac{e^{- a R} e^{i a y}}{1 + e^{-R} e^{iy}} dy \right| & = \lim_{R \rightarrow \infty} e^{- a R} \left| \int_0^Q \frac{e^{i a y}}{1 + e^{-R} e^{iy}} dy \right|
\nonumber \\
& \leq \lim_{R \rightarrow \infty} e^{- a R} \int_0^Q \frac{1}{\sqrt{e^{-2R} + 2 e^{-R} \cos y + 1}} dy
\nonumber \\
& \leq \lim_{R \rightarrow \infty} e^{- a R} \int_0^Q \frac{1}{\sqrt{e^{-2R} - 2 e^{-R} + 1}} dy
\nonumber \\
& \leq \lim_{R \rightarrow \infty} e^{- a R} \frac{1}{1-e^{-R}} Q
= 0 \quad (**)
\end{align*}

because ##a > 0##.We define a rectangular contour integral:

\begin{align*}
& x \text{ integrated from } -R \text{ to } R, \; y =0
\nonumber \\
& x = R , \; y \text{ integrated from } 0 \text{ to } 2 \pi
\nonumber \\
& x \text{ integrated from } R \text{ to } -R, \; y = 2 \pi
\nonumber \\
& x = -R , \; y \text{ integrated from } 2 \pi \text{ to } 0
\end{align*}

This encloses a pole at ##z_0 = i \pi##. Near ##z_0 = i \pi## we have ##1 + e^z = 1 + e^{z_0} e^{z - z_0} = 1 - e^{z - z_0} = - (z - z_0)##, so that

$$
\oint \frac{e^{a z}}{1 + e^z} dz = - 2 \pi i e^{i a \pi}
$$

Then

\begin{align*}
- 2 \pi i e^{i a \pi} & = \oint \frac{e^{a z}}{1 + e^z} dz
\nonumber \\
& = \int_{-R}^R \frac{e^{a x}}{1 + e^x} dx + \int_0^{2 \pi} \frac{e^{a R} e^{i a y}}{1 + e^R e^{iy}} idy - e^{i 2 \pi a} \int_{-R}^R \frac{e^{a x}}{1 + e^x} dx - \int_0^{2 \pi} \frac{e^{- a R} e^{i a y}}{1 + e^{-R} e^{iy}} idy
\end{align*}

In the limit ##R \rightarrow \infty## we have:

$$
- 2 \pi i e^{i a \pi} = (1 - e^{i 2 \pi a}) \int_{-\infty}^\infty \frac{e^{a x}}{1 + e^x} dx
$$

where we have used ##(*)## and ##(**)##. So that

$$
\int_{-\infty}^\infty \frac{e^{a x}}{1 + e^x} dx = \frac{\pi}{\sin \pi a} .
$$
 
Last edited:
  • #4
The correct answer to the integral

$$
\int_{-\infty}^\infty \frac{e^{a x}}{1 + e^x} dx
$$

can be obtained by closing the contour by a large semi circle in the upper or lower HP. But you need to give an argument for why the semi circle part of the integral can be taken to be zero. I'm not sure how rigorous an argument you can give.

Or you can follow @pasmith's suggestion.

Or you can evaluate the integral that @pasmith alludes to by a contour integral.
 
  • #5
julian said:
The correct answer to the integral

$$
\int_{-\infty}^\infty \frac{e^{a x}}{1 + e^x} dx
$$

can be obtained by closing the contour by a large semi circle in the upper or lower HP. But you need to give an argument for why the semi circle part of the integral can be taken to be zero. I'm not sure how rigorous an argument you can give.

The idea would be that [tex]
\begin{split}
\left| \int_0^\pi \frac{e^{aR\cos \theta+iaR\sin\theta}}{1 + e^{R\cos\theta + iR\sin\theta}}iRe^{i\theta}\,d\theta \right|
&\leq \pi R \sup_{\theta \in [0,\pi]} \left| \frac{e^{-(1-a)R\cos \theta}}{e^{-R\cos\theta} + e^{iR\sin\theta}} \right|
\end{split}
[/tex] tends to zero as [itex]R \to \infty[/itex].
 
  • #6
pasmith said:
The idea would be that [tex]
\begin{split}
\left| \int_0^\pi \frac{e^{aR\cos \theta+iaR\sin\theta}}{1 + e^{R\cos\theta + iR\sin\theta}}iRe^{i\theta}\,d\theta \right|
&\leq \pi R \sup_{\theta \in [0,\pi]} \left| \frac{e^{-(1-a)R\cos \theta}}{e^{-R\cos\theta} + e^{iR\sin\theta}} \right|
\end{split}
[/tex] tends to zero as [itex]R \to \infty[/itex].
Do you have to be concerned about ##\theta = \pi/2##?
 
  • #7
julian said:
Do you have to be concerned about ##\theta = \pi/2##?

Indeed.

Multiplying numerator and denominator by [itex]e^{-R\cos\theta} + e^{-iR\sin\theta}[/itex] gives [tex]
\left| \frac{e^{-(1-a)R\cos\theta}}{e^{-R\cos\theta} + e^{iR\sin\theta}}\right| = \frac{e^{-(1-a)R\cos\theta}}{\sqrt{e^{-2R\cos\theta} + 2e^{-R\cos\theta}\cos(R\sin\theta) + 1}}.[/tex] For [itex]\theta = \frac12\pi[/itex] the denominator vanishes for [itex]R = (2n + 1)\pi[/itex], but in this case the pole at [itex]x =( 2n+1)i\pi[/itex] actually lies on the contour, so I think we exclude these values of [itex]R[/itex].

Fortunately, substituting [itex]u = e^x[/itex] and proceeding as in my original post avoids this difficulty.
 
  • #8
Putting ##\theta = \pi/2## into:

$$
\left| \frac{R e^{-(1-a)R\cos \theta} d \theta}{e^{-R\cos\theta} + e^{iR\sin\theta}} \right|
$$

we get

$$
\frac{R}{\sqrt{2 + 2 \cos (R)}} \qquad (*)
$$

Which doesn't have a decaying term: ##e^{-C R}##. Also, as ##R## gets closer to ##(2n+1) \pi## the value of ##(*)## will increase without bound.

We could write

\begin{align*}
\left| \int_0^\pi \frac{e^{aR\cos \theta+iaR\sin\theta}}{1 + e^{R\cos\theta + iR\sin\theta}}iRe^{i\theta}\,d\theta \right|
& \leq \int_0^{\frac{\pi}{2} - \frac{\epsilon (R)}{2}} \left| \frac{R e^{-(1-a)R\cos \theta} d \theta}{e^{-R\cos\theta} + e^{iR\sin\theta}} \right| + \int_{\frac{\pi}{2} + \frac{\epsilon (R)}{2}}^\pi \left| \frac{R e^{aR\cos \theta} d \theta}{1 + e^{R\cos\theta + iR\sin\theta}} \right|
\nonumber \\
& + \left| \int_{\frac{\pi}{2} - \frac{\epsilon (R)}{2}}^{\frac{\pi}{2} + \frac{\epsilon (R)}{2}} \frac{e^{aR\cos \theta+iaR\sin\theta}}{1 + e^{R\cos\theta + iR\sin\theta}} Re^{i\theta}\,d\theta \right|
\end{align*}

for some appropriate ##\epsilon (R)##. The first two integrals on the RHS trivially tend to zero as ##R \rightarrow \infty##. The third integral is zero because the integrand oscillates wildly as ##R \rightarrow \infty##?
 
Last edited:
  • #9
The problem is that we have a sequence of poles [itex](2n + 1)i\pi[/itex]. So instread of taking a continuous liimit [itex]R \to \infty[/itex], we must take a discrete sequence [itex]R_n = 2n\pi[/itex] (say) and let [itex]n \to \infty[/itex].

Let [tex]
f_n(\theta) = \left| \frac{n\pi e^{-(1-a)n\pi \cos\ \theta}}{e^{-n\pi\cos\theta} + e^{in\pi\sin\theta}}\right| =
\frac{n\pi e^{-(1-a)n\pi\cos\theta}}{\sqrt{e^{-2n\pi\cos\theta} + 2e^{-n\pi\cos\theta}\cos(n\pi\sin\theta) + 1}}.[/tex] Then [itex]f(\frac\pi 2) = \frac {n\pi}2[/itex], so [itex]\sup f_n(\theta)[/itex] will not tend to zero; we need a tighter bound.

For [itex]\theta < \pi/2[/itex], the denominator tends to 1 as [itex]n \to \infty[/itex] and the numerator tends to zero. For [itex]\theta > \pi/2[/itex], the denominator tends to [itex]e^{n\pi|\cos\theta|}[/itex] so [itex]f_n (\theta) \sim n\pi e^{-an\pi|\cos\theta|} \to 0[/itex]. So we should be able to find a bound of the form [tex]
f_n(\theta) \leq \begin{cases} m_n & 0 \leq \theta < \frac 12(\pi - \epsilon_n) \\
\frac{n\pi}2 & \frac12 (\pi - \epsilon_n) \leq \theta \leq \frac 12(\pi + \epsilon_n) \\
m_n' & \frac12(\pi + \epsilon_n) < \theta \leq \pi\end{cases}[/tex] where [itex]m_n[/itex], [itex]m_n'[/itex] and [itex]n\pi\epsilon_n[/itex] all tend to zero as [itex]n \to \infty[/itex]. Then [tex]
\int_0^\pi f_n(\theta)\,d\theta \leq \frac{(m_n + m_n')(\pi - \epsilon_n)}{2} + \frac{n\pi}2 \epsilon_n \to 0.[/tex]
 
  • #10
Hi, @pasmith. Did you make slight typo? You said you where going to use ##R_n = 2n \pi##, but you seem to have plugged ##n \pi## into your formula instead.

Assuming you meant to plug in ## R_n = 2n \pi##, I think I have found a bound of the form you postulated...

Starting with

$$
f_n (\theta) = \dfrac{2 n \pi e^{-(1-a) 2n \pi cos \theta}}{\sqrt{e^{-4n \pi cos \theta} + 2 e^{-2n \pi cos \theta} \cos \left( 2 n \pi \sin \theta \right) + 1}}
$$

and using

\begin{align*}
\sin (\frac{\pi}{2} + \phi) = \sin \frac{\pi}{2} \cos \phi + \cos \frac{\pi}{2} \sin \phi = \cos \phi
\nonumber \\
\cos (\frac{\pi}{2} + \phi) = \cos \frac{\pi}{2} \cos \phi - \sin \frac{\pi}{2} \sin \phi = - \sin \phi
\end{align*}

we have

$$
\frac{f_n (\frac{\pi}{2} + \phi)}{n \pi} = \dfrac{2 e^{(1-a) 2n \pi \sin \phi}}{\sqrt{e^{4n \pi \sin \phi} + 2 e^{2n \pi \sin \phi} \cos \left( 2 n \pi \cos \phi \right) + 1}} \qquad (*)
$$

First assume ##\phi \geq 0##. For ##0 \leq \phi \leq \frac{\pi}{2}## we have ##\sin \phi \leq \phi## and ##\cos \phi \geq 1 - \frac{2}{\pi} \phi##. Note for ##\phi## close enough to ##0##,

\begin{align*}
\cos \left( 2 n \pi \cos \phi \right) & \geq \cos \left( 2 n \pi \left(1 - \frac{2}{\pi} \phi \right) \right)
\nonumber \\
& = \cos \left( 2 n \pi - 4 n \phi \right)
\nonumber \\
& = \cos \left( 4 n \phi \right)
\end{align*}

This inequality obviously holds on the interval ##[0, \frac{\pi}{8n}]## (at ##\phi = \frac{\pi}{8n}##, ##\cos \left( 4 n \phi \right)## becomes zero).

Using these inequalities in ##(*)##, we have

\begin{align*}
\frac{f_n (\frac{\pi}{2} + \phi)}{n \pi} & \leq \dfrac{2 e^{(1-a) 2n \pi \phi}}{\sqrt{1 + 2 \cos \left( 2 n \pi (1 - \frac{2}{\pi} \phi) \right) + 1}}
\nonumber \\
& = \dfrac{2 e^{(1-a) 2n \pi \phi}}{\sqrt{2 + 2 \cos \left( 4n \phi \right)}}
\end{align*}

while ##0 \leq \phi \leq \phi_{max} = \frac{\pi}{8n}##, with the RHS of this inequality monotonically increasing on ##[0,\frac{\pi}{8n}]##. And so

$$
\frac{f_n (\frac{\pi}{2} + \phi)}{n \pi} \leq \sqrt{2} e^{(1-a) \pi^2 / 4}
$$

while ##0 \leq \phi \leq \phi_{max} = \frac{\pi}{8n}##.

Now assume ##\phi \leq 0##. Write

\begin{align*}
\frac{f_n (\frac{\pi}{2} + \phi)}{n \pi} & = \dfrac{2 e^{-a 2n \pi \sin \phi}}{\sqrt{e^{-4n \pi \sin \phi} + 2 e^{-2n \pi \sin \phi} \cos \left( 2 n \pi \cos \phi \right) + 1}}
\nonumber \\
& = \dfrac{2 e^{a 2n \pi |\sin \phi|}}{\sqrt{e^{4n \pi |\sin \phi|} + 2 e^{2n \pi |\sin \phi|} \cos \left( 2 n \pi \cos \phi \right) + 1}}
\end{align*}

By similar reasoning, we have

$$
\frac{f_n (\frac{\pi}{2} + \phi)}{n \pi} \leq \sqrt{2} e^{a \pi^2 / 4}
$$

while ##- \frac{\pi}{8n} = - \phi_{max} \leq \phi \leq 0##.

Put

$$
M = max \left( \sqrt{2} e^{(1-a) \pi^2 / 4} , \sqrt{2} e^{a \pi^2 / 4} \right)
$$

we have for ##\frac{1}{2} \epsilon_n \leq \phi_{max} = \frac{\pi}{8n}##,

$$
f_n (\theta) \leq \left\{
\begin{matrix}
m_n & 0 \leq \theta < \frac{1}{2} \left( \pi - \epsilon_n \right) \\
n \pi M & \frac{1}{2} \left( \pi - \epsilon_n \right) \leq \theta \leq \frac{1}{2} \left( \pi + \epsilon_n \right) \\
m_n' & \frac{1}{2} \left( \pi + \epsilon_n \right) < \theta \leq \pi
\end{matrix}
\right.
$$

We choose:

$$
\epsilon_n = \frac{\pi}{4 n^2}
$$

so that

$$
n \pi M \epsilon_n \rightarrow 0 .
$$
 
Last edited:
  • #11
Hang on, I don't think it works. Because then we couldn't have ##m_n, m_n' \rightarrow 0##.

Generally, not sure it can work. If we want

$$
n \pi M \epsilon_n \rightarrow 0
$$

then we need ##\epsilon_n = \frac{C}{n^\beta}## with ##\beta > 1##. But then

\begin{align*}
f_n ( \frac{\pi}{2} - \frac{C}{2n^\beta}) & =
\dfrac{2 n \pi e^{-(1-a) 2n \pi \cos( \frac{\pi}{2} - \frac{C}{2n^\beta})}}{\sqrt{e^{-4n \pi \cos (\frac{\pi}{2} - \frac{C}{2n^\beta})} + 2 e^{-2n \pi \cos (\frac{\pi}{2} - \frac{C}{2n^\beta})} \cos \left( 2 n \pi \sin (\frac{\pi}{2} - \frac{C}{2n^\beta}) \right) + 1}}
\nonumber \\
& = \dfrac{2 n \pi e^{-(1-a) 2n \pi \sin (\frac{C}{2n^\beta})}}{\sqrt{e^{-4n \pi \sin (\frac{C}{2n^\beta})} + 2 e^{-2n \pi \sin (\frac{C}{2n^\beta})} \cos \left( 2 n \pi \cos (\frac{C}{2n^\beta}) \right) + 1}}
\nonumber \\
& \rightarrow \dfrac{2 n \pi e^{-(1-a) \pi \frac{C}{n^{\beta - 1}}}}{\sqrt{e^{-2 \pi \frac{C}{n^{\beta - 1}}} + 2 e^{- \pi \frac{C}{n^{\beta} - 1}} + 1}}
\nonumber \\
& \rightarrow \infty
\end{align*}
 
Last edited:
  • #12
Seeing as I had already given an answer in post #3, I gave a suggestion as how to evaluate the integral via another method in post #4 (the third suggestion). Seeing as the two weeks are nearly up, I'll post how to do it this other way:

First use the substitution ##u = e^x## and obtain

$$
I = \int_{-\infty}^\infty \frac{e^{a x}}{1 + e^x} dx = \int_0^\infty \frac{u^{a-1}}{1+u} du
$$

This ##u## integral can be evaluated via a contour integral. First note that the function

$$
f(z) = \frac{z^{a-1}}{1+z}
$$

has a branch point at ##z=0##. We take the branch line to be the positive real axis. We employ the contour integral shown in the figure. We integrate along the positive real axis just above the branch cut, then around a large circle, then along the positive real axis in the negative direction just below the cut, but on the same branch, and then around a small circle.

contour.jpg


Just below the cut, but on the same branch, we have

$$
\frac{(e^{i 2 \pi} u)^{a-1} }{1+e^{i 2 \pi} u} = e^{i 2 \pi (a-1)} \frac{u^{a-1}}{1+u} = e^{i 2 \pi a} \frac{u^{a-1}}{1+u}
$$

For the integration around a small circle around the origin put ##z = \epsilon e^{i \theta}##, then

\begin{align*}
\left| \int_0^{2 \pi} \frac{\epsilon^a e^{i \theta a} i}{1 + \epsilon e^{i \theta}} d \theta \right| & \leq \epsilon^a \int_0^{2 \pi} \frac{d \theta}{\sqrt{\epsilon^2 + 2 \epsilon \cos \theta + 1}}
\nonumber \\
& \leq \epsilon^a \int_0^{2 \pi} \frac{d \theta}{\sqrt{\epsilon^2 - 2 \epsilon + 1}}
\nonumber \\
& = 2 \pi \epsilon^a \frac{1}{1 - \epsilon} \rightarrow 0 \quad \text{as} \; \epsilon \rightarrow 0
\end{align*}

For integration around circle of large radius put ##z = R e^{i \theta}##, then

\begin{align*}
\left| \int_0^{2 \pi} \frac{R^a e^{i \theta a} i}{1 + R e^{i \theta}} d \theta \right| & \leq \int_0^{2 \pi} \frac{R^{a-1} d \theta}{\sqrt{R^{-2} + 2 R^{-1} \cos \theta + 1}}
\nonumber \\
& \leq \int_0^{2 \pi} \frac{R^{a-1} d \theta}{\sqrt{R^{-2} - 2 R^{-1} + 1}}
\nonumber \\
& = 2 \pi \frac{R^{a-1}}{1 - R^{-1}} \rightarrow 0 \quad \text{as} \; R \rightarrow \infty
\end{align*}

Therefore,

$$
\oint \frac{z^{a-1}}{1+z} = \int_0^\infty \frac{u^{a-1}}{1+u} du - e^{i 2 \pi a} \int_0^\infty \frac{u^{a-1}}{1+u} du
$$

The function ##f(z)## has a simple pole at ##z=-1##, so that

$$
\oint \frac{z^{a-1}}{1+z} =2 \pi i (e^{i \pi})^{(a-1)} = - 2 \pi i e^{i \pi a} .
$$

Combining the above gives

$$
(1 - e^{i 2 \pi a}) \int_0^\infty \frac{u^{a-1}}{1+u} du = - 2 \pi i e^{i \pi a}
$$

or

$$
\int_0^\infty \frac{u^{a-1}}{1+u} du = \frac{\pi}{\sin a \pi} .
$$
 
  • Like
Likes Euge

Related to Exponential-Type Integrals

1. What are exponential-type integrals?

Exponential-type integrals are mathematical expressions that involve an exponential function raised to a power. They are commonly used in many areas of science and engineering, such as in the study of growth and decay processes.

2. How are exponential-type integrals different from other types of integrals?

Exponential-type integrals are different from other types of integrals because they involve an exponential function, which can make the integration process more challenging. They also have unique properties and techniques for solving them.

3. What are some common applications of exponential-type integrals?

Exponential-type integrals have many applications in science and engineering, such as in the study of radioactive decay, population growth, and signal processing. They are also used in fields like physics, chemistry, and economics.

4. How are exponential-type integrals used in real-world problems?

Exponential-type integrals are used in real-world problems to model and predict various phenomena, such as the growth or decay of populations, the spread of diseases, and the behavior of electrical circuits. They are also used in statistical analysis and data fitting.

5. What are some techniques for solving exponential-type integrals?

There are several techniques for solving exponential-type integrals, including substitution, integration by parts, and the use of special functions such as the gamma function. These techniques can be combined and applied depending on the specific form of the integral.

Similar threads

  • Math POTW for University Students
Replies
3
Views
776
  • Math POTW for University Students
Replies
6
Views
929
  • Math POTW for University Students
Replies
8
Views
691
  • Math POTW for University Students
Replies
1
Views
2K
  • Math POTW for University Students
Replies
3
Views
660
  • Math POTW for University Students
Replies
3
Views
614
  • Math POTW for University Students
Replies
15
Views
598
  • Math POTW for University Students
Replies
5
Views
936
Replies
1
Views
2K
  • Calculus and Beyond Homework Help
Replies
7
Views
960
Back
Top