摘要
沥青的老化是沥青路面使用过程中最受关注的问题之一。然而，沥青通常以薄膜的形式存在于混合料中。氧气扩散是影响沥青老化的重要因素，为了探究不同工况下氧气在沥青膜中的扩散行为，该研究建立了氧气‒沥青‒集料 （OAA）模型，以模拟在实际情况下集料表面沥青膜状态。在不同时间、温度、氧气分子浓度和集料类型条件下，对 OAA 模型中氧气扩散和沥青组分迁移规律进行研究。通过计算分子的相对分子浓度、径向分布函数、均方位移和扩散系数，比较各因素对模型分子行为的影响，进一步分析各参数之间的联系。从分子模拟结果可知，氧气分子浓度随扩散时间的增加而增加，且集料与沥青界面处的氧气分子浓度明显上升。随着温度的升高，氧气向沥青膜内的扩散能力增强，氧气分子浓度随之线性增加。增加氧分子数量，沥青膜内的氧气相对分子浓度整体上升，CaCO3集料更利于氧气向沥青膜内扩散。
Abstract
Asphalt aging is one of the most important concerns during the use of asphalt pavement. However， asphalt usually exists in the form of a film in the mixture. Oxygen diffusion is an important factor affecting asphalt aging. To investigate the oxygen diffusion behavior in asphalt film under different working conditions， an oxygen-asphalt-aggregate （OAA） model was developed in this study to simulate the state of the asphalt film on the aggregate surface under real conditions. Oxygen diffusion and asphalt component migration in the OAA model were investigated under different conditions of time， temperature， oxygen concentration， and aggregate type. The influence of each factor on the molecular behavior of the model was compared by calculating the relative concentration of molecules， radial distribution function， mean square displacement， and diffusion coefficient， and the connection between the parameters was further analyzed. From the molecular simulation results， it is clear that the oxygen concentration increases with the increase in diffusion time， and the oxygen concentration at the interface between aggregate and asphalt increases significantly. Higher temperature promotes the diffusion of oxygen into the asphalt film， and oxygen concentration varies nearly linearly with temperature. Increasing the number of oxygen molecules can make the relative concentration of oxygen molecules within the asphalt film increase overall， and the CaCO3 aggregate is more favorable for oxygen diffusion into the asphalt film.
0 引言
沥青路面具有优越的服务性能和相对较低的成本，现已经被认可，并广泛应用[1]。虽然沥青黏结剂的质量在沥青混合料中最多占 7%，但对路面抵抗车辆荷载和环境影响有很大的贡献[2-3]。沥青混合料的老化和随之而来的性能恶化是路面研究领域的关键问题[4]。对于沥青路面来说，老化是指沥青混合料在生产过程中不可避免的化学、物理和机械变化（短期老化），以及在使用过程中持续暴露在环境条件下产生的老化（长期老化）[5]。此外，道路表面的光化学老化也起着关键作用[6]。沥青材料的老化主要分为 4 种机制：挥发、氧化、紫外线辐射和物理硬化，这些机制共同导致沥青使用性能下降[7-8]。一般来说，随着沥青混合料的老化，其模量上升，水分减少、开裂增加和抗疲劳性能下降，最终导致其使用寿命缩短[9-10]。因此，研究沥青材料的老化行为对沥青路面的设计、施工和养护具有重要意义。
众所周知，热氧化是导致沥青材料劣化的重要因素[11]。目前，大多研究集中在氧化的作用机制及其对沥青性能的影响[12-13]。在实际的沥青路面中，黏结剂是以薄薄的沥青膜形式存在于集料表面或集料之间[14]。通常情况下，沥青薄膜上的热氧化过程是一个从薄膜表面逐渐深入的渐进现象[15]，这导致了不同厚度的沥青膜在老化行为上存在差异。沥青膜的梯度老化行为直接影响其内部的均匀性，使得在车辆荷载和外部环境的共同作用下，沥青膜内部极易出现微裂缝[16]。这些微裂缝逐渐扩展，最终会导致沥青路面的开裂、松散及剥落等[17]。
沥青的氧化程度受氧气扩散行为的影响，随着时间的推移，氧气会扩散到沥青中，为沥青的氧化反应提供氧气[18]。因此，活性沥青成分可以与氧气反应，导致沥青硬化和脆化[19]。关于氧化后沥青的化学、物理化学和机械特性，已经在进行相关研究。沥青的化学成分随原油来源的不同而有很大差异，不同类型的沥青的测试数据可能显示出不同的性能和特点[20]。沥青中的氧化反应主要取决于沥青的微观结构和化学成分，但也取决于温度、压力和紫外线[21]。沥青中来自氧化反应的极性基团和强关联基团的比例随着老化时间的推移而变化，随之而来的是分子团聚和结构重排。氧化反应不断消耗沥青中扩散的氧气，导致氧气的进一步扩散[22]。当沥青中没有更多的反应物，且没有氧的浓度梯度时，就可以达到终结状态[23]。沥青中的氧气扩散和反应过程的耦合可以显示动态的氧气浓度和反应动力学，还可以用于寿命预测。扰动的微观结构和沥青的极性增加会通过更密集的分子分布阻碍氧气的进一步扩散。氧气是路面恶化的重要因素，导致沥青的持续降解。因此，研究不同因素下氧气在沥青中的行为尤为重要。
1 分子模型
从分子的角度来看，沥青的热氧老化主要是由于轻质组分分子的损失，以及沥青质和树脂分子与沥青中的氧气反应生成羰基、亚砜基等极性基团，最终导致沥青变硬、变脆。氧气分子浓度的大小决定了氧化反应进行的程度。而温度为氧气扩散提供能量，能量越高氧分子越活跃。除温度的影响外，现有老化试验通过增加压强的方式来加速氧气在沥青中的扩散。本研究采用 Materials Studio （MS）2020 软件，建立了四组分十二分子组成的 AAA-1 沥青、集料和氧气模型。为探究在热氧老化条件下集料表面沥青材料的氧气扩散和组分迁移特性，最终组建了氧气‒沥青‒集料（OAA）三层模型。
1.1 沥青
为建立沥青模型，将各组分分子按照 AAA-1 沥青模型的分子数量添加到分子动力学计算软件 MS 中 Modules-Amorphous Cell模块，AAA-1 沥青模型的分子组成如表1 所示。模型初始密度设定为 0.1 g/cm3，以消除沥青组分之间的聚集效应。然后，采用智能优化算法对初始结构进行 5 000 次迭代优化，之后在恒定压力和温度（NPT）系综下进行 500 ps（1 ps= 10-12 s）的 MD 模拟，进一步获得平衡的分子结构，时间步长为1 fs（1 fs=10-15 s）。其中，压力为101.325 kPa，温度设置为 298.15 K，最后在正则系综（NVT）下扩散 500 ps 得到沥青模型如图1 所示。AAA-1 沥青模型的相关性质参数列于表2 中。
表1 AAA-1 沥青模型的分子组成
Table1 Composition of AAA-1 asphalt model
图1 AAA-1 沥青模型
Figure1 AAA-1 asphalt model
表2 AAA-1 沥青模型性质参数
Table2 Properties of AAA-1 asphalt model
1.2 集料
采用 XRD 测试确定了本研究选用的石灰岩矿料的主要成分为 CaCO3。SiO2为沥青混合料常用集料的主要氧化物成分之一，较多研究将其作为集料的分子模型。因此，本研究分别建立了 CaCO3 和 SiO2 晶体模型，用于模拟集料对氧气和沥青分子行为的影响。以 CaCO3 模型为例，先从 MS 软件中导出 CaCO3晶胞模型，通过 Morphology tools 工具确定切面，选取暴露可能性最大的（1-1 2）晶面进行切割，厚度设置为 25.22 Å，截断半径 15.50 Å。最后，根据沥青模型尺寸进行扩胞得到 CaCO3 超晶胞分子模型，具体过程如图2 所示。
图2 集料分子模型建立过程示意图
Figure2 Aggregate molecular modelling process
1.3 氧气
先将 300 个优化后的氧分子导入 Amorphous Cell 模块，并根据沥青模型尺寸设定氧气模型的初始密度，使得氧气模型的尺寸与沥青模型接近，得到氧气模型。然后，对氧气模型进行 5 000 步的几何优化，最终得到如图3 所示的氧气模型。
图3 氧气模型
Figure3 Oxygen model
2 分子动力学模拟过程
2.1 组建氧气-沥青-集料三层分子模型
为了从分子尺度模拟集料表面沥青膜在热氧条件下的氧气及沥青组分的扩散和迁移行为，分别建立氧气、沥青和集料模型，然后通过 MS 软件中的 Build layers 功能依次导入集料、沥青和氧气分子模型，并在氧气模型上方增加 50 Å 真空层得到氧气‒沥青‒集料（OAA）三层分子模型，如图4 所示。
图4 氧气-沥青-集料三层分子模型
Figure4 OAA three-layer model
2.2 OAA 模型分子动力学模拟过程
探究集料表面沥青膜中氧气扩散和沥青组分迁移过程，对研究沥青膜梯度老化的分子机制具有重要意义。本研究借助分子动力学计算软件 MS，对建立的 OAA 三层分子模型进行分子模拟。采用 COMPASS 力场，首先将 OAA 模型（SiO2 和 CaCO3）进行 5 000 步的几何优化，在 NVT 系综（压强为 101.325 kPa，温度分别为 436.15 K、353.15 K、298.15 K、253.15 K） 下，进行了 1 000 ps、2 000 ps 和 3 000 ps 的分子扩散模拟。氧分子数量分别为 200个、300个，用于模拟不同氧气分子浓度。最后，利用 Forcite 模块中的 Analysis 计算不同模型的沥青组分和氧分子的相对分子浓度、径向分布函数和均方位移，并进行数据分析。
3 分子动力学参数
3.1 相对分子浓度
相对分子浓度是指某一分子相对于整个分子体系的浓度，通过计算体系中某一组分在坐标轴方向上的相对分子浓度来研究组分的分布情况。根据公式（1）、（2）计算相对分子浓度。
(1)
(2)
式中：CRC，i为 i组分的相对分子浓度；Ci为 i组分的分子浓度；Cbulk为体积分数；Ni为 i组分的颗粒数；Vi为组分 i 的体积。
3.2 径向分布函数
利用径向分布函数（RDF），表征沥青系统中分子在一定距离内的分布，反映系统中存在的分子的吸附行为。因此，RDF 可以用来评估沥青成分、橡胶和废油分子的分布。RDF 由式（3）计算：
(3)
式中：r 为粒子间的距离；N 为体系的粒子个数；ρ 为体系的平均密度。
3.3 均方位移和扩散系数
均方位移（SMSD）和扩散系数 D 之间存在对应关系，如式（4）、（5）所示。SMSD和扩散系数被用来研究不同因素对分子扩散行为的影响。
(4)
(5)
式中：a 为常数，与模型维度和移动方向有关；ri（t）为粒子 i 在时间 t 的位置向量；ri（0）为原子 i 在初始时刻的位置向量。
4 分子模拟结果与分析
4.1 沥青组分和氧气相对浓度影响因素分析
SARA 法将沥青分为沥青质（Asphaltene）、芳香分（Aromatic）、胶质（Resin）和饱和分（Saturate）。沥青四组分在沥青中具有不同的作用，沥青质是硬质的分散相，为沥青材料提供黏附性和强度；饱和分和芳香分作为轻质组分是沥青质的分散介质，轻质组分含量多会增加热沥青的流动性，同时常温时沥青较软；胶质与沥青质形成胶团，均匀地分散在轻质组分中，形成稳定的胶体结构，使得沥青成为均质的胶结材料。四组分的分布情况决定了沥青膜的性质，沥青膜中组分的迁移会导致其不同厚度位置的材料结构和特性发生变化。通常用于道路铺筑的沥青材料为均质状态。然而，在沥青混合料高温拌和、与集料接触以及沥青路面的实际使用过程中，受温度、集料和环境因素等影响，沥青膜内的组分会发生迁移。因此，探究组分的迁移规律对研究沥青梯度老化行为具有重要意义。
沥青组分在热氧状态下发生氧化反应是沥青老化的重要原因之一，沥青内部氧气分子浓度决定了沥青的氧化程度。然而，沥青内部的氧气从表面逐渐向内部扩散，其分布规律直接影响沥青膜的老化速度和程度。本研究利用 MS 软件中的分子浓度分析功能，对 OAA 模型中沥青内部的氧气分子浓度进行了计算。
4.1.1 模拟时间
老化时间是影响沥青老化程度的重要因素之一，沥青组分的迁移和氧气的扩散行为随着时间的增加逐渐发展。本研究在 NVT 系综下对 OAA 模型分别进行了 1 000 ps、2 000 ps 和 3 000 ps 的扩散模拟，并计算了体系中沥青各组分和氧分子的相对分子浓度，结果如图5 所示。图5（a）为沥青质的结果，图中有两处明显的峰值，分别位于 25 Å 和 40 Å 附近。在 25 Å 处，3 种模拟时间下，沥青质的浓度较为接近。而在 40 Å 处三者的差异较大，表现为随着模拟时间的增加相对分子浓度增大。40 Å 接近于沥青的中部，该结果表明沥青质分子从沥青膜表面向中间迁移。从图5（b）中可以看出：随着模拟时间的增加，胶质的相对分子浓度变化明显。在 25.26 Å、 31.26 Å 和 50.26 Å 出现了 3 处峰值，在 25.26 Å 处的相对分子浓度：3 000 ps>2 000 ps>1 000 ps，与 50.26 Å 处的浓度变化规律相反。而在 31.26 Å 处： 2 000 ps>3 000 ps，1 000 ps 的结果中没有峰值。这充分说明随着时间的增加胶质从沥青表面向集料界面迁移。图5（c）为芳香分相对分子浓度的结果，随着模拟时间的增加，在集料界面处饱和分分子浓度先增加后减小。而在氧气界面处的变化规律为先减小后增加。出现这种现象的主要原因是集料对轻质组分的吸附，而随着模拟时间的增加轻质组分挥发。饱和分的相对分子浓度如图5（d）所示，随着时间的增加，饱和分的变化规律与芳香分类似，会同时受到集料吸附和轻质组分挥发的影响。
图5（e）为不同模拟时间沥青膜内氧气分子浓度的结果，从图5（e）中可以看出：随着扩散时间和距离的变化，沥青中氧气分子浓度差异较为明显。50 Å 处 3 种扩散时间沥青内的氧气分子浓度呈现出显著的变化。从沥青内部（40~50 Å）氧气分子浓度可以看出：在常压 436.15 K 条件下，氧气分子浓度随扩散时间的增加而增加。集料与沥青界面处的氧气分子浓度均明显上升，氧气分子浓度的大小决定了沥青热氧老化的程度，氧气分子浓度的纵向差异促使沥青膜梯度老化特性的产生。
图5 不同模拟时间OAA模型的沥青组分和氧气相对分子浓度
Figure5 Relative concentrations of asphalt components and oxygen molecules in OAA model at different simulation time
4.1.2 温度
温度也是影响沥青老化程度的重要因素之一，温度越高分子活性越高，从而提升分子间反应的概率和速率。目前，通过分子动力学模拟还无法实现分子间的化学反应。然而，不同温度下分子浓度同样能显示反应发生的概率。为了模拟低温、常温和高温 3 种实际工作温度下氧气在沥青膜中的扩散行为，本研究将 OAA 模型在-20℃、25℃和 80℃条件下进行扩散模拟。图6 为 3 种不同温度下沥青组分和氧分子浓度的模拟结果。从图6（a）可以看出：在 25℃和 80℃下，沥青质的相对分子浓度曲线较为接近。当温度为-20℃时，30 Å 附近的峰值明显小于另外两种温度。这说明低温条件沥青质分子不易发生迁移。另外，随着温度的升高，沥青质分子表现为从中间向两界面处迁移的现象；而从图6（b）中胶质相对分子浓度曲线峰值分布的状况来看，当温度增加时，胶质分子有整体向集料表面运动的趋势；从图6（c）中可以看出：随着温度的升高，集料界面和沥青膜表面的饱和分分子浓度分别表现为先减小后增大和先增大后减小的变化趋势。这表明温度升高导致饱和分分子活性增加迁移规律动态变化；图6（d）中，在 60 Å 附近饱和分分子浓度曲线出现了明显的峰值，分子浓度的变化规律为：-20℃>25℃>80℃，且 25℃与-20℃的分子浓度差大于其与 80℃的分子浓度差。这表明温度是影响饱和分迁移的重要因素，低温条件下分子惰性显著；从图6（e）中可以看出：在氧气与沥青的界面处，随着温度的升高氧气分子浓度逐渐下降，3 种温度之间的浓度差较为接近。这说明温度升高可以促进氧气向沥青膜内的扩散，且氧气分子浓度随温度接近线性变化。通常，增加温度能加剧沥青老化的程度，其中一部分原因为氧气扩散行为的改变。
4.1.3 氧分子数量
热氧化是沥青材料老化的主要原因，氧气分子浓度会影响沥青老化的程度。不同海拔地区的沥青路面的使用环境中氧气的含量不同，因此有必要研究不同氧气分子浓度下沥青组分和氧分子的扩散行为。图7 为不同氧分子数量（200 个，300 个）的 OOA 模型 x、y 和 z 轴 3 个方向上相对分子浓度的结果。从图7 可以看出：沥青组分和氧气相对分子浓度在 x 和y 轴方向上的波动较小，且不同氧气分子浓度间的差异性较小，而 z 轴方向上存在显著不同。从图7（a）可以看出：当氧分子数量为 200 时，沥青质相对分子浓度最大值出现在 30 Å 附近，而 Z-300 曲线的最大值出现在 40 Å 附近。对比两曲线发现，氧气分子浓度增加，沥青质呈现向沥青膜表面移动的趋势。这是因为沥青质分子易于与氧分子发生反应；在图7（b）中，胶质的相对分子浓度变化与沥青质表现出相反的规律，当氧分子数量增加时，接近集料表面的胶质分子相对分子浓度会随之增加；从图7（c）和（d）可知：随着氧分子数量的增加，芳香分浓度的变化规律与胶质一致，饱和分则与沥青质类似；图7（e）为氧分子相对分子浓度曲线。从图7（e）中可以看出：增加氧分子数量，沥青膜内的氧气相对分子浓度整体上升。Z-300 曲线在集料‒沥青界面处出现明显峰值，说明氧气会在界面处聚集，这将提升界面处沥青发生氧化的概率。
图6 不同温度 OAA 模型的沥青组分和氧气相对分子浓度
Figure6 Relative concentrations of asphalt components and oxygen molecules in OAA model at different temperatures
图7 不同氧分子数OAA模型的沥青组分和氧气相对分子浓度
Figure7 Relative concentrations of asphalt components and oxygen molecules in OAA model with different numbers of oxygen molecules
4.1.4 集料类型
目前，沥青混合料生产使用较多的集料类型为石灰岩和玄武岩，两者的主要成分为碳酸钙和二氧化硅。为研究不同集料类型对沥青组分和氧分子扩散行为的影响，本研究选用 CaCO3和 SiO2两种集料的 OAA 模型进行模拟扩散。图8 为沥青组分和氧气相对分子浓度曲线。
从图8 可以看出：集料类型对沥青质和胶质相对分子浓度的影响较为明显，对芳香分和饱和分的影响相对较小。具体来说，从图8（a）可以看出：SiO2曲线在 35 Å 附近出现明显峰值，并且远大于 CaCO3对应的相对分子浓度值。这说明二氧化硅对沥青质的吸附能力较弱，沥青质分子向氧气界面移动；图8（b） 中，CaCO3的胶质相对分子浓度曲线出现最大峰值，但两者均表现为在沥青‒氧气界面处的分子浓度最高；芳香分和饱和分均为沥青的轻质组分，从图8（c） 和（d）可以看出：两者的分子浓度变化规律类似，均表现为在集料界面附近的分子浓度较高。这主要是由于集料的吸附作用所致，从分子浓度值的变化可以看出：SiO2对芳香分的吸附能力更强，而对饱和分的吸附能力弱于 CaCO3；图8（e）为两种不同集料模型的氧气相对分子浓度曲线图。从图中可以看出：在 40~70 Å 范围内，CaCO3 模型的氧气相对分子浓度高于 SiO2，这表明在 CaCO3 模型中更利于氧气向沥青膜内扩散。并且在更接近集料界面的位置，出现了大于 SiO2模型的氧气相对分子浓度峰值，这进一步为前述结论提供了有力的佐证。
图8 不同集料类型OAA模型的沥青组分和氧气相对分子浓度
Figure8 Relative concentrations of asphalt components and oxygen molecules in OAA model with different aggregate types
4.2 沥青组分与氧分子聚集态分析
径向分布函数是指在某个粒子半径为 r 的球壳空间内另外一种粒子出现的概率，可用于表征沥青内分子的聚集和分布状态。本研究分别以集料模型和沥青组分分子质心为中心，计算了沥青四组分和氧分子的 RDF，以此来表征集料表面组分的分布特点和四组分与氧分子间的聚集状态。
图9 和图10 为径向分布函数的结果。
图9 OAA 模型中沥青组分与集料间径向分布函数
Figure9 Radial distribution function between asphalt components and aggregates in OAA model
图10 OAA 模型中沥青组分与氧分子间径向分布函数
Figure10 Radial distribution function between asphalt components and oxygen molecules in OAA model
从图9 可以看出：随着 r 的增大，g（r）值呈现明显的波动，且峰值数量增加，并且不同集料类型沥青组分的径向分布函数结果也存在显著差异。在图9（a）中，半径约 15 Å 处出现了明显峰值分别为 11.66 和 6.09。这表明沥青质组分向集料的位置移动，与沥青质相对分子浓度的结果一致。CaCO3的径向分布函数值约为 SiO2的 2 倍，说明 CaCO3对沥青质的吸附能力更强；从图9（b）可以看出：CaCO3的曲线出现了峰值，而 SiO2 的曲线虽有波动但没有出现明显峰值。这说明胶质分子的移动受 SiO2的影响较小；根据图9 （c）和（d）可知，约在 20 Å 附近，CaCO3模型对应的芳香分和饱和分 g（r）均出现了最大峰值，而 SiO2模型在较大半径波动，说明无聚集现象。综上所述，石灰岩集料对沥青组分的吸附能力较强，这与组分相对分子浓度的结果存在对应关系。从图10 可以看出：除沥青质之外，另外 3 种沥青组分分子的 g（r）均无明显差别。说明胶质、饱和分和芳香分与氧分子之间的相互作用受集料类型的影响较小。比较图10（a） 两种集料类型对应的沥青质径向分布函数结果可知，当半径为 10 Å 时，CaCO3对应的曲线有出现峰值的趋势。因此，认为沥青质分子与氧气分子在石灰岩表面更容易发生聚集，增大了两者发生反应的概率。在图10（b）、（c）和（d）中，两种集料对应的胶质、饱和分和芳香分径向分布函数的结果差异较小，且均未出现明显峰值。
4.3 沥青组分与氧分子的扩散
为了研究沥青组分和氧分子的扩散行为计算均方位移和扩散系数。图11、12 分别为两种集料类型模型的沥青组分和氧分子的均方位移曲线，扩散系数的结果列于表3 中。
从图11 可以看出：不同沥青组分的均方位移曲线存在明显差异，两种集料类型对应均方位移曲线的变化趋势不同。曲线从上至下依次为饱和分、胶质、芳香分、沥青质，这主要受分子质量和数量的影响。在图12 中，随着模拟时间的增加，两种模型中氧分子的均方位移曲线差异逐渐增大。结合表3 中扩散系数值，SiO2模型对应的沥青质和氧分子扩散系数大于 CaCO3模型，而其余分子的扩散系数值均小于 CaCO3模型。由扩散系数的结果说明，沥青质和氧分子更易在 SiO2集料表面扩散，与相对分子浓度结果并不一致。计算扩散系数时，统计了模型内所有的氧分子，根据相对分子浓度的结果发现，在 CaCO3模型中沥青膜内的氧分子数量较多，而这部分氧分子的移动受沥青分子的影响会导致整体氧分子的扩散系数下降。
图11 OAA 模型中沥青组分的均方位移
Figure11 Mean square displacement of asphalt components in OAA model
图12 OAA 模型中氧分子均方位移
Figure12 Mean square displacement of oxygen molecules in OAA model
表3 沥青组分和氧分子的扩散系数
Table3 Diffusion coefficients of asphalt components and oxygen molecules
5 结论
（1）随着模拟时间的增加，沥青质分子从沥青膜表面向中心迁移，而胶质分子从沥青表面向集料界面迁移。由于集料对轻质组分的吸附，在集料界面处，饱和分分子浓度先增加后减小，而在氧气界面处的变化规律为先减小后增加，饱和分的变化规律与芳香分类似，均受集料吸附和轻质组分挥发的双重影响。同时，氧气分子浓度随扩散时间的增加而上升。特别是在集料与沥青界面处，氧气分子浓度上升明显。
（2）低温条件下沥青质分子不易发生迁移。然而，随着温度的升高，沥青质分子表现为从中间向两界面处迁移的现象。胶质分子有整体向集料表面运动的趋势。集料界面和沥青膜表面的饱和分分子浓度变化各异，分别表现为先减小后增大和先增大后减小的变化趋势。温度是影响饱和分迁移的重要因素，低温条件下分子惰性显著，而高温可以促进氧气向沥青膜内的扩散，且氧气分子浓度随温度接近线性增长。
（3）当氧气分子浓度增加时，沥青质分子有向沥青膜表面移动的趋势。同时，接近集料表面的胶质分子相对分子浓度也随之增加。芳香分分子浓度的变化规律与胶质相似，而饱和分则与沥青质的变化趋势相近。氧气在界面处的聚集现象，提高了界面处沥青发生氧化的概率。
（4）集料类型对沥青质和胶质分子相对分子浓度的影响较为明显，但对芳香分和饱和分的影响相对较小。在 CaCO3模型中，氧分子的相对分子浓度高于 SiO2模型，表明 CaCO3模型更有利于氧气向沥青膜内扩散。
（5）沥青质组分向集料位置的移动与沥青质相对分子浓度的变化结果相吻合，进一步证实了碳酸钙对沥青质的吸附能力更强。相比之下，胶质分子的移动受 SiO2 的影响较小。石灰岩集料对沥青组分的吸附能力较强，这与组分相对分子浓度的结果存在对应关系。尽管胶质、饱和分和芳香分与氧分子之间的相互作用受集料类型的影响较小，但沥青质分子与氧气分子在石灰岩表面更容易发生聚集。
