How can I properly compute this contour integral with the residue theorem?

Click For Summary
SUMMARY

The discussion focuses on computing the integral \(\int^{\infty}_{-\infty}\frac{e^{ix}dx}{\sqrt{a - x^2 + ib}}\) using the residue theorem and Cauchy's theorem. Participants suggest transforming the integral by letting \(a + ib = c^2\) to identify branch points at \(\pm c\). The integral is evaluated along a semicircular contour that loops around these branch points, leading to expressions involving branch cuts extending towards infinity. Numerical evaluations in Mathematica reveal discrepancies likely due to branch selection.

PREREQUISITES
  • Complex analysis fundamentals, specifically contour integration.
  • Understanding of the residue theorem and Cauchy's theorem.
  • Familiarity with branch cuts and multi-valued functions in complex analysis.
  • Experience with numerical integration tools, particularly Mathematica.
NEXT STEPS
  • Study the application of Cauchy's theorem in complex integrals.
  • Learn about branch cuts and their implications in complex analysis.
  • Explore the residue theorem in greater depth, focusing on its applications in evaluating integrals.
  • Investigate numerical integration techniques in Mathematica, particularly handling multi-valued functions.
USEFUL FOR

Mathematicians, physicists, and students engaged in complex analysis, particularly those working on integrals involving branch points and contour integration techniques.

freechus9
Messages
8
Reaction score
0
For my research I am trying to solve a particular integral. I know that I have to transform it properly to use the residue theorem, but I am having difficulty attaining this form. Here is the integral:

\int^{\infty}_{-\infty}\frac{e^{ix}dx}{\sqrt{a - x^2 + ib}}

where a and b are real constants. If you have any suggestions I would be more than grateful.

Thanks in advance.
 
Physics news on Phys.org
Let a+ib =c^2, so that the integrand has poles at +-c.Integrate along a semicircle which contains a part of the real axis and loops around +-c.
 
Then \pm c would be branch points, not poles. In this case, the integral could be evaluated using Cauchy's Theorem but not the Residue Theorem. How about doing this first: Just let a+ib=1+i. Then in order to capture the entire real axis, let the branch cuts extend outward towards infinity from each branch-point. So the one in the upper half-plane would be 2^{1/4}e^{\pi i/8}. Now, let the branch-cut from that one extend along the line re^{\pi i/8}. Now, close the contour in a semi-circular arc which traverses over that branch-cut and just compute the two terms along the branch-cut. Does it agree with direct numerical integration of the integral?
 
I understand what you mean by which lines the cuts are on, but I am having trouble explicitly evaluating the contributions from either side of the cut in the upper half plane.

For example, along the cut on the right side, x would have the value re^{i*\pi/8}. Therefore x^{2} is r^{2}e^{i*\pi/4}?

Along the cut on the left side, x would have the value re^{-i*15\pi/8}. Therefore x^{2} is r^{2}e^{-i*15\pi/4}?

Is this correct?
 
The value on one side of a square-root branch-cut is just the negative of the other side. Remember, the multifunctioned square root is \pm. So on the lower side, the integral would be:

\int_{\infty}^{2^{1/4}} \frac{e^{iz}}{\sqrt{z_0-z^2}} e^{\pi i/8}dr,\quad z\to re^{\pi i/8}

and on the other side, it would be:

\int_{2^{1/4}}^{\infty} \frac{e^{iz}}{-\sqrt{z_0-z^2}} e^{\pi i/8}dr,\quad z\to re^{\pi i/8}

and after analyzing all the other parts of the contour and simplifying, we can write:

\int_{-\infty}^{\infty} \frac{e^{iz}}{\sqrt{z_0-z^2}}dz=2\int_{\infty}^{2^{1/4}} \frac{e^{iz}}{-\sqrt{z_0-z^2}} e^{\pi i/8}dr,\quad z\to re^{\pi i/8}

at least for this problem. May have other issues if z0 is on the real axis or may have to modify the analysis somewhat.

Not sure if that's any easier to evaluate analytically than the original problem. May be or probably is an easier way to evaluate it but it's good practice to do it any way possible to help understand and perhaps lead to a better way of doing it.
 
Last edited:
Thank you very much!
 
jackmell said:
and after analyzing all the other parts of the contour and simplifying, we can write:

\int_{-\infty}^{\infty} \frac{e^{iz}}{\sqrt{z_0-z^2}}dz=2\int_{\infty}^{2^{1/4}} \frac{e^{iz}}{-\sqrt{z_0-z^2}} e^{\pi i/8}dr,\quad z\to re^{\pi i/8}

That's not correct. Remember it's the sum over all the contours and by Cauchy's Theorem, that's zero so I should have written:

\int_{-\infty}^{\infty} \frac{e^{iz}}{\sqrt{z_0-z^2}}dz+\left(2\int_{\infty}^{2^{1/4}} \frac{e^{iz}}{-\sqrt{z_0-z^2}} e^{\pi i/8}dr,\quad z\to re^{\pi i/8}\right)=0

so that:

\int_{-\infty}^{\infty} \frac{e^{iz}}{\sqrt{z_0-z^2}}dz=\left(-2\int_{\infty}^{2^{1/4}} \frac{e^{iz}}{-\sqrt{z_0-z^2}} e^{\pi i/8}dr,\quad z\to re^{\pi i/8}\right)

However when I compute them numerically in Mathematica, there is a sign difference. I suspect that's because Mathematica is integrating over two different branches but I'm not sure.
 

Similar threads

  • · Replies 8 ·
Replies
8
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 2 ·
Replies
2
Views
6K
  • · Replies 9 ·
Replies
9
Views
2K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 19 ·
Replies
19
Views
4K
  • · Replies 3 ·
Replies
3
Views
4K