очередные мудрые мысли от VTur

вопросы строения молекул и квантовой химии
Ответить
Аватара пользователя
Droog_Andrey
Сообщения: 2700
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: очередные мудрые мысли от VTur

Сообщение Droog_Andrey » Чт мар 14, 2013 3:02 pm

sanya1024 писал(а):А насчет максимума в перпендикулярном сечении -- я уже объяснила, откуда там сигмы. Это SGG 1s(A)+1s(B) и SGG 2s2px(A)+2s2px(B). Они (в виде квадратов) дают ненулевые вклады.
Они вместе с SGU дают тупо две несвязывающие s-орбитали. Два чуть поляризованных шара. В плоскости, равноудалённой от ядер, они дают очень пологий максимум. На фоне его пи-минимум должен проявляться. Но, как видим, этого нет.

Так что, снова повторюсь, все эти узлы и т.п. - условности. Атомные орбитали s, p, d, f и т.д. - просто сферические гармоники, математические ортогональные функции, которыми удобно описывать распределение электронной плотности вокруг ядра. Точно теми же функциями описывают, например, распределение реликтового фона Вселенной. Это ничуть не лучше описания профиля пилы, висящей на стене сарая, с помощью синусов разных частот.

Все представления об орбиталях и переходах между ними возникают благодаря удобству и простоте описания многоэлектронных состояний как совокупности одноэлектронных, т.к. набор и тех и иных дискретен. Иногда корреляция между этими наборами хорошая, и удаётся отследить какие-то закономерности, делать предсказания и т.п., но иногда она плохая (посмотрите на орбитали атома церия в основном состоянии 1G° на разных уровнях теории), поэтому, строго говоря, весь этот аппарат - "a rule of a thumb", не более того. Думаю, VTur именно это хотел донести.
amge писал(а):but also one weak ‘inverted’ bond...
Вот и начинаются всякие надстройки над моделью, чтобы залатать её слабые места. Вспоминаю про аналогичный случай с методом ВС и молекулой NOF3 :)
2^74207281-1 is prime!

Аватара пользователя
sanya1024
Сообщения: 1672
Зарегистрирован: Чт янв 20, 2011 3:24 pm

Re: очередные мудрые мысли от VTur

Сообщение sanya1024 » Пт мар 15, 2013 1:16 am

Droog_Andrey писал(а):Они вместе с SGU дают тупо две несвязывающие s-орбитали. Два чуть поляризованных шара. В плоскости, равноудалённой от ядер, они дают очень пологий максимум. На фоне его пи-минимум должен проявляться. Но, как видим, этого нет.
Мне все-таки так не кажется... Пи орбитали не вычитают эл. плотность на оси. Просто дают нулевой вклад. Так что если там что-то было (а оно было, Вы не отрицаете), никуда оно не денется.

Если б я могла открыть чем-то Ваш *.fch файл :( Почему-то ChemCraft на нем падает. Чем их еще открывают? Вы можете по отдельности показать тут (или мне в личку, если не хотите засорять тему) картинки отдельных квадратов 1SGG, 1SGU, 2SGG и 2SGU орбиталей? просто посмотрим, что это за сигмы.
Droog_Andrey писал(а): Так что, снова повторюсь, все эти узлы и т.п. - условности. Атомные орбитали s, p, d, f и т.д. - просто сферические гармоники, математические ортогональные функции, которыми удобно описывать распределение электронной плотности вокруг ядра.


А что будет, если вместо базиса атомно-центрированных функций взять функции двухатомных фрагментов? были такие работы где-то в 70-е или в 80-е, но быстро не найду. А если плоские волны? Плотность никуда не денется, а то, из чего она будет строиться, будет сохранять симметрию (в широком, топологическом смысле) исходного базиса.
Droog_Andrey писал(а): Все представления об орбиталях и переходах между ними возникают благодаря удобству и простоте описания многоэлектронных состояний как совокупности одноэлектронных, т.к. набор и тех и иных дискретен. Иногда корреляция между этими наборами хорошая, и удаётся отследить какие-то закономерности, делать предсказания и т.п., но иногда она плохая (посмотрите на орбитали атома церия в основном состоянии 1G° на разных уровнях теории), поэтому, строго говоря, весь этот аппарат - "a rule of a thumb", не более того. Думаю, VTur именно это хотел донести.
Так, да не так. Визуализировать одноэлектронные функции -- вполне реально, это наглядная конструкция. А вот визуализировать компоненты плотности в конструкции из диатомиков или плоских волн... в общем, если мне кто-нибудь покажет, как это выглядит, скажу спасибо :)

Насчет орбиталей тяжелого атома -- да, наверное, только серьезное спин-орбитальное взаимодействие способно поломать стройную музыку сфер картинку сферических гармоник (построенную в предположении сохранения орбитального углового момента). Ну так для тяжелых атомов и полноэлектронные базисы строят соответствующие, и гамильтонианы берут совсем другие. Кстати, а картинки для этого церия не покажете? что с ним не так? неужели перемешались s,p,d,f функции в одной оболочке?

Что же хотел нам сказать VTur -- трудно угадать, его мысль блуждает довольно хаотически, постоянно возвращаясь к необходимости разрушить Карфаген "неправильности" орбитальных представлений (там, наверное, странный аттрактор).
Вот и Татевский, покойничек, тоже: как, бывало, услышит слово орбиталь, так взовьется весь, кричит, аж трясется: "Никаких орбиталей нет!" Сама видела, как на защитах люди по полчаса оправдывались за сорвавшееся с языка запретное слово.
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)

VTur
Сообщения: 7357
Зарегистрирован: Пт авг 31, 2007 1:36 pm

Re: очередные мудрые мысли от VTur

Сообщение VTur » Пт мар 15, 2013 8:03 am

Droog_Andrey, Вы правильно поимаете. Орбитали - это математическая абстракция. Это нормальный термин. Нормальное математическое описание. Они не являются наблюдаемыми и относится к ним надо, как к математическим абстрактам. Об этом говрится во многих книгах, например, см. в руководстве для чайников
"Foresman J.B., Frisch A. Exploring chemistry with electronic structure methods" на стр. 19 справа внизу (на красном фоне).

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

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

Аватара пользователя
Droog_Andrey
Сообщения: 2700
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: очередные мудрые мысли от VTur

Сообщение Droog_Andrey » Пт мар 15, 2013 8:24 am

sanya1024 писал(а):Мне все-таки так не кажется... Пи орбитали не вычитают эл. плотность на оси. Просто дают нулевой вклад. Так что если там что-то было (а оно было, Вы не отрицаете), никуда оно не денется.
Вы как будто разговариваете с кем-то другим, а не со мной. Максимумы на ядрах, нулевая электронная плотность, вычитание... откуда это всё? Я говорил всего лишь о минимуме.
sanya1024 писал(а):можете по отдельности показать тут (или мне в личку, если не хотите засорять тему) картинки отдельных квадратов 1SGG, 1SGU, 2SGG и 2SGU орбиталей? просто посмотрим, что это за сигмы.
ОК, вечерком, как вернусь с работы.
sanya1024 писал(а):А что будет, если вместо базиса атомно-центрированных функций взять функции двухатомных фрагментов?
Отлично будет работать. Можно, например, взять молекулу азота и для каждого атома указать только одну базисную s-функцию, а всё остальное описать функциями, центрированными на середине связи между атомами.
sanya1024 писал(а):Так, да не так. Визуализировать одноэлектронные функции -- вполне реально, это наглядная конструкция. А вот визуализировать компоненты плотности в конструкции из диатомиков или плоских волн... в общем, если мне кто-нибудь покажет, как это выглядит, скажу спасибо :)
Это можно точно так же визуализировать, и в этом будет примерно столько же смысла :)
sanya1024 писал(а):картинки для этого церия не покажете? что с ним не так? неужели перемешались s,p,d,f функции в одной оболочке?
Не покажу, давно дело было. А функции эти почти всегда перемешиваются :)
2^74207281-1 is prime!

Аватара пользователя
sanya1024
Сообщения: 1672
Зарегистрирован: Чт янв 20, 2011 3:24 pm

Re: очередные мудрые мысли от VTur

Сообщение sanya1024 » Пт мар 15, 2013 3:27 pm

