diff --git a/ho_basis_3body.jl b/ho_basis_3body.jl index 0eb3d4e..38c84c5 100644 --- a/ho_basis_3body.jl +++ b/ho_basis_3body.jl @@ -7,10 +7,8 @@ E_max = 20 Λ = 0 m = 1.0 -Va = -5 -Ra = sqrt(3) -Vb = 2 -Rb = sqrt(10) +Va = -2 +Ra = 2 μ1 = m * 1/2 μ2 = m * 2/3 @@ -27,7 +25,7 @@ println("Constructing KE matrices") @time "T2" T2 = sp_T_matrix(n2s, l2s; ω=ω, μ=c2^2 * μ2) 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 "V relative" V_relative = V1 + sp_V_matrix(V_l, n2s, l2s) @time "Moshinsky brackets" U = Moshinsky_transform(Es, n1s, l1s, n2s, l2s, Λ)