|
This is technically a homework thread, but I'm not asking for someone to do all my work for me and solve it completely, all I need is a nudge in the right direction because I'm fairly lost and can't find help anywhere else. If even this is not allowed, then I'm sorry 
The question:
+ Show Spoiler +Find the one dimensional particle motion in the Morse Potential
What I have so far:
+ Show Spoiler +I now need to change the variables to get something integrable, but I'm unsure what to use. For example: doesn't work because I need the du expression to be free of the variable x. I think I need to divide by something but I just can't work out what that is. Any tips or hints would be greatly appreciated, thanks.
Oh, and http://thornahawk.unitedti.org/equationeditor/equationeditor.php is what I use for the equations.
|
is awesome32277 Posts
nvm this doens't work ;_;
|
use the lagrangefunction.
L=T-U , T beeing kinetic and U potential energy.
this gives you directly the differential equation for the motion.
q beeing x in your case.
The solution of the resulting motion equation depends on your initial values ofc. depending on them you could solve at least the linear part.
You could also expand tp the 2nd order your inital potential around its minimum and approximate the tail to its right. use the solution of one as boundary condition for the other.
|
On October 31 2009 09:13 aqui wrote:use the lagrangefunction. L=T-U , T beeing kinetic and U potential energy. this gives you directly the differential equation for the motion. q beeing x in your case.
Thats pretty much what I used to get to the first step I've written down, its from there that the problem I have is.
I know how to get the EoM, its just that specific integral I'm stuck on. Thanks though. 
|
United States24752 Posts
I'm so glad I don't have to do this anymore lol
What course is this exactly and what program is it for? When I learned lagrangians in mechanics I don't think I ever was given this potential to use.
|
On October 31 2009 09:38 micronesia wrote:
I'm so glad I don't have to do this anymore lol
What course is this exactly and what program is it for? When I learned lagrangians in mechanics I don't think I ever was given this potential to use.
Its just called advanced mechanics, its part of 2nd year in my degree. (theoretical physics) This question is part of a section on integrating the equations of motion.
|
You can just use computer software to evaluate that integral, or look it up in a table of integrals. If you really have to do it by hand, I think you would have to substitute some sort of trig function.
|
On October 31 2009 09:30 deconduo wrote:Show nested quote +On October 31 2009 09:13 aqui wrote:use the lagrangefunction. L=T-U , T beeing kinetic and U potential energy. this gives you directly the differential equation for the motion. q beeing x in your case. Thats pretty much what I used to get to the first step I've written down, its from there that the problem I have is. I know how to get the EoM, its just that specific integral I'm stuck on. Thanks though.
thought you just solved t=x/v. because thats what stands there^^ the integral looks like a bitch though^^ no idea how to solve it, sry:p
|
I dont know how to solve that integral... but:
http://www93.wolframalpha.com/input/?i=integrate [1%2F%28%28k-%281-e^%28-ax%29%29%29%29^0.5%2C+x]%2C+a%3D{-infinite%2Cinfinite}%2Ck%3D{-infinite%2Cinfinite}
dont know if everything is set up correctly, so give a look before follow it blindly 
edit: sorry.. thats the best way to put the link.. copy all of it
|
|
|
If you have some idea that 'a' is very small or that 'x' will be very small, then you can use the binomial expansion or the Taylor expansion to simplify the expression to something integrable. Without any paper immediately on hand, the obvious substitution looks like u = ax. Then you can probably do something like partial fractions.
Other than that, just try the Integrator.
http://integrals.wolfram.com/index.jsp?expr=1/sqrt(1-(1-exp(-3x))^2)&random=false
|
Thanks for the integrator thing, but I have to do it by hand and show the work unfortunately. I'll see if I can work backwards from the result though.
|
Now that I look at it more closely, I don't see where you got that integral at all.
If you use the Lagrangian (or in this case, just say F = -grad(u), which I think is faster), you should get the expression
mx'' = -2Aa(exp(-ax)-exp(-2ax))
From there you should be able to separate variables and integrate twice. There is a substitution needed (for the x side) and some algebraic manipulation, then you will end up with two integrals for the x side which should solvable. Then just backtrack your substitution carefully after you've solved.
|
after the substitution
we have
the next substitution is
where
if E>A (infinite motion) then the equation in the denominator has 2 different real roots, if E<1 this equation has no real roots. These two cases should be considered separatly.
Let us consider the first case.
where
where
The second case is similar but there will be arctan. It looks messy, so maybe some simplifications are available.
|
impressive... most impressive.
|
|
|
|
|
|
EPT![[image loading]](http://www.codecogs.com/eq.latex?U(x)%20=%20A(1-e^{-ax})^2)
![[image loading]](http://upload.wikimedia.org/math/a/e/9/ae9c64c6d51a4dac2002030f436866fd.png)
![[image loading]](http://www.texify.com/img/%5CLARGE%5C%21%5Cbegin%7Bequation%7D%5Cint%20%5Cfrac%7Bdx%7D%7B%5Csqrt%7BE-A%281-e%5E%7B-ax%7D%29%5E2%7D%7D%5Cend%7Bequation%7D.gif)
![[image loading]](http://www.texify.com/img/%5CLARGE%5C%21%5Cbegin%7Bequation%7Dt%3D1-e%5E%7B-ax%7D%5Cend%7Bequation%7D.gif)
![[image loading]](http://www.texify.com/img/%5CLARGE%5C%21%5Cbegin%7Bequation%7D%5Cint%20%5Cfrac%7Bdt%7D%7B%281-t%29%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D-t%5E2%7D%7D%5Cend%7Bequation%7D.gif)
![[image loading]](http://www.texify.com/img/%5CLARGE%5C%21%5Cbegin%7Bequation%7Dt%3D%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D%7D%5Csin%20%5Cphi%5Cend%7Bequation%7D.gif)
![[image loading]](http://www.texify.com/img/%5CLARGE%5C%21%5Cbegin%7Beqnarray%7D%26%26%5Cfrac%7B1%7D%7B%5Csqrt%7BA%7Da%7D%5Cint%20%5Cfrac%7Bd%5Cphi%7D%7B1-%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D%7D%5Csin%20%5Cphi%7D%3D%5Cfrac%7B1%7D%7B%5Csqrt%7BA%7Da%7D%5Cint%20%5Cfrac%7Bd%5Cphi%7D%7B%5Ccos%5E2%5Cfrac%7B%5Cphi%7D%7B2%7D%2B%5Csin%5E2%5Cfrac%7B%5Cphi%7D%7B2%7D-2%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D%7D%5Ccos%5Cfrac%7B%5Cphi%7D%7B2%7D%5Csin%5Cfrac%7B%5Cphi%7D%7B2%7D%7D%3D%5C%5C%26%26%3D%5Cfrac%7B1%7D%7B%5Csqrt%7BA%7Da%7D%5Cint%20%5Cfrac%7B%5Cfrac%7B1%7D%7B%5Ccos%5E2%5Cfrac%7B%5Cphi%7D%7B2%7D%7Dd%5Cphi%7D%7B1%2B%5Ctan%5E2%5Cfrac%7B%5Cphi%7D%7B2%7D-2%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D%7D%5Ctan%5Cfrac%7B%5Cphi%7D%7B2%7D%7D%3D%5Cfrac%7B2%7D%7B%5Csqrt%7BA%7Da%7D%5Cint%20%5Cfrac%7Bdy%7D%7B1%2By%5E2-2%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D%7Dy%7D%5Cend%7Beqnarray%7D.gif)
![[image loading]](http://www.texify.com/img/%5CLARGE%5C%21%5Cbegin%7Bequation%7Dy%3D%5Ctan%5Cfrac%7B%5Cphi%7D%7B2%7D%5Cend%7Bequation%7D.gif)
![[image loading]](http://www.texify.com/img/%5CLARGE%5C%21%5Cbegin%7Beqnarray%7D%5Cfrac%7B2%7D%7B%5Csqrt%7BA%7Da%7D%5Cint%20%5Cfrac%7Bdy%7D%7B1%2By%5E2-2%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D%7Dy%7D%3D%5Cfrac%7B2%7D%7B%5Csqrt%7BA%7Da%7D%5Cint%20%5Cfrac%7Bdy%7D%7B%28y-y_1%29%28y-y_2%29%7D%5Cend%7Beqnarray%7D.gif)
![[image loading]](http://www.texify.com/img/%5CLARGE%5C%21%5Cbegin%7Beqnarray%7D%26%26y_1%3D%5Csqrt%5Cfrac%7BE%7D%7BA%7D%2B%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D-1%7D%5C%5C%26%26y_2%3D%5Csqrt%5Cfrac%7BE%7D%7BA%7D-%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D-1%7D%5Cend%7Beqnarray%7D.gif)
![[image loading]](http://www.texify.com/img/%5CLARGE%5C%21%5Cbegin%7Beqnarray%7D%26%26%5Cfrac%7B2%7D%7B%5Csqrt%7BA%7Da%7D%5Cint%20%5Cfrac%7Bdy%7D%7B%28y-y_1%29%28y-y_2%29%7D%3D%5Cfrac%7B1%7D%7B%5Csqrt%7BA%7Da%7D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D-1%7D%7D%20%5Cint%20%5Cleft%28%20%5Cfrac%7B1%7D%7By-y_1%7D-%5Cfrac%7B1%7D%7By-y_2%7D%20%5Cright%29dy%3D%5Cnonumber%5C%5C%26%26%3D%5Cfrac%7B1%7D%7B%5Csqrt%7BA%7Da%7D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D-1%7D%7D%20%5Cleft%28%20%5Cln%7Cy-y_1%7C-%5Cln%7Cy-y_2%7C%20%20%5Cright%29%3D%5Cnonumber%5C%5C%26%26%3D%5Cfrac%7B2%7D%7B%5Csqrt%7BA%7Da%7D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D-1%7D%7D%20%5Cleft%28%20%5Cln%7C%5Ctan%5Cfrac%7B%5Cphi%7D%7B2%7D-y_1%7C-%5Cln%7C%5Ctan%5Cfrac%7B%5Cphi%7D%7B2%7D-y_2%7C%20%20%5Cright%29%3D%5Cnonumber%5C%5C%26%26%3D%5Cfrac%7B1%7D%7B%5Csqrt%7BA%7Da%7D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D-1%7D%7D%20%5Cleft%28%20%5Cln%7C%5Cfrac%7B%5Csin%20%5Cphi%7D%7B1%2B%5Ccos%20%5Cphi%7D-y_1%7C-%5Cln%7C%5Cfrac%7B%5Csin%20%5Cphi%7D%7B1%2B%5Ccos%20%5Cphi%7D-y_2%7C%20%20%5Cright%29%3D%5Cnonumber%5C%5C%26%26%3D%5Cfrac%7B1%7D%7B%5Csqrt%7BA%7Da%7D%20%5Cfrac%7B1%7D%7B%5Csqrt%7B%5Cfrac%7BE%7D%7BA%7D-1%7D%7D%20%5Cleft%28%20%5Cln%7C%5Cfrac%7Bt%5Csqrt%7BA%7D%7D%7B%5Csqrt%7BE%7D%2B%5Csqrt%7BE-t%5E2A%7D%7D-y_1%7C-%5Cln%7C%5Cfrac%7Bt%5Csqrt%7BA%7D%7D%7B%5Csqrt%7BE%7D%2B%5Csqrt%7BE-t%5E2A%7D%7D-y_2%7C%20%20%5Cright%29%5Cnonumber%5Cend%7Beqnarray%7D.gif)
