Difference between the average position and the most likely position of a particle

  • Thread starter TFM
  • Start date
  • #1
TFM
1,026
0

Homework Statement



part (a):

What is the average position of a particle

part (b)

What is the most likely position of the particle

Homework Equations



N/A

The Attempt at a Solution



What is the difference between these two questions, because on first glance, I would have thought they were the same?

TFM
 

Answers and Replies

  • #2
malawi_glenn
Science Advisor
Homework Helper
4,786
22


They are not a priori the same, consider as an analogy the maxwell veolcity distribution:

http://www.ac.wwu.edu/~vawter/PhysicsNet/Topics/Thermal/gifs/KinTheoryGas05.gif [Broken]

The peak is most probable-speed and the location on the x-axis where you have equal amount of probability to the left as to the right of that point is the mean-speed.

This is general for all probability distributions, e.g. wave functions in quantum mechanics. The most probable value and mean value only coincide if the distribution is symmetric around its maximum.

As an example, consider my Cheat - Dice:

It has 6sides, 4 of them has the value 1, and the rest of them has the value 6. What is the most probable outcome? And what is the mean-value of my dice? (a normal dice has 3.5)
 
Last edited by a moderator:
  • #3
TFM
1,026
0


Okay, that makes sense, since the two questions are:

(a)

Calculate the average position of the particle on the x axis, as a function of β. Mark the average position in your sketch.

(b)

Mark the most likely position of the particle in your sketch and calculate it as a function of β.

So for this, you will have to work out (a), using [tex] \int^\infty_\infty x P(x) dx [/tex], and for (b), locate its position on the graph. But how would you work out it's function?

TFM
 
  • #4
malawi_glenn
Science Advisor
Homework Helper
4,786
22


you have still not said anything about what course it is, and what probability density function you have. Also, show attempt to solution ....
 
  • #5
TFM
1,026
0


Ops, Sorry,

Okay so, I have the P(x) to be [tex] B^2 x e^{-\beta x} [/tex]

Where B is a function of beta

We also know from the wave function that the integrals are actually between 0 and infinity, because the wave function is 0, thus:

[tex] <x> = \int^{\infty}_{0} B^2 x e^{-\beta x} [/tex]

Using integration by parts gives me:

[tex] <x> = [-\beta x e^{-\beta x} - \beta^2 e^{-\beta x}]^{\infty}_0 [/tex]
 
  • #6
malawi_glenn
Science Advisor
Homework Helper
4,786
22


ok, now you are not quite correct, if the probability function [itex]
B^2 x e^{-\beta x}
[/itex] is the wave function modulus square, then this is the mean value:

[tex] <x>
= \int^{\infty}_{0} x P(x) \, dx
[/tex]

So please make sure you have done this, you have written something else in post #3 ...

Then most probable, it is just to find where maximum of P(x) is located.
 
Last edited:
  • #7
TFM
1,026
0


Isn't it normally though:

[tex] <x> = \int^\infty_\infty x P(x) dx [/tex]

But changes sometimes due to limits impsoed on the system, eg with this specific function it goes to:

[tex] <x> = \int^{\infty}_{0} x P(x) dx [/tex]
 
  • #8
malawi_glenn
Science Advisor
Homework Helper
4,786
22


Yes, but you said that [tex]
P(x) = B^2 x e^{-\beta x}
[/tex]

and then you performed this integral:

[tex]
<x> = \int^{\infty}_{0} B^2 x e^{-\beta x} dx
[/tex]

Which is just [tex]
\int^{\infty}_{0} P(x) dx \neq <x>
[/tex]

So be careful.
 
  • #9
TFM
1,026
0


Okay fair, enough, I can also see what I have done...

P(x) actually is
[tex] B^2 \sqrt{x} e^{-\beta x} [/tex]

But the version I copied from my workings I had already inserted the extra x, making

[tex] B^2 x e^{-\beta x} [/tex]

So does my integral look right now:

[tex] <x> = [-\beta x e^{-\beta x} - \beta^2 e^{-\beta x}]^{\infty}_0 [/tex]

