YMO

GNU/Linux 環境

2009-11-21T01:42:00.005+09:00

現状で使いこなせる（と思い込んでいる）ツール・言語・ライブラリetc

C++（0x対応コンパイラを待ちわびている信者なので）
関連諸ツール（ platex とか dvipdfmx とか、 Word も PowerPoint も、ましてや OpenOffice なんて使う気がしないので必須）
gnuplot （usingとかsetを使いこなせるようになると飛躍的に便利になる）

そこそこ使えるが、まだまだ勉強が足りないツール・言語・ライブラリetc

AWK （たぶん文法はほとんど知ってるけど毎回忘れる、これさえあれば数値実験のデータ解析は大体用が足りて、 Excel はほとんど不要）
Ruby （まだまだ勉強中）
PHP （研究には絶対使わないし、大規模システムは絶対 PHP では組みたくないけど、このサイトとかちょっとしたウェブアプリには便利）
XHTML （文法はほとんど忘れてる、別に勉強したところで得るものは W3C に準拠しているという満足感だけだろうけど）
CSS （同上）
boost （C++信者には欠かせないライブラリ）
make （これも研究で必要になることはないだろうけど）
gdb （効率的に使えるようになるとデバッグが進んでうれしいな）
ssh （自分の知らない機能が大量にあるはず）
apache （研究で Web サーバ立てることなどまずない）
Coreutils （コマンド使いこなせるようになるとちょっとスキルアップした気がする）
Emacs （完璧に gEdit っぽいキーバインドになるのなら便利なエディタだけど）
vim （こっちはキーバインドどころではないが、使いこなせたら便利そう）
zsh （「漢のシェル」だが、機能が多すぎてぜんぜん全体像までたどり着けない）

ほとんど使えないが勉強してみたいツール・言語・ライブラリetc

R （統計関数が山ほど入っているのがうれしい、どうでもいいが一語の命名は検索にかからないからやめて）
subversion （数値実験のソースをバージョン管理した方がいい気がしてきた）
gprof （性能解析も今後必要になるだろうな）
JavaScript （その気になればすぐ習得できるだろうけど面倒なので労力を割かない）
SWIG （ Ruby と C++ を一緒に使えるのはかなり魅力的）
autoconf （まあ、研究で使う機会は皆無だろうね）
bison, flex （同上）
as （コンピュータを使うものの教養程度には……　アセンブラ使った方が早くなる計算処理とかあるのかなあ）
Postfix （まあ、サーバ立てる機会があれば）
GSL （数値計算やるものとしては魅力的だが如何せん言語が C だし）
GTK++ （ X プログラミングで、計算結果を画像で見れるとうれしい）
OpenGL （3次元の数値計算では、やっぱりこういうので絵で見られればいいな）
ImageMagick （ convertしか使えません）
ghostscript （ TeX やってると、関連ソフトにはお世話にはなるけど）
XSL （昔真面目に使ったことはあるけど、もはや完全に忘れてしまった）

そこそこユーザはいるものの今更使う気にならないツール・言語

FORTRAN （この業界未だにいるよね、自分で使ってる分にはいいけど何も学生に教えなくてもいいのにと思う）
C （例えクラスもテンプレートも使わなくても、 C++ を better C として使った方が圧倒的によい）
Java （中間言語って、こういうことやってるとあらゆる意味で中途半端）
Perl （読みにくさは犯罪的、 zsh 以上 ruby 以下の作業があればもしかしたら使うかもしれないけど）
Python （一面的な価値観の押し付けは嫌い、そういう意味で Ruby も同等だけど、あっちの方がだいぶマシ）
bash その他シェル（ zsh に勝てるところってあったっけ？　軽さ？）

なんか XHTML とか JavaScropt とかが場違いな気がするけど、まあいいや。

卒研進捗状況

2009-11-21T00:21:00.005+09:00

相変わらず卒研は遅々として進まない。化学反応は相変わらず基本的なことが分かっていないし、超流動（・超伝導）に関してはもはや忘れ去られてしまった。並行して Goldenfeld を読み進めていくという話が、未だに Introduction すら読み終わっておらず、研究に費やす集中が切れたところで思い出したように今日開いてみると、そもそもどこまで読み終わったかを忘れているくらいなので前の節に何が書いてあるかはとうに記憶の彼方（そもそも、集中が切れた状態で開いたので英語も数式も全然頭に入ってこずに、読むのは結局諦める）。

ただ、しばしば人から聞くように研究が辛いかというと、目下のところ、そのようなことは全くない。今のところは特に期限を意識していないというのが非常に大きく働いていて、「○日までにとにかく結果を出さなければいけない」というのが差し当たってないので（もちろん最終的には卒研発表があるのだけれど）、純粋に趣味のようなものとして、上手くいかない試行錯誤の過程も含めた上で研究を楽しめている（これには、師匠がいつぞやか言っていた「研究は上手くいかないのがデフォルト」という言葉も大きい）。ただ、裏を返せば、差し当たっては修論の提出期限が迫って、研究成果を近いうちに出せなければ自分の人生設計が大きく狂うとかいった状況に追い込まれれば非常に苦しくなるだろうが。

そういう、長期的な期限に追われているわけではない状況において、毎週1回の師匠との面会はペースを保つ上で非常に有難い。さすがに貴重な時間を割いてもらってアポをとっている以上、単なる個人的な趣味のように「今週は何もやってません、すいません」というわけにもいかないので、結果が出ようと出まいと、とにかく試行錯誤をして、うまくいったならいったでめでたいし、詰まったなら詰まった箇所を説明できる程度には手を動かす。将来的にはこの定期面談に甘え切っているようではいけないが、とりあえず今は研究の右も左も分からないので、ちょっとくらいはペースメーカーに頼ってもよいだろう。

