Оптимизация геометрий депротонированных молекул в растворах

[Drg]

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

[Гесс]

Вы про сольватацию анионов? Explicit vs implicit solvation?
Ну да, молекула "живой" воды рядом с местом специфичной адсорбции полезна, и потомвсе это в континиумник.

[Drg]

Вы про сольватацию анионов? Explicit vs implicit solvation?
Ну да, молекула "живой" воды рядом с местом специфичной адсорбции полезна, и потомвсе это в континиумник.

Да, вопрос в частности про это и про ряд сопутствующих вопросов.
1) Правильно ли я понимаю, что в Gaussian придется в геометрию для каждого атома ввести опции типа (Fragment=1) и (Fragment=2), где для отрицательных ионов кислорода придется выставить свой заряд (-1) и мультиплетность (2)? Будет ли это корректно считать возбужденные синглетные состояния для всей молекулы в целом в Gaussian?
2) Достаточно ли по одной молекуле воды около O^{-1}, чтобы это правдоподобно локализовало (экранировало) воздействие O^{-1} на скелет молекулы (без больших отклонений расчета от реалий раствора)?

[Гесс]

1) нет, необязательно. Можно все считать вместе. Да, спектры от всей системы будет корректно. Нет, присваивать минус одному атому кислорода [может быть] некорректно, скажем карбоксилат анион симметричен.
2) кто ж его знает. "Реалии раствора" и без отрицательных ионов могут нехило отличаться от дфт-минимума.

[Drg]

1) нет, необязательно. Можно все считать вместе. Да, спектры от всей системы будет корректно. Нет, присваивать минус одному атому кислорода [может быть] некорректно, скажем карбоксилат анион симметричен.
2) кто ж его знает. "Реалии раствора" и без отрицательных ионов могут нехило отличаться от дфт-минимума.

Спасибо за ответы. Но я не очень понял как иначе, кроме как с опцией Fragment и явным описанием его спецификации (знака заряда и мультиплетности) еще другим способом явно добавить лишний электрон при нужном атоме. Если я просто сотру водород в геометрии, то он уберется вместе с электроном (это будет не депротонирование). Вариант с заданием некоего эффективного заряда при атоме тоже сильно про другое (электронная конфигурация не та будет все-равно).
Ну и, вроде, пишется, что опция Fragment не значит раздельный расчет (за исключением некоторых методов, не DFT), а просто как раз позволяет указывать в молекуле локальные отличия электронной структуры от default.

[Гесс]

Вы просто даете суммарный заряд системы -1, мультиплетность с уходом протона не изменится.
ДФТ-расчет должен сам распределить электроны по системе, не надо ему навязывать где плюс где минус.