研究関連の日誌については別システムにしようかと思い、hnsとは別のシステムを捜してみました。
簡単な日誌システムを構築しようということで、fdusを導入してみる。Perlなので、いろいろ改造してみようかと。
ようやく気力が充実したので, いよいよ本丸. 以前書いた準備稿に沿いつつ, 追加データやらプロットしなおしやらが入るので, なかなか進まない.
銀河団やγ線背景輻射の評価に対して, 銀河団の生成過程を加味することが重要であるというのがこの論文の主旨. で, その生成過程において, 宇宙線を閉じこめることができていたかどうかが決め手となる, と言いたい. そのために, まず拡散で近似をして, 散逸の効果をとりこむことを主張する.
で, そのためのプロットを4枚.
- 1枚目は「標準質量」の進化. 現在ある質量の天体は, その昔は(平均的には)どれくらいの質量だったのかを示した図. 宇宙論パラメータの違いに対する質量進化の違いもあわせてプロットしてある.
- 2枚目は恒星生成率(SFR)の進化. 昔ほど星の生成活動が活発だったことを, 解析的な式にしてあるので(Madau et al.), それを利用することにし, 式の紹介が面倒なので, 図示. ただし, その式を本格的に使うのは後ほど.
- 3枚目は拡散係数の変化. 拡散を効かせるほど宇宙線は散逸する. 当然といえば当然の結果が出る. 特に高赤方偏移で顕著な効果. これは特にγ線背景を考えるときに効くということを意味する.
- 最後に天体質量に対する依存性. 重くなるほど閉じこめが効くことを図示. これは重くなるほど半径が増大することで理解できる.
と, ここまで書いたら時間および気力が尽きたので, 明日継続.
昨日はサッカー観戦で執筆が全然進まなかったので, 巻き直し. 時系列に沿って記録してみる. (16:07)
まずZiniに関する依存性を示すためのデータ作成. 実はこれはすでに先週木曜に作ってあるので, それをそのまま利用する. (16:09)
先週の木曜の執筆が進まなかったのは, このデータが原因. timestep をdynamical timeの何倍に設定するかというパラメータを用意したので, そのパラメータに対する数値の依存性を調べようと思ったら, 1/100に設定したときの計算が12h/runになってしまった(大学のDEC/Alphaを用いて. Pen4あたりだと1.5倍くらいになってもおかしくない. やる気もおきないので検証してないけれど). 木曜がこれで潰れたのが敗因... 他の部分から書けばいいのに, と後になって思ったりもしたけれど, 後の祭り. で, 1/100 と 1/10 は大学のコンピュータが頑張ってくれたので, 今日は1倍と10倍についての計算. (16:12)
と, 書いているうちに計算が終わるはずだったのに, 一向に終わらないので, Makefileを見てみたら, 1/100 の設定のままだった. これでは終わらない. ということで, 書き直して再度挑戦... 終了. (16:14)
で, Makefile があまりにも汚いので, Makefile の分割を試みる. Makefileに include 機能があればいいなと思いつつ捜すも, Web 上になかなか情報が見付からず. 仕方がないから英語なページを捜したら, 結局ここに辿りつく... って, 結局 GNU じゃねぇか. ってなわけで, 無事include を発見. 分割してみた. Good! (16:43)
プロットするために, Coma cluster内の宇宙線量を求めるプログラムを噛ませる. プロットデータを作成している間, コンピュータがハングったり, 1/100のせいでやはり時間がかかったりで, 意外に大変. (17:40)
データはできたので, 5枚目のグラフを作成... と. また調子が悪くなったのでリブートしたり, TV観賞が入ったりしたが, とりあえずサンプルプロット. 10倍はあまりにもデータが汚いので削除. 1/10と1/100はほとんど差がないことが示される. しかし1倍とはそれなりに差があるのと, 従来使ってきた1/2も少し離れていることが見えている. どうするのが妥当なのだろうか...(18:27)
1/100をグレー, 1/10を実線, 1/2を破線, 1倍を点線で描いておくことにしよう. 凡例をつけて, とりあえず見栄えを良くしてみた. (18:47)
父親の迎え, 食事その他を経て, 再開. (23:11)
ほーっとしていたら, いつの間にやら25:25(.hack観賞終わり). で, 書きだしていたらグラフが1枚足りないことに気付く. Ziniを変えた図を描いていなかったので, これを描いてみる. こっちはもっと味気ない. 結局どの赤方偏移で計算しても結果は変わらないことを図にするだけ. 3本の線がほぼ重なってしまっているので... ま, とりあえず載せておこう. (25:55)
で, プロット完了. (26:11)
いきなりハング. 再起動したらMeな起動画面からちっとも動かなくなる. でもHDDは何やら一生懸命働いてるし. 何やってたんだろう. 何はともあれ, ようやく復旧. (26:52)
父親が起きてきたのでその相手をしたりもしつつ, 本文を書きあげる. あとは図のキャプションを入れなきゃ. (27:31)
計算時にtimestepを決めるパラメータ, および計算開始の時刻を決めるパラメータが存在する. それらパラメータには物理的根拠がないため, これらに大きく依存するような計算であってはならない. その影響を図により示す.
- まずtimestepの影響. 基準となるのは力学的時間(重力によるboundの時間スケール). 1/10と1/100はほとんど差がない. ところが計算時間は10倍じゃきかない. ということで, 1/10あたりが妥当かも. でも通常の計算は今のところ1/2でやっている. この数字に根拠があるわけじゃないんだけど.
- つぎに赤方偏移の初期値に対する依存性. これは完全に無視できる. ただし閉じこめ効果が強いときは, 密度が十分上がるまでは差異が生じているのが見える.
とりあえずこんなところでよかろう. (28:07)
プロットしてみたら, Makefileのデータ生成にミスがあったことが発覚. 書きなおしてデータを作り直す. 計算を待つ間にII.を書きあげるものの, まだ終わらず... そういえば flat universe の時の計算て結構時間がかかるんだった... (28:25)
どうやら, 10min/run かかる模様. とりあえずデータ生成が成功したかどうかを見極めて, 今日は終了ということになりそう... 本当はγ線データの生成までやりたかったんだけどなぁ... 見積りで4時間くらいかかる仕事量なのに... (28:26)
とりあえずプロット完了. 12min/runでした. 凡例はEdS(灰)とLambda(黒). あとはグループごとに1pc, 100pc, 10kpcとコメントを付加. まずまず, かな. ということで, 本日はここまで. (28:45)
昨日の続きを再開. まずは昨日描いたグラフのチェック. 他のグラフにもちょこちょこ手直しをしてみた. あとはキャプションと本文. (15:24)
7枚目のグラフ用のMakefileの作成開始. (15:35)
ようやく終了... これからcaptionを直す作業. (16:26)
縦軸のミスを直しつつ, 本文が完成. ということで, しばしblack out. (16:55)
black outの間, 9枚目のグラフ用にデータを生成. 起きてみたものの, まだ終わってない. (18:00)
キャプション, 本文ともに完成. 以前のコードを流用するだけなのでさほどでもなかった. (19:17)
で, 9枚目を描きたいのだが, まだ計算が終わらない. とりあえずEdSだけでもプロットしてみる. (19:18)