and how do you insert the limit of infinity?
 
  • #10
malawi_glenn
Science Advisor
Homework Helper
4,786
22


if [tex]

P(x) = B^2 \sqrt{x} e^{-\beta x}

[/tex]

And if [tex]P(x) = \Psi (x)^* \, \Psi (x) [/tex]

Then

[tex] <x> = \int _0^{\infty} x\,P(x) \, dx = B^2\int _0^{\infty} x\, \sqrt{x} e^{-\beta x} \, dx[/tex]

Now I can give you a hint: [itex]xe^{-x} \rightarrow 0[/itex] when [itex]x\rightarrow \infty [/itex]
 
  • #11
TFM
1,026
0


I feel silly, I should have know that [tex] e^{-\infty} \approx 0 [/tex]

[tex] = \int _0^{\infty} x\,P(x) \, dx = B^2\int _0^{\infty} x\, \sqrt{x} e^{-\beta x} [/tex]

Okay, so now:

[tex] B^2\int _0^{\infty} x\, x^{1/2} e^{-\beta x} = B^2\int _0^{\infty} x^{3/2} e^{-\beta x} [/tex]

Integration by parts:

[tex] u = x^{3/2}, du = (x/2)^{3/2} dv = e^{-\beta x}, v = -\beta e^{-\beta x} [/tex]

This gives:

[tex] -\beta x^{3/2} e^{-\beta x} - \int {(x/2)^{3/2} e^{-\beta x}}[/tex]

Hmm, I need to use the integration by parts again, but if I do, I will still be left with another integral that needs int by parts again... Have I done something wrong?
 
  • #12
malawi_glenn
Science Advisor
Homework Helper
4,786
22


you might want to do the integration starting all over again...

