Вопрос по выходным файлам Gaussian.

[Дмитрий Викторович]

Уважаемые коллеги прощу вашей помощи. Меня интересует информация о так называемых Archive в выходных файлах Gaussian. Вот нашел на одном сайте:

At the very end of each Gaussian calculation, an archive entry in a very compact format is printed to summarize the results. This archive entry is frequently used as supplemental material in publications of theoretical results.
(Enter /scr1/g03/l9999.exe)
1\1\GINC-TERMINUS\SP\RHF\6-31G(d)\C1H2O1\ZIPSE\18-Oct-2004\0\\#P HF/6-
31G(D) SCF=TIGHT\\test1 HF/6-31G(d) sp formaldehyde\\0,1\C\O,1,1.2\H,1
,1.,2,120.\H,1,1.,2,120.,3,180.,0\\Version=x86-Linux-G03RevB.03\State=
1-A1\HF=-113.8529674\RMSD=2.649e-09\Dipole=0.,0.,-1.1094497\PG=C02V [C
2(C1O1),SGV(H2)]\\@

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

[VTur]

Лично мне необходимо знать каким образом представляется в этом разделе информация о производных дипольного момента и тензорах поляризуемости после расчетов молекулярных колебаний.

По координатам атомов? Имеются в виду GAPT заряды?

[Дмитрий Викторович]

По координатам атомов?

Судя по всему именно так. Вот как это выглядит после расчетов колебаний молекулы аммиака:

1|1|UNPC-UNK|Freq|RB3LYP|6-311++G(3df,3pd)|H3N1|PCUSER|28-Apr-2009|0||
#P GEOM=ALLCHECK GUESS=READ SCRF=CHECK GENCHK RB3LYP/6-311++G(3DF,3PD)
FREQ||Title Card Required||0,1|N,0.,0.,-0.1123097988|H,0.9417386299,0
.,0.2620561972|H,-0.470869315,-0.8155695772,0.2620561972|H,-0.47086931
5,0.8155695772,0.2620561972||Version=x86-Win32-G03RevB.01|State=1-A1|H
F=-56.5875478|RMSD=1.219e-009|RMSF=3.483e-007|Dipole=0.,0.,0.5959892|D
ipoleDeriv=-0.3622183,0.,-0.000002,0.,-0.3622171,0.,-0.0000019,0.,-0.5
496972,0.0892838,0.,-0.1089679,0.,0.1521944,0.,-0.1272669,0.,0.1832317
,0.136467,-0.027241,0.0544844,-0.0272419,0.1050114,0.09437,0.0636339,0
.110217,0.1832322,0.136467,0.027241,0.0544844,0.0272419,0.1050114,-0.0
9437,0.0636339,-0.110217,0.1832322|Polar=13.3354812,0.,13.33548,0.0001
727,0.,15.1035902|PolarDeriv=-3.2618769,0.,3.2610298,-3.7759934,0.,0.0
000463,0.,3.2608872,0.,0.,-3.7759486,0.,-4.7676436,0.,-4.7676176,-0.00
00682,0.,-0.0306124,7.1751718,0.,1.0838799,2.3983612,0.,2.5136819,0.,0
.871377,0.,0.,0.1189224,0.,2.4531254,0.,0.7249534,0.3388392,0.,0.01012
9,-1.9566475,-1.6959669,-2.1724549,0.6888161,0.9868639,-1.2568641,-1.8
799441,-2.0661321,-5.2721347,0.9869252,1.8285131,-2.1769951,1.1572591,
0.7483085,2.0213321,-0.1693855,-0.2934533,0.0102417,-1.9566475,1.69596
69,-2.1724549,0.6888161,-0.9868639,-1.2568641,1.8799441,-2.0661321,5.2
721347,-0.9869252,1.8285131,2.1769951,1.1572591,-0.7483085,2.0213321,-
0.1693855,0.2934533,0.0102417|HyperPolar=12.2234683,0.,-12.2162348,0.,
8.6212597,0.,8.621742,-0.0048827,0.,28.3703796|PG=C03V [C3(N1),3SGV(H1
)]|NImag=0||0.64191649,0.,0.64191673,0.00000006,0.,0.20858588,-0.36824
441,0.,-0.11692517,0.40226785,0.,-0.05969993,0.,0.,0.05815830,-0.17011
789,0.,-0.06952867,0.13656538,0.,0.06826099,-0.13683604,-0.13360373,0.
05846257,-0.01701174,-0.03313270,0.01677623,0.14418568,-0.13360365,-0.
29110835,0.10126017,0.00233237,0.00077084,0.00165354,0.14900381,0.3162
4046,0.08505893,0.14732644,-0.06952864,-0.00982012,-0.01370188,0.00063
384,-0.06828269,-0.11826909,0.06826099,-0.13683604,0.13360373,0.058462
57,-0.01701174,0.03313270,0.01677623,0.00966213,-0.01773254,-0.0069561
1,0.14418568,0.13360365,-0.29110835,-0.10126017,-0.00233237,0.00077084
,-0.00165354,0.01773254,-0.02590303,-0.01535541,-0.14900381,0.31624046
,0.08505893,-0.14732644,-0.06952864,-0.00982012,0.01370188,0.00063384,
-0.00695611,0.01535541,0.00063384,-0.06828269,0.11826909,0.06826099||0
.,0.,0.00000067,-0.00000053,0.,-0.00000022,0.00000027,0.00000046,-0.00
000022,0.00000027,-0.00000046,-0.00000022|||@

Соответствующие фрагменты выделил жирным шрифтом. Меня интересуют размерности этих величин и порядок их следования. В справке по Gaussian наткнулся на упоминание "Browse Quantum Chemistry Database System User’s Guide" для желающих получить более детальную информацию. Дальнейшие мои поиски успехом не увенчались.

[VTur]

Если честно, разбираться не очень хочется. По GAPT зарядам и этим производным я ссылки уже приводил viewtopic.php?f=11&t=24079&p=186270&hilit=GAPT#p186270

[Дмитрий Викторович]

И все таки может быть кто-нибудь знает ответы на мои вопросы?

[Yu/2]

И все таки может быть кто-нибудь знает ответы на мои вопросы? :(

По поводу производных дипольного момента:

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

http://www.ccl.net/chemistry/resources/messages/2006/01/20.005-dir/index.html

[Дмитрий Викторович]

По поводу производных дипольного момента:

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

http://www.ccl.net/chemistry/resources/messages/2006/01/20.005-dir/index.html

Спасибо, попробую. Судя по всему, не только у меня есть вопросы.

[VTur]

К тому же, непонятно, как их использовать. Компоненты этих тезоров есть относительные величины. Они зависят от вида используемых координат (декартовы, внутренние), положения начала координат и ориентации системы координат. При небольшом повороте их значения меняются. Скорее всего нужны собственные числа или след, или 1/3 следа - GAPT заряды.