VTur писал(а):Орбитали - это математическая абстракция. Это нормальный термин. Нормальное математическое описание. Они не
являются наблюдаемыми и относится к ним надо, как к математическим абстрактам.
Вот тут 100% согласна, что возводить орбитали в абсолют и строить работу целиком на орбитальном анализе -- не нужно. Просто не нужно так нервно реагировать на то, что кто-то для удобства пользуется орбиталями.
VTur писал(а): Некоторые не знают квантовой механики и смешивают принципиально разные модели в одну кучу.
При этом они еще и не догадываются, что орбитали неинвариантны относительно унитарных преобразований. Инвариатна только сумма их квадратов, поэтому говорить о свойства отдельных орбиталей странно.
И тут кто б спорил, чего нервничать-то. Тем не менее о топологических свойствах отдельных орбиталей говорить можно -- с ними ничего не делается ни при каких преобразованиях, потому что кое-какую полезную информацию из них извлечь можно, и пренебрегать этим не нужно. Но не хотите -- и не надо, Ваши проблемы.
VTur писал(а): sanya1024, я не собираюсь отвечать на Ваши реплики. Вы уже достаточно сказали. Вы используете многоконфигурационные методы, пожалуйста, объясните, что такое электронная плотность многодетерминатной волновой функции.
Диагональная часть матрицы плотности, а что? Вы теперь и многоконфигурационные методы хотите выкинуть на помойку?
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)

Аватара пользователя
sanya1024
Сообщения: 1672
Зарегистрирован: Чт янв 20, 2011 3:24 pm

Re: очередные мудрые мысли от VTur

Сообщение sanya1024 » Пт мар 15, 2013 3:50 pm

Droog_Andrey писал(а):Максимумы на ядрах, нулевая электронная плотность, вычитание... откуда это всё? Я говорил всего лишь о минимуме.
Ну, максимумы на ядрах я на картинке увидела (разве нет?). Вычитаться там ничего не вычитается, я именно это и говорю: при построении эл. плотности все складывается. А минимум -- ну хоть от руки нарисуйте, откуда он там должен взяться.
Droog_Andrey писал(а): ОК, вечерком, как вернусь с работы.
Спасибо, жду.
Droog_Andrey писал(а): Отлично будет работать. Можно, например, взять молекулу азота и для каждого атома указать только одну базисную s-функцию, а всё остальное описать функциями, центрированными на середине связи между атомами.
Да можно... только если взять только функции s симметрии, кое-какие важные эффекты будут потеряны. Именно поэтому в базис добавляют поляризационные функции, к-рые на первый взгляд не нужны -- они помогают описывать отклонения эл. плотности в стороны от линии связи. Наверняка абсолютную энергию молекулы N2 можно описать, набрав в базис достаточное кол-во одних s функций. Но как только Вы захотите описывать взаимодействие этой молекулы с чем-то еще, без p функций не обойтись, а может быть, даже и d. Так что выбор базиса диктуется задачей.
Droog_Andrey писал(а): Это можно точно так же визуализировать, и в этом будет примерно столько же смысла :)
Ну, диатомиками в чистом виде действительно никто не пользуется, а плоскими волнами -- еще как, даже здесь, на форуме есть люди, работающие с плосковолновыми пакетами. Вот только с визуализацией там как-то тяжело.
Droog_Andrey писал(а): Не покажу, давно дело было. А функции эти почти всегда перемешиваются :)
Спин-орбита их перемешала. В тяжелых элементах уже не LS-связь (т.е., сохраняются по отдельности угловой и спиновый момент), а jj (сохраняется суммарный момент), поэтому сферические гармоники запросто могут превратиться в кашу.
А сколько времени займет повторить эксперимент?
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)

Аватара пользователя
sanya1024
Сообщения: 1672
Зарегистрирован: Чт янв 20, 2011 3:24 pm

Re: очередные мудрые мысли от VTur

Сообщение sanya1024 » Пт мар 15, 2013 5:03 pm

Кстати, вдогонку о многоконфигурационных методах. Если это просто КВ без оптимизации самих МО, то какие орбитали взяли, такие и будут. Можно брать хартри-фоковские (обычно так и делают), но можно и кон-шемовские. Хорошо это или плохо -- решается сравнением целевых величин с экспериментом. Если это МКССП, то орбитали оптимизируют вместе с коэффициентами КВ, и тогда без разницы, какие стартовые орбитали брать (разница только в скорости сходимости МКССП) -- все равно к моменту окончания итераций они будут совсем не такими, как исходные. Тем не менее, пи орбиталь никогда не станет сигмой, разве что размазанная по всей ароматике пи съедет на одно кольцо, или некоторые пи и сигмы могут поменяться местами.
Для анализа электронных распределений в многоконфигурационных методах используют натуральные орбитали, получаемые из матрицы плотности. Они четко показывают, куда и откуда переместилась эл. плотность при электронном переходе.
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)

