Register to reply

Brownian Particle bound by a Spring / internal Energy

Share this thread:
Abigale
#1
Mar25-13, 11:13 AM
P: 44
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):

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

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

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


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!
Phys.Org News Partner Physics news on Phys.org
Vibrational motion of a single molecule measured in real time
Researchers demonstrate ultra low-field nuclear magnetic resonance using Earth's magnetic field
Bubbling down: Discovery suggests surprising uses for common bubbles
klawlor419
#2
Mar25-13, 11:54 AM
P: 117
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
#3
Mar25-13, 04:16 PM
P: 44
Hey thx,
i continued calculating the Expected Value for the internal Energy of this System.

Might this way be okay?
I am not sure.

[itex]
\begin{equation}
\bar{U}=\int^{+\infty}_{-\infty}~U(x)P_{eq}dx
\end{equation}\\
~~~=(\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
[/itex]

Than i used a trick:

[itex]
\lbrace
I=\int dx \exp{-\alpha x^2}\\
~~\Rightarrow
\frac{dI}{d\alpha}=\frac{d}{d\alpha}\int dx \exp{-\alpha x^2}\\
~~~~~~~~~~~~=\int dx \frac{d}{d\alpha} \exp{-\alpha x^2}\\
~~~~~~~~~~~~=\int dx~-x^2 \exp{-\alpha x^2}
\rbrace
[/itex]


So it ensues:
[itex]
\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\\
~~~~=-(\frac{\beta}{2\pi})^{1/2}\int^{+\infty}_{-\infty}~\frac{d}{d\beta}
\exp[{-\frac{\beta}{2}(x-x_{0})^2}]~~dx \\

~~~~=-(\frac{\beta}{2\pi})^{1/2}\frac{d}{d\beta}\int^{+\infty}_{-\infty}~
\exp[{-\frac{\beta}{2}(x-x_{0})^2}]~~dx\\
\\
\\
[/itex]

And with the Gaussian Function:

[itex]
\\
\bar{U}=-(\frac{\beta}{2\pi})^{1/2}\frac{d}{d\beta}\sqrt{\frac{2\pi}{\beta}}\\
~~~=\frac{1}{2}\frac{1}{\beta}


[/itex]



[itex]
\beta
[/itex] 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
#4
Mar25-13, 07:54 PM
P: 117
Brownian Particle bound by a Spring / internal Energy

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
#5
Mar25-13, 08:09 PM
P: 117
I cannot see the solution to the equation yet.. will let you know if I get it.
klawlor419
#6
Mar25-13, 08:44 PM
P: 117
This could be a useful resource http://www.ks.uiuc.edu/~kosztin/PHYC...Notes/chp7.pdf
Abigale
#7
Mar26-13, 01:07 PM
P: 44
Hallo physicist ,

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

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

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

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?):

[itex]
\begin{split}
\bar{U} &=\int\limits_{-\infty}^{\infty} \frac{1}{2} (x-x_{0}^{'})^2 P(x,t_{0}) ~dx\\
\end{split}
[/itex]

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

[itex]
\Rightarrow\begin{split}
\bar{U} &=\int\limits_{-\infty}^{\infty} \frac{1}{2} (x-x_{0}^{'})^2
\cdot
(\frac{\beta}{2\pi})^{1/2}\exp{[-\frac{\beta}{2}(x-x_{0})^{2}]}


~dx\\
\end{split}
[/itex]






Second way (maybe even wrong?):

I regard the "Diffusion-Equation" and include the Potential, with the changed rest-position [itex]
x_{0}^{'}
[/itex] :

[itex]
U(x)=\frac{1}{2}(x-x_{0}^{'})^{2}
[/itex]

So i get the equation:

[itex]


\partial_{t} P(x,t_{0}) = \frac{D}{2} \partial_{x}^{2} P(x,t_{0}) - \Gamma \partial_{x} [
-\frac{d}{dx}(\frac{1}{2}(x-x_{0}^{'})^{2}) \cdot P(x,t_{0})
]
[/itex]


After some partial derrivations i get:

[itex]
\begin{split}
\partial_{t} P(x,t_{0})& = \underbrace{~\sqrt{\frac{\beta}{2}} \exp{[-\frac{\beta}{2}(x-x_{0})^{2}]}}
~~~
\underbrace{
\lbrace~ ~
\frac{D}{2}(-\beta+(\beta(x-x_{0}))^{2})
+ \Gamma (1-\beta(x-x_{0})(x-x_{0}^{'}))~~
\rbrace}\\
&=~~~~~~P(x,t_{0})~~~~~~~~~~~~~~~~~~~~~~~\cdot~~~~~~~~~~~~~~~~~~~~~~~~~ ~~~~~~~~~~~~~~~f
\end{split}

[/itex]


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

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

And after that i do the integration:

[itex]
\begin{split}

\bar{U} &=\int\limits_{-\infty}^{\infty} \frac{1}{2} (x-x_{0}^{'})^2
\cdot P(x,t_{1}) dx\\
&=\int\limits_{-\infty}^{\infty} \frac{1}{2} (x-x_{0}^{'})^2
\cdot P(x,t_{0})\exp{[\delta t~ f]} ~~dx\\

\end{split}
[/itex]



Please help me
and thx 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
#8
Mar26-13, 02:30 PM
P: 117
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.
Attached Files
File Type: pdf brownianOscillator.pdf (61.6 KB, 6 views)
klawlor419
#9
Apr6-13, 12:22 AM
P: 117
Any further luck with this problem?
Abigale
#10
Apr6-13, 05:33 AM
P: 44
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 thx for your letter.
P.s.
Latex can be used by writing two times the dollar-sign:
dollardollar
/xxx
dollardollar


Register to reply

Related Discussions
Free particle -> bound particle Advanced Physics Homework 1
Emissions of a particle bound in a cavity Quantum Physics 0
Bound state of a spin 3/2 particle in a potential Advanced Physics Homework 0
Energy from Brownian motion Classical Physics 16
Particle bound by quadratic potential Advanced Physics Homework 9