From e42ee2fae6861d505e46eae559c1d08e8654bc74 Mon Sep 17 00:00:00 2001
From: Nuwan Yapa <yapanuwan@gmail.com>
Date: Sat, 26 Sep 2026 00:22:10 -0400
Subject: [PATCH] Nucleon spectrum plot

---
 .claude/worktrees/naughty-bassi | 1 +
 test/Pb208_nucleon_spectrum.jl  | 9 +++++++--
 2 files changed, 8 insertions(+), 2 deletions(-)
 create mode 160000 .claude/worktrees/naughty-bassi

diff --git a/.claude/worktrees/naughty-bassi b/.claude/worktrees/naughty-bassi
new file mode 160000
index 0000000..4211227
--- /dev/null
+++ b/.claude/worktrees/naughty-bassi
@@ -0,0 +1 @@
+Subproject commit 4211227d6a6e437accd33a2c1ebccebc0a8b8db2
diff --git a/test/Pb208_nucleon_spectrum.jl b/test/Pb208_nucleon_spectrum.jl
index eded7cc..1e185f4 100644
--- a/test/Pb208_nucleon_spectrum.jl
+++ b/test/Pb208_nucleon_spectrum.jl
@@ -40,6 +40,11 @@ annotate!(κs .+ 0.3, Ebinds .+ 1, [(o, 9) for o in spec.occ])
 scatter!(bench_κs, bench_Ebinds, label="benchmark"; markershape=:cross, markercolor="red")
 annotate!(bench_κs .+ 0.3, bench_Ebinds .- 1, [(o, 9, :red) for o in bench_occ])
 
-xlabel!("κ")
+# label each κ in spectroscopic notation l_j (GR renders $...$ as LaTeX)
+spectroscopic(κ::Int) = "\$\\mathrm{" * "spdfghiklmno"[l_κ(κ) + 1] * "}_{" * string(Int(2 * j_κ(κ))) * "/2}\$"
+tick_κs = sort(unique(Int.([κs; bench_κs])))
+xticks!(tick_κs, spectroscopic.(tick_κs))
+xlabel!("Orbital")
 ylabel!("Binding energy (MeV)")
-yflip!()
\ No newline at end of file
+yflip!()
+savefig("temp/Pb208_nucleon_spectrum.pdf")