とか、そういう内情を書き散らかそうとブログを作ったつもりだったが、研究しているうちに気になったことを書いているうちに次第にアクセスが増えてきた。「熱力学ポテンシャル」だとか「等重率の原理」だとかで検索するとこのブログがかなり上位にヒットしてしまうが、そんなのでいいのかなあ。一応、度を越して恥ずかしいこと（えらく見当違いの記事だったり、身の程を越えて背伸びした記事だったり）は書いてないつもりだけど。まあ、ちょっと検索すると出てくる今の感じだと、例えば研究室の人なんかには読まれてると思った方がいい（というか少なくとも田崎さんはコメント残していったし）。まあ、このブログは研究以上に義務感はないので、ぼちぼちやっていきます。

(int)

2009-11-19T23:59:00.005+09:00

例えば、1/2が0として（意に反して）認識されるというバグはよく埋め込まれる。今日ちょうどこれで数週間悩んでいたのが発覚して、あまりの呆気なさに愕然とした。指摘されればすぐ気づくバグではあるが、なかなか自分では気づかないのが悲しい。

ところが、これも同じく今日気づいたのだが（これで悩んだのは数十分）、こちらはどう考えても言語設計が悪いだろうという罠。


container<double> x; // 0 <= x[i] < 1
std::cout << std::accumulate(
 x.begin(),
 x.end(),
 0,
 std::plus<double>()
) << std::endl;
// prints 0.

期待した結果が欲しければ、以下のようにせねばならない。


std::accumulate(
 x.begin(),
 x.end(),
 0.0,
 std::plus<double>()
)

そう、std::accumulateは第3引数で型を推論しているので、第4引数にどんな関数オブジェクトを指定しても出力は第3引数の型となる。これは非常に厄介だ。まあ、実装する側になって考えれば、この罠を回避するうまい術は思いつかないのだが。

科学技術関係予算の確実な確保について

2009-11-19T23:46:00.003+09:00

例の記事で紹介した以外のところでも署名は動いているらしく、その話がうちの研究室の先輩に届いた所なのだが、科学者から反対意見が相次いでも民主党は予算削減を断行できるのだろうか。

科学技術関係予算の確実な確保について（緊急提言）

先日の事業仕分けに対して、このような文書が総合科学技術政策会議の有識者議員からさえも出ている。東大のような組織がこのような声明を出したら民主党は打撃になると思っていたが、まさか政府の公的会議から来るとは。

しかし、民主党もこうなると、撤回即ち自分の非を認めることになるので、却って泥沼化しそうな気もしなくはない。だが、まあ少なくとも現段階で科学者側からの「脅し」をかけておかないと、今後政権交代が頻繁に起こるような世の中において、それとは無関係に科学技術予算を恒常的に確保するのは難しいだろう。

それにしても、幸か不幸か、今まで、科学技術に伴う利権を得る代わりに科学技術の予算を守ってくれるような「科学技術族」の議員はあまりいないように思う。科学者というのは概して面倒な政治を嫌うものだから、そこのところの発信が抜けてきたツケが回ってきたというところもあるだろうが。

科学研究費補助金の一部の執行停止に対する反対署名

2009-11-17T00:26:00.007+09:00

民主党の政策で、予算の配分を色々変えようといった問題がある。筆者は確かに、子供手当てや後期中等・高等教育に関する負担の軽減には賛成である（高速道路無料化はメリットが感じられず反対だが）。少子高齢化が進むと、福祉政策が成立しないわけで、そこから逆算すると、子供を産むインセンティヴ、産んだところできちんと育てられる保証が国から出れば、ちょっとは子供の数も増えるはず（ついでに、シングルマザーをもっと保護するだとか、特別養子制度を推進するだとか、もっとすればいいと思うけど）。

だがしかし、民主党の科学技術政策には、基礎科学の世界に生きる者としては決して賛同しがたいものがある。

科学研究費補助金の一部の執行停止に対する反対署名

基礎科学の研究は広く人類の文化の問題であり、経済的に余裕のある国は、そこに相応の負担をする義務があると考える。加工貿易云々という論点はいささか狭いと感じるが、大きな趣旨に賛同するので署名する。（田崎晴明氏）

呼びかけ人や、すでに署名をされた方々への応援の意味をこめて、この署名をさせていただきます。呼びかけ文のすべてにわたって賛同するわけではありませんが、このような声を挙げること自体は意味のある行動だと思います。私自身の意見は、以下のweb siteに載せていますので、一読いただければ幸いです。　吉川研一（京大理・物理・教授）
http://www.sci.kyoto-u.ac.jp/modules/tinycontent1/index.php?id=5

自分の分野に近いところだけでも、このような方達が署名に賛同している。自分も、科学の世界に生きるものとしてはまだまだひよっこの中のひよっこだが、何らかの意味を持てばと考え署名してきた。

自分自身の意見としては、田崎さんの意見に近いものがある。というか、基本的に理学にいるものはそういう実感ではなかろうか。自分の研究成果はまず世界中の同僚に共有されるわけで、基礎科学が振興したところで基本的にはある種排他的な「科学技術立国」という枠組みには結びつかない。だからこそ、予算が付きにくいというのもあるのだが、そこに予算をつけるのが先進国の役割であると考える。あと、そういう高尚な動機に加え、若手研究者に対する予算が減るというのはそのまま自分の経済状況に跳ね返るため、黙してはいられないというのもあるが。

ともあれ、このサイトのアクセス数も幾分増えてきたことであるので、ここで署名の広報をすることは一定の意味があると考え、掲載する。

追記。直筆署名に移行するらしい。しかし、当分ネット署名は残っていそうなので、告知は残しておきます。自分自身は直筆署名をする機会があればしたいとは思うが。

追記。署名が一旦打ちきられたようです。トップのリンクも消しておきます。何か続報があればまた。

事業仕分け

2009-11-17T00:11:00.005+09:00

事業仕分けが着々と進んでいるらしい。まあ、従来の常識から考えるならば、これははっきり言って小泉劇場のやり口その物で、「抵抗勢力」と改革側がチャンチャンバラバラのやりとりをしている様を敢えて国民の目前でやっているという、パフォーマンス。

