24062009 q-chem
面白そう。
読んでみよう
24062009 q-chem
複合体結晶構造から重要な残基とのインタラクションを見積もるのに、リガンドと特定の残基のシンプルなモデルを切り出してきて、きちんと計算すればよかろうと。
とか書いておいて、FMOの結果から相互作用の強いアミノ酸残基とリガンドを切り出してきて、morokuma-analysisしてみたけど、ちっとも収束しなかった。
うーん、これは困ったと、切り出したモデルでオプティマイズかけて再度計算してみたけどやっぱり収束しなかった。
結晶構造の座標みたいにあまり安定でないデータからエナジー出して、その成分をELECTROSTATIC,EXCHANGE REPULSION,POLARIZATION,CHARGE TRANSFERに分解する技術はきちんと押さえておく必要があるんだけど、どうやればいいんだろうか?
やっぱPIEDAの論文読んでおく必要ある?
43 Pair interaction energy decomposition analysis, D. G. Fedorov, K. Kitaura, J. Comp. Chem. 28 (2007) 222-237.
23062009 life
娘がDSで遊びたがるので買ってみた。
本当は中古で1000円くらいになっているのを期待していたのだけど、値崩れしないのね。娘は楽しそうに遊んでいるので良しとしたい。
あと、ちょっと触ってみた感じでは、ひらがな書き取り判定は結構シビアな気がする。僕が書き取りしても、書き順があってても、気を抜いて線からはみ出て「ぶーーーー」とか言われてちょっとむかついた。
【タイトルデザイン】修悦文字:東京新宿駅や日暮里駅などでガムテープ案内人として話題の佐藤修悦氏による制作。今回、はじめてこども向けの商品タイトルを手がける。
正直ガムテフォントが幼児の脳トレにとってどういう意味があるのかわからん。
MOROKUMA ANALYSISのサンプル(EXAM28)を。
水とアンモニアのダイマー計算
$contrl scftyp=rhf runtyp=morokuma coord=zmt $end
$basis gbasis=n31 ngauss=4 $end
$guess guess=huckel $end
$morokm iatm(1)=3 $end
$data
water-ammonia dimer
Cs
H
O 1 rOH
H 2 rOH 1 aHOH
N 2 R 1 aHOH 3 0.0
H 4 rNH 3 aHNaxis 1 180.0
H 4 rNH 3 aHNaxis 5 +120.0
H 4 rNH 3 aHNaxis 5 -120.0
rOH=0.956
aHOH=105.2
rNH=1.0124
aHNaxis=112.1451 ! makes HNH=106.67
R=2.93
$end
結果
HARTREE KCAL/MOLE
ELECTROSTATIC ENERGY ES= -0.022336 -14.02
EXCHANGE REPULSION ENERGY EX= 0.014308 8.98
POLARIZATION ENERGY PL= -0.001786 -1.12
CHARGE TRANSFER ENERGY CT= -0.003780 -2.37
HIGH ORDER COUPLING ENERGY MIX= -0.000687 -0.43
TOTAL INTERACTION ENERGY, DELTA-E= -0.014282 -8.96
DECOMPOSITION OF CT
CHARGE TRANSFER ENERGY, MON= 1 CT= -0.000494 -0.31
CHARGE TRANSFER ENERGY, MON= 2 CT= -0.003286 -2.06
DECOMPOSITION OF PL
EPL, MON= 1 PL= -0.001068 -0.67
EPL, MON= 2 PL= -0.000584 -0.37
HIGH ORDER COUPLING FOR PL, PMIX= -0.000133 -0.08
複合体結晶構造から重要な残基とのインタラクションを見積もるのに、リガンドと特定の 残基のシンプルなモデルを切り出してきて、きちんと計算すればよかろうと。
21062009 life
シールでチューンされた。

ボードの板みたいだ。
Apple MacBook 2.0GHz Core 2 Duo/13.3"/2G/160G/8xSuperDrive DL/Gigabit/802.11n/BT/Mini DisplayPort MB466J/A21062009 Python
整数で割る際にfloatの数値が欲しいときに2.0と書いてたけど0は必要ないらしい
>>> 1/2
0
>>> 1/2.
0.5
>>> 1./2
0.5
でもちょっと読みにくいので、今まで通りきちんと書くかな。
18062009 life
夏までにどれか一冊消化したい。
Programming Clojure (Pragmatic Programmers)