[论文解读] Effects of the tilted and wavy current sheet on the solar modulation of galactic cosmic rays

Shoko Miyake, S. Yanagita|arXiv (Cornell University)|Oct 26, 2006

Solar and Space Plasma Dynamics被引用 3

一句话总结

本研究开发了一个三维、时变的随机微分方程（SDE）模型，用于模拟日球层中银河宇宙射线（GCR）的调制，该模型结合了动态演化、倾斜且具有波纹状的日球层电流片（HCS）。该模型成功再现了22年的太阳调制周期，并定性地捕捉了正负电荷粒子的不对称强度调制，但未能再现质子在 qA > 0 时的观测能谱，这可能是由于HCS结构过于简化，且在极区缺失了横向磁场分量。

ABSTRACT

Transport equation of the galactic cosmic ray (GCR) is numerically solved for qA>0 and qA<0 based on the stochastic differential equation (SDE) method. We have developed a fully time-dependent and three-dimensional code adapted for the wavy heliospheric current sheet (HCS). Results anticipated by the drift pattern are obtained for sample trajectories and distributions of arrival points at the heliospheric boundary for GCR protons. Our simulation reproduced a 22-year cycle of solar modulation which is qualitatively consistent with observations. Energy spectra of protons at 1 AU are calculated and compared with the observation by BESS.

研究动机与目标

• 研究日球层电流片（HCS）的倾斜度与波纹性如何影响银河宇宙射线（GCR）的太阳调制。
• 基于随机微分方程（SDE）开发一个完全时变的三维数值模型，以在真实的HCS几何结构下模拟GCR输运。
• 再现观测到的GCR强度22年调制周期，并将模拟的质子能谱与BESS观测结果进行比较。
• 识别当前模型的局限性，特别是未能再现正电荷粒子（qA > 0）的观测能谱。

提出的方法

• 本研究采用与Parker输运方程等价的随机微分方程（SDE）公式，包含太阳风流、扩散张量以及梯度-曲率漂移速度。
• 日球层磁场（HMF）被建模为Parker螺旋状，HCS通过一个随时间与径向变化的倾斜函数来表示，该函数随11年太阳周期演化。
• HCS结构由一个动态方程定义，包含太阳自转和随时间变化的倾斜角，使信号能以太阳风速度（约400 km/s）从太阳传播至80 AU。
• 粒子从1 AU反向追踪至日球层边界（80 AU），初始条件设定在1 AU的赤道面上，记录到达时间、能量损失及轨迹分布。
• 采用漂移速度场方法计算粒子穿越HCS的漂移，其行为具有电荷符号依赖性（qA > 0 与 qA < 0）。
• 将模拟的质子能谱归一化至1998年和2000年BESS观测数据，并进行对比以评估模型准确性。

实验结果

研究问题

• RQ1 HCS的时变、径向传播的倾斜度与波纹性如何影响银河宇宙射线的调制？
• RQ2 三维、时变的SDE模型能否再现观测到的GCR强度22年调制周期？
• RQ3 尽管捕捉到了调制趋势，为何该模型未能再现qA > 0质子的观测能谱？
• RQ4 由于HCS漂移效应，qA > 0与qA < 0粒子的轨迹与到达分布有何不同？

主要发现

• 该模型成功再现了银河宇宙射线强度的22年调制周期，但强度变化存在时间延迟，这是由于信号从太阳传播至日球层边界存在有限传播时间。
• 对于qA < 0，模拟的质子能谱与2000年BESS观测结果高度一致，尤其在倾斜角达70°且太阳极性为负的时期。
• 对于qA > 0，模型低估了观测到的质子强度，表明调制不足，这可能是由于极区缺失横向磁场分量，而这些分量本应抑制漂移运动。
• qA > 0粒子平均仅需4.5天即可到达1 AU，表明调制作用微弱；而qA < 0粒子则需80天，与更强的约束及沿HCS的漂移行为一致。
• 在80 AU处的到达点纬度分布反映了HCS结构：qA < 0粒子聚集在HCS附近，而qA > 0粒子则分布在极区附近，具体取决于观测者相对于HCS的位置。
• 该模型未能再现qA > 0能谱的事实凸显了未来模型中需引入更真实的HCS动力学，如随机横向磁场与非定常演化过程。