Аватара пользователя
Yu/2
Сообщения: 541
Зарегистрирован: Сб ноя 15, 2003 12:33 pm

Re: очередные мудрые мысли от VTur

Сообщение Yu/2 » Пт мар 15, 2013 6:30 pm

sanya1024 писал(а): Ну, диатомиками в чистом виде действительно никто не пользуется, а плоскими волнами -- еще как, даже здесь, на форуме есть люди, работающие с плосковолновыми пакетами. Вот только с визуализацией там как-то тяжело.
Хм, а что в этом тяжелого? Вот только что визуализировал орбитали из "плосковолнового" DFT.
У вас нет необходимых прав для просмотра вложений в этом сообщении.

Аватара пользователя
sanya1024
Сообщения: 1672
Зарегистрирован: Чт янв 20, 2011 3:24 pm

Re: очередные мудрые мысли от VTur

Сообщение sanya1024 » Пт мар 15, 2013 8:17 pm

О, спасибо! это был именно плосковолновой базис, а не гауссовский с периодическими граничными условиями? А чем считали и чем визуализировали?
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)

Аватара пользователя
Droog_Andrey
Сообщения: 2700
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: очередные мудрые мысли от VTur

Сообщение Droog_Andrey » Пт мар 15, 2013 10:07 pm

sanya1024 писал(а):Вычитаться там ничего не вычитается, я именно это и говорю
А зачем говорить очевидные вещи, не относящиеся к теме? :)
sanya1024 писал(а):ну хоть от руки нарисуйте, откуда он там должен взяться.
Ну есть же узел у пи-орбитали. Значит, при их складывании на оси должен быть минимум. Рисовать влом :)
sanya1024 писал(а):Спасибо, жду.
ОК, скоро сделаю, вот только с работы приполз :)
sanya1024 писал(а):Да можно... только если взять только функции s симметрии, кое-какие важные эффекты будут потеряны.
И снова Вы говорите очевидные истины, которые не относятся к теме. Ясен пень, что только s-функциями адекватно описать ничего не выйдет. Нужно брать p, d, f, g, h, i, k и так далее, чем более "высокочастотные" компоненты распределения плотности хотим учесть - тем дальше.
sanya1024 писал(а):к-рые на первый взгляд не нужны
Это ж какие вещества надо употреблять, чтобы такое почудилось? :)
sanya1024 писал(а):Наверняка абсолютную энергию молекулы N2 можно описать, набрав в базис достаточное кол-во одних s функций.
Нет, нельзя. Но можно обойтись вообще без центрированных на ядрах функций, если повышать угловой момент.
sanya1024 писал(а):Спин-орбита их перемешала. В тяжелых элементах уже не LS-связь (т.е., сохраняются по отдельности угловой и спиновый момент), а jj (сохраняется суммарный момент), поэтому сферические гармоники запросто могут превратиться в кашу.
Это, опять же, с точки зрения модели. И не только в спин-орбите там дело. Реально они перемешиваются везде и всегда, вопрос только, насколько оправдан учёт этого.
sanya1024 писал(а):А сколько времени займет повторить эксперимент?
Честно - не знаю, у меня тут сейчас много задач висит. Попробуйте сами, на bse недавно опубликовали хорошие базисы, вот для того же церия например:

Код: Выделить всё

Ce     0 
S   5   1.00
  66920.6810000              0.0000050        
   7142.4190000              0.0000620        
   1149.2279000              0.0004080        
    626.0474000              0.0000800        
    137.2813000              0.0035590        
S   1   1.00
     36.6434000              1.0000000        
S   1   1.00
     25.9742000              1.0000000        
S   1   1.00
     11.8859000              1.0000000        
S   1   1.00
      3.0284000              1.0000000        
S   1   1.00
      1.5664000              1.0000000        
S   1   1.00
      0.5937000              1.0000000        
S   1   1.00
      0.2630000              1.0000000        
S   1   1.00
      0.0490000              1.0000000        
S   1   1.00
      0.0207000              1.0000000        
P   6   1.00
   3813.8026000              0.0000050        
   1216.9447000              0.0000510        
    496.6641000              0.0000530        
    212.5431000              0.0006710        
     27.6306000              0.0088050        
     19.6040000              0.1508610        
P   1   1.00
     13.8918000              1.0000000        
P   1   1.00
      4.4389000              1.0000000        
P   1   1.00
      2.3374000              1.0000000        
P   1   1.00
      1.1067000              1.0000000        
P   1   1.00
      0.5287000              1.0000000        
P   1   1.00
      0.2305000              1.0000000        
P   1   1.00
      0.0800000              1.0000000        
D   6   1.00
    367.7157000              0.0001200        
    109.8798000              0.0009910        
     36.0211000              0.0077780        
     14.7637000             -0.0629580        
      7.3281000              0.1803420        
      3.9441000              0.4325290        
D   1   1.00
      2.0202000              1.0000000        
D   1   1.00
      0.9649000              1.0000000        
D   1   1.00
      0.3273000              1.0000000        
D   1   1.00
      0.1032000              1.0000000        
F   5   1.00
    123.4821000              0.0015660        
     43.9881000              0.0181010        
     19.4518000              0.0761570        
      8.6013000              0.1926830        
      3.8049000              0.3243320        
F   1   1.00
      1.6176000              1.0000000        
F   1   1.00
      0.6364000              1.0000000        
F   1   1.00
      0.2164000              1.0000000        
G   4   1.00
     19.4518000              0.0021990        
      8.6013000              0.0374280        
      3.8049000              0.0303780        
      1.6176000              0.3556640        
G   1   1.00
      0.6364000              1.0000000        
G   1   1.00
      0.2164000              1.0000000        
****

CE     0
CE-ECP     5     28
h-ul potential
  1
2      1.000000000            0.000000000      
s-ul potential
  1
2     20.137829000          580.083457000      
p-ul potential
  1
2     15.998482000          310.302833000      
d-ul potential
  1
2     14.974187000          167.813944000      
f-ul potential
  1
2     23.402455000          -49.390229000      
g-ul potential
  1
2     16.570553000          -21.331879000      
2^74207281-1 is prime!

Аватара пользователя
Droog_Andrey
Сообщения: 2700
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: очередные мудрые мысли от VTur

Сообщение Droog_Andrey » Пт мар 15, 2013 10:58 pm

sanya1024 писал(а):Если это МКССП, то орбитали оптимизируют вместе с коэффициентами КВ, и тогда без разницы, какие стартовые орбитали брать (разница только в скорости сходимости МКССП)
Не только в скорости. Может сойтись не туда, если исходные орбитали были слишком кривыми. Особенно в случае анионов с диффузными функциями.
2^74207281-1 is prime!

Аватара пользователя
Yu/2
Сообщения: 541
Зарегистрирован: Сб ноя 15, 2003 12:33 pm

Re: очередные мудрые мысли от VTur

Сообщение Yu/2 » Пт мар 15, 2013 11:40 pm

sanya1024 писал(а):О, спасибо! это был именно плосковолновой базис, а не гауссовский с периодическими граничными условиями? А чем считали и чем визуализировали?
Да, это точно посчитано с базисом на плоских волнах в Quantum Espresso. Визуализацию делал в XCrysDen.

Аватара пользователя
Droog_Andrey
Сообщения: 2700
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: очередные мудрые мысли от VTur

Сообщение Droog_Andrey » Пт мар 15, 2013 11:44 pm

Droog_Andrey писал(а):ОК, скоро сделаю, вот только с работы приполз :)
http://www.primefan.ru/stuff/chem/C2_1sg.png
http://www.primefan.ru/stuff/chem/C2_1su.png
http://www.primefan.ru/stuff/chem/C2_2sg.png
http://www.primefan.ru/stuff/chem/C2_2su.png
Как видим, в плоскости, равноудалённой от ядер, максимум на оси очччень пологий, как я и говорил :)

http://www.primefan.ru/stuff/chem/C2_piu.png
А вот у пи-орбитали здесь явный минимум, и он намного резче. Его должно быть видно на фоне сигм, но - увы :)

http://www.primefan.ru/stuff/chem/C2_3sg.png
http://www.primefan.ru/stuff/chem/C2_3su.png
http://www.primefan.ru/stuff/chem/C2_pig.png
Ну и вдогонку НСМО заодно.
2^74207281-1 is prime!

Аватара пользователя
sanya1024
Сообщения: 1672
Зарегистрирован: Чт янв 20, 2011 3:24 pm

Re: очередные мудрые мысли от VTur