小泉は郵政民営化は成し遂げたが、道路公団改革は骨抜きとなってしまった。まあ、そこまでいくこともないだろうが、こういう改革で大概「総論賛成、各論反対」とならないものは想像できないわけで、今回の事業仕分けも、「あれも削減、これも削減」と言って国民の注目を集めつつ、結局のところ削減することはするものの、その幅は適当なところで手を打って、落ち着くところに落ち着くんではなかろうか。そういう希望的観測を持っているのだが。

科学技術予算でも、色々削減となっている。まあ、スパコン世界1位は要らないと思うけど、個人的には1位は要らないが、浮いた予算で世界10位以内を複数欲しいと思っているわけで、おそらく科学者側から反対意見を唱えている人の多くはそれに賛同するところだと思う。スパコンの予算を削ってしまえと言っている科学者は、（削って自分の分野によこせと言っているわけでもない限り）少ないのではなかろうか。

この件に関連して、署名を集めているようだが、それに関しては別の記事で。

Stratonovich Rule

2009-11-16T02:24:00.012+09:00

未だに Langevin 方程式を考えている。質量は1で、ポテンシャルは簡単のため調和型。

これについて、運動エネルギーの変化の平均を考えると、となる。ここで、因果律よりは t<t₀ の物理量にしか依存しないはずだから ξ とは独立であり、第三項は0である。従って、 γ→0 (γT=const) なら平均的にエネルギーは保存するはず。

と、師匠と面会の時に言ったら笑われる。まあ、少し考えればその通りで、質量0の極限だと速度が「 Brown 運動」するはずだから、どんどん速度の二乗平均は増加していくはずだ。しかし、先の議論が間違っているようにも思えない。

試しに先の極限の差分化を行って、先の極限で運動エネルギーを見てみると、やっぱり増えているが、は誤差の範囲で0だ。

お手上げ状態で師匠と面会して、「差分化した E の増分を真面目に計算しろ、まずはポテンシャル0の場合。」と言われてすぐに気づく。ランダム力の2乗の項が E の増分に入っている。ポテンシャルがある場合には、となる。そうか、これが Stratonovich Rule ってやつか！

最初の議論では軽々しくなんて書いたが、そもそもこれは揺らぐ量同士の積だから定義されていない。こういう積を議論したい時は適切な定義を選ばなければいけないが、 Itō Rule を定義として用いても物理的描像とは合致せず、物理的に意味のある結果を出したければ Stratonovich Rule を採用しなければならない。 v_i+1 は当然 ξ と相関がある。しかもランダム力の平方の項が出るから常に正の相関だ。こうやって考えればエネルギーが増えるのも自然だ。

こういう罠に引っかかって、改めて確率微分方程式は怖いなあと思う。まあ、こうやって自分で色々失敗して身につけた知識こそが身になるものなのだろう。

本題とは関係ないが、 xmodmap でデッドキーの設定はどうやるのかなあ。 Itō の “ō” の入力が面倒くさい。フランス語・ドイツ語に登場する文字は AltGr + アルファベットで出るようにしてるんだけど（だから Poincaré とか Schrödinger はすぐに入力できる）、 ō なんて滅多に入力する機会ないしなあ。

加重線型最小二乗法

2009-11-16T00:11:00.009+09:00

誤差付き観測データに直線を当てはめる場合を考える。ここで、 x についての誤差は考えない（自分が扱っているのが数値実験だから、基本的に x はパラメタであり誤差はないので）。また、直線以外の当てはめは基本的に精度が悪くなるので、適当な変形で直線に直すことをまず考える。

ここで a, b は残差の平方和 e²、を最小にするように定めればよい。これは a, b についての連立一次方程式を解けばよいだけで、簡単のための略記を用いればとなる。

さらに、と仮定すると（ N は正規分布、二つの引数は平均と分散）であることがいえるので、これを用いてフィッティングの信頼性を判定できる（が、そもそも誤差の分布が Cauchy 分布に従うとかいうとんでもない状況であればこんなことにはならない）。

これに加えて、各測定の誤差を独立でかつ相対誤差は十分小さいと仮定すると、誤差の伝播則を用いることができ、真面目に計算するとが得られる。

加重最小二乗法についてあまり触れているサイトもないし、自分の計算結果をメモ代わりに。

大地の息吹水の唱声宇宙とのハーモニー

2009-11-15T21:59:00.004+09:00

一部の方から好評を博しているトンデモ Adsense 集だけど、ついにここまできたかという広告。クリックして爆笑してしまった。皆さんもぜひクリックしてみてください。

恩寵宇宙エネルギーを水に受けとり高波動水をつくるアンテナ

まさかこんなのを科学的と信じる人はいないだろうから、トンデモに含めるのは違うかもしれないけど、面白かったので紹介。

Pileon

2009-11-09T02:01:00.004+09:00

先週の日曜、UbuntuをJauntyからKarmicにアップデートした（検索に引っかかるようにUbuntu 9.10とも書いておこう）。GCCは4.4.0が標準なようなので、C++-0xのautoを使いたいがために野良ビルドしたGCC4.4.0を削除する（でも4.5が出たらLambda使いたいがためにまた野良ビルドすることになると思う）。zshも4.3.10にアップデートされているな（RPROMPTに`zsh --version`を表示しているのですぐに気づく）。Open Officeもアップデートされているようだが、ほとんど使わないので気にしない。GCCの他に野良ビルドしているアプリケーションといえばpteTeXなのだが、LiveTeXの方はUTF-8のソースコードに対応しているのだろうか。対応しているようならそちらもCanonicalの方のパッケージに戻してもよいんだが（誰か情報ください）。先週の日曜は音が出ない問題等々に忙殺されてしまった。

