Molar Susceptibility Plot

General forum for EasySpin: questions, how to's, etc.
Post Reply
Satpathi7
Newbie
Posts: 1
Joined: Wed Oct 30, 2024 2:39 am

Molar Susceptibility Plot

Post by Satpathi7 »

Hi!
I have two queries; first one may be very naive, but can anyone please let me know how easyspin is calculating the molar susceptibility of a molecule without the knowledge of molar mass and mass density! Any reference on this would be really helpful.

Second query is; why the chimolT vs T plot is changing drastically when full 3x3 g tensor given as input and when principal values are given with Euler gFrame. I believe the values should come same in this two cases. Please let me know where I am going wrong. My code is given below:

clear, clf, clc
cm = 100clight/1e6; % Conversion constant from cm-1 to MHz
B = 100; % mT
T = 0.01:0.05:300; % K
Cr.S = [3/2, 3/2];
Cr.L = [1 1];
Cr.soc = [-17.7; 232.5]
cm;
%Cr.D = [0,-24.474]cm;
Cr.g = [1.992813; 1.793090];
Cr.gFrame = [0 0 0;78.8607 19.8848 -89.2224];
%Cr.g = [0.643345 1.850015 0.020856; -1.772939 0.623834 -0.656308;-0.604806 0.187924 1.927663; 0.118414 -0.305717 -0.078879; 1.461836 0.511324 -0.315881;-3.786474 -1.551041 1.782057];
Cr.J = 77
cm; % cm^-1 -> MHz
Exp.Field = B;
Exp.Temperature = T;

Opt.Output = 'chimolT';
Opt.Units = 'CGS';
chimolT = curry(Cr,Exp,Opt);

plot(T,chimolT);
xlabel('temperature (K)');
ylabel('molar magnetic susceptibility (K cm3 mol^{-1})');
grid on;

thanasis
Local Expert
Posts: 253
Joined: Thu Jan 21, 2016 6:28 am
Location: Strasbourg

Re: Molar Susceptibility Plot

Post by thanasis »

  1. Models calculating susceptibilities, do so for the molar values of these quantities. The experimentalist must previously bring the measurements to these quantities by knowledge of sample mass, FW and experimental magnetic field, and also diamagnetic and sample holder corrections. The model is agnostic towards these values.

  2. I won't get into why you define a SOC for Cr(III) (which if octahedral should have L=0). I will just focus on the fact that you define either:
    -two isotropic g-factors: Cr.g = [1.992813; 1.793090]
    -or two anisotropic ones whose principal values have nothing to do with the previous isotropic ones. Try diagonalizing the anisotropic g-tensors:

    Code: Select all

    [V1,D1] = eig(Cr.g(1:3,1:3))
    [V2,D2] = eig(Cr.g(4:6,1:3))

    and you will see that D1 and D2 are all over the place.

Stefan Stoll
EasySpin Creator
Posts: 1108
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

Re: Molar Susceptibility Plot

Post by Stefan Stoll »

Other answer to 1.: EasySpin's curry calculates the magnetic susceptibility of an individual paramagnetic center and then multiplies the result by Avogadro's number (avogadro in EasySpin) to obtain the molar susceptibility.

Post Reply