How to define a spin system of lanthanide ions in easyspin?

General forum for EasySpin: questions, how to's, etc.
Post Reply
lnxj99
Newbie
Posts: 1
Joined: Thu Jun 07, 2018 8:51 am

How to define a spin system of lanthanide ions in easyspin?

Post by lnxj99 »

Hello, I am a new user of easyspin. I met problems in treating lanthanide ions using easypin.
As the document "EPR Reference> Orbital angular momentum and spin-orbit interaction" said to treat lanthanide ions, the soc Hamiltonian usually include higher orders terms (lambda n*(sigma*L*S)^n).
And in the user guide "The spin system", in the soc part an example was given as below:
---
% three electron spins, with linear soc parameter of 170, 250 cm-1 and 0:
Sys.soc = [170; 250; 0]*clight*1e-4;
% as above, plus the first spin has a cubic soc parameter of 10 cm-1:
Sys.soc = [170 0, 10; 250, 0, 0; 0, 0, 0]*clight*1e-4;
------
Is there a mistake that a comma between 170 and 0 is lost? Hence the right format is Sys.soc = [1st order, 2nd order, 3rd order, 4th order, ....], corresponding to lambda n?

If I use this definition with lambda n, the soc of excited states will be also considered and the system will be very large. Is there any method/approximation to accelerate the calculations?
Stefan Stoll
EasySpin Creator
Posts: 1073
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

Re: How to define a spin system of lanthanide ions in easysp

Post by Stefan Stoll »

You can add a comma, but Matlab treats [1 2 3] and [1, 2, 3] identically.

Regarding performance, it is not possible to speed this up if spin-orbit coupling is present.
Post Reply