Fourier transform of Langevin equation (integral cancellation problem)

Click For Summary

Discussion Overview

The discussion revolves around the application of Fourier transforms to the Langevin equation in the context of modeling the motion of a tethered bead. Participants explore the mathematical steps involved in transforming the equation and the implications for deriving the power spectrum.

Discussion Character

  • Technical explanation
  • Mathematical reasoning
  • Debate/contested

Main Points Raised

  • One participant, Stephen, questions whether it is valid to cancel the integral parts in Eq. 6 to simplify to Eq. 7, expressing uncertainty about the Fourier transform process.
  • Another participant suggests an alternative approach by taking the Fourier transform of the original equation directly, proposing different forms for the Fourier transforms of the variables involved.
  • A third participant challenges the correctness of the power spectrum equation provided by Stephen, suggesting it should be expressed as the product of the function and its complex conjugate.

Areas of Agreement / Disagreement

Participants do not reach a consensus on the validity of canceling the integrals in Eq. 6, and there are differing opinions on the correct formulation of the power spectrum. The discussion remains unresolved regarding these points.

Contextual Notes

There are unresolved assumptions regarding the application of Fourier transforms and the conditions under which the integrals can be manipulated. The discussion also highlights potential discrepancies in the formulation of the power spectrum.

stephenx_86
Messages
14
Reaction score
0
Hi,

(To cut a long story short, can I cancel the integrals in Eq. 6 to leave me with Eq. 7?)

I am trying to follow the method for modelling the motion of a tethered bead from a couple of papers ("Te Velthuis, A. J. W. et al. (2010) Biophys. J. 99 1292–1302" and "Lansdorp, B. M., & Saleh, O. A. (2012) Rev. Sci. Instrum. 83 025115"), but am getting stuck going from the Langevin equation (where γ is the friction coefficient and k is the tether stiffness):

[1] [itex]\gamma \dot{x} (t)+kx(t)=F_{therm}[/itex]

to the power spectrum:

[2] [itex]P(\omega )\equiv \left |x^{2}(\omega )\right| =\frac{2\gamma k_{B}T}{\gamma ^{2}\omega ^{2}+k^{2}}[/itex]

I've followed a few lecture handouts I've found online and think I understand the majority of the problem. From what I understand, I need to take the three components of the first equation and individually apply the Fourier transforms (all in the range ∞→-∞):

[3] [itex]\dot{x}(t)=\frac{d}{dt}\left|\frac{1}{2\pi }\int x(\omega )e^{i\omega t}d\omega \right|=\frac{1}{2\pi }\int x(\omega ) i\omega e^{i\omega t}d\omega[/itex]

[4] [itex]x(t)=\frac{1}{2\pi }\int x(\omega )e^{i\omega t}d\omega[/itex]

[5] [itex]F(t)=\frac{1}{2\pi }\int F(\omega )e^{i\omega t}d\omega[/itex]

If I substitute all these back into the original equation I get:

[6] [itex]\gamma \frac{1}{2\pi }\int x(\omega ) i\omega e^{i\omega t}d\omega +k\frac{1}{2\pi }\int x(\omega )e^{i\omega t}d\omega = \frac{1}{2\pi }\int F(\omega )e^{i\omega t}d\omega[/itex]

From the lecture handouts I've found, I should be able to rearrange this to get the following (or something similar):

[7] [itex]\gamma i\omega x(\omega) +kx(\omega) = F(\omega)[/itex]

Now, here's my question. Is it possible for me to cancel the parts of the integral that appear in all three components of in Eq. 6 (this is [itex]\frac{1}{2\pi}\int e^{i\omega t} d\omega[/itex]) to leave me with Eq. 7?

Any help on this would be greatly appreciated, since I've been stuck on this for a couple of days now!

Thanks
Stephen

p.s. I also seem to end up with a rogue i, which will become -1 when I square P(x). Have I done my Fourier transform equations wrong?
 
Physics news on Phys.org
Yes, but it is not the way I would do it.

Instead, I would just take the Fourier transform of the original equation.

instead of

stephenx_86 said:
[itex]x(t)=\frac{1}{2\pi }\int x(\omega )e^{i\omega t}d\omega[/itex]

use this
[itex]x(\omega )=\int x(t )e^{-i\omega t}dt[/itex]
[itex]i\omega x(\omega )=\int x'(t )e^{-i\omega t}dt[/itex]
[itex]F_{therm}(\omega )=\int F_{therm}(t )e^{-i\omega t}dt[/itex]

the FT is linear:
[itex]y(\omega ) + z(\omega )=\int (y(t) + z(t))e^{-i\omega t}dt[/itex]

Just apply the FT to both sides original equation.
 
As for the question about power spectrum, I think the equation is wrong. It should be
[itex]P(\omega) = \left|x(\omega) \right|^2 = (x(\omega))^* x(\omega)[/itex]
 
That's brilliant, thanks for your help.
 

Similar threads

  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 4 ·
Replies
4
Views
2K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 5 ·
Replies
5
Views
2K
Replies
1
Views
3K
  • · Replies 2 ·
Replies
2
Views
1K
  • · Replies 5 ·
Replies
5
Views
2K
  • · Replies 11 ·
Replies
11
Views
3K
  • · Replies 3 ·
Replies
3
Views
2K
  • · Replies 3 ·
Replies
3
Views
4K