One integral (square root and gaussian peak)

[jostpuur]

An integral
<br /> \int\limits_{-\infty}^{\infty} \sqrt{1+x^2} e^{-Ax^2+Bx} dx<br />
is the problem. Anyone knowing the solution?

 

[AiRAVATA]

Try completing squares...

 

[jostpuur]

A variable change that simplifies the exponent only makes the expression under the square root more complicated, like this
<br /> \int\limits_{-\infty}^{\infty} \sqrt{x^2+Ax+B}e^{-x^2} dx<br />
which doesn't seem any easier. (A and B are different constants than in the previous equation). Or was this what you meant by completing the squares?

 

[CompuChip]

I think he meant this (also check out the example on that page), although I must admit I don't immediately see how that helps you... think about it :)

 

[jostpuur]

An equation
<br /> -Ax^2 + Bx = -\big(\sqrt{A}x - \frac{B}{2\sqrt{A}}\big)^2 + \frac{B^2}{4A}<br />
would naturally have preceded my variable change
<br /> x=\frac{1}{\sqrt{A}}y + \frac{B}{2A}<br />
and renaming y\mapsto x, in simplifying the exponent.

 

[Gib Z]

Sorry to tell you but nothing you've done helps you very much. In fact the only thing that might help in the slightest was post 3, where one can at least conclude that the integral even exists.

The integrand has no elementary anti derivative, or even one in known non-elementary functions. You have to use more powerful tools.

Define:

I(a)= \int^{a}_{-a} \sqrt{1+x^2} e^{-Ax^2 + Bx} dx = \int^{a}_{-a} \sqrt{1+y^2} e^{-Ay^2+By} dy

I^2(a) = \int^{a}_{-a} \int^{a}_{-a} \sqrt{(1+x^2)(1+y^2)} e^{ -A(x^2+y^2) +B(x+y)} dx dy.

Now applying Fubini's Theorem we can see this double integral is the same as the integral:

\int_C \sqrt{(1+x^2)(1+y^2)} e^{ -A(x^2+y^2) +B(x+y)} d(x,y), where C is the square with vertices {(−a, a), (a, a), (a, −a), (−a, −a)} on the xy-plane.

Now switch to polar co-ordinates, sox=r\,\cos \theta, y= r\,\sin\theta, d(x,y) = r\, d(r,\theta).

Hopefully it will work out.

 

[tookan]

Gaussian integral with square root problem

Hi - I don't suppose that anyone figured this out yet? Been giving me a hard time for a while now - I can't see how the transformation to polar helps - still left with a gaussian and square root product in the integrand.

 

[jostpuur]

Hi - I don't suppose that anyone figured this out yet? Been giving me a hard time for a while now - I can't see how the transformation to polar helps - still left with a gaussian and square root product in the integrand.

Do you mean that you took this as a challenge now when you saw this post, or that you had been trying to solve this same integral earlier, and found this post when you were trying to find solution?

The GibZ's hint doesn't work, I checked it. Looks like the circles (in shperical coordinates) should be replaced with some other shapes. In fact, with those shapes, that satisfy

<br /> (1+x^2)(1+y^2)=c^2<br />

for some constant c, but I couldn't make that stuff work.

Other idea I had, was to define

<br /> I(A,B) = \int\limits_{-\infty}^{\infty} dx\;\sqrt{1+x^2}e^{-Ax^2+Bx}<br />

Accept I(A_0,0) as some special constant, where A_0 is some constant, and try Taylor series

<br /> I(A,B)=I(A_0,0) + (A-A_0,B)\cdot(\partial_A, \partial_B) I(A_0,0) + \frac{1}{2}\big( (A-A_0,B)\cdot(\partial_A,\partial_B)\big)^2 I(A_0,0) +\cdots<br />

but the partial derivatives produced just more difficult integrals, that I couldn't make return to the original one.

I haven't been trying much after these attempts.

 

[CompuChip]

Is this integral even analytically solvable?

 

[jostpuur]

Is this integral even analytically solvable?

I don't know how to prove such things, but the Wolfram Integrator at least doesn't give integral function for that. That of course doesn't mean, that it would be impossible to carry out the integration over some given domain.

 

[Gib Z]

When one runs out of ideas, crazy things come to mind. Perhaps the squeeze theorem on the exponent will turn out fine... Give me some time on this ill actually bother to do this now.

 

[jostpuur]

When one runs out of ideas, crazy things come to mind. Perhaps the squeeze theorem on the exponent will turn out fine... Give me some time on this ill actually bother to do this now.

GibZ, are you interested in physics? We had some interesting discussion about quantum theory in threads

https://www.physicsforums.com/showthread.php?t=175155

https://www.physicsforums.com/showthread.php?t=176514&page=3

and about the operator \sqrt{-\nabla^2+m^2}. I though, that being able to calculate this integral, would bring some understanding to how this operator behaves. That was where I got the original motivation.

 

[tookan]

integral fun...

Hey people - I actually came across this integral whilst doing a functional integral calculation for the truncated (in real lattice space) Hubbard model. If you don't know then the Hubbard model is a lattice model for stongly correlated electron systems in condensed matter. If I can get this integral then some serious physics should drop out analytically. The problem may be very similar in form to some different field theories. jostpuur - nice one for your help - I had come to the conclusion as well that some sort of expansion is the way to proceed but the details of my specific problem (where x is a complicated function of multivalued auzilliary fields) suggest that the square root remains in the answer! maybe can expand integral, do individual integrals, the gaussian integrals with higher powers of x as a product in the integrand are gamma integrals and can be done easily, and then regroup...?

 

[jostpuur]

One way would be to write \sqrt{1+x^2} as Taylor series. Unfortunately those series don't converge everywhere. An ugly way would be to press on by integrating small pieces somewhere near the origo always using some series that converge there. This would not give precise solution, but an approximation at least.

One might also ask, that can it be written as a different kind of series, that converge everywhere, and whose terms can be integrated.