hi,
to your first question:
I'm not particularly sure about the k-point sampling, kindly take a look and direct me
the k-sampling is indeed odd. you probably want a k-grid of (nkx, nky, nkz) = (5, 1, 5) as opposed to (5, 5, 1) as in your script.
apply the following change in mos2structure.py:
numerical_accuracy_parameters = NumericalAccuracyParameters(
k_point_sampling=(5, 1, 5),
)
or use the gui.
With that change applied, i can run it in roughly 20 minutes on a single core.
hope that helps you to get started.
best regards, carsten