フリーセルをしようと思ったらない。一瞬焦ったが、AisleRiotソリティアを初めて開いてみたら、そちらに収録されていた。山ほどゲームがあるので色々やってみた。個人的には完全情報ゲームの方が好き（なのでWindowsを使っていた時代にもソリティアはほとんどやらなかった）だが、その中でも特に面白いと思ったのがPileon。日本語ではピレノンと書かれているが、どう考えても音写ミスだろうな。意味から考えても"Pile on"ととるべき。

最初は全然クリアできなかったが、だんだんコツを掴んできた。ほとんどの場合、最初の数手、つまり空きスペースにどれを置くかでそのゲームをクリアできるか否かが決まる。始めてすぐは、**Xの形の山と**XXの形の山ができてしまうせいで手詰まりになるケースが多かったので、それを避けることくらいは学んだ。

それにしても、RhythmBoxとかPidginとか、画面右上に表示されるポップアップ（って呼び名でいいの？）のフォントって変えられないの？　なんかKarmicにアップデートしたらフォントがセリフ体に変わってしまった。これに関しても情報持ってる方がいたらお願いします。

Windowsは良くできてますよ。

自分が（多少Windowsより不便なところがあっても）Linuxを使うのは、研究やる上で空気のようにUNIX環境を使いたい（使う必要がある）から。その辺、Mac OS Xでもまあいいんだけど、WindowsでCygwinはあまりにストレスがたまる。まあ、でもDTMやるときはWindows使ってるけどね。一応WindowsとUbuntuの両方で使う可能性のもっとも高いFirefoxはプロファイル共有してる。

Dr. Giovanni Jona-Lasinio

2009-11-09T01:31:00.002+09:00

先日の水曜、Jona-Lasinnioさん（既に敬意を表して敬称を省くべき先生かもしれない）の講演を聴いてきた。しかし、圧倒的に英語力が足りない。Readingは専門用語を知らないくらいで、辞書引いて時間かければ何とかなるし、Writingもたどたどしいが何とかなる。ただ、Listeningはゆっくりはっきり話せば聴き取れはするものの、ちょっとマイクに入ってなかったり、聴衆の咳払いがあるだけですぐ分からなくなる。あと、今回は関係ないが、Speakingは圧倒的にダメ。

印象に残った点：

「対称性の自発的破れのような面白い現象が主に物性で見られるのはなぜか」という質問に対し、「物性は温度とか圧力とかいろいろパラメタいじれるけど、素粒子のパラメタは自然が決めちゃってるからね」との答え
Noticeable asymmetries still awaiting a convincing explanation in which spontaneous symmetry breaking may play a role:
1. matter dominance over antimatter in the near universe
2. dominance of left-handed molecules over right-handed in living matter
（講演会のスライド）
この2つを並列にしている点

特に後者はなかなか普通の素粒子屋さんからは出ない発想。さすが南部先生の共同研究者だ。素粒子屋さんとか宇宙屋さん（はそうではないかもしれないけど、特にそれを目指している学生）が陥りがちな行き過ぎた還元主義についてはそのうち別の記事を書くと思う。

どうでもいいけど<p>タグはブロック要素を包含できないのね。今日始めて気づいた（or昔は知っていたが完璧に忘却していた）。今までのは面倒なので放置するけど、これからは気をつけよう。

新たに文章を書くのであればDTDに沿った形にマーク付けしていけば良いだけだが、既存の文書をマーク付けするときには、「段落の中のリスト」のように、マーク付けしたい文書の構造がどう考えてもDTDに合わないので困るという問題が生じる。用途を詳しく分析した上で DTDを設計すれば、そのようなことはあまり起きない (かもしれない)。しかし、HTMLのように一つのDTDを様々なタイプの文書に用いるのでは、この種の不都合は避け難い。

（矢野啓介「HTMLにおける「段落」をめぐって」）

ちなみに、これは1998年の文書だが、XHTMLもかなり使われるようになった今日となっては、HTMLといえども終了タグの省略はあまり推奨されたものではないとも思う。

VIVA GIANTS 2009

2009-11-08T02:53:00.005+09:00

今年の日ハムは強かった。セ・リーグの他チーム（特に中日以外）が弱かっただけに余計にそう思える。ファイターズ選手、コーチ、スタッフ、そしてファンは誇りに思っていい。そして、その素晴らしいチームをさらに上回って優勝したジャイアンツ、おめでとう！

打線は小笠原、ラミレス、谷、阿部の主軸が2つのクリーンナップに別れ、そこに坂本・松本コンビや繋ぎの5番、亀井がいい味を出していた。投手陣も、グライシンガーやクルーンだけに留まらず、思わぬ活躍のゴンザレス、オビスポに加え、内海、尚成、東野のローテーションに、豊田、山口、越智のセットアッパーの堅さ。もちろん、キムタク、大道、脇谷、古城、寺内、鈴木、工藤の「脇役」陣も忘れてはならない。これだけの選手が躍動した戦い様は、まさに横綱相撲といってよいものだった。

コーチ陣の活躍も極めて大きかった。これだけの選手を活躍させたのに、育成コーチの手腕が影響していないはずがない。そして、伊原、尾花の両コーチが参謀となり攻撃・守備の戦術を立て、ルンバと揶揄されながらも人望の厚い原が彼らを一丸にまとめ上げた。

来年は選手・コーチとも着々と補強が進んでいる。キムタクと尾花の離脱は極めて大きな打撃であるが、来年もぜひ連覇を目指してもらいたい。それ以上に、来年はもう少し他チーム（特に横浜）も、育成を頑張ってほしいところ。

それにしても、今年は野球ファンにとって寂しい年になった。巨人に縁のあった選手だけでも、キムタクに加え（それにしても、生え抜きではないのにこんなにも生え抜き同様に巨人ファンに愛される選手もなかなかいない）、主な所では仁志、清水、江藤、工藤、斉藤、木田。セ・リーグV3へのキムタクの貢献は、小笠原やラミレスと同等に大きいものであったし、引退は他チームで迎えることになったが、生え抜きの仁志や清水、斉藤、そして外様ではあるが江藤や工藤の巨人軍への貢献も極めて大きかった。球界全体では立浪、緒方、ジェフ、小宮山、磯部とか。本当にお疲れ様でした。

