Brownian Particle bound by a Spring / internal Energyby Abigale Tags: brownian motion, nonequilibrium, statistical mechanic, statistical physics, thermo dynamics 

#1
Mar2513, 11:13 AM

P: 44

Hi,
i regard a Brownian Particle connectet to a Spring and there is a heatreservoir. The distribution of the xcoordinate of the particle follows the DiffusionEquation (FokkerPlanckEquation): [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}(xx_{0})^{2} [/itex] is a Potential. Also i know that the equilibriumequation is a GaussianFunction. [itex] P_{eq}=(\frac{\beta}{2\pi})^{1/2}\exp[{\frac{\beta}{2}(xx_{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! 



#2
Mar2513, 11:54 AM

P: 112

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.




#3
Mar2513, 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}(xx_{0})^2 \exp[{\frac{\beta}{2}(xx_{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}(xx_{0})^2 \exp[{\frac{\beta}{2}(xx_{0})^2}]~~dx\\ ~~~~=(\frac{\beta}{2\pi})^{1/2}\int^{+\infty}_{\infty}~\frac{d}{d\beta} \exp[{\frac{\beta}{2}(xx_{0})^2}]~~dx \\ ~~~~=(\frac{\beta}{2\pi})^{1/2}\frac{d}{d\beta}\int^{+\infty}_{\infty}~ \exp[{\frac{\beta}{2}(xx_{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 



#4
Mar2513, 07:54 PM

P: 112

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 onehalf 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 FokkerPlanck equation, for which you can extract some timedependence on the average potential energy. Of course for the larget limit it will go to the equilibrium distribution, but the behavior in the mean time is interesting 



#5
Mar2513, 08:09 PM

P: 112

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




#6
Mar2513, 08:44 PM

P: 112

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




#7
Mar2613, 01:07 PM

P: 44

Hallo physicist ,
Let us consider a change of the restposition of the spring from [itex]x_0[/itex] to [itex]x_{o}^{'}[/itex]. This happened within the short timeintervall [[itex]t_0,t_1[/itex]]. The change of the restposition is very fast. So the assumption can be done, that the particlecoordinate [itex]x[/itex] and the distribution [itex]P(x,t)[/itex] is not changed. Now i want to callculate the Expexted Value of the InternalEnergy 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} (xx_{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} (xx_{0}^{'})^2 \cdot (\frac{\beta}{2\pi})^{1/2}\exp{[\frac{\beta}{2}(xx_{0})^{2}]} ~dx\\ \end{split} [/itex] Second way (maybe even wrong?): I regard the "DiffusionEquation" and include the Potential, with the changed restposition [itex] x_{0}^{'} [/itex] : [itex] U(x)=\frac{1}{2}(xx_{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}(xx_{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}(xx_{0})^{2}]}} ~~~ \underbrace{ \lbrace~ ~ \frac{D}{2}(\beta+(\beta(xx_{0}))^{2}) + \Gamma (1\beta(xx_{0})(xx_{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} (xx_{0}^{'})^2 \cdot P(x,t_{1}) dx\\ &=\int\limits_{\infty}^{\infty} \frac{1}{2} (xx_{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 restposition. But this later guys. Bye Abby 



#8
Mar2613, 02:30 PM

P: 112

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. 



#9
Apr613, 12:22 AM

P: 112

Any further luck with this problem?




#10
Apr613, 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 dollarsign: 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 