Again, I have a question regarding exchange coupled system when simulating hyperfine splittings.
The example case is for biradicals known with strong exchange coupling interactions. I wanted to simulate the hyperfine splittings compare to the spectrum.
The normal code for an isotropic simulation for the single radical is like below:
Code: Select all
Sys.S=1/2
Sys.g=2
Sys.Nucs='14N'
Sys.lw=0.2
Sys.A=20
Exp.Range=[330,350]
Exp.mwFreq=9.5
Code: Select all
Sys.S=[1/2 1/2]
g=[2 2]
Sys.A=[20 0;0 20]
Sys.lw=0.2
Sys.Nucs='14N,14N'
Sys.J=-3.3e7
Here is the problem, I am working on a more complicated system where I have a lot more Nuclei contributing to the hfi, what I thought might be reasonable to reduce the computing time is an example shown below for the same biradicals:
Code: Select all
Sys.S=1
g=2
Sys.A=[20 0;0 20]
Sys.lw=0.2
Sys.Nucs='14N,14N'
Code: Select all
Sys.S=1
g=2
Sys.A=[20 20]
Sys.lw=0.2
Sys.Nucs='14N,14N'
The first code gave me a mixing feature of hyperfine splittings, those on the edge will be half of hfi parameters, and the central ones are not.
The second code will just give a same looking from the single radical which is reasonable, the program just treats the S=1 spin as a single spin system.
When I define the spin system as a triplet, I thought I told the program that only triplet will populate and hence strong exchange coupling, while the result didnot like exactly what I wanted.
This might be no difference for a small system like the codes I showed here, while this become a huge difference when the spin system is much larger. The actual case I am working on have at least 4 nuclei, the A matrix will be very large and difficult to compute when I introduce an interacting spin system with exchange coupling value.
Another question related to this topic is about zero-field splitting term.
An example is shown below:
Code: Select all
S: [0.5000 0.5000]
g: [2 2 2; 2 2 2]
D: [550 0]
lw: 2
J: -7900000
So my question is, in order to do so, what should I do?
Thanks a bunch!