What distinguishes the average and most likely position of a particle?

Click For Summary

Homework Help Overview

The discussion revolves around distinguishing between the average position and the most likely position of a particle, particularly in the context of probability distributions and wave functions in quantum mechanics.

Discussion Character

  • Conceptual clarification, Mathematical reasoning, Assumption checking

Approaches and Questions Raised

  • Participants explore the definitions and implications of average position versus most likely position, referencing probability distributions and specific examples such as the Maxwell velocity distribution and a biased die. There are attempts to clarify the mathematical expressions for these concepts and how they relate to the given probability density function.

Discussion Status

The discussion is ongoing, with participants providing insights into the mathematical relationships and potential errors in calculations. Some guidance has been offered regarding integration techniques and the proper formulation of the probability function, but no consensus has been reached on the correct approach or solution.

Contextual Notes

Participants note the importance of the specific probability density function and the constraints of the problem, including the limits of integration and the nature of the wave function provided. There is an acknowledgment of possible misunderstandings regarding the relationship between the wave function and the probability function.

TFM
Messages
1,016
Reaction score
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
 
Physics news on Phys.org


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

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:


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
 


you have still not said anything about what course it is, and what probability density function you have. Also, show attempt to solution ...
 


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]
 


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><br /> = \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:


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]
 


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.
 


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


if [tex] <br /> P(x) = B^2 \sqrt{x} e^{-\beta x} <br /> [/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


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


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


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


f 'prime' means derivative of f

that is standard notation

are you sure you have done integration by parts before?
 
Last edited:
  • #15


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


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


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


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


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


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


Do we not need to use integration by parts then?
 
  • #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


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


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


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...
 
  • #26


malawi_glenn said:
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

The integral goes from 0 to infinity correct?...If so a simple substitution [itex]u=\beta\sqrt{x}[/itex] gives

[tex]\int_0^{\infty} x^{-1/2}e^{-\beta x} dx=\frac{2}{\beta} \int_0^{\infty} e^{-u^2} du[/tex]

For which the solution is well known and easily derivable.
 
  • #27


but oh my godness... wave function is NOT probability function, look at post #10.

why didn't you confirm this to me?
 
  • #28


gabbagabbahey said:
The integral goes from 0 to infinity correct?...If so a simple substitution [itex]u=\beta\sqrt{x}[/itex] gives

[tex]\int_0^{\infty} x^{-1/2}e^{-\beta x} dx=\frac{2}{\beta} \int_0^{\infty} e^{-\beta x} du[/tex]

For which the solution is well known and easily derivable.

yeah, that is true ;-)

but you still have x, but you integrate over u, that is strange
 
  • #29


malawi_glenn said:
yeah, that is true ;-)

but you still have x, but you integrate over u, that is strange

Take a look at the edited version of my post :wink:
 
  • #30


gabbagabbahey said:
Take a look at the edited version of my post :wink:

I saw, just want to be a bit rude since I did't remember that trick!

Anyway, OP do not need this integral anymore since we found a mistake on the way
 

Similar threads

Replies
14
Views
3K
  • · Replies 12 ·
Replies
12
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K
  • · Replies 1 ·
Replies
1
Views
4K
Replies
6
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 8 ·
Replies
8
Views
3K
Replies
1
Views
2K
  • · Replies 2 ·
Replies
2
Views
2K
  • · Replies 1 ·
Replies
1
Views
2K