From 49fcd1315a05a54d34d1f229dfd0c65c3ba1fc56 Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Thu, 7 May 2026 14:02:14 -0400
Subject: [PATCH] MT-V benchmark (incorrect)

---
 calculations/MT-V.jl | 33 +++++++++++++++++++++++++++++++++
 1 file changed, 33 insertions(+)
 create mode 100644 calculations/MT-V.jl

diff --git a/calculations/MT-V.jl b/calculations/MT-V.jl
new file mode 100644
index 0000000..13a130b
--- /dev/null
+++ b/calculations/MT-V.jl
@@ -0,0 +1,33 @@
+# Malfliet–Tjon V (MT-V) — the standard few-body benchmark
+# Reference 2-body bound state: E₂ ≈ −0.4107 MeV
+
+include("../Hamiltonian.jl")
+mode = cpu_tensor
+T = Float32 # single-precision mode
+
+V_R = +1438.72 # MeV·fm   
+μ_R = 3.11 # fm⁻¹
+V_A = −626.885 # MeV·fm 
+μ_A = 1.55 # fm⁻¹
+
+r2_reg = 0.000001 # fm², regularization parameter to avoid singularity at r=0
+V_MTV(r2) = V_R * exp(−μ_R * sqrt(r2))/sqrt(r2 + r2_reg)  +  V_A * exp(−μ_A * sqrt(r2))/sqrt(r2 + r2_reg)
+
+# ℏ²/m_N = 41.47 MeV·fm²   (m_N c² = 938.92 MeV)
+μ = (1 / 41.47) / 2
+
+ϕ = 0
+d = 3
+n = 2
+N = 80
+L = 20
+n_imag = 1
+
+s = system{T}(d, n, N, L, μ)
+H = Hamiltonian{T}(s, V_MTV, ϕ, n_imag, mode)
+@time evals, _, info = eig(H, 10, resonances = false)
+
+sort!(evals)
+display(evals[1:3])
+
+# Result incorrect and does not converge, most probably due to singularity at r=0
\ No newline at end of file