Hi, I am using Quantica for finding the band structure of twisted bilayer graphene. Here is the part of my code that finds the hamiltonian and plots the band structure:
b = hamiltonian(tbg, model) |>
bands(subdiv((0, 2, 3, 4), 20);
mapping=(0, 2, 3, 4) => (:Γ, :K, :M, :Γ),
solver=ES.ShiftInvert(ES.ArnoldiMethod(nev=10), 0.0))
qplot(b, axis = (; xticks = ([0, 2, 3, 4], ["Γ", "K", "M", "Γ"]), xlabel="k-space",ylabel = "ϵ", title="θ = $angle ° (n=$n)"))
I want to be able to extract the raw data from the plot.
Thank you.