追記。工藤は西武に再就職決まりましたね。おめでとう。実力での一勝をぜひ期待します。

誤差論

2009-11-03T00:09:00.003+09:00

誤差論では（唯一の）「真の値」の存在を仮定する
我々が得たいのは真の値である
計測値は様々な理由により真の値からずれている
このずれを誤差(error)という
誤差は系統誤差と偶然誤差に分類される
系統誤差は発覚しにくいが、発覚したら除去は容易である
偶然誤差は正規分布にしたがって発生すると仮定する
系統誤差を除去した後、複数の標本をとることにより、真の値は区間推定に従ってエラーバー付きで推定する

というのが、一般的な誤差論の肝。だから、例えば量子力学に従う物理量を何回も測定して、その平均をプロットして、標準偏差をエラーバーにするという行為にはあまり意味がない。その測定の過程には「唯一の真の値」は存在しないから。

今考えてみれば全く当然のことだけど、自分で計画を立てて数値実験をしてみると全然分かってなくて、確率変数の標準偏差をエラーバーにして「あれ、エラーバーえらいでかいな」とか悩んでいた。まあこういう試行錯誤を通してみんな分かっていくんだろうけど。

水は答えを知っている

2009-11-02T20:25:00.005+09:00

伊藤園の水素水の記事のAdsenseが大変なことになってます。

電解水素水
波動水
水素豊富水
還元水
飲料水

誰かこの電波飛び交う状況から助けてください。自分がネタにするからそんな広告集まってくるんだろって話だけど。それにしても波動水ってなんだよ波動水って。

小学生

2009-11-01T23:45:00.005+09:00

【韓国ブログ】なぜ韓国人は日本の目を気にしすぎるのか

この記事に、「日本人は韓国人と違って韓国（朝鮮半島）に興味ないよ！」という趣旨のコメントが大量についている。本当に興味がないのであればわざわざコメントなど残さないわけで（これが韓国人ではなくトルクメニスタン人であったと仮に考えてみればよい）、しっかり日本人（というかネトウヨ？）が韓国に興味津々である様子が表れている。

好きな子が気になって気になって仕方ないあまり、ついつい悪戯をしてしまう小学生のようで、見ていて微笑ましい。

どうでもよいが、昨日友人と食ったキムチ鍋は非常に旨かった。

酸化還元電位

2009-11-01T17:42:00.004+09:00

再びAdSenseの広告（もちろん広告はエネルギーの項目から）。今度はメジャー企業、伊藤園。

伊藤園の"還元性水素水"
鹿鳴越山系の伏流水使用。弱アルカリ性の軟水が体にやさしく浸透します
www.kenkotai.jp

おいおいおい。ページを開いてみたら酸化還元電位がマイナスとか書いてあって、そりゃ確かに基準をpH 7の水に設定すればそうなるでしょうが、あんまり書いてる人が酸化還元電位を分かってるとは思えない。まあ、水素は優秀な還元剤だけどね。

まあ、トンデモと言い切っちゃうのはちょっとはばかられるけど、だからって効能があると宣言するにはまだまだ早い代物だと思う。

差別化

2009-10-31T00:45:00.002+09:00

民主党は「国民の生活が第一」と唱えている（本当に実践できているかは別として）。製造業派遣の禁止にはっきり現れるように、多少の経済成長は犠牲にしても国民1人1人が多少は豊かに生活できるようになるという方向だろう。一方、自民党は従来通りの均霑論というか、（多少痛みを伴っても）経済成長をすることでやがては国民全体が豊かになるという、（少なくとも短期的には）「持てる階層」にとって好ましい政策をとるよう、民主党と差別化を図っている印象を代表質問から受けた。

（個人的にどちらを支持するかは別として）これは日本の政治にとって非常によい傾向だと思う。今までは自民党内の政権交代で、この2つの間を適当に行き来してきたわけだが、国政選挙の結果が直接反映されない以上、必ずしもそれは民意を反映しているとは言い難かった。一方、今回の総選挙も無論そうだが、その前の郵政選挙も、民意を受けてこのうちどちらかの政策が支持されたと考えることができる。そして、その政策が「やり過ぎ」だと思われれば、次回の選挙で政権交代が起こるだろう。

外交・安保に関してもうだが、これと並んで経済財政政策は非常に大きな論点である。このような論点に、（それほど大きくは違わないが）方向性の異なる2つの選択肢が国民に提示され、それを選挙によって選択するというのは、今までより格段に成熟した民主主義のあり方だと思う。

Enjoy programming

2009-10-30T21:04:00.004+09:00

Enjoy Programmingとかいう割には、Rubyは細かいところで手が届かない。特に、イテレータは大好きで、よく使っているんだけど（C++にもさっさとLambda入らないかなあ……0xは10進法じゃなさそうだし）、ちょっとやはり物足りない。each_with_indexが入ったのも1.9でようやくだし、てかmap!とかにも実装してほしいんですけど。この類のメソッドはeach_with_indexの挙動をデフォルトにして、|item,i|のiをオプションにすればいいのに。

あと、


std::transform(
 InputIteratorBegin,
 InputIteratorEnd,
 OutputIteratorBegin,
 UnaryFunctor
);

は


OutputArray = InputArray.map { |x| Sentence }

とか書けるけど、


std::transform(
 InputIterator1Begin,
 InputIterator1End,
 InputIterator2Begin,
 OutputIteratorBegin,
 BinaryFunctor
);

はどうすんの？　全般的にもうちょっとこの辺関数型プログラミング的にならないかなあ。

全然話は変わるが、AWKは使えるときには超便利。しかし、たまにしか使わないから書く度に文法忘れてるので、毎回文法調べて書いてる。そしてそんなAWKをRubyから（別にperlでもPHPでもいいけど）簡単に呼び出せる``演算子はもっと超便利。

ちなみに、研究室のサーバに自動認証かけられなかったのは自分の手違いでした。

アグネスタキオン

2009-10-28T20:42:00.010+09:00

Adsenseの広告は自分でも面白いと思ったのをクリックすることがある（必要以上に多数回クリックするのはダメでも、自分で「超伝導入門」買いたいと思ってクリックすることは別に規約違反じゃないよね）。右のところにエネルギーという項目があるので開いてみる。すると、

タキオン製品カタログ
タキオンエネルギーを利用するための製品販売
www2u.biglobe.ne.jp/~ocean/

ん？

好奇心に負けてクリックしてしまう。

量子物理学によると、この物質的な宇宙はエネルギーの非常に凝縮された形状（フォーム）以外の何者でもないと理解されています。この宇宙の中に存在する全てのもの（最も微妙で、活性化した組織帯の純化した領域から、全体、物資の最も凝縮された領域まで）は、全てエネルギーの連続体の中で一列に並んでいます。

エネルギー連続体を通って流れているエネルギーは、1つの源泉から来ます。インドでは、それは神母と呼ばれます。キリスト教はそれを聖霊と呼び、多くの現代的な、ニューエイジの教えにおいて、それは宇宙エネルギーと呼ばれます。現代の物理学はそれをゼロポイントエネルギーまたは自由エネルギーと呼びます。

（出典：「タキオンとは何ですか？」http://www2u.biglobe.ne.jp/~ocean/tachyon/productnote/whattachyon.htm）

えええええ！！！　それどこの物理学だよ！　一応物理学では零点エネルギー(Zero-point energy)や自由エネルギーというタームがあって、ちゃんとした定義が与えられてるけど、そもそもこの2つは全然違う意味だし、ましてや「エネルギー連続体を通って流れているエネルギー」の源泉とか、何がなんだか。

まあこのサイト読んでる人はこんなのに騙されるわけないけど、こういう商売が罷り通る世の中ってどうなのよ。タキオンは馬だけで結構。

等重率の原理とエルゴード仮説

2009-10-25T02:49:00.010+09:00

田崎さんの『統計力学』について、こんなレスを見つける。

田崎さんは「系がマクロになると、許される状態のほとんど全てがほとんど区別できない」を公理にして「マクロな系が平衡状態に緩和する」を導出してたが、そこで暗に等重率の原理のようなものを使っている気がするんだけど・・

「あれ、田崎さんの本に緩和なんてブッ飛んだこと書いてあったっけ？」と思って久しぶりにIの方を開いてみる……確かに書いてあるな。

考えてみても、この指摘の本質は全く正しいように思う。筆者は平衡状態について

平衡状態についての基本的な仮定：ある系での（熱力学でいうところの）平衡状態の様々な性質は，対応する「許される量子状態」の中の「典型的な状態」が共通にもっている性質に他ならない。

（田崎晴明『統計力学I』培風館初版p85）

と特徴づけているが、これに基づき非平衡からの緩和を

次に平衡状態への緩和の問題を考えよう。（中略）こういった非平衡状態も，エネルギーがほぼUである限りは，やはり「許される量子状態」の中に含まれている（脚注：このような非平衡状態に対応するのは，ほとんどの場合，エネルギー固有状態ではなく，複数のエネルギー固有状態の重ね合わせになると考えられる）。しかし，これは典型的な状態ではあり得ないから，ごく少数の例外的な状態の一つである。（中略）そこから出発して，系が時間発展したとき，特別な事情がないかぎりは，系の状態はごくごく小さな例外的な領域から外に出て，典型的な領域に入っていくと期待される（脚注：例外的な状態を作るには，複数のエネルギー固有状態を絶妙の係数をかけて足し合わせる必要があると考えられる。系がシュレディンガー方程式に従って時間発展すると，各々のエネルギー固有状態の位相が変化し，絶妙の重ね合わせは失われていくはずだ。）。

（引用元同上pp86-87、強調はこのサイトの筆者）

と議論している。「期待される」「はずだ」ってのは、あくまで仮定だ。もちろん、実験事実と照合することで正当化されている仮定ではあるけど、著者がこの前に置いている唯一の統計力学的仮定（もう一つ仮定があるが、それは前述の「平衡状態についての基本的な仮定」であくまで平衡系に限った話）

マクロな系の基本的性質：マクロな量子系では，ある平衡状態に対応する「許される量子状態」のほとんど全てが（マクロな物理量の測定値で比較されるかぎり）ほとんど区別できない。

（引用元同上p85）

から導かれるものでは、全くないように思える。

この「マクロな系の基本的性質」は、これを認めることで等重率の原理がエルゴード仮説を仮定することなく自然に導入され、それに立脚して統計力学の議論を進めることができたという点で、（著者はこの様な論理展開は明文化されているにせよいないにせよ真新しいことではないとしているが、少なくともはっきり第一原理として明文化したことについては）この書籍の極めて非凡なところだと思う。しかし、これによって説明されるのは、理論の中に「時間」という変数が全く入ってこない平衡統計力学だけであって、それを平衡統計力学の範疇外である緩和の説明にも使おうとしたのは、ちょっと勇み足ではなかろうか。

「例外的な量子状態」から出発して時間発展したときに、未来永劫「例外」に留まりつづけるかもしれないし（「許される量子状態」全体の空間から「例外的な量子状態」全体の部分空間を抜き出したとき、その部分空間は時間発展に対して不変かもしれない！）、そこまでしなくても、「典型的な量子状態」から出発した時より極めて高頻度に「例外的な量子状態」を訪れ続けるかもしれない。こういう事態を防ぐためには、新たな仮定

「任意の（あるいは、「例外的な量子状態」の中のさらに例外を除いたほとんどの）初期条件を選んで系を時間発展させることで得られる状態のほとんどは『典型的』である。」

を加えないといけない。清水さんの方の草稿には、「強いエルゴード仮説」として、これに似た（もっと強い）仮定が紹介されている（一般的なエルゴード仮説とは異なり、対象を巨視的変数に限定しているが）。しかし、清水さんもその草稿で言及しているが、これじゃあ「非平衡状態は平衡状態に緩和する」という仮説を、大して変わらない別の仮説を使って言い換えただけに見えてきてしまって、ほとんど旨味はない。

こういうことに関してなかなか深く考える機会はないが、平衡統計力学の基礎付けすら難しいのであるから、況んや本質的に非平衡である緩和過程をやと思った。そして、田崎さんの本はそれでもやっぱりいい本だなあと思った（熱力学の方はHelmholtzエネルギーから入るところがエントロピーから入る佐々さん・清水さんと比べてあんまり好きではないんだけど……でもあっちも網羅的だし、いい本）。

追記。田崎さん本人のコメント（これ以下）で、この批判は誤読だと指摘されています。ぜひご一読ください。

Log

2009-10-25T02:24:00.003+09:00

このサイトのログを見るのが楽しい。アクセス解析された結果を見るのもいいが、やっぱり一番面白いのは生ログ。こんなこと書いてるとストーカーじみてるように思われるかもしれないけど。

おそらく知人と思しきログや、熱力学・統計力学の記事に関して検索してきたログ（うちを含んだ大学からのアクセスが多い）が多い。中には教科書とか書いている偉い先生かもしれないログも（でも、他のページ見ずにすぐ帰ってるから別人かもしれない）。

OSはやっぱりWindowsが多いけど、Macもそこそこいる。Linux使ってるのは自分くらいかと思ってたら、どうやら他にもLinuxの閲覧があった。ブラウザはWindowsはIEが多い（あんなブラウザ、よく使う気になれるなと思うけど）が、2番手はChrome、意外。その後にFirefox。なんと懐かしきネスケのログもあった。MacユーザはSaafriとFirefoxが大体同数。

ちなみに、Adsenseのクリック数は0。まあ、別に期待なんか最初からしていないが、それにしてももう少し物理関係の広告だしてほしいな。

確率微分方程式

2009-10-19T03:08:00.008+09:00

化学反応のモデルとして次の師匠との面会までにLangevin方程式を解くことになっている。そんなの余裕だろと思っていたら、先日の面会では「罠があるよ」と言われた。いざコードを書こうと思ったら、4次のRunge-Kuttaではランダム力を単純に外力として入れるわけにいかないことくらいはすぐ気付く。

4次のRunge-Kuttaは関数を4回使うからランダム力をどうやって入れたらいいのか分からないのであって、単純なEuler法なら行ける気がするけどなあ……と思いながら、でも確率微分方程式の数値解法は今後も頻繁にお世話になると思われるので、Runge-Kuttaの単純な拡張があるのかとりあえずGoogleで検索かけてみる。そうすると、どうも既に「終わった」話題ではなさそうで、不安定性とかそういうキーワードがすぐに目につく。そりゃそうか、Fokker-Planckが不安定性生じる可能性あるなら、当然Langevinでも不安定性が生じていいよな。

普通の常微分方程式なんてのは、4次のRunge-Kuttaなる万能解法があるわけで、保存量を離散化させても保存させたいからシンプレクティック解法使うとか、まあそういう特殊な事情がなければこれ以外を使う必要はない。一方、偏微分方程式はちょっと気を抜くとすぐに解が発散するから、方程式に応じて適切な解法を考える必要がある。Fokker-Planckも、そのものを解いたことはないが、過去に反応拡散解いたときはえらい苦労させられた。確率微分方程式は、形の上では常微分方程式なので前者のクラスに属するかと思ったら、偏微分方程式同様、方程式ごとに解法考えないといかんという理解でいいのかなあ。まあ、LangevinとFokker-Planckは等価だから、多分そうなんだろうなあ。

でもとりあえず、まずは深いことを考えずに（とりあえず誤差評価だけは予めしておくけど）Euler法を回してみよう。対応するFokker-Planckも回して、結果を比較した法がいいかもしれない。

今まで特に何も考えず、Fokker-Planckと例のドリフトのある拡散方程式のことを言っていたけど、このPlanckって、かのMax Planckのことか。今日まで気付かなかった。やっぱり物理学会の大巨人だなあ。化学反応の話を師匠としたとき、師匠は物理の観点からArrheniusを導出したのはKramersと言っていたが、彼もKramers-KronigとかWentzel-Kramers-Brillouinとか、大はつかないまでも巨人。BrillouinもWKBだけじゃなくBrillouinゾーンとかに名を残している。とか、こんなことを考えていると、自分もすっかり物理オタクになったもんだなあと感慨深い。太田先生には足元にも及ばないが。

研究テーマ

2009-10-16T23:32:00.008+09:00

昨日は、研究室セミナー（自己紹介）と歓迎会があった。自己紹介（研究の概要紹介）はさすがによく分からない。歓迎会ではとある大先生の武勇伝（研究室でキレて醤油ブチまけたとか）が聞けた。ほかにも「あの人は一時期株にはまってた」とか「この人は偏見が強いからなあ」とか。M1の先輩によると熊本ラーメン（今年度秋の物理学会は物性は熊本）はまずかったらしい。

今日は師匠と面談して、研究テーマ決め。5つ考えてこいとの難題で、2つは内容を詰められたが3つは苦し紛れ。結局苦し紛れな方の3つは捨てて、残り2つの、反応速度論を非平衡から真面目にやるか、超流動（超伝導）の輸送係数の消失を非平衡で真面目に出すかというところに決定。まずはこの2本立てでやってみて、あとで絞ろうという話に。

