yapanuwan
/
BergEC-jl
1
0
Fork 0
Code Issues Pull Requests Packages Projects Releases Wiki Activity

New 3 body system

This commit is contained in:
Nuwan Yapa 2024-03-12 20:00:36 -04:00
parent 073972d0eb
commit ead7380919
1 changed files with 3 additions and 5 deletions
Show Stats Download Patch File Download Diff File Expand all files Collapse all files

8
ho_basis_3body.jl
View File

@ -7,10 +7,8 @@ E_max = 20
Λ = 0 Λ = 0
m = 1.0 m = 1.0
Va = -5 Va = -2
Ra = sqrt(3) Ra = 2
Vb = 2
Rb = sqrt(10)
μ1 = m * 1/2 μ1 = m * 1/2
μ2 = m * 2/3 μ2 = m * 2/3
@ -27,7 +25,7 @@ println("Constructing KE matrices")
@time "T2" T2 = sp_T_matrix(n2s, l2s; ω=ω, μ=c2^2 * μ2) @time "T2" T2 = sp_T_matrix(n2s, l2s; ω=ω, μ=c2^2 * μ2)
println("Constructing PE matrices") println("Constructing PE matrices")
V_l(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω) + Vb * V_Gaussian(Rb, l, n1, n2; ω=ω) V_l(l, n1, n2) = Va * V_Gaussian(Ra, l, n1, n2; ω=ω)
@time "V1" V1 = sp_V_matrix(V_l, n1s, l1s) @time "V1" V1 = sp_V_matrix(V_l, n1s, l1s)
@time "V relative" V_relative = V1 + sp_V_matrix(V_l, n2s, l2s) @time "V relative" V_relative = V1 + sp_V_matrix(V_l, n2s, l2s)
@time "Moshinsky brackets" U = Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ) @time "Moshinsky brackets" U = Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ)