Particle simulation of cosmic ray acceleration/transport in global heliosphere

NASAのOMNIデータを用いて日ごとの太陽風変動の情報を入力として取り入れた3次元MHD計算を行い、太陽圏の大域的構造を再現した。Voyager 1号、2号の軌道に沿った太陽圏界面の位置が120天文単位付近にあることを確認した。太陽風中の乱流の効果で、終端衝撃波の局所的な衝撃波角が50度程度になる場合があることに着目し、この場合の終端衝撃波の2次元フル粒子計算を行った。このときピックアップイオンの存在を陽に考慮して、背景電子、陽子、およびピックアップイオン（相対密度25%の陽子）から成るプラズマを想定した。アルフベンマッハ数は5.1、上流電子ベータ値は0.25、上流電子のプラズマ周波数とサイクロトロン周波数の比は4、イオン・電子間質量比を100とした。斜め衝撃波特有の現象であるイオンの反射が見られたが、反射されるのはピックアップイオンのみで、背景イオンは反射されないことがわかった。一部のピックアップイオンはイオンジャイロ周波数の逆数の100倍程度の時間をかけて非熱的エネルギーにまで加速された。加速機構として衝撃波ドリフト加速が確認されたが、詳しい粒子軌道解析によって、衝撃波での反射直後には異なる加速機構がはたらくことがわかった。反射イオンが励起する上流波動がピックアップイオンの初期加速に影響している可能性について、詳細を解析中である。定常太陽風条件下での太陽圏の3次元MHD計算結果を用いて、銀河宇宙線の太陽圏内への進入過程をテスト粒子計算した。従来よりも桁違いに多くの粒子軌道を追うことで統計精度を上げることに成功しており、現在データを精査中である。

NASAのOMNIデータを用いて日ごとの太陽風変動の情報を入力として取り入れた3次元MHD計算を行い、太陽圏の大域的構造を再現した。Voyager 1号、2号の軌道に沿った太陽圏界面の位置が120天文単位付近にあることを確認した。太陽風中の乱流の効果で、終端衝撃波の局所的な衝撃波角が50度程度になる場合があることに着目し、この場合の終端衝撃波の2次元フル粒子計算を行った。このときピックアップイオンの存在を陽に考慮して、背景電子、陽子、およびピックアップイオン（相対密度25%の陽子）から成るプラズマを想定した。アルフベンマッハ数は5.1、上流電子ベータ値は0.25、上流電子のプラズマ周波数とサイクロトロン周波数の比は4、イオン・電子間質量比を100とした。斜め衝撃波特有の現象であるイオンの反射が見られたが、反射されるのはピックアップイオンのみで、背景イオンは反射されないことがわかった。一部のピックアップイオンはイオンジャイロ周波数の逆数の100倍程度の時間をかけて非熱的エネルギーにまで加速された。加速機構として衝撃波ドリフト加速が確認されたが、詳しい粒子軌道解析によって、衝撃波での反射直後には異なる加速機構がはたらくことがわかった。反射イオンが励起する上流波動がピックアップイオンの初期加速に影響している可能性について、詳細を解析中である。定常太陽風条件下での太陽圏の3次元MHD計算結果を用いて、銀河宇宙線の太陽圏内への進入過程をテスト粒子計算した。従来よりも桁違いに多くの粒子軌道を追うことで統計精度を上げることに成功しており、現在データを精査中である。