Analytic Approximation for an Oscillatory Integral

  • Context: Graduate 
  • Thread starter Thread starter csmallw
  • Start date Start date
  • Tags Tags
    Approximation Integral
Click For Summary
SUMMARY

The forum discussion focuses on deriving an analytic approximation for the oscillatory integral defined as \int_0^\infty \frac{k \sin(kr)}{\sqrt{1+v^2(k-k_F)^2}}dk, particularly under the assumption that v k_F >> 1. The user reports unsuccessful attempts with Mathematica but has generated numerical solutions. A closed-form upper bound is provided as (2 k_F ArcSinh[k_F v])/v, while the lower bound is suggested to be the negation of the upper bound. The discussion also explores the potential of integrating by parts to refine the approximation, especially for large values of r.

PREREQUISITES
  • Understanding of oscillatory integrals and their properties
  • Familiarity with Fourier transforms and their applications
  • Proficiency in using Mathematica for symbolic computation
  • Knowledge of asymptotic analysis techniques
NEXT STEPS
  • Explore advanced techniques in asymptotic analysis for oscillatory integrals
  • Learn about the properties of Fourier transforms in the context of physical applications
  • Investigate the use of integration by parts in deriving approximations for integrals
  • Experiment with Mathematica's numerical integration capabilities for complex integrands
USEFUL FOR

Mathematicians, physicists, and engineers interested in analytic approximations for oscillatory integrals, particularly in the context of Fourier transforms and numerical analysis.

csmallw
Messages
25
Reaction score
0
I'm looking for a way to write down an analytic approximation for the following integral:

[tex]\int_0^\infty \frac{k \sin(kr)}{\sqrt{1+v^2(k-k_F)^2}}dk[/tex]

Let's assume that v kF >> 1, so that the the oscillating piece at large k doesn't contribute much uncertainty. Ideas? Thus far, Mathematica has failed me, though I have been able to generate some numerical solutions. Is there some way to take advantage of the fact that the integrand is peaked at kF?
 
Physics news on Phys.org
Code:
In[1]:= $Assumptions = v > 0 && kf > 0 && r > 0;
 Integrate[k /Sqrt[1 + v^2 (k - kf)^2], {k, 0, 2 kf}]

Out[1]= (2 kf ArcSinh[kf v])/v

Your integral will be strictly less than that, but at least it gives a closed form upper bound and it doesn't depend on substituting in some arbitrary values for constants. The only lower bound I can see from this would be negating the upper bound. That also doesn't integrate out to infinity, but from a few numerical examples it looks like by the time you are out to 2kf or some other appropriate multiple of kf or v kf that, as you have mentioned, you are far from the peak and the additional contribution may be modest.

You can substitute a variety of reasonable constants and plot your expression and the expression without the Sin and see how these behave.
 
The integral I am trying to solve is a version of a Fourier transform, so it would be better if the approximation were r-dependent.

I was thinking that integrating by parts twice would sharpen the 1/√(...) piece and make it look like a delta-function. Then the integration would be easy if I constrained "r" to be sufficiently small (and if the boundary terms don't diverge). Still, I think I'm stuck in terms of providing an approximation at large "r," which is where I might be most interested in the integral's value...
 

Similar threads

  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
1K
  • · Replies 7 ·
Replies
7
Views
4K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 9 ·
Replies
9
Views
3K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 1 ·
Replies
1
Views
1K