6月19日(土)のつぶやき

09:02 from web
@statabo SWoPPやっぱりでようと思ってます。予定はまだつまってないはず
09:03 from web
やったあああCCSD(T)に勝ってた!!!
09:20 from web
@yasuhito 飛行機で寝てたら、隣に座ってたおばあちゃんがいきなり手のひらに飴を押し付けてきたことがあります。お金いらなかったら僕にください。
09:23 from web
@sdpaninf 内点法ってたしかXY=\mu I といてますよね。\muのさがりかたは普通ですよ。
14:15 from web
@yasuhito いっぱいもらってきてください。おすそ分けお願いします。中学の頃近鉄奈良線で、おばあさんに「あんた大丈夫か」といわれて席を譲ってもらったことがあります。
14:16 from web
@ tota_rtfm FLOSS関連企業に転職しないの?
18:23 from web
宮崎のマンゴーおいしいよね! 実家に帰ったときちょっとだけ食わせてもらった。@shouts_out_to ここ数日、マンゴー食べたくて。
18:26 from web
@anti_tank_rifle ありがとうございます。というか君の立候補も待ってたけどね。本なんかの翻訳、コーディネーターに了解とってやってみるのもいいかも。出版もできるかもね。
18:29 from web
これってPRで返事すべき問題では?? RT @tota_rtfm ports/139448 は 8.1-RELEASE 後に対応します。
19:43 from web
@anti_tank_rifle こんど京都いくから会いましょう。
by NakataMaho on Twitter

6月18日(金)のつぶやき

07:58 from web
寺田さんja.oo.o top pageのfunding requestの件。一応問題意識を投げられたように思える。いけそうな気もするがわからん。スレッドhttp://tinyurl.com/27tz5da #oooja
09:02 from web
能力のある人はたくさんいる。色々やってもらって、出来もよい。指導力を持った人を見つける、となるとこれはちょっとやそっとではうまくゆかない。#oooja
09:13 from web
ドイツの教育の状況 1/6がいわゆるNEET。社会が高度化すると、こんなふうになるのかね。皆が労働にむいているとは思えん。http://tinyurl.com/28qvq2l
10:00 from web
DFTも波動関数使わんよね...
10:01 from web
@anti_tank_rifle みんながそんなにあくせく働かなくていいんじゃない? あとみんなが大学にいかなくてもいいんじゃない? とかは思うけどなぁ。
10:01 from web
縮約密度行列でgoogleかけたら2chのスレッドが出てきて吹いた。
10:02 from web
ALDETを使うとFullCIができそう。ただ、FreeBSDでは動かずLinuxでは動いた。微妙な計算方法の違いとかあるんだろうなぁ...
10:06 from web
@anti_tank_rifle ギムナジウムとハウプトシューレだったっけ。大学なんて専門的過ぎてあまり一般の人が学ものではないし、中学生くらいで、ある程度興味が分かれる気もする。
11:29 from web
くそーCCSD(T)に負けた!!!
13:17 from web
right. jaでは大っぴらに言わないけど、自分の中でそうやってます。@anti_tank_rifle 「達成可能な小さな目標をいくつも立てよう」
13:18 from web
@anti_tank_rifle OOoの戦略考えてみぃひんか?
13:41 from web
だよなぁ...@anti_tank_rifle OOoの仮想敵にMSOを考えると、... 細かい地道な成果を重ねることが大切だと思います。
13:41 from web
わかるけど、これ誰がどうやってやるかが問題。@anti_tank_rifle jaオリジナルの特製テンプレート的なものでも充分だと思います。
14:01 from web
@anti_tank_rifle http://templates.services.openoffice.org/ja
みたいなのはある。
16:35 from web
document projectで、書いてみたら? 意外と皆さんノウハウに飢えてます。どういう使い方しているのかは個人的にもとても興味あります。@anti_tank_rifle 「Draw＋Writer＋Calcの連携等」
16:46 from web
@anti_tank_rifle まぁ、今日、OOoのために何ができたか? こうtweetsしてても、ほとんど進まないよね...こちらはDEV300_m83チェックアウトしたぞ #oooja
17:14 from web
@anti_tank_rifle ちょっとしたことで、毎日のようにできる、メンテナンス的なことがいいかもね。一ヶ月離れると頭が忘れるんで。
18:12 from web
@anti_tank_rifle 月曜までになんか一つ超簡単なドキュメント[ja-doc]に投稿してください。
by NakataMaho on Twitter

