How Do You Compute the Complex Gaussian Integral I = ∫ e^(-ax^2 + ibx) dx?

homer
Messages
46
Reaction score
0

Homework Statement


Let a,b be real with a > 0. Compute the integral
<br /> I = \int_{-\infty}^{\infty} e^{-ax^2 + ibx}\,dx.<br />

Homework Equations


Equation (1):
\int_{-\infty}^{\infty} e^{-x^2}\,dx = \sqrt{\pi}

Equation (2):
-ax^2 + ibx = -a\Big(x - \frac{ib}{2a}\Big)^2 - \frac{b^2}{4a}

The Attempt at a Solution


Completing the square in -ax^2 + ibx gives me Equation (2), so that my integral is now
<br /> I = e^{-b^2/4a}\int_{-\infty}^{\infty} e^{-a(x-ib/2a)^2}\,dx.<br />
Making the substitution u = \sqrt{a}(x-ib/2a) I get du = \sqrt{a}\,dx so tht my integral becomes
<br /> I = \frac{1}{\sqrt{a}}\,e^{-b^2/4a}\int_{-\infty - ib/2\sqrt{a}}^{\infty - ib/2\sqrt{a}} e^{-u^2}\,du.<br />

But this doesn't seem right.
 
Physics news on Phys.org
You've got it right. Just note that ##-\infty - ib/2\sqrt{a}=-\infty## and ##\infty - ib/2\sqrt{a}=\infty##, and you're basically done.
 
  • Like
Likes homer
Thanks Zeta. Huge brain fart on my part in making it rigorous. The integral I was trying to compute is the limit of

I(R_1, R_2) = \int_{-R_1}^{R_2} e^{-ax^2 + bx}\,dx

as R_1, R_2 \to \infty. Then I can make a subsitution z = \sqrt{a}(x-ib/2a) to get the integral

I(R_1,R_2) = \frac{1}{\sqrt{a}}\,e^{-b^2/4a}\int_{-\sqrt{a}R_1-ib/2\sqrt{a}}^{\sqrt{a}R_2 - ib/2\sqrt{a}}e^{-z^2}\,dz.

The integrand e^{-z^2} is analytic on the entire complex plane, so the integral is path independent. So in particular I can take it on a contour consisting of:

(1) A straight line up from z = -\sqrt{a}R_1 - ib/\sqrt{a} up to the real axis at point z = -\sqrt{a}R_1.

(2) A straight line on the real axis from z = -\sqrt{a}R_1 to z = \sqrt{a}R_2.

(3) A straight line down from the real axis at point z = \sqrt{a}R_2 to z = \sqrt{a}R_2 - ib/2\sqrt{a}.

Taking the limit as R_1, R_2 \to \infty the integrals on contour sections (1) and (3) vanish since \lvert z\rvert \to \infty and thus e^{-z^2} \to 0 on these two vertical lines. The integral on contour (2) then becomes

\lim_{R_1, R_2 \to \infty}\int_2 e^{-z^2}\,dz = \int_{-\infty}^{\infty} e^{-x^2}\,dx = \sqrt{\pi}.

Thus we have

\begin{align*}
\int_{-\infty}^{\infty} e^{-ax^2 + bx}\,dx
& = \lim_{R_1,R_2 \to \infty} I(R_1,R_2) \\
& =
\lim_{R_1,R_2 \to \infty} \frac{1}{\sqrt{a}}\,e^{-b^2/4a}\int_{-\sqrt{a}R_1-ib/2\sqrt{a}}^{\sqrt{a}R_2 - ib/2\sqrt{a}}e^{-z^2}\,dz \\
& = \frac{1}{\sqrt{a}}\,e^{-b^2/4a}\sqrt{\pi}.
\end{align*}
 
Oops, contour (1) should be from z = -\sqrt{a}R_1 - ib/2\sqrt{a} to z = -\sqrt{a}R_1.
 
To solve this, I first used the units to work out that a= m* a/m, i.e. t=z/λ. This would allow you to determine the time duration within an interval section by section and then add this to the previous ones to obtain the age of the respective layer. However, this would require a constant thickness per year for each interval. However, since this is most likely not the case, my next consideration was that the age must be the integral of a 1/λ(z) function, which I cannot model.
Back
Top