Drexel dragonThe Math ForumDonate to the Math Forum



Search All of the Math Forum:

Views expressed in these public forums are not endorsed by Drexel University or The Math Forum.


Math Forum » Discussions » Software » comp.soft-sys.math.mathematica

Topic: Exact Solution to the double well potential problem in mathematica
Replies: 1   Last Post: Mar 19, 2014 4:24 AM

Advanced Search

Back to Topic List Back to Topic List Jump to Tree View Jump to Tree View  
Bob Hanlon

Posts: 906
Registered: 10/29/11
Re: Exact Solution to the double well potential problem in mathematica
Posted: Mar 19, 2014 4:24 AM
  Click to see the message monospaced in plain text Plain Text   Click to reply to this topic Reply


hbar = 105*10^-36;
m = 96*911*10^-36;
a = 5*10^-9;
b = 5*10^-10;
V = 3*16*10^-21;
k = Sqrt[(2*m*e)/((hbar)^2)];
l = Sqrt[(2*m*(V - e))/(hbar)^2];


f[e_] = k*Cot[(-k)*a] - l*Tanh[(-l)*b/2] // Simplify;


Plot[f[e], {e, 0, 10^-18}]


FindRoot[f[e], {e, #}, WorkingPrecision -> 25] & /@
Table[est, {est, {6, 15, 35, 57, 85}*10^-18}]


{{e -> 5.795052681985343157761426*10^-18}, {e ->
1.522389834765808661449499*10^-17}, {e ->
3.474550534236380053861731*10^-17}, {e ->
5.707872905176405519130226*10^-17}, {e ->
8.492682781188889164140891*10^-17}}



Bob Hanlon




On Sat, Mar 15, 2014 at 3:47 AM, Ragavendran Nagarajan <
ragavendran.nagarajan@gmail.com> wrote:

> A little background. This is a Schrodinger equation problem to find the
> energy stages e in the double well where there is a step potential inside a
> infinite rectangular well. After solving the boundary conditions
> analytically, I tried to solve the following program in mathematica
>
> hbar := 1.05*10^-34;
> m := 0.096*9.11*10^-31;
> a := 5*10^-9;
> b := 0.5*10^-9;
> V := 0.3*1.6*10^-19;
> k = Sqrt[(2*m*e)/((hbar)^2)];
> l = Sqrt[(2*m*(V - e))/(hbar)^2];
> Block[{e}, e /. First@Solve[k*Cot[(-k)*a] - l*Tanh[(-l)*b/2] == 0, e]]
>
> As you can see since the values of the constants are very small, I am
> getting a lot of errors and I am not sure how to get the solution.
>
> Please let me know if there is a way to solve this equation using
> Mathematica. I have to submit this assignment in a couple of days, so I
> wold be really really grateful if you guys can show me the way to solve
> this.
>
> Thanks for your time
>
> Raga
>
>






Point your RSS reader here for a feed of the latest messages in this topic.

[Privacy Policy] [Terms of Use]

© Drexel University 1994-2014. All Rights Reserved.
The Math Forum is a research and educational enterprise of the Drexel University School of Education.