VDD charges

[surius]

Коллеги, прошу помочь советом в вопросе расчета зарядов по Вороному, а именно: в чем можно это сделать?
Я сталкивался с реализацией VDD только в ADF, но вот он мне как раз и не подходит. А проще говоря нет желания заморачиваться с орбиталями
Слейтеровского типа. В общем буду рад любой наводке.

[alxyppv]

Ну, если Вам необязательно аналитически их считать, то можно делать куб с плотность а потом считать при помощи программы Bader (http://theory.cm.utexas.edu/bader). Правда при таком методе есть проблемы с численным интегрированием вблизи ядра, поэтому куб надо делать с частой сеткой

[surius]

Ну, если Вам необязательно аналитически их считать, то можно делать куб с плотность а потом считать при помощи программы Bader (http://theory.cm.utexas.edu/bader). Правда при таком методе есть проблемы с численным интегрированием вблизи ядра, поэтому куб надо делать с частой сеткой

интересно, можно попробовать. Правда для Бадера "Bader" мне нравиться существенно меньше чем @IM@LL, но это другая песня. А вы про математику ничего не знаете? Я в плане реализации его там и там? Он в АДФ реализован аналитически?

[alxyppv]

так там же совсем по разному все сделано. Понятно что @IM@LL гораздо лучше по ряду причин. Хотя бы потому, что там расчет аналитический, и поэтому много всего посчитать можно. "Bader" это код именно для численного интегрирования плотности на сетке (ну типа cube файла). Т.е. он берет этот куб, ищет там границы атомов, а потом складывает начения на сетке в пределах найденных границ. Для средних и больших молекул это может оказаться на порядки быстрее честного аналитического расчета @IM@LL. Но проблема, как я уже говорил выше, в интегрировании плотности возле ядра, где значения меняются очень резко, и поэтому для качественного расчета, особенно для тяжелых атомов, нужна ОЧЕНЬ частая сетка. Про адф ничего сказать не могу, не работал я с ним. Но по моему аналитического варианта в нем нет...

[alxyppv]

вот еще можно посмотреть код от китайских товарищей. Я когда-то пробовал, что-то мне не понравилось, но это давно было, а автор его постоянно улучшает. В общем, я бы посмотрел.
http://multiwfn.codeplex.com/