Ошибка #2070 в Gaussian 09W

[VTur]

И не могли бы вы подсказать какой лучше метод\базис использовать для расчета dG dH

Для полной энергии молекул из атомов 1-3 строки т. М. я использую B3LYP/aug-cc-pvQZ. Ошибка где-то до 10 кДж/моль. Но считаю в ангармоническом приближении. Ангармонические частоты и постоянные ангармоничности считаю в B3LYP/6-311++g(3df,3pd)

и для расчета смещения электронной плотности и величины полузаряда на атомах???

Термодинамика здесь неприменима. Термодинамические свойства относятся к веществу, а не к молекуле. Молекулы не характеризуются т/д функциями.
Что Вы имеете в виду под

смещения электронной плотности и величины полузаряда на атомах

[chronos377]

Спасибо Вам большое за ответы. Под смещением электронной плотности и величиной полузаряда я понимаю вероятность нахождения электрона к в области около ядра атома и отрицательным зарядом, который он создает в результате. Если глупости несу извините.

[VTur]

Под смещением электронной плотности и величиной полузаряда я понимаю вероятность нахождения электрона к в области около ядра атома и отрицательным зарядом, который он создает в результате. Если глупости несу извините.

Все равно непонятно

Разность электронной плотности промолекулы и молекулы, что считает программа Цирельсона?

Вероятность нахождения электрона в какой-то области есть интеграл от электронной плотности по этому бассейну, деленный на число электронов (безразмерный интеграл).

[Prokop]

Здравствуйте! Не получается оптимизировать молекулы имидазолидов (пример в архиве). Нарисовал при помощи GaussView 5.0. При запуске расчета выскакивает ошибка 2070. Опыта работы с программой не имею. Помогите разобраться, в чем ошибка. За помощь буду очень благодарен. Спасибо!

AdCOIm.rar

[VTur]

Попробуйте убрать путь к чек файлу. Оставьте только
%chk=AdCOIm.chk

[Prokop]

Спасибо! Расчет длился немного дольше. Однако опять выскакивает та же ошибка, но в этот раз в output файле указано:
Subroutine NAOANL could not find a S-type valence orbital on atom C 1.
IVAL : 1 1 0 0 M : 1 LA : 1
Error
Error termination via Lnk1e in C:\G09W\l607.exe at Fri Mar 15 18:12:11 2013.
Job cpu time: 0 days 0 hours 0 minutes 2.0 seconds.
File lengths (MBytes): RWF= 9 Int= 0 D2E= 0 Chk= 1 Scr= 1

Подскажите, пожалуйста, что это значит и как это исправить?

[Rusia]

Да там метод PM6, потянет.

Тогда надо пользоваться мопаком.

[o-oxhem]

нужно удалить pop=nbo (если не удалить, то расчет заселенностей будет выполняться на каждом цикле оптимизации геометрии). pop=nbo не совместимо с pm3 (крайней мере nbo 3.1). Кстати, для полуэмипирических методов применяются манимальные базисные наборы (STO-3G) не думаю, что с такими наборами имеет смысл выполнять анализ NBO.
Файл приложен. Комментарии: для большинства квантово-химических методов нет необходимости указывать geom=connectivity и соответствующую матрицу связности. Рекомендую указывать %mem=выделенная оперативная память

[Prokop]

Спасибо! Теперь нормально удается считать. Вы мне очень помогли!

[Nbveh]

Помогите пожалуйста разобраться. Выскакивает ошибка #2070. Не идут расчеты йодпроизводных бензола во вложении input и uotput

Расчеты.rar

Вроде все правильно, молекула предварительно рассчитана полуэмперикой в hyperhchem (для длины связи I-C). Аналогичные соединения (Br, Cl ...) считаются без проблем.

[Гесс]

Помогите пожалуйста разобраться. Выскакивает ошибка #2070. Не идут расчеты йодпроизводных бензола во вложении input и uotput

Расчеты.rar

Вроде все правильно, молекула предварительно рассчитана полуэмперикой в hyperhchem (для длины связи I-C). Аналогичные соединения (Br, Cl ...) считаются без проблем.

Там в последних строках написано:
Standard basis: 6-31G(d,p) (6D, 7F)
Atomic number out of range for 6-31G basis set.

И бром в этом базисе тоже считаться не должен.
С помощью Gen можно дозадать базис для Йода, например 6-311, LANLDZ, ну или на выбор с https://bse.pnl.gov/bse/portal

[Nbveh]

Про базис то забыл...
Спасибо!

[Гесс]

Только учтите что задав для всех атомов 6-311 вы получите то же самое, 6-311 для йода есть, но в Гауссиан он не вписан, его надо добавлять именно через Gen.

[Любитель_Манниха]

Если гауссиан выдаёт ошибку 2070 и последняя строка в аутпуте:
open-new-file
- тупо памяти не хватает, или возможны варианты?

[Гесс]

Если гауссиан выдаёт ошибку 2070 и последняя строка в аутпуте:
open-new-file
- тупо памяти не хватает, или возможны варианты?

а можно чуть больше последних строк?

[Любитель_Манниха]

бла-бла копирайты-ссылки......
******************************************
%chk=E:\filename.chk
%mem=1GB
fname=E:\filename.chk fd = -1

open-new-file

Усё
Инпут генереровал в гауссвью на основе мол-файла

[Гесс]

Ну методом научного тыка:
1) а с меньшей задачей ранится?
2) а место на диске есть?

[Любитель_Манниха]

1) То же самое
2) > 100 Gb

[Гесс]

Шапку инпута покажите....

[Любитель_Манниха]

Инпут:
[cut]

%chk=E:\2-Ac-pyrrole.chk
# opt b3lyp/3-21g geom=connectivity

Title Card Required

0 1
C 17.62080000 -14.52000000 0.00000000
N 18.40260000 -13.44400000 0.00000000
C 18.40260000 -15.59600000 0.00000000
C 19.66750000 -13.85500000 0.00000000
C 19.66750000 -15.18500000 0.00000000
C 16.29080000 -14.52000000 0.00000000
O 15.62580000 -13.36820000 0.00000000
C 15.62580000 -15.67180000 0.00000000
H 18.24809887 -12.96846935 0.86602540
H 18.07196757 -16.61363559 0.00000000
H 20.53315022 -13.22607259 0.00000000
H 20.53315022 -15.81392741 0.00000000
H 14.57381634 -15.47627410 0.00000000
H 15.88428943 -16.23288515 0.87365134
H 15.88428943 -16.23288515 -0.87365134

1 2 1.0 3 2.0 6 1.0
2 4 1.0 9 1.0
3 5 1.0 10 1.0
4 5 2.0 11 1.0
5 12 1.0
6 7 2.0 8 1.0
7
8 13 1.0 14 1.0 15 1.0
9
10
11
12
13
14
15

[/cut]