Deriving the expression for the diffusion constant

Click For Summary
SUMMARY

The discussion centers on deriving the expression for the diffusion constant, D, as presented in R. Kubo's paper on the fluctuation-dissipation theorem. The key equation is D=lim_{t\rightarrow\infty}\frac{1}{2t}<\{x(t)-x(0)\}^2>, which is transformed through variable substitution and integration of velocity correlations. The participant seeks clarification on the transition to D=lim_{t\rightarrow\infty}\frac{1}{t}\int{dt_{1}}\int{dt', specifically regarding the upper limit of integration and the application of Fubini's theorem to evaluate the double integral.

PREREQUISITES
  • Understanding of the fluctuation-dissipation theorem
  • Familiarity with integral calculus and double integrals
  • Knowledge of statistical mechanics and diffusion processes
  • Experience with variable substitution in mathematical expressions
NEXT STEPS
  • Study the fluctuation-dissipation theorem in detail, focusing on its applications in statistical mechanics
  • Learn about Fubini's theorem and its implications for evaluating double integrals
  • Explore the mathematical derivation of diffusion constants in various physical systems
  • Investigate the role of velocity correlation functions in diffusion theory
USEFUL FOR

Physicists, mathematicians, and researchers in statistical mechanics who are interested in understanding the mathematical foundations of diffusion processes and the fluctuation-dissipation theorem.

csk08
Messages
1
Reaction score
0
I have a question of the form, "how did they go from this step to that step?" I am following a paper from R. Kubo, called "The fluctuation-dissipation theorem". I suspect variable substitution and a simple assumption are used, but I don't know how. Here is the basic expression from which everything else follows.

[tex]D=lim_{t\rightarrow\infty}\frac{1}{2t}<\{x(t)-x(0)\}^2>[/tex]​

Since we have

[tex]x(t)-x(0)=\int{v(t')dt'}[/tex]​

where the limits of the integration are 0 to t, we substitute the second equation into the first equation to get

[tex]D=lim_{t\rightarrow\infty}\frac{1}{2t}\int{\int{dt_{1}dt_{2}<v(t_{1})v(t_{2})}}>[/tex]​

again, the limits of both integrals are from 0 to t. (Here comes the step I don't understand)

[tex]D=lim_{t\rightarrow\infty}\frac{1}{t}\int{dt_{1}}\int{dt'<v(t_{1})v(t_{1}+t')}>[/tex]​

where the limits of the outer integral are from 0 to t, and the limits of the inner integral are from 0 to [tex](t-t_{1})[/tex]

It looks like the substitution [tex]t'=t_{2}-t_{1}[/tex] was made. But what happened to the upper limit of integration?
 
Last edited:
Physics news on Phys.org
I think it is not that difficult. let's call [tex]t_1=x[/tex] and [tex]t_2=y[/tex] and:

[tex]<v(x) v(y)> = f(x,y)[/tex]​
The integral

[tex]I = \int \int dx dy f(x,y)[/tex]​
is a 'volume integral' it should be called

[tex]\int_{R(t)} f(x,y)[/tex]​
Now notice that [tex]f(x,y) = f(y,x)[/tex], and since we want to integrate this over the box (bound by x=0, y=0, x=t and y=t) it is equal to integrate this twice over the triangle bound by x = 0, y=x, y=t. Then by Fubini's theorem

[tex]I = 2 \int_{T(t)} f(x,y) = 2 \int_0^t dx \int_{x}^{t} dy f(x,y)[/tex]​
Now substitute [tex]y=x+t'[/tex] then you should get the result.
 

Similar threads

  • · Replies 6 ·
Replies
6
Views
4K
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 0 ·
Replies
0
Views
466
  • · Replies 1 ·
Replies
1
Views
3K
  • · Replies 4 ·
Replies
4
Views
2K
Replies
5
Views
3K
  • · Replies 3 ·
Replies
3
Views
3K
  • · Replies 19 ·
Replies
19
Views
5K
  • · Replies 38 ·
2
Replies
38
Views
3K
  • · Replies 4 ·
Replies
4
Views
3K