Brownian Particle bound by a Spring / internal Energy

[Abigale]

Hi,
i regard a Brownian Particle connectet to a Spring and there is a heat-reservoir.
The distribution of the x-coordinate of the particle follows the Diffusion-Equation (Fokker-Planck-Equation):

\partial_{t}P(x,t)=\frac{D}{2}<br /> \partial_{x}^{2}P(x,t)-<br /> \Gamma\partial_{x}[f(x)P(x,t)]<br />

A deterministic Force is given by f(x)=-\frac{d}{dx}U(x).
Whereby
<br /> U(x)=<br /> \frac{1}{2}(x-x_{0})^{2}<br /> is a Potential.

Also i know that the equilibrium-equation is a Gaussian-Function.
<br /> P_{eq}=(\frac{\beta}{2\pi})^{1/2}\exp[{-\frac{\beta}{2}(x-x_{0})^2}]<br />


I want to determine the Expected Value of the (potential) internal Energy.
But i don't know how I can get it. ;-(

Please help me ;) and
thank you a lot!
Bye Abigale

Sorry for my bad english!

 

[klawlor419]

For any function say f(x), to find its average value you take the product of it with the probability of finding a value x (your probability distribution in this case) and sum it up over all possible values of x which is in your case the displacement of the particle.

 

[Abigale]

Expected Value of the Brownian Particle/ continued

Hey thx,
i continued calculating the Expected Value for the internal Energy of this System.

Might this way be okay?
I am not sure.

<br /> \begin{equation}<br /> \bar{U}=\int^{+\infty}_{-\infty}~U(x)P_{eq}dx<br /> \end{equation}\\<br /> ~~~=(\frac{\beta}{2\pi})^{1/2}\int^{+\infty}_{-\infty}~\frac{1}{2}(x-x_{0})^2 \exp[{-\frac{\beta}{2}(x-x_{0})^2}]~~dx<br />

Than i used a trick:

<br /> \lbrace<br /> I=\int dx \exp{-\alpha x^2}\\<br /> ~~\Rightarrow <br /> \frac{dI}{d\alpha}=\frac{d}{d\alpha}\int dx \exp{-\alpha x^2}\\<br /> ~~~~~~~~~~~~=\int dx \frac{d}{d\alpha} \exp{-\alpha x^2}\\<br /> ~~~~~~~~~~~~=\int dx~-x^2 \exp{-\alpha x^2}<br /> \rbrace<br />


So it ensues:
<br /> \bar{U}=-(\frac{\beta}{2\pi})^{1/2}\int^{+\infty}_{-\infty}~-\frac{1}{2}(x-x_{0})^2 \exp[{-\frac{\beta}{2}(x-x_{0})^2}]~~dx\\<br /> ~~~~=-(\frac{\beta}{2\pi})^{1/2}\int^{+\infty}_{-\infty}~\frac{d}{d\beta}<br /> \exp[{-\frac{\beta}{2}(x-x_{0})^2}]~~dx \\<br /> <br /> ~~~~=-(\frac{\beta}{2\pi})^{1/2}\frac{d}{d\beta}\int^{+\infty}_{-\infty}~<br /> \exp[{-\frac{\beta}{2}(x-x_{0})^2}]~~dx\\<br /> \\<br /> \\<br />

And with the Gaussian Function:

<br /> \\<br /> \bar{U}=-(\frac{\beta}{2\pi})^{1/2}\frac{d}{d\beta}\sqrt{\frac{2\pi}{\beta}}\\<br /> ~~~=\frac{1}{2}\frac{1}{\beta}<br /> <br /> <br />



<br /> \beta <br /> is the inverse temperature.
So the Expected Value of the internal Energy should be proportional to the Temperature, this seems good.
But i am not sure if the result/callculations i have done are correct.
Could someone help me?

Thanks a lot.
Bye Abigale

 

[klawlor419]

That looks good and makes sense. In fact, its easily confirmed by the equipartition theorem which states that any quadratic degree of freedom will contribute one-half kT to the average energy.

The next thing I think you need to consider is that the probability distribution may get slightly altered by the oscillator potential. I recognize your P(eq) as the equilibrium probability distribution for purely an oscillator potential. If you notice though it doesn't satisfy the diffusion equation with time present.

Point being that there is a more general solution to this adjusted Fokker-Planck equation, for which you can extract some time-dependence on the average potential energy. Of course for the large-t limit it will go to the equilibrium distribution, but the behavior in the mean time is interesting

 

[klawlor419]

I cannot see the solution to the equation yet.. will let you know if I get it.

 

[klawlor419]

