Exchange interaction for non parallel spin

General forum for EasySpin: questions, how to's, etc.
Post Reply
Qing Zhang
Newbie
Posts: 2
Joined: Wed Dec 06, 2017 6:54 am

Exchange interaction for non parallel spin

Post by Qing Zhang »

Hello,

My name is Qing Zhang from Prof. Myriam P. Sarachik's group. Could you please help me about the question of the Dy3 paper (Phys. Rev. Lett. 100, 247205):

I am interested in anisotropic exchange in molecular magnets. I am able to rotate the local tensor frame for anisotropic g and zero field splitting terms. However I do not know how to tilt tensors associated with anisotropic exchange. An important example would be the ising model applied for the interpretation single crystal of magnetization of Dy3 (Phys. Rev. Lett. 100, 247205). I make an attempt to reproduce the characteristic Zeeman splitting diagram (fig 3)

% code before exchange:
Sys.S = [1/2 1/2 1/2];
Sys.g = [0 0 1; 0 0 1; 0 0 1]; % anisotropic g tensors
Sys.gFrame = [0 0 0; 0 120 0; 0 240 0]*pi/180;
Sys.ee = [0;0;0]*(-1e5); % coupling, in MHz

% code after adding exchange:
Sys.ee = [0 0 1;0 0 1;0 0 1]*(-1e5); % coupling, in MHz
Sys.eeFrame = [0 0 0; 0 120 0; 0 240 0]*pi/180;

Could you please give me any suggestion about how to write the code about the exchange or how to calculate the Zeeman splitting diagram? Thank you very much!


Best regards,

Qing Zhang
Attachments
Zeeman energy of Dy3 molecule
Zeeman energy of Dy3 molecule
image.png (134.81 KiB) Viewed 2796 times
katarkon
Local Expert
Posts: 182
Joined: Mon Jan 12, 2015 4:01 am

Re: Exchange interaction for non parallel spin

Post by katarkon »

Use function levelsplot.
Qing Zhang
Newbie
Posts: 2
Joined: Wed Dec 06, 2017 6:54 am

Re: Exchange interaction for non parallel spin

Post by Qing Zhang »

Thank you for your comment. I am familiar with levelsplot. I should be clearer in my question; the code does not match the desired field dependence (i.e. a double degenerate singlet ground state at low field, like a) of the attached figure). This is most likely due to the required rotation of the anisotropic exchange. Could you please tell me how to show the following expression by the easyspin spin system: j_zz * S_zi * S_zk (the S_zi is equal 1/2, but S_zi is not parallel with S_zk)? Thank you.
Matt Krzyaniak
EasySpin Guru
Posts: 153
Joined: Tue Jul 22, 2014 11:01 am
Location: Northwestern University

Re: Exchange interaction for non parallel spin

Post by Matt Krzyaniak »

Here's a little snippet of code that might help you in thinking about how the rotations work and how things sum together to generate the effective interactions(you can do the same with the exchange coupling)

Code: Select all

Sys.g = [0 0 1; 0 0 1; 0 0 1];
Sys.gFrame = [0 0 0; 0 120 0; 0 240 0]*pi/180;

g = diag(Sys.g(1,:));
R_m2g = erot(Sys.gFrame(1,:)); % mol frame -> g frame
R_g2m = R_m2g.'; % g frame -> mol frame
g1 = diag(R_g2m*g*R_g2m.')

g = diag(Sys.g(2,:));
R_m2g = erot(Sys.gFrame(2,:)); % mol frame -> g frame
R_g2m = R_m2g.'; % g frame -> mol frame
g2 = diag(R_g2m*g*R_g2m.')

g = diag(Sys.g(3,:));
R_m2g = erot(Sys.gFrame(3,:)); % mol frame -> g frame
R_g2m = R_m2g.'; % g frame -> mol frame
g3 = diag(R_g2m*g*R_g2m.')


gavg=(g1+g2+g3)/3
thanasis
Local Expert
Posts: 236
Joined: Thu Jan 21, 2016 6:28 am
Location: Strasbourg

Re: Exchange interaction for non parallel spin

Post by thanasis »

If what I have so far understood is correct, the gFrame that turns logal g-tensors like the green arrows indicate is:

Code: Select all

Sys.gFrame = [0 90 0; -30 -90 0; 30 90 0]*pi/180;
and the way the red arrows indicate it is:

Code: Select all

Sys.gFrame = [90 90 0; 60 -90 0; -60 90 0]*pi/180;
If I am correct, I hope this helps.
If I am wrong, could someone correct me?
Post Reply