Re: Schroedinger EQ
- To: mathgroup at smc.vnet.net
- Subject: [mg122055] Re: Schroedinger EQ
- From: "Kevin J. McCann" <Kevin.McCann at umbc.edu>
- Date: Tue, 11 Oct 2011 04:23:21 -0400 (EDT)
- Delivered-to: l-mathgroup@mail-archive0.wolfram.com
- References: <j6q61b$p1g$1@smc.vnet.net> <j6uafl$f2c$1@smc.vnet.net>
There are a few problems. First, I would load the Notation package and then Symbolize a_so, V_so, etc. or else just use aso, Vso, ... Next, in one of your equations you have simply R instead of a subscripted one, I assume you left this off. In your function eq[V_,x_,x0_] you have ... == -e = y[x]. I am not sure what you mean by this, but the syntax is not good. Also the "e" you use should be either "E" or esc-ee-esc, i.e. 2.71828 I didn't get any further, because there are others I think. Kevin On 10/10/2011 4:30 AM, raj kumar wrote: > sorry, A =82 > > > On Oct 9, 2:50 am, raj kumar<rajesh779... at gmail.com> wrote: >> dear esteemed experts, >> >> i wonder if somebody can help me >> >> i have been trying to find a a certain value of a parameter V that >> will "match" the logarithmic derivative (of solutions to the time >> independent Schroedinger eq ) at both sides of a particular point >> called the matching point. But cannot seem to find the correct V value >> that will make >> bc1[V_] = bc2[V_].. ...mathematica keeps giving an error message .See >> below for the code. >> >> any help will be most appreciated >> >> a=0.63; >> A8; >> j=9/2; >> L=4; >> mu=(931.5 (208 1.))/(208+1.); >> Z=82; >> Subscript[a, so]=0.5; >> R=1.25 A^(1/3); >> Subscript[V, so]=7; >> Subscript[R, c]=1.25 A^(1/3); >> Subscript[R, so]=1.1 A^(1/3); >> V1[x_,V_]:=-(V/(E^((x-R)/a)+1)); V2[x_]:=-((2 ((j+1) j-L (L+1)-3/4= > ) >> Subscript[V, so] E^((x-Subscript[R, so])/Subscript[a, so]))/ >> (Subscript[a, so] (E^((x-Subscript[R, so])/Subscript[a, so])+1)^2)); >> pott[x_,V_]=V1[x,V]+V2[x]; >> >> emin=-55; >> emax= -5; >> xmax=10; >> xmin=0.1; >> xmatch=4.5; >> e=3.94; >> >> eq[V_, x_, x0_] = {-( >> \!\(\*SuperscriptBox["y", "\[Prime]\[Prime]", >> MultilineFunction->None]\)[x]/( >> 2 mu)) + (pott[x, V] + L (L + 1)/(2 mu (x^2))) y[x] == -e= > y[x], >> y[x0] == 0, >> \!\(\*SuperscriptBox["y", "\[Prime]", >> MultilineFunction->None]\)[x0] == 1/10^6}; >> y1[V_, x_] := y[x] /. NDSolve[eq[V, x, xmin], y, {x, xmin, xmatch}]; >> bc1[V_] := \!\( >> \*SubscriptBox[\(\[PartialD]\), \(x\)]\(y1[x]\)\)/y1[x] /. x -> >> xmatch; >> y2[V_, x_] := y[x] /. NDSolve[eq[V, x, xmax], y, {x, xmax, xmatch}]; >> bc2[V_] := \!\( >> \*SubscriptBox[\(\[PartialD]\), \(x\)]\(y2[x]\)\)/y2[x] /. x -> >> xmatch; >> bc[V_?NumericQ] := bc1[V] - bc2[V]; >> Vvalue = V /. >> If[emax == emin, V, >> FindRoot[bc[V], {V, emin, emax}, AccuracyGoal -> 10, >> WorkingPrecision -> 20]]; >> Print["the value of V is =" , Vvalue] > >