EasySpin forum

Dear all,

I would like to know how to include the nuclear Zeeman interaction into the Spin Hamiltonian.

Sys = struct('S',1/2,'g',2,'lw',0.001);
Sys = nucspinadd(Sys,'9Be',[1 2 3])

Would this already include the nuclear Zeeman term?
If yes, how could this be shown?
If no, what would have to be done?

The nuclear Zeeman interaction is automatically included whenever you add a magnetic nucleus, just like the electron Zeeman interaction is automatically included if you add an electron spin.

Thank you. I was only confused because the nuclear Zeeman (and quadrupole) interaction are neglected for a purely perturbative approach (e.g. Opt.Method = 'perturb2';).

I currently have to deal with a system of one electron spin and several quadrupolar nuclei.
Is there any possibility except Opt.Method = 'matrix' to include the nuclear quadrupole interaction?

Please post a self-contained script that shows what you want to do. Maybe start a new thread, since this different from your original question.

Understood. I consider this thread closed.