縮約密度行列の変分法(IV): おそらく唯一の多項式程度の電子構造計算の方法

縮約密度行列の変分法は、おそらく唯一の多項式程度の電子構造計算の方法だろう。SCFは? 実はglobal minimumを得るのは多項式程度ではできない。

我々の方法は、一粒子基底の選び方に依存しない。つまり直交化したatomic orbitalを用いてもまったく同一の結果を得る。ここらへんはfullCIと同じである。

なぜ多項式程度のリソースで計算できるか。これは凸計画なので最適化としてはやりやすい。もう少しわかりやすく言えば、円を変分領域とし、線形関数を最小化するというのを考えよう。この場合は、円の端っこで最小値となる。

半正定値計画法は、行列を半正定値に保ちながら変分するというのではあるが、意外にも計算負荷が軽い。最小値は数学的に検証できる(primal or dual feasibility, gapなどから)。

Post Hartree-Fockはすべて難しい。Hatree-Fockの難しさを引きずっているからだ。DFTはいうにおよばず。

他に多項式程度のアルゴリズムで解ける、電子相関理論にはどんなものがあろうか。私は知らない。

縮約密度行列の変分法(III): 大域的な最小値がもとまり局所解が存在しない

縮約密度行列の変分法で、N-representability条件を、必要条件のみとし、それらが凸であれば、大域的な最小解(global optimum)が求まる。

この条件を満たすのは比較的簡単である。P, Q, G, T1, T2' などはすべて必要条件かつ、凸であることがすぐわかる。

このglobal optimumが求まる理由は、線形計画の双対定理を理解していればなんとなく解るかもしれない。PQGT1T2'などの必要上件を持ってくれば、半正定値計画問題となるが、これは双対定理、相補性定理などから、最小解では、gapがゼロとなる、ということからその大域性が解るのである。ただ、我々の定式化では、内点がないため必ずしも主双対内点法の枠組みの範囲に収まらず、計算上のgapが存在するかもしれないが、少なくともgapがゼロの場合は正しく、物理的にもgapはゼロであることがわかる。つまり安心してよいということである。

SCFのように、求まった解がglobalかsub optimalか、悩む必要はない。縮約密度行列の変分法ならば、主、双対のgapを見れば自ずと大域的な最小化と解る。

6月17日(木)のつぶやき

07:35 from web
O2+の8-CIが300Gのdiskを使ってコケてた。もうやる気なし。
07:44 from web
もうこの大きさだとfullCIは出来んなぁ。ここから先は超絶技巧? どうするかなぁ。CCSD(T)かなあ。DMRGかなぁ。いずれにせよ厳密な結果の限界が見えた。
07:48 from web
@BuchikoYossyTK お元気ですか?
07:59 from web
@shigetay いや～時間とHDDかけるのもいいですが、たった20軌道です。DMRGならもう少し大きくできると思いますが正しくもとまったかの判定は非常に困難であります。
08:00 from web
どなたか密度行列の研究やってくれないかなぁ...
08:20 from web
@shigetay 密度行列の研究やりません?
08:22 from web
今時のCASってどれくらいとるのが普通なんかな? 20軌道20電子ってできるの?
09:02 from web
@shigetay frozen coreずっとまえにやりました...さらっと
09:44 from web
@shigetay やってないですね～。
09:45 from web
@shigetay やっと最近18軌道できるようになりました。
17:56 from web
@anti_tank_rifle 聡明なところと、tweetがマニアックなことかな?
19:39 from web
gamessが手持ちのFreeBSD 8/amd64でうごいた。他のFreeBSD8/amd64では動くのだが...原因はkern.ipc.shmemが尽きていたことらしい。うーん、なかなかわからんのでpkg-messageにもいれとくか.
by NakataMaho on Twitter