Pcs2chi.m

From Spinach Documentation Wiki
Jump to: navigation, search

Uses linear least squares procedure on top of Equation 10 from http://dx.doi.org/10.1039/c4cp03106g to extract the rank 2 com- ponent of the susceptibility tensor from DFT hyperfine coupling tensors and experimentally observed pseudocontact shifts.

                chi=pcs2chi(hfcs,shifts,isotopes)

where hfcs is a cell array of 3x3 hyperfine coupling tensors in MHz, shifts is a vector of the observed pseudocontact shifts, ex- cluding the diamagnetic contribution, and isotopes is a cell ar- ray of character strings specifying isotopes that exhibit each of the chemical shifts supplied, for example {'1H','13C'}.