This could be a useful resource http://www.ks.uiuc.edu/~kosztin/PHYCS498NSM/LectureNotes/chp7.pdf

 

[Abigale]

Hallo physicist ,

Let us consider a change of the rest-position of the spring from x_0 to x_{o}^{&#039;}.
This happened within the short time-intervall [t_0,t_1].

The change of the rest-position is very fast.
So the assumption can be done, that the particle-coordinate x and the distribution P(x,t) is not changed.

Now i want to callculate the Expexted Value of the Internal-Energy immediately after the change, so \bar U (t_1).

I am thinking for many hours and got maybe an idea, but i am not sure if it is right.
Could you help me to find what is wrong, or do you have a better idea?





First way, but seems to easy (What is wrong?):

<br /> \begin{split}<br /> \bar{U} &amp;=\int\limits_{-\infty}^{\infty} \frac{1}{2} (x-x_{0}^{&#039;})^2 P(x,t_{0}) ~dx\\<br /> \end{split}<br />

So here i think P(x,t_{0}) must then be the P_{eq} from my first post.

<br /> \Rightarrow\begin{split}<br /> \bar{U} &amp;=\int\limits_{-\infty}^{\infty} \frac{1}{2} (x-x_{0}^{&#039;})^2 <br /> \cdot<br /> (\frac{\beta}{2\pi})^{1/2}\exp{[-\frac{\beta}{2}(x-x_{0})^{2}]}<br /> <br /> <br /> ~dx\\<br /> \end{split}<br />






Second way (maybe even wrong?):

I regard the "Diffusion-Equation" and include the Potential, with the changed rest-position <br /> x_{0}^{&#039;}<br /> :

<br /> U(x)=\frac{1}{2}(x-x_{0}^{&#039;})^{2}<br />

So i get the equation:

<br /> <br /> <br /> \partial_{t} P(x,t_{0}) = \frac{D}{2} \partial_{x}^{2} P(x,t_{0}) - \Gamma \partial_{x} [<br /> -\frac{d}{dx}(\frac{1}{2}(x-x_{0}^{&#039;})^{2}) \cdot P(x,t_{0}) <br /> ]<br />


After some partial derrivations i get:

<br /> \begin{split}<br /> \partial_{t} P(x,t_{0})&amp; = \underbrace{~\sqrt{\frac{\beta}{2}} \exp{[-\frac{\beta}{2}(x-x_{0})^{2}]}}<br /> ~~~<br /> \underbrace{<br /> \lbrace~ ~<br /> \frac{D}{2}(-\beta+(\beta(x-x_{0}))^{2})<br /> + \Gamma (1-\beta(x-x_{0})(x-x_{0}^{&#039;}))~~<br /> \rbrace}\\<br /> &amp;=~~~~~~P(x,t_{0})~~~~~~~~~~~~~~~~~~~~~~~\cdot~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~~f<br /> \end{split}<br /> <br />


Now, i regard the small time intervall [t_{0},t_{1}]=\delta t.
So i claim that :

<br /> \begin{split}<br /> P(x,t_{1})&amp;=P(x,t_{0}+\delta t)\overset{\text{taylor}}=P(x,t_{0})+\delta t \partial_{t}P(x,t_{0})\\ <br /> &amp;=P(x,t_{0})\lbrace 1+ \delta t f \rbrace\\<br /> &amp;\overset{\exp{\delta t f \approx 1 + \delta t f}}=P(x,t_{0})\exp{[\delta t~ f]}<br /> \end{split}<br />

And after that i do the integration:

<br /> \begin{split}<br /> <br /> \bar{U} &amp;=\int\limits_{-\infty}^{\infty} \frac{1}{2} (x-x_{0}^{&#039;})^2 <br /> \cdot P(x,t_{1}) dx\\<br /> &amp;=\int\limits_{-\infty}^{\infty} \frac{1}{2} (x-x_{0}^{&#039;})^2 <br /> \cdot P(x,t_{0})\exp{[\delta t~ f]} ~~dx\\<br /> <br /> \end{split}<br />



Please help me
and thanks for the help.

Also i am interessted in the work i have to do for the change of the rest-position.
But this later guys.


Bye Abby

 

[klawlor419]

Abigale

I wrote up a LaTeX file and will attach it. I tried to address your questions and I asked a few as well. Still have not figured out how to use LaTeX on the forum.

 

[klawlor419]

Any further luck with this problem?

 

[Abigale]

Not yet, I regard an similar problem with other boundary conditions for better understanding the fpe. But I slowly understand it. I study physics and I hear the first time statistical mechanics. I will write you later more.
And thanks for your letter.
P.s.
Latex can be used by writing two times the dollar-sign:
dollardollar
/xxx
dollardollar