natural isotopes of equivalent nuclei

General forum for EasySpin: questions, how to's, etc.
burkhard
Newbie
Posts: 3
Joined: Fri Apr 07, 2017 4:47 am

natural isotopes of equivalent nuclei

Post by burkhard »

Dear Stephan,
I tried the discussed issue of an radical with one Cl couplings (of 12 all equivalent nuclei) :

Code: Select all

Sys.g = 2;
Sys.Nucs = 'Cl';
Sys.n = [12];
Sys.A = [3];
Sys.lwpp = [0 0.01];
Exp.mwFreq = 9.7;
Exp.CenterSweep = [346.5 28];
garlic(Sys,Exp);
I do get the Error:
Error using isotopologues
Cannot compute isotope patterns for equivalent nuclei.
Error in expandcomponents
Error in garlic
Error in isotest (line 8)
garlic(Sys,Exp);
So I do something wrong, it looked so easy when you showed it to me. :oops:
Burkhard.
katarkon
Local Expert
Posts: 182
Joined: Mon Jan 12, 2015 4:01 am

Re: natural isotopes of equivalent nuclei

Post by katarkon »

EasySpin does not support equivalent nuclei for tme mixtures of isotopes. You should define each nucleus separeately. For your case it should be defined like this:

Code: Select all

Sys.Nucs = 'Cl,Cl,Cl,Cl,Cl,Cl,Cl,Cl,Cl,Cl,Cl,Cl';
Sys.A = [3 3 3 3 3 3 3 3 3 3 3 3];
Stefan Stoll
EasySpin Creator
Posts: 1041
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

Re: natural isotopes of equivalent nuclei

Post by Stefan Stoll »

Burkhard, katarkon has it right. In the current ES version, for multiple equivalent nuclei, you have to spell them out if you want natural abundance mixtures. I know it would be nice if this wasn't necessary and one could just specify Sys.Nucs='Cl' and Sys.n=12. We'll put this on our to-do list.
burkhard
Newbie
Posts: 3
Joined: Fri Apr 07, 2017 4:47 am

Re: natural isotopes of equivalent nuclei

Post by burkhard »

Thank you Stepan and katarkon, with this method I was able to get the whole hyperfine pattern. As you guessed right, the direct way to simulate natural abundance equivalent nuclei would be a nice feature to have.
Stefan Stoll
EasySpin Creator
Posts: 1041
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

Re: natural isotopes of equivalent nuclei

Post by Stefan Stoll »

Thinking about this a little more, this is actually an important feature we should include in the next EasySpin version. The savings can be enormous.

In your case of N=12 Cl sites, there are 2^N = 4096 isotopologues (since we have one of 2 Cl isotopes, 35Cl and 37Cl, per site). If the sites were not magnetically equivalent, we would have to simulate all 4096 spectra and add them up with the appropriate abundance weights. And that's what the ES code with Sys.Nucs='Cl,Cl,Cl,Cl,Cl,Cl,Cl,Cl,Cl,Cl,Cl,Cl' does.

Now, however, the sites are all equivalent. This means, for example, that the isotopologue with 37Cl in sites 1 and 5 gives the same spectrum as the isotopologues with 37Cl in site 6 and 12, or any other pair of sites. Therefore, there are only a total of N+1=13 distinct spectra. Each spectrum is due to a number of equivalent isotopologues: [1 12 66 220 495 792 924 792 495 220 66 12 1], where the first entry refers to the 35Cl-only isotopologue, the second to the isotopologues with one 35Cl, etc. The simulation would therefore run much faster: 13/4096 = 0.003.

It is nontrivial to implement this for the general case of an element with more than 2 isotopes (e.g. Sn), but we will attempt it.
Stefan Stoll
EasySpin Creator
Posts: 1041
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

Re: natural isotopes of equivalent nuclei

Post by Stefan Stoll »

After some initial work, we are optimistic that this feature will be implemented in the next version.

Meanwhile, here is a list of isotopologue spin systems and their abundances that you need to run separately in the current version (5.1.10). Make sure the hyperfine couplings are properly scaled between 35Cl and 37Cl. The last two or three isotopologues can probably be skipped, since they are of very low relative abundance, and will probably not show above the noise in your experimental spectrum.

Code: Select all

 abundance  Sys.Nucs     Sys.n
 0.035864   35Cl         12 
 0.137549   35Cl,37Cl    11  1 
 0.241790   35Cl,37Cl    10  2 
 0.257594   35Cl,37Cl     9  3 
 0.185242   35Cl,37Cl     8  4 
 0.094728   35Cl,37Cl     7  5 
 0.035322   35Cl,37Cl     6  6 
 0.009676   35Cl,37Cl     5  7 
 0.001933   35Cl,37Cl     4  8 
 0.000275   35Cl,37Cl     3  9 
 0.000026   35Cl,37Cl     2 10 
 0.000002   35Cl,37Cl     1 11 
 0.000000   37Cl         12 
Stefan Stoll
EasySpin Creator
Posts: 1041
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

Re: natural isotopes of equivalent nuclei

Post by Stefan Stoll »

This is now implemented in version 5.2. The script you posted works now as expected.
burkhard
Newbie
Posts: 3
Joined: Fri Apr 07, 2017 4:47 am

Re: natural isotopes of equivalent nuclei

Post by burkhard »

Dear Stefan,
thank you :D
Without simulation all isotope contributions of 12 Chlorine and 4 Silicon separably, I could get the same with one garlic call :lol:
Burkhard.
Stefan Stoll
EasySpin Creator
Posts: 1041
Joined: Mon Jul 21, 2014 10:11 pm
Location: University of Washington

Re: natural isotopes of equivalent nuclei

Post by Stefan Stoll »

Great! Let us know once you've published on this. This could be a great case to include in the EasySpin examples.
katarkon
Local Expert
Posts: 182
Joined: Mon Jan 12, 2015 4:01 am

Re: natural isotopes of equivalent nuclei

Post by katarkon »

I see that the output of the function isotopologues() was changed. New version does not provide data about hyperfine splitting more. Is any way to rid it out?
Post Reply