Уже дошёл до такого:
Код: Выделить всё
OPT=(gdiis,verytight) Int=UltraFine B3LYP/6-311++G(d,p) EmpiricalDispersion=GD3BJ GFInput NMR SCRF=(PCM, Solvent=Chloroform) Freq pop(FULL,nboread) gfoldprin
Код: Выделить всё
OPT B3LYP/6-311++G(d,p).....
Код: Выделить всё
OPT=(gdiis,verytight) Int=UltraFine B3LYP/6-311++G(d,p) EmpiricalDispersion=GD3BJ GFInput NMR SCRF=(PCM, Solvent=Chloroform) Freq pop(FULL,nboread) gfoldprin
Код: Выделить всё
OPT B3LYP/6-311++G(d,p).....
Гесс писал(а): ↑Ср окт 11, 2017 4:43 am1) а может это реально не минимум? Как вообще выглядит колебание - вращение метильной группы?
2) попробуйте сделать IRC, правда если это артефакт то расчет умрет на первых же итерациях сколько не играй с шагом.
3) очень может быть эффектом от континиумного растворителся. Если проблема не воспроизводится в газофазном расчете надо либо пересчитывать все с сольватацией несколькими молекулами настоящего растворителя либо долго играть с параметрами SCRF
4) дальнейшее увеличение базиса для столь малой молекулы мне кажется бесполезным и вредным.
5) ежели че - то могу найти до десятка архивных блоков гауссиана для тетраметилсилана с частотами на B3LYP, разные меньшие базисы, все минимумы, кажется растворителя нет.
Код: Выделить всё
OPT B3LYP/6-311++G(d,p) EmpiricalDispersion=GD3BJ...
Когда просто нужна геометрия "для себя", небольшие значения можно почти смело игнорировать. Но если людям показывать - не солидно.
Спасибо! ПопробуюYuraM писал(а): ↑Ср окт 11, 2017 10:15 amjust a few quick hints/comments:
a) when increasing the grid size (Integral(Grid=ultrafine)) did you try to add "CPH(Grid=ultrafine)"? The problem is that for some strange reasons Gaussian 09 reduce the grid size significantly upon calculation of the energy second derivative, no matter which grid size is initially specified. This was pointed out by the reviewer of one of our paper, I checked it out and it turned out to be true.
b) The problem might be also coming from the dispersion correction and/or SCRF. If the option above does not help, then could you perform two separate calculations: one with dispersion/ no SCRF and another with SCRF and no dispersion.
Код: Выделить всё
CPH(Grid=ultrafine)
Код: Выделить всё
CPHF(Grid=ultrafine)
Ну "крайний" это зависит от настройки амплитуды колебаний, с тем размахом как на гифке думаю будет перебор. Но да, с несколько сдвинутой геометрии или IRC. Опять же можно подрезать базис для скорости. Если это не анион (желтый это сера? это тиоанион?) то можно было бы с чистой совестью отказаться скажем от второго плюса.
Гесс писал(а): ↑Ср окт 11, 2017 7:14 pmНу "крайний" это зависит от настройки амплитуды колебаний, с тем размахом как на гифке думаю будет перебор. Но да, с несколько сдвинутой геометрии или IRC. Опять же можно подрезать базис для скорости. Если это не анион (желтый это сера? это тиоанион?) то можно было бы с чистой совестью отказаться скажем от второго плюса.
Код: Выделить всё
CPH(Grid=ultrafine)
Насколько я понял, для TMS это ещё может и сработает - все остальные соединения не симметричные.
Сейчас этот форум просматривают: нет зарегистрированных пользователей и 20 гостей