Hint: consider [tex]\int f g^{'} dx = fg - \int f^{'} g dx[/tex]
with
[tex] f(x) = x^{3/2} [/tex]
[tex] g^{'} (x) = e^{-\beta x} [/tex]

you have done the derivative of x^{3/2} wrong...
 
  • #13
TFM
1,026
0


Okay, so:

[tex] f(x) = x^{3/2}, f'(x) = \frac{2}{5}x^{5/2} [/tex]

[tex] g'(x) = e^{-\beta x}, g(x) = -\beta e^{-\beta x} [/tex]

Thus:

[tex] -x^{3/2}\beta e^{-\beta x} - \int{-\frac{2}{5}x^{5/2}\beta e^{-\beta x} } [/tex]

Okay so far?
 
  • #14
malawi_glenn
Science Advisor
Homework Helper
4,786
22


f 'prime' means derivative of f

that is standard notation

are you sure you have done integration by parts before?
 
Last edited:
  • #15
TFM
1,026
0


Yeah, I did it for A Level a few years back, although, we used u/v and du/dv. I can do most of them okay...

So:

[tex] f(x) = x^{3/2}, f'(x) = \frac{1}{2}x^{1/2} [/tex]

[tex] g'(x) = e^{-\beta x}, g(x) = -\frac{1}{\beta} e^{-\beta x} [/tex]

So:

[tex] = -x^{3/2}\frac{1}{\beta} e^{-\beta x} - \int{-\frac{1}{2}x^{1/2}\frac{1}{\beta} e^{-\beta x} } [/tex]

[tex] = -x^{3/2}\frac{1}{\beta} e^{-\beta x} - -\frac{1}{2\beta} \int{x^{1/2}e^{-\beta x} } [/tex]

Is this okay now?
 
  • #16
malawi_glenn
Science Advisor
Homework Helper
4,786
22


Nope, what is derivative of x^{3/2} ?
 
  • #17
TFM
1,026
0


Okay:

[tex] f(x) = x^{3/2}, f'(x) = \frac{3}{2}x^{1/2} [/tex]

Thus:

[tex] = -x^{3/2}\frac{1}{\beta} e^{-\beta x} + \frac{3}{2\beta} \int{x^{1/2}e^{-\beta x} } [/tex]

Better?
 
  • #18
malawi_glenn
Science Advisor
Homework Helper
4,786
22


yes, now the second integral you have to do integration by parts once more.
 
  • #19
TFM
1,026
0


Okay, so:

[tex] = -x^{3/2}\frac{1}{\beta} e^{-\beta x} + \frac{3}{2\beta} \int{x^{1/2}e^{-\beta x} } [/tex]

So:

[tex] f(x) = x^{3/2}, f'(x) = \frac{1}{2}x^{-1/2} [/tex]

[tex] g'(x) = e^{-\beta x}, g(x) = -\frac{1}{\beta}e^{-\beta} [/tex]

Giving:

[tex] -x^{3/2}\frac{1}{\beta}e^{-\beta} - \int{-\frac{1}{2}x^{-1/2}\frac{1}{\beta}e^{-\beta} } [/tex]

[tex] -x^{3/2}\frac{1}{\beta}e^{-\beta} + \frac{1}{2\beta} \int{x^{-1/2}e^{-\beta} } [/tex]

I still seem to get another integral that needs another int. by parts
 
  • #20
malawi_glenn
Science Advisor
Homework Helper
4,786
22


What are you doing???

you have to evaluate integral

[tex]\int x^{1/2}e^{-\beta x} \, dx [/tex]

Why can't you do it?
 
  • #21
TFM
1,026
0


Do we not need to use integration by parts then?
 
  • #22
malawi_glenn
Science Advisor
Homework Helper
4,786
22


your last post is non-sense.

you had

[tex]
-x^{3/2}\frac{1}{\beta} e^{-\beta x} + \frac{3}{2\beta} \int x^{1/2}e^{-\beta x} \, dx
[/tex]

now you have an integral

[tex]
\int x^{1/2}e^{-\beta x} \, dx
[/tex]

which you have to integrate by parts
 
  • #23
TFM
1,026
0


That's what you get when you copy stuff and don't bother to check properly.

So:

[tex] \int x^{1/2}e^{-\beta x} dx [/tex]

Thus:

[tex] f(x) = x^{1/2}, f(x) = \frac{1}{2}x^{-1/2} [/tex]

[tex] g'(x) = e^{-\beta x}, g(x) = -\frac{1}{\beta}e^{-\beta x} [/tex]

Giving:

[tex] = -x^{1/2}\frac{1}{\beta}e^{-\beta x} - \int {\frac{1}{2}x^{-1/2} \frac{1}{\beta}e^{-\beta x}} [/tex]

[tex] = -x^{1/2}\frac{1}{\beta}e^{-\beta x} + \frac{1}{2\beta}\int {x^{-1/2}e^{-\beta x}} [/tex]

Is this better?
 
  • #24
malawi_glenn
Science Advisor
Homework Helper
4,786
22


oh crap, you can't solve this since these are not standard integrals. Sorry, I only recognize them when there is x^1/2.. I must see it as Sqrt(x) ;-)

The answer is

[tex]
\int x^{1/2}e^{-\beta x} dx = \frac{\sqrt{\pi}\text{erf} (\sqrt{\beta x})}{2\beta ^{3/2}} - \dfrac{1}{\beta}\sqrt{x}e^{-\beta x}
[/tex]

who gave you this assignment?

http://en.wikipedia.org/wiki/Error_function
http://mathworld.wolfram.com/Erf.html
 
  • #25
TFM
1,026
0


I am wondering if I made a mistake with the probability function, We was given the wave function:

[tex] \psi = B \sqrt{x}e^{-\beta x} for x \geq 0 [/tex]

Now I looked through my calculations, and I had a graph from this, and I assumed that this was also the probabiltiy function, so I assume I may have made an error somewhere...
 

Related Threads on Difference between the average position and the most likely position of a particle

Replies
6
Views
2K
Replies
9
Views
4K
Replies
1
Views
2K
  • Last Post
Replies
0
Views
941
Replies
1
Views
1K
Replies
3
Views
642
Replies
15
Views
3K
Top