## Rotation patterns

General forum for EasySpin: questions, how to's, etc.

### Rotation patterns

Hi All,
I succeeded to simulate the rotation patterns of a spin ladder system imbedded in triclinic single crystal using the simplest Hamiltonian constituting of Zeeman term with anisotropic g-tensor and ZFS D and E. The rotations were around 3 perpendicular axes in the respect of the ladder which is aligned along the [011] crystal axis ‘or the ladder legs direction’. But in the simulation, I used [100] and [010] and [210] as rotN which are almost perpendicular regarding to the crystal symmetry parameters.
The point is that for each rotation the main g-tensor and D-tensor are the same but the g-frame and D-frame were different. The difference between the g and D Frame in each rotation pattern experiment can be seen in term of non-collinearity of the principal axes which is a common effect in such system. But the frames of both tensors are not the same from rotation experiment to another. In first glance, I can say that is the effect of non-perpendicular axes that I chose to make simulation, but as you can see below the difference if g frame from rotation to another was so great and I could not make a connection between the frame changing in respect of the rotation.
Below you may find the tables of the Hamiltonian parameters and the cod used to simulate those rotation pattern experiments.
My question is how it could be possible physically to explain this disparity of g and D frames as function of rotation axes and the difference between g frame from rotation to another where the fit was almost perfect?
Best regards.
Khalif.
Orientation of g- and D-tensors within the crystal frame
Rotation axis [100] [010] [210]
Euler angle α β γ α β γ α β γ
g-Frame (°) 0 46 60 0 40 32 95 70 85
D-Frame (°) 0 0 0 37 42 12 71 83 24
cod:
Sys1 = struct ('S',1, 'g',[2.041 0 0;0 2.105 0;0 0 2.001], 'D', mt2mhz([12 0 0;0 26 0;0 0 -38]));
Sys1.gFrame=[0 42 60]*pi/180;
Sys1.DFrame=[0 0 0]*pi/180;
Exp.mwFreq = 9.8;
Exp.Range = [250 450];
Exp.CrystalSymmetry = 'P1';
rotN = [1 0 0]; N = 31;
[phi,theta] = rotplane(rotN,[0 pi],N);
chi = zeros(N,1);
Exp.CrystalOrientation = [phi(:) theta(:) chi];
Opt.Output = 'separate';
Opt.Method='perturb2';
Bres1 = resfields(Sys1,Exp,Opt);
plot(Bres1,(theta*180/pi))
Khalif
Newbie

Posts: 5
Joined: Wed Nov 01, 2017 5:54 pm

### Re: Rotation patterns

Be careful how you specify the rotation axis. EasySpin needs it represented in a coordinate frame with orthogonal axes. In a triclinic system, the angles alpha, beta, gamma might deviate from 90 degrees. In that case you can't use [1 0 0], [0 1 0], [2 1 0], which are defined in the oblique abc frame, but you need to transform these vectors to an orthogonal frame (e.g. ab*c* IIRC).
Stefan Stoll
EasySpin Creator

Posts: 528
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

### Re: Rotation patterns

Dear Stefan,
I am agreeing, the rotations must be around orthogonal axes. In fact, the ladder is aligned along [011], where the experimental rotation was usually done around this later and the axes constituting the perpendicular plan of [011]. But EasySpin does not simulate the entire rotation pattern for rotation axes [hkl] with k≠0. For example, for a rotation around [011], the simulation could be show only the part between the magic angles;
Would like to explain to me why please ?
Best regards.
Khalif.
Khalif
Newbie

Posts: 5
Joined: Wed Nov 01, 2017 5:54 pm

### Re: Rotation patterns

The two lines
Code: Select all
`rotN = [1 0 0]; N = 31;[phi,theta] = rotplane(rotN,[0 pi],N);`
in your script define the rotation completely. `rotN` is the axis, `N` is the number of points, and the chosen angle range is `[0 pi]`. If you want a full rotation, then use `[0 2*pi]` as angle range.
Stefan Stoll
EasySpin Creator

Posts: 528
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

### Re: Rotation patterns

Dear Stephan,
It does not work even for [0 2*pi], the figure was the same without any changing (see the cod and the figure).
Best regards
Khalif.

Code: Select all
`close allclear allclcSys1 = struct ('S',1, 'g',[2.041 0 0;0 2.105 0;0 0 2.001], 'D', mt2mhz([12 0 0;0 26 0;0 0 -38]));Sys2 = struct ('S',1, 'g',[2.041 0 0;0 2.105 0;0 0 2.001], 'D', mt2mhz([10 0 0;0 16 0;0 0 -26]));Sys3 = struct ('S',1/2, 'g', [2.041 0 0; 0 2.105  0;0 0 2.001]);% Experimental parametersExp.mwFreq = 9.48; Exp.Range = [250 470];Exp.CrystalSymmetry = 'P1';rotN = [0 1 1];  N = 31;[phi,theta] = rotplane(rotN,[0 2*pi],N);chi = zeros(N,1);Exp.CrystalOrientation = [phi(:) theta(:) chi];% Simulate spectraOpt.Output = 'separate';  Opt.Method='perturb2';Bres1 = resfields(Sys1,Exp,Opt);Bres2 = resfields(Sys2,Exp,Opt);Bres3 = resfields(Sys3,Exp,Opt);% plottinghold allplot(Bres1,(theta*180/pi),'r')plot(Bres2,(theta*180/pi))plot(Bres3,(theta*180/pi))%legendxlabel('magnetic field (mT)');ylabel('theta (°)');`
Attachments
Rot Pat.jpg (23.25 KiB) Viewed 536 times
Khalif
Newbie

Posts: 5
Joined: Wed Nov 01, 2017 5:54 pm

### Re: Rotation patterns

You need to plot the resonance fields against the actual rotation angle, and not against theta:
Code: Select all
`rho = linspace(0,2*pi,N)plot(Bres1,rho*180/pi)`
Stefan Stoll
EasySpin Creator

Posts: 528
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

### Re: Rotation patterns

Hello,
It’s work very well with rho, after defining a non-collinear adequate g- and D-frame, but what angle is rho? between which axes is it?
‘resfields’ calculates the resonance field as function of theta (angle between anisotropic axes and the magnetic field) right…., but rho seems to be independent angle!!!!!
I need to understand those points, I know that they are naïve questions but is very important to me to understand them ‘very well’ to make good point.
Khalif
Newbie

Posts: 5
Joined: Wed Nov 01, 2017 5:54 pm

### Re: Rotation patterns

rho is the rotation angle around your axis of rotation, corresponding to the angle on the goniometer in the experimental setup.
Stefan Stoll
EasySpin Creator

Posts: 528
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

### Re: Rotation patterns

Many Thanks
Khalif
Newbie

Posts: 5
Joined: Wed Nov 01, 2017 5:54 pm