Сообщение sanya1024 » Сб мар 16, 2013 4:26 am

Droog_Andrey писал(а):Ну есть же узел у пи-орбитали. Значит, при их складывании на оси должен быть минимум. Рисовать влом :)
И вовсе не обязательно. Разве что максимум сигмы совсем уж пологий. А он, судя по картинке, не такой уж и пологий. По крайней мере мне показалось, что плотность белого тона в центре картинки C2_2sg гораздо больше, чем плотность белого тона в любом месте на картинке C2_piu. Т.е., пи орбитали дают небольшую добавку в общую плотность, а плотность в центре действильно происходит главным образом от довольно резкого (а вовсе не пологого) пика 2s-сигма орбитали. В общем, если попробовать нарисовать графики квадратов орбиталей вдоль оси, проходящей через центр молекулы перпендикулярно C-C связи, наш спор разрешится: пологий максимум 2s-сигма орбитали или не очень, и что будет в сумме, если все квадраты сложить
Droog_Andrey писал(а): Это ж какие вещества надо употреблять, чтобы такое почудилось? :)
А почему тогда, случается, используют базисы без поляризационных функций? Я-то не сомневаюсь, что без них никуда, но случаи разные бывают...
Droog_Andrey писал(а): Нет, нельзя. Но можно обойтись вообще без центрированных на ядрах функций, если повышать угловой момент.
Если повышать угловой момент, а центрировать функции где угодно кроме ядер -- Вы и межмол. взаимодействия вполне прилично сможете описать. Может даже лучше, чем с помощью функций, центрированных только на ядрах.
Droog_Andrey писал(а): Это, опять же, с точки зрения модели. И не только в спин-орбите там дело. Реально они перемешиваются везде и всегда, вопрос только, насколько оправдан учёт этого.
Модельность тут состоит только в том, пренебрегаем мы спин-орбитой или нет и в какой степени мы ею пренебрегаем. Это когда речь идет о свободном атоме. У отдельных атомов в молекуле, конечно, угловой момент не сохраняется, но сохранение углового момента в молекуле как целом как раз и приводит к разделению орбиталей по узловой структуре. Ну, или по симметрии -- если повезет.
Droog_Andrey писал(а): Честно - не знаю, у меня тут сейчас много задач висит. Попробуйте сами, на bse недавно опубликовали хорошие базисы, вот для того же церия например:
Да у меня тоже сейчас работы выше ушей и все машины заняты. Может, руки дойдут -- погляжу. Но чего мелочиться, надо тогда уж брать полноэлектронный базис, и гамильтониан насколько возможно релятивистский. Я подумаю.
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)

Аватара пользователя
sanya1024
Сообщения: 1672
Зарегистрирован: Чт янв 20, 2011 3:24 pm

Re: очередные мудрые мысли от VTur

Сообщение sanya1024 » Сб мар 16, 2013 4:38 am

Yu/2 писал(а): Да, это точно посчитано с базисом на плоских волнах в Quantum Espresso. Визуализацию делал в XCrysDen.
Ага, спасибо. Оказывается, и от такой флудильной темы есть польза :)
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)

Аватара пользователя
Droog_Andrey
Сообщения: 2700
Зарегистрирован: Сб сен 29, 2007 8:29 pm
Контактная информация:

Re: очередные мудрые мысли от VTur

Сообщение Droog_Andrey » Сб мар 16, 2013 9:14 am

sanya1024 писал(а):А почему тогда, случается, используют базисы без поляризационных функций?
Потому что в простейших случаях и для очень грубого приближения этого иногда оказывается достаточно :)
sanya1024 писал(а):надо тогда уж брать полноэлектронный базис, и гамильтониан насколько возможно релятивистский. Я подумаю.
Ну тогда ANO-RCC только.
2^74207281-1 is prime!

VTur
Сообщения: 7357
Зарегистрирован: Пт авг 31, 2007 1:36 pm

Re: очередные мудрые мысли от VTur

Сообщение VTur » Сб мар 16, 2013 10:02 am

sanya1024 писал(а):
VTur писал(а): Вы используете многоконфигурационные методы, пожалуйста, объясните, что такое электронная плотность многодетерминатной волновой функции.
Диагональная часть матрицы плотности, а что?
Скажите, а к какому детерминанту относится используемая Вами диагональная часть матрицы плотности или эл. плотность в многоконфигурационном методе? Для все сразу? Для смеси всех состояний? Или только первый детерминант?
sanya1024 писал(а):Вы теперь и многоконфигурационные методы хотите выкинуть на помойку?
Мнгоконфигурационные методы я использую для поиска строения, энергии, частот. Для анализа эл. плотности я использую DFT.
После отстоя требуйте долива

