Commit 98cc64b4 authored by stah32's avatar stah32
Browse files

Changed all outer spline radii to be <1.8 to (hopefully) remove discrepancies...

Changed all outer spline radii to be <1.8 to (hopefully) remove discrepancies between spline and non-splined tables
parent 34e11687
......@@ -9,7 +9,7 @@ bks_buck = potentialforms.buck(5131.179, 0.3040, 248.0)
coul = potentialforms.coul(-2,-2)
buck_coul = atsim.potentials.plus(bks_buck,coul)
buck_coul_zbl = atsim.potentials.plus(buck_coul,zbl)
spline = atsim.potentials.SplinePotential(zbl,buck_coul,0.5,2.0)
spline = atsim.potentials.SplinePotential(zbl,buck_coul,0.5,1.7)
#Buckingham
atsim.potentials.plot( 'buck.dat', 0.1,16,bks_buck,5000)
......
This diff is collapsed.
This diff is collapsed.
This diff is collapsed.
......@@ -7,7 +7,7 @@ bks_buck = potentialforms.buck(698.547, 0.3224, 0)
coul = potentialforms.coul(2,-2)
buck_coul = atsim.potentials.plus(bks_buck,coul)
buck_coul_zbl = atsim.potentials.plus(buck_coul,zbl)
spline = atsim.potentials.SplinePotential(zbl,buck_coul,0.9,1.9)
spline = atsim.potentials.SplinePotential(zbl,buck_coul,0.9,1.7)
atsim.potentials.plot( 'buck_coul.dat', 0.1, 16.0, buck_coul, 5000)
atsim.potentials.plot( 'zbl.dat', 0.1, 16.0, zbl, 5000)
......
This diff is collapsed.
This diff is collapsed.
......@@ -7,7 +7,7 @@ bks_buck = potentialforms.buck(5134.176, 0.3140, 283.8)
coul = potentialforms.coul(2,2)
buck_coul = atsim.potentials.plus(bks_buck,coul)
buck_coul_zbl = atsim.potentials.plus(buck_coul,zbl)
spline = atsim.potentials.SplinePotential(zbl,buck_coul,0.3,2.0)
spline = atsim.potentials.SplinePotential(zbl,buck_coul,0.3,1.7)
atsim.potentials.plot( 'buck.dat', 0.1,16,bks_buck,5000)
atsim.potentials.plot( 'buck_coul.dat', 0.1, 16.0, buck_coul, 5000)
......
This diff is collapsed.
Markdown is supported
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment