Я бы в целом рекомендовал поизучать тот же ReaxFF, чтобы примерно представлять, насколько замороченные потенциалы надо брать, чтобы получалось "нормально". Там как раз тоже потенциалы отдаленно похожие на Морзе, но не настолько примитивно и в лоб, как вы предлагаете. Ну и еще много всего типа электростатики, которую вы почему-то решили выкинуть.Vit Nhoc писал(а): ↑Сб авг 19, 2023 1:51 pmЯ могу сформулировать свою идею так. Когда мы считаем силовое поле молекулы и решаем колебательную задачу в приближении гармонического осциллятора, мы как бы тоже используем простую модель молекулярной механики, крайне грубую – для всех параметров энергия описывается параболой. Если же использовать более физичную модель ММ, результаты должны выйти лучше.
Предположим, мы посчитали геометрию и частоты молекулы BH3:
У молекулы 4 атома и 6 степеней свободы. Расчёт дал одну энергию, шесть градиентов энергии (все нулевые) и кажется 6*5/2=15 силовых постоянных. Прошу меня поправить если что-то путаю. Суммарно имеем 22 параметра, по которым мы тренируем ММ функционал (из-за симметрии реальных параметров в несколько раз меньше).
В обычных методах ММ, как я понимаю, мы имеем квадратичную функцию для вкладов по связям, т.е. это два параметра – положение минимума и изгиб. Таких же два параметра имеем для валентных углов. Т.е. суммарно 4 параметра.
Если же использовать более физичную модель, мы можем описать вклад в энергию по связям через потенциал Морзе:
U(R)=D(1-exp(-a(R-Re))^2
Тут три параметра, предположим один мы оцениваем эмпирически, из типа связи. И ещё мы описываем взаимодействия между несвязанными протонами через потенциал Леннарда-Джонса:
U(R)=4*E*(sigma/R^12-sigma/R^6)
Тут ещё два параметра. Суммарно таким образом мы имеем снова четыре параметра, но более физичных; и по идее, когда мы их посчитаем, мы сможем решить ангармоническую задачу и найти ангармонические частоты/обертоны.
Вот ссылка https://pubs.acs.org/doi/abs/10.1021/jp709896w, там в дополнительных материалах расписаны формулы.