とりあえず、今日は疲れたから自分の手は動かさずにまずは勉強。現象論としてGinzburg-Landau理論を調べてみる。巨視的波動関数って何じゃ一体。まずは現象論を抑えて、次に先学期聴いた院の講義（最後フェードアウトしてしまったけど）にも出てきたBCS理論、Bogoliubov理論で平衡統計力学の範疇の理解をしなきゃ。あと、線型応答も極めて怪しいので、応答係数と時間相関関数の関係とかKubo公式の導出あたりは土日にやっといた方がいいかもしれない。

そんなこと考えてると気が滅入ってくるので、研究室のサーバの設定もしてみようと思い、まずはSSHでログイン。authorized_key2をSCPで送ってみる。しかし、一向にパスワードの入力がサボれない。これはひょっとするとRSA認証を許可してないのか？　今度管理してる先輩に会ったら聞いてみよう。それにしても、SFTPが使えないようだ。これは、テキストファイルのちょっとした編集をするのにはEmacs（かvi）を使うしかなさそう。Emacsはあのキーバインドにどうしても慣れないので、今まではgEdit使ってたけど、さすがにそろそろ宗旨変えが必要かなあ。気が乗らないなあ。

今日の講義の名言。「Kramers-Kronig関係式の導出には因果律しか使ってないので、どれだけ強い非平衡であっても、例えば生物をどついてもその応答はKramers-Kronig関係式が成立する」

平均場近似と擬似自由エネルギー

2009-10-12T03:02:00.011+09:00

例えばIsing模型の平均場近似では（実際にやったのは格子気体だけど、簡単のためIsingで記載）となるわけだけど（zは隣接サイト数、Jは結合強度）、これの対数とって自由エネルギーG=F-HMを求めて、それをHで微分してmに関するself-consistent と持ち込もうと思ったら(M=mN)、師匠から「<s>はβとHの関数でしょ、それもちゃんと微分しないとダメじゃん」という趣旨のツッコミを受け、ホワイトボードの前で大混乱する。

結局その時は、師匠の助け船で、「分配関数（この場合は平均場近似により事実上一体問題に帰着できたので、「一粒子」分配関数で良い）でBoltzmann因子を除したものが確率だから、それを使って直接平均値を出せば良い」ということで事なきを得た。ついでに「下手に熱力学を使うと混乱するから使うな」とも言われ、「平均場が分かっていない」と叱責（語調はそうではなかったけど、後々のためにそういうこととして捉えておくことにする）を受けた。

しかし、例えば長距離相互作用模型では平均場と全く結果は同じになるはずだから、統計力学的確率分布が使えて熱力学が使えないはずがない気がする。そういうわけで、相転移・平均場近似のあたりを復習してみると、ここで求めたGは自由エネルギーではなく、「擬似自由エネルギー」というやつかと気付く（というか田崎さんの本に書いてある）。これはG(N,β,H)と捉えるべきではなく、G'(N,M,β,H)と捉えるべき関数なのか。真の自由エネルギーはLegendre変換によって求められ、self-consistent方程式を解くことがGを最小にするMを求めている操作に他ならないことに気付いた。

似たようなところで混乱する人はきっといるはずなので、記録を残しておくことにする。（といっても、好き放題書き散らかした文章なので、下手したらこれ読んで余計混乱する人がいるかも……あらかじめ、不十分な理解により誤りもしくはミスリーディングな部分が含まれている可能性があると予防線は張っておく）

無明の酒の酔い心

2009-10-10T05:02:00.003+09:00

大学で薪能をやるというので、見てきた。演目は狂言「萩大名」能楽「紅葉狩」。狂言はともかく、能楽はほとんど詞章が聞き取れないので、外国語のオペラみたいなものだと思って舞と囃子を楽しむ。

最初に教授の解説があって、「紅葉狩は世阿弥のような幽玄を追求したのみの作品ではなく、エンターテイメント性も重視しており、その点歌舞伎に通ずるところがある」とあったが、納得。女の裳束を羽織った鬼が一気にそれを脱ぐ演出とか、確かに圧巻で、こういう演出が歌舞伎の外連に繋がっていくのだろうと実感。能楽と歌舞伎は対照的な芸能だとか見たことない人に説明される風潮があるけど。

それにしても、ちょっと音響が微妙だったな。台詞、地謡、囃子全部をピンマイクか何かで拾ってスピーカーから流してたけど、臨場感が激減。まあ、そうでもしないと聞こえないんだろうけど。ミュージカルとかでも、ソロになるとマイクを使うときとかあるけど、ああいうの大嫌い。役者のいるところから声が聞こえてくってのが、一つ生の舞台の臨場感を演出してると思うんだが、どこに誰がいようと一切関係なく全部ミキシングして同じスピーカーから流すってのは、はっきり言って興ざめ。

以下、全くどうでもいい感想だけど、普通現代の芝居（無論ルーツは西洋演劇）って、例えばスタッフなんかが見切れるのを極端に嫌うけど、例えば能楽の後見とか、歌舞伎の黒子なんてのは、彼らの衣裳を「舞台上に見えないものとする」という記号として、堂々と舞台上に上がってくるんだよね。現代演劇の常識で考えたら新鮮だったけど、実は歴史的に演劇一般とかだとそっちの方が主流だったりするんだろうか。

YMO

GNU/Linux 環境

卒研進捗状況

(int)

科学技術関係予算の確実な確保について

科学研究費補助金の一部の執行停止に対する反対署名

事業仕分け

Stratonovich Rule

加重線型最小二乗法

大地の息吹 水の唱声 宇宙とのハーモニー

Pileon

Dr. Giovanni Jona-Lasinio

VIVA GIANTS 2009

誤差論

水は答えを知っている

小学生

酸化還元電位

差別化

Enjoy programming

アグネスタキオン

等重率の原理とエルゴード仮説

Log

確率微分方程式

研究テーマ

平均場近似と擬似自由エネルギー

無明の酒の酔い心

大地の息吹水の唱声宇宙とのハーモニー