Аватара пользователя
sanya1024
Сообщения: 1672
Зарегистрирован: Чт янв 20, 2011 3:24 pm

Re: очередные мудрые мысли от VTur

Сообщение sanya1024 » Сб мар 16, 2013 2:22 pm

VTur писал(а):Скажите, а к какому детерминанту относится используемая Вами диагональная часть матрицы плотности или эл. плотность в многоконфигурационном методе? Для все сразу? Для смеси всех состояний? Или только первый детерминант?
Либо к конкретному многодетерминантному состоянию (а не детерминанту), либо -- в state-averaged случае -- усредняется для нескольких состояний.
VTur писал(а):Мнгоконфигурационные методы я использую для поиска строения, энергии, частот. Для анализа эл. плотности я использую DFT.
А Вы не обращали внимания, что иногда результаты DFT и CASSCF расчета оказываются качественно разными?
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)

Аватара пользователя
sanya1024
Сообщения: 1672
Зарегистрирован: Чт янв 20, 2011 3:24 pm

Re: очередные мудрые мысли от VTur

Сообщение sanya1024 » Сб мар 16, 2013 3:24 pm

Droog_Andrey писал(а):Ну тогда ANO-RCC только.
Я подумаю.
Я считала ион Eu3+ в CASSCF/XMCQDPT (FireFly) с обычным гамильтонианом с последующим спин-орбитальным расчетом в GAMESS-US. Но у меня был базис со Штуттгартским small-core псевдопотенциалом, а орбитали визуализировались только на CASSCF стадии -- и все было вполне штатно, как в учебнике, хотя при генерации стартовых орбиталей в ХФ незанятая f орбиталь уехала фиг знает куда, и пришлось ее искать, чтобы вернуть в акт. пространство. На стадии учета спин-орбиты делалось спин-орбитальное КВ на CAS-овских орбиталях, так что увидеть влияние спин-орбиты на вид АО не получилось. А вот в полноэлектронном базисе в Природе с квазирелятивистским гамильтонианом все орбитали перемешались в кашу: вперемешку по энергии шли f, d и s (p вроде не попали). Но все-таки f по форме оставались f, d -- d, а s -- s.
Может быть, ион Eu3+ сильно отличается от нейтрального атома Ce в сторону облегчения расчетов? Церий когда-то давно пыталась считать Природой, но быстро забросила, получив нечто неинтерпретируемое, а вернуться -- не сложилось.
А вообще, наверное, тут нужен бескомпромиссный DIRAC, с к-рым я не умею обращаться.
Вот и вся моя работа. Стеречь ребят над пропастью во ржи. (Дж. Д. Сэлинджер)

VTur
Сообщения: 7357
Зарегистрирован: Пт авг 31, 2007 1:36 pm

Re: очередные мудрые мысли от VTur

Сообщение VTur » Сб мар 16, 2013 6:54 pm

sanya1024 писал(а):
VTur писал(а):Скажите, а к какому детерминанту относится используемая Вами диагональная часть матрицы плотности или эл. плотность в многоконфигурационном методе? Для все сразу? Для смеси всех состояний? Или только первый детерминант?
Либо к конкретному многодетерминантному состоянию (а не детерминанту), либо -- в state-averaged случае -- усредняется для нескольких состояний.
1. То есть Вас не смущает, что в "электронную плотность" основного состояния входят орбитали возбужденных состояний? И это уже друга плотность.
2. Но тогда, как Вы находите эл. плотность конкретного возбужденного состояния - это же несколько детерминантов одной спиновой конфигурации.
sanya1024 писал(а):
VTur писал(а):Мнгоконфигурационные методы я использую для поиска строения, энергии, частот. Для анализа эл. плотности я использую DFT.
А Вы не обращали внимания, что иногда результаты DFT и CASSCF расчета оказываются качественно разными?
Что имеется в виду?
После отстоя требуйте долива

Ответить

Вернуться в «квантовая химия и моделирование»

Кто сейчас на конференции

Сейчас этот форум просматривают: нет зарегистрированных пользователей и 12 гостей