Расчет парамагнитных комплексов в EasySpin

вопросы строения молекул и квантовой химии
Ответить
Няшка
Сообщения: 601
Зарегистрирован: Вс фев 24, 2013 8:44 am

Расчет парамагнитных комплексов в EasySpin

Сообщение Няшка » Пн мар 04, 2013 11:42 am

Всем доброго времени суток!

В связи с интенсивным изучением EasySpin возник вопрос. В documentation, опубликованном на сайте http://www.easyspin.org/documentation/u ... manynuclei, приводится пример расчета спектра Cu2+ porphyrin (расчет на нескольких ядрах):

Sys.S = 1/2;
Sys.g = [2 2.2];
Sys = nucspinadd(Sys,'63Cu',[50 500]);
A = [20 30];
Sys = nucspinadd(Sys,'14N',A);
Sys = nucspinadd(Sys,'14N',A);
Sys = nucspinadd(Sys,'14N',A);
Sys = nucspinadd(Sys,'14N',A);
Sys.lwpp = 0.5;
Exp.mwFreq = 9.5;
Exp.Range = [260 380];
Opt.Method = 'matrix';
[x,y1] = pepper(Sys,Exp,Opt);
Opt.Method = 'perturb2';
[x,y2] = pepper(Sys,Exp,Opt);
plot(x,y1,'r',x,y2,'b');.

Когда я ввожу все это в MatLab, то на экране появляется следующее:

??? Error using ==> pepper at 146
Sys.lw and Sys.lwpp cannot be used simultaneously.
Error in ==> C:\EasySpin\easyspin\pepper.p>pepper at 105.

Что это означает?

Заранее спасибо за ответы и прошу прощения за идиотский вопрос.

alxyppv
Сообщения: 560
Зарегистрирован: Сб апр 07, 2007 11:23 am

Re: Расчет парамагнитных комплексов в EasySpin

Сообщение alxyppv » Пн мар 04, 2013 2:51 pm

Sys.lw и Sys.lwpp это параметры, задающие ширину линии
Под рукой сейчас программы нет (вечером смогу поиграть), но что если попробовать только один метод использовать (Opt.Method = 'perturb2';)
А более простые примеры у Вас уже получаются?
А.П.

Няшка
Сообщения: 601
Зарегистрирован: Вс фев 24, 2013 8:44 am

Re: Расчет парамагнитных комплексов в EasySpin

Сообщение Няшка » Пн мар 04, 2013 4:30 pm

alxyppv, спасибо!

Я пробовала только один метод (по очереди), но программа выдавала тоже самое, я, просто ради баловства, убирала Sys.lw и Sys.lwpp, результат тот же.
Да, более простые примеры у меня уже получаются, около десятка комплексов с Ni и Co (предварительно я сняла для них ЭПР спектры, проанализировала их и определенные мной параметры совпали с заявленными в литературе), правда, наверное, из-за того, что я учитываю только один атом спектры получаются как бы нарисованные, хотя и форма и ширина и .... хорошо воспроизводятся.

Няшка
Сообщения: 601
Зарегистрирован: Вс фев 24, 2013 8:44 am

Re: Расчет парамагнитных комплексов в EasySpin

Сообщение Няшка » Вт мар 05, 2013 6:55 am

Если это куму-нибудь понадобится, то у меня все посчиталось, но вот в таком виде:

Sys.S = 1/2;
Sys.g = [2 2.2];
Sys.lwpp = 0.5;
Sys = nucspinadd(Sys,'63Cu',[50 500]);
A = [20 30];
Sys = nucspinadd(Sys,'14N',A);
Sys = nucspinadd(Sys,'14N',A);
Sys = nucspinadd(Sys,'14N',A);
Sys = nucspinadd(Sys,'14N',A);
Exp.mwFreq = 9.5;
Exp.Range = [260 380];
Opt.Method = 'matrix';
[x,y1] = pepper(Sys,Exp,Opt);
Opt.Method = 'perturb2';
[x,y2] = pepper(Sys,Exp,Opt);
plot(x,y1,'r',x,y2,'b');.

Ответить

Вернуться в «квантовая химия и моделирование»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 48 гостей