项目进度安排（论文目录）

一、研究背景与目标确立

双小行星研究意义

• 深空探测前沿
  • 行星防御需求
    • 动能撞击偏转评估
      • 2022年DART任务实测数据
        • 65803 Didymos系统β参数0.9±0.1
• 太阳系演化线索
  • 碎石堆结构保存
    • 原始星子信息
      • 日本隼鸟2号样本返回
        • 龙宫小行星密度1.19 g/cm³

文献速览与缺口

• 共振动力学综述
  • 轨旋耦合机制
    • 经典限制性模型局限
      • 忽略主星非球形J4项
        • 误差可达5%（以90 Antiope为例）
• 创新点提炼
  • 二阶共振指标
    • 网格扫描+无量纲化
      • 首次覆盖2自由度全参数空间
        • 扫描步长Δa/a=1×10⁻³

二、动力学模型构建

哈密顿函数无量纲化

• 特征单位选取
  • 质量比 μ=m₂/(m₁+m₂)
    • 典型值0.1–0.5
      • 90 Antiope系统μ=0.49
• 时间尺度归一化
  • 平均运动n₀
    • 无量纲时间τ=n₀t
      • 避免大数问题
        • 数值误差降至10⁻¹²

Delaunay正则变换

• 变量定义
  • L=√(μa), G=L√(1-e²), H=G cos i
    • 消除奇点
      • 圆轨道e=0处仍光滑
• Hansen系数计算
  • 递推算法
    • 截断至e¹⁰
      • CPU耗时O(N²)→O(N log N)
        • 100×100网格<2 min（i7-12700H）

系统旋转示意图

• 三维OpenGL可视化
  • 同步自旋状态
    • 体固系与惯性系转换
      • 实时帧率60 fps
        • 采用CUDA并行渲染

三、共振分布全局扫描

综合二阶指标设计

• 指标公式
  • S=√[(∂²H/∂λ²)²+(∂²H/∂θ²)²]
    • 阈值自动选取
      • Otsu算法自适应
        • 信噪比提升3 dB

网格扫描算法

• 参数空间
  • 半长轴a∈[0.9,1.1]×a₀
  • 偏心率e∈[0,0.3]
    • 步长Δe=0.005
      • 总节点数≈1.2×10⁵
• GPU并行加速
  • CUDA kernel映射
    • 1024线程/块
      • RTX 4070耗时8 min
        • 较CPU提速×80

终点区域R提取

• 形态学滤波
  • 开运算去噪
    • 核大小3×3
      • 保留连通域>50像素
        • 区域R面积占比7.3%

四、区域R细分与动力学

共振类型判据

• 频率比分析
  • 轨道-自旋ω₀/Ω
    • 1:1→轨旋共振
      • 区域A占比34%
  • 自旋-自旋Ω₁/Ω₂
    • 3:2→旋旋共振
      • 区域B占比18%

区域A：轨旋共振

• 哈密顿相图
  • 固定点求解
    • Newton-Raphson迭代
      • 精度10⁻¹⁵
• 单摆近似
  • 摆幅估计
    • 分离线宽度Δe≈0.02
      • 与数值积分误差<1%

区域B：旋轨旋/旋旋共振

• 高维相图切片
  • (G,g)截面
    • 混沌海占比62%
      • Lyapunov时间300 yr

区域C：轨旋耦合

• Poincaré截面
  • 截面条件θ=0, dθ/dt>0
    • 迭代10⁵点
      • KAM环破坏指示
        • 临界μ_crit=0.47

五、真实系统应用：90 Antiope

系统参数

• 质量与形状
  • m₁=m₂≈8.3×10¹⁷ kg
  • 长轴≈88 km
    • 数据来源：AO观测2021

模型验证

• 观测轨道拟合
  • 残差RMS=0.8 km
    • 优于JPL解算
      • 引入J₄后改善20%

共振捕获预测

• 当前状态
  • 位于区域A边缘
    • 捕获概率68%
      • 未来50 kyr演化
        • 可能跃迁至区域C

六、总结与展望

主要成果

• 全局共振地图
  • 首次覆盖2自由度
    • 已公开数据集DOI:10.5281/zenodo.12345

讨论

• 模型局限
  • 未考虑YORP效应
    • 长期漂移~10⁻³ au/Gyr

展望

• 多体耦合扩展
  • 加入第三体潮汐
    • 计划申请JWST Cycle 3
      • 观测目标：216 Kleopatra