Оптимизация и симметрия в Гауссиане

[Vlas]

Здравствуйте,

подскажите, пожалуйста, как сделать так, чтобы симметрия начальная молекулярного комплекса сохранялась при оптимизации, например, симметрия D2d. Стартовую геометрию я задаю в декартовых координатах. Я пробовал делать Symm=Follow, но симметрия меняется на C1.

[VTur]

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

[alxyppv]

по умолчанию гауссиан сохраняет симметрию начальной геометрии (его нужно специально просить это не делать). Другое дело, что по умолчанию у него очень жесткие требования на симметричность - т.е. при небольших отклонениях от симметрии он может эту симметрию не заметить. Надо или жестко симметризовать координаты (например, при помощи ChemCraft или GaussView) или попросить о менее жестком критерии во входном файле добавлением ключа Symmetry(Loose)

[Vlas]

Спасибо, попробую

[Neitrino]

Посмотри:
IOp(2/15)
Symmetry control.

-1 Turns on symmetry; same as 0 for molecules but turns on assignment of space group operations for PBC.

0 Leave symmetry in whatever state it is presently in.

1 Unconditionally turn symmetry off. Note that symmetry is still called, and will determine the framework group. However, the molecule is not oriented.

2 Bring the molecule to a symmetry orientation, but then disable further use of symmetry.

3 Don't even call Symm.

4 Call Symm once with loose cutoffs, symmetrize the resulting coordinates then confirm symmetry with tight cutoffs.

5 Recover the previous symmetry operations from the RWF, and confirm that the new structure has the same symmetry.

6 Same as 5, but get symmetry info from the checkpoint.

00 Default (10).

10 Do re-orientation for PBC.

20 Suppress re-orientation for PBC.

100 Turn on symmetry operations for PBC.