M. Karakoc and I. Boztosun, IJMPE Vol. 15, No. 6 (2006) 1253–1262.
https://www.worldscientific.com/doi/abs/10.1142/S0218301306004806
In [1]:
# Python program to use AIM tools
from asymptotic import *
In [2]:
# symengine (symbolic) variables for lambda_0 and s_0
En, m, hbar, L, r, r0 = se.symbols("En, m, hbar, L, r, r0")
beta, alpha, A = se.symbols("beta, alpha, A")
In [3]:
# lambda_0 and s_0
l0 = 2*beta-2/r
s0 = -beta**2 - (2*En*m)/hbar**2 + L/r**2 + L**2/r**2 + (2*beta)/r - (2*A*se.exp(-alpha*r)*m)/(hbar**2*r)
In [4]:
## Case: A = 4, L=0
nbeta = 3
nr0 = o* 1/nbeta
# parameters of lambda_0 (pl0) and s_0 (ps0)
pl0 = {beta: nbeta}
ps0 = {beta: nbeta,
alpha: o* 2/10,
A: 4,
m: o* 1/2,
L: 0,
hbar: 1,
r0: nr0}
In [5]:
%%time
# pass lambda_0, s_0 and variable values to aim class
yukawa_A4L0 = aim(l0, s0, pl0, ps0)
yukawa_A4L0.display_parameters()
yukawa_A4L0.display_l0s0(0)
yukawa_A4L0.parameters(En, r, nr0, nmax=201, nstep=10, dprec=500, tol=1e-101)
In [6]:
%%time
# create coefficients for improved AIM
yukawa_A4L0.c0()
yukawa_A4L0.d0()
yukawa_A4L0.cndn()
In [7]:
%%time
yukawa_A4L0.get_arb_roots(showRoots='-r', printFormat="{:25.20f}")
In [8]:
%%time
## Case: A = 8, L=0
nbeta = 3
nr0 = o* 1/nbeta
# parameters of lambda_0 (pl0) and s_0 (ps0)
pl0 = {beta: nbeta}
ps0 = {beta: nbeta,
alpha: o* 2/10,
A: 8,
m: o* 1/2,
L: 0,
hbar: 1,
r0: nr0}
# pass lambda_0, s_0 and variable values to aim class
yukawa_A8L0 = aim(l0, s0, pl0, ps0)
yukawa_A8L0.display_parameters()
yukawa_A8L0.display_l0s0(2)
yukawa_A8L0.parameters(En, r, nr0, nmax=201, nstep=10, dprec=500, tol=1e-101)
# create coefficients for improved AIM
yukawa_A8L0.c0()
yukawa_A8L0.d0()
yukawa_A8L0.cndn()
In [9]:
# the solution
yukawa_A8L0.get_arb_roots(showRoots='-r', printFormat="{:25.20f}")
In [10]:
%%time
## Case: A = 8, L=1
nbeta = 3
nr0 = o* 1/nbeta
# parameters of lambda_0 (pl0) and s_0 (ps0)
pl0 = {beta: nbeta}
ps0 = {beta: nbeta,
alpha: o* 2/10,
A: o* 8,
m: o* 1/2,
L: o* 1,
hbar: o* 1,
r0: nr0}
# pass lambda_0, s_0 and variable values to aim class
yukawa_A8L1 = aim(l0, s0, pl0, ps0)
yukawa_A8L1.display_parameters()
yukawa_A8L1.display_l0s0(2)
yukawa_A8L1.parameters(En, r, nr0, nmax=201, nstep=10, dprec=500, tol=1e-101)
# create coefficients for improved AIM
yukawa_A8L1.c0()
yukawa_A8L1.d0()
yukawa_A8L1.cndn()
In [11]:
# the solution
yukawa_A8L1.get_arb_roots(showRoots='-r', printFormat="{:25.20f}")
In [12]:
%%time
## Case: A = 16, L=0
nbeta = 3
nr0 = o* 1/nbeta
# parameters of lambda_0 (pl0) and s_0 (ps0)
pl0 = {beta: nbeta}
ps0 = {beta: nbeta,
alpha: o* 2/10,
A: o* 16,
m: o* 1/2,
L: o* 0,
hbar: o* 1,
r0: nr0}
# pass lambda_0, s_0 and variable values to aim class
yukawa_A16L0 = aim(l0, s0, pl0, ps0)
yukawa_A16L0.display_parameters()
yukawa_A16L0.display_l0s0(2)
yukawa_A16L0.parameters(En, r, nr0, nmax=201, nstep=10, dprec=500, tol=1e-101)
# create coefficients for improved AIM
yukawa_A16L0.c0()
yukawa_A16L0.d0()
yukawa_A16L0.cndn()
In [13]:
# the solution
yukawa_A16L0.get_arb_roots(showRoots='-r', printFormat="{:25.20f}")
In [14]:
%%time
## Case: A = 16, L=1
nbeta = 3
nr0 = o* 1/nbeta
# parameters of lambda_0 (pl0) and s_0 (ps0)
pl0 = {beta: nbeta}
ps0 = {beta: nbeta,
alpha: o* 2/10,
A: o* 16,
m: o* 1/2,
L: o* 1,
hbar: o* 1,
r0: nr0}
# pass lambda_0, s_0 and variable values to aim class
yukawa_A16L1 = aim(l0, s0, pl0, ps0)
yukawa_A16L1.display_parameters()
yukawa_A16L1.display_l0s0(2)
yukawa_A16L1.parameters(En, r, nr0, nmax=201, nstep=10, dprec=500, tol=1e-101)
# create coefficients for improved AIM
yukawa_A16L1.c0()
yukawa_A16L1.d0()
yukawa_A16L1.cndn()
In [15]:
# the solution
yukawa_A16L1.get_arb_roots(showRoots='-r', printFormat="{:25.20f}")
In [16]:
%%time
## Case: A = 24, L=0
nbeta = 3
nr0 = o* 1/nbeta
# parameters of lambda_0 (pl0) and s_0 (ps0)
pl0 = {beta: nbeta}
ps0 = {beta: nbeta,
alpha: o* 2/10,
A: o* 24,
m: o* 1/2,
L: o* 0,
hbar: o* 1,
r0: nr0}
# pass lambda_0, s_0 and variable values to aim class
yukawa_A24L0 = aim(l0, s0, pl0, ps0)
yukawa_A24L0.display_parameters()
yukawa_A24L0.display_l0s0(2)
yukawa_A24L0.parameters(En, r, nr0, nmax=201, nstep=10, dprec=500, tol=1e-101)
# create coefficients for improved AIM
yukawa_A24L0.c0()
yukawa_A24L0.d0()
yukawa_A24L0.cndn()
In [17]:
%%time
# the solution
yukawa_A24L0.get_arb_roots(showRoots='-r', printFormat="{:25.20f}")
In [18]:
%%time
## Case: A = 24, L=1
nbeta = 3
nr0 = o* 1/nbeta
# parameters of lambda_0 (pl0) and s_0 (ps0)
pl0 = {beta: nbeta}
ps0 = {beta: nbeta,
alpha: o* 2/10,
A: o* 24,
m: o* 1/2,
L: o* 1,
hbar: o* 1,
r0: nr0}
# pass lambda_0, s_0 and variable values to aim class
yukawa_A24L1 = aim(l0, s0, pl0, ps0)
yukawa_A24L1.display_parameters()
yukawa_A24L1.display_l0s0(2)
yukawa_A24L1.parameters(En, r, nr0, nmax=201, nstep=10, dprec=500, tol=1e-101)
# create coefficients for improved AIM
yukawa_A24L1.c0()
yukawa_A24L1.d0()
yukawa_A24L1.cndn()
In [19]:
%%time
# the solution
yukawa_A24L1.get_arb_roots(showRoots='-r', printFormat="{:25.20f}")
In [20]:
%%time
## Case: A = 24, L=2
nbeta = 3
nr0 = o* 1/nbeta
# parameters of lambda_0 (pl0) and s_0 (ps0)
pl0 = {beta: nbeta}
ps0 = {beta: nbeta,
alpha: o* 2/10,
A: o* 24,
m: o* 1/2,
L: o* 2,
hbar: o* 1,
r0: nr0}
# pass lambda_0, s_0 and variable values to aim class
yukawa_A24L2 = aim(l0, s0, pl0, ps0)
yukawa_A24L2.display_parameters()
yukawa_A24L2.display_l0s0(2)
yukawa_A24L2.parameters(En, r, nr0, nmax=201, nstep=10, dprec=500, tol=1e-101)
# create coefficients for improved AIM
yukawa_A24L2.c0()
yukawa_A24L2.d0()
yukawa_A24L2.cndn()
In [21]:
%%time
# the solution
yukawa_A24L2.get_arb_roots(showRoots='-r', printFormat="{:25.20f}")
In [ ]: