一颗宇宙线粒子从被加速的那一刻起，就踏上了一条远非“直线飞行”的路。它在银河系的湍流磁场中反复偏转、散射，进入日球层后又遭遇太阳风与行星际磁场的动态挤压与筛选。真正抵达探测器的粒子，携带的不只是自身起源的信息，更是一整条传播路径上磁场结构与等离子体活动的叠加印记。理解宇宙线传播，因此不是一项“背景修正”工作，而是解读宇宙线观测数据的核心环节本身。
银河系中的随机漫步：扩散不是近似，而是物理现实
宇宙线在银河系尺度上的传播，通常被描述为“扩散”。这个描述容易被误解为一种数学上的简化，仿佛粒子实际上走的是某条确定路线，只是我们不知道罢了。真实情况更为彻底：星际磁场中充满了尺度不一的湍流结构，带电粒子沿磁场运动时不断被这些不规则结构散射，行进方向持续改变，其轨迹本质上是一个随机过程。
扩散的具体效率，取决于粒子回旋半径与磁场湍流特征尺度之间的比值。当回旋半径小于湍流尺度时，粒子的运动受磁场的随机性主导，扩散系数表现出一种能量依赖关系；当回旋半径超过湍流尺度后，传播则更接近经典的平均自由程图像。这意味着，不同能量的宇宙线在同一个银河系中所“感受”到的传播环境并不相同——高能粒子更容易穿透局部湍流，路径相对直接；低能粒子则被更频繁地散射，在银河系中滞留更久。
近年来的传播模型进一步表明，单一扩散系数不足以同时解释质子、氦核的能谱和大尺度各向异性。一个包含银河系盘面成分、银心持续注入成分以及近邻源（如与Geminga脉冲星相关的超新星遗迹）的多尺度扩散框架，能够在合理参数下同时再现这些观测量。传播不是均匀、各向同性的“背景板”，它的空间结构本身就在塑造我们看到的能谱形状和各向异性图案。
进入日球层：太阳调制作为动态过滤器
如果说银河系扩散决定了宇宙线“大致从哪里来、用多久来”，那么太阳调制则决定了“此时此刻有多少能通过”。日球层是太阳风撑开的等离子体泡，其中的行星际磁场携带着太阳风的流速与不规则结构。宇宙线进入这一区域后，会经历磁场的散射、太阳风的对流携带、以及因向外膨胀导致的绝热能量损失。
太阳调制最容易被忽视的特征是它的电荷依赖性。行星际磁场的漂移效应——包括大尺度梯度和电流片结构——使正电荷粒子与负电荷粒子在日球层中的传播路径出现系统性差异。CALET在国际空间站上对质子和电子计数率的同步观测显示，在同一平均刚度下，电子计数率随电流片倾角的变化幅度显著大于质子，这一差异可以用包含漂移项的数值传播模型再现。AMS-02在2011至2021年间对电子和质子日流量的精确测量，则将这种电荷符号相关的调制从定性图像推进到了定量检验的层面。宇宙线传播在日球层内不是一个“对称过滤器”，它对不同电荷的粒子执行着不同的通过规则。
更剧烈的证据来自太阳爆发活动。LHAASO在2021年11月一次日冕物质抛射事件中，于TeV能段观测到宇宙线大尺度各向异性的瞬时变化：在激波抵达后、磁通量绳结构到达前的鞘层阶段，来自太阳方向的宇宙线通量明显高于背向太阳方向，各向异性幅度约1%。这一观测打破了“TeV宇宙线各向异性在数十年尺度上恒定”的传统认知，说明宇宙线传播对行星际磁场结构的响应可以发生在小时量级上。传播不是静态参数，而是一个与空间天气耦合的过程。
用碎片反推路径：化学组成作为传播的“路标”
宇宙线传播的另一个可观测窗口，藏在粒子的化学组成里。碳核和氧核主要在恒星核合成过程中产生，属于原初粒子；硼核则极少在源区直接产生，主要是碳、氧核在传播途中与星际气体碰撞后生成的次级粒子。硼/碳流量比因此成为一个“传播深度计”：比值越高，意味着原初粒子在抵达探测器之前穿过的物质越多，碰撞产生次级粒子的概率越大。
这一比值随能量的变化方式，直接编码了传播路径的几何特征。高能宇宙线受磁场偏转较弱，倾向于走更“直接”的路径，穿越的物质较少，硼/碳比随之下降；低能粒子则被磁场反复弯曲，在银河系中逗留更久，碰撞更频繁。经典传播理论据此预言，在Kolmogorov湍流谱主导的扩散 regime 下，硼/碳比应近似按能量的−1/3次幂下降。AMS-02的高精度数据在百GeV/n以下确实与这一预言高度吻合，为扩散图像提供了有力的化学证据。
然而，“悟空”号将精确测量延伸到了TeV/n能区，发现了在约100 GeV/n处硼/碳比和硼/氧比偏离单一幂律、出现变硬的新结构。这一偏离提示，经典传播模型至少在高能端需要修正：要么银河系中的扩散系数并非简单的幂律形式，要么源的离散性、局域磁场的结构在更高能量上开始以不同方式介入传播过程。传播的“路标”并没有指向一个完结的答案，而是指向更复杂的、能量依赖的传播物理。
各向异性：传播结构在到达方向上的投影
如果宇宙线传播是完美的各向同性扩散，那么抵达地球的粒子应该来自天空各个方向、强度几乎均匀。实际观测却显示，到达方向分布存在微小但系统性的不均匀，相对强度变化在 (10^{-4}) 到 (10^{-2}) 量级，这就是宇宙线各向异性。它的存在本身并不反常——扩散传播并不要求完美各向同性，局部磁场的结构和源分布的不均匀都会在到达方向上留下投影。
各向异性的价值在于它充当了传播模型的“反向探针”。在垂直于扩散系数远小于平行扩散系数的情形下，偶极各向异性的方向倾向于指向局域规则磁场的方向，100 TeV以下的观测相位可以在这一框架下得到再现；而100 TeV以上的偶极相位则强烈依赖于局域扩散随能量的演化方式。这意味着，各向异性不只是传播的“副产品”，它可以反过来约束传播系数在局域星际介质中的行为。传播研究的困难也在这里集中体现：能谱、各向异性、化学组成从不同侧面触及同一组传播参数，任何单一观测量都不足以独立锁定它们。
传播不是旅途的“背景”，而是信号的一部分
回看整条链条，宇宙线传播的面貌已经不同于一个世纪前“粒子从A点飞到B点”的朴素图像。在银河系尺度上，扩散是湍流磁场与带电粒子相互作用的物理结果，不同能量、不同电荷的粒子经历不同的有效路径。在日球层内，太阳调制以电荷相关的方式动态调节通过率，太阳风暴甚至能在TeV能段留下小时级的各向异性瞬变。化学组成则从另一个维度证明，传播路径的长度和物质穿越量随能量系统性地变化，而高能端的新结构又在要求既有传播模型作出修正。
对天文网站而言，这意味着一条观测记录的含义从来不止于“探测器记录到了一个粒子”。它同时携带了源的信息和路径的信息，而这两者在数据中混叠在一起。宇宙线传播研究的工作，正是用多信使、多观测量、多能量段的交叉约束，把路径的效应从源的效应中剥离出来。每一次对传播模型的修正，都是对银河系磁场结构、星际介质湍动属性和日球层等离子体过程的一次更精细的读取。
问：宇宙线从产生到抵达地球，典型需要多长时间？
答：没有单一答案。低能粒子在银河系中扩散滞留的时间可达千万年量级，而高能粒子路径更直接，滞留时间显著缩短。进入日球层后，传播时间还取决于太阳活动水平与粒子电荷符号。传播时间不是一个固定常数，而是能量、电荷和太阳调制状态的函数。
问：太阳调制只影响低能宇宙线吗？
答：传统上太阳调制对低能粒子的影响最为显著，但LHAASO在TeV能段观测到的瞬态各向异性表明，太阳风暴驱动的大尺度磁场结构可以在很高能量上留下可测量的传播效应。太阳调制的影响边界比过去认为的更宽。
问：硼/碳比值为什么能用来研究传播？
答：硼几乎完全由碳、氧等原初核在传播途中与星际气体碰撞产生，属于次级粒子。硼/碳比越高，说明原初粒子在路径上穿过的物质越多。这一比值随能量的变化方式，直接反映了不同能量宇宙线在银河系磁场中的扩散路径长度和滞留时间。
问：“悟空”号发现的硼/碳比高能变硬意味着什么？
答：经典传播模型预期硼/碳比按能量的−1/3次幂平滑下降。DAMPE发现约100 GeV/n处出现变硬，说明实际传播过程在该能段偏离了单一幂律扩散的预言。这可能指向扩散系数的能量依赖形式需要修正，或源离散性、局域磁场结构在高能端的作用增强，但目前尚无定论。
问：宇宙线各向异性是否说明传播模型出了问题？
答：恰恰相反。扩散传播并不要求到达方向完全均匀；局部磁场方向、源分布和扩散系数的各向异性都会产生可观测的偶极或更高阶各向异性。各向异性是传播模型的重要检验量，而非对扩散图像的否定。
问：为什么宇宙线传播研究需要多种观测量联合？
答：能谱、各向异性和次级/原初比值对传播参数的敏感方式不同，且存在参数简并。例如，能谱的变化可能由源注入谱形或扩散系数共同导致，各向异性则对局域磁场方向和扩散各向异性敏感。只有联合拟合多类观测，才能有效约束传播模型，避免把路径效应误判为源的性质。
问：银河系扩散和日球层调制，哪个对抵达地球的宇宙线影响更大？
答：两者作用于不同能段和不同层面，不可简单比较“大小”。银河系扩散决定了粒子在银河系中的滞留时间、能谱形状和化学组成的基本框架；太阳调制则在此基础上对进入日球层的粒子进行电荷相关的动态筛选和能量损失。对低能粒子，太阳调制可以造成通量数量级的变化；对高能粒子，银河系扩散的结构效应更为根本。二者是串联关系，而非竞争关系。
问：宇宙线传播研究与暗物质探测有什么关联？
答：次级粒子（如反质子、正电子）的传播计算是暗物质间接探测中天体物理背景的核心组成部分。对硼/碳比等传播参数的精确约束，直接影响到对正电子过量或反质子过量是否源于暗物质的判断。传播模型的不确定性，在很大程度上也是暗物质间接探测背景不确定性的来源。