下载的是科学复算包，先认清记录来源
压缩包包含分析脚本、两个输入CSV、出处JSON和四份期望输出。主记录的探测频率是合成天体轨道频率的二分之一、一倍和两倍；每行列出重力Love数的实部、虚部及各自独立标准差。这些是本站生成的合成数据，不是Cassini跟踪残差，不是土卫六Love数时间序列，也不是从论文提取的数据集。
生成真值采用五十千米同密度壳、三点五吉帕弹性模量和十的十四点五次方帕秒黏度。高斯扰动使用固定种子12095，实部标准差零点零零二，虚部零点零零一六。保存数值是复算依据；分析时不重新抽噪声，也不暗中改变输入。
使用Python 3.11或更新版本，不需要第三方库。分发目录内有数据和期望输出子目录，脚本不导入网站代码。结果会写出输入的SHA-256哈希；遇到差异时，先与出处记录比较，再判断物理结果。原论文链接提供近似背景，公开包只分发原创教学材料。
对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 · Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025)
先运行，再逐个看残差
下载materials.zip并解压，在解压目录运行下面的命令。脚本会打印有效刚度、松弛时间、残差平方和以及条件剖面摘要。拟合分量、每项残差和条件剖面都保存为小型文本表格。把重跑结果写到另一个目录，方便保留原始期望输出。
打开fitted_response.csv。每行用观测减预测，再除以列出的标准差；把六个标准化残差平方并相加，期望结果约为四点零零四零二九五。拟合使用两个有效参数，在指定独立高斯模型下留下四个残差自由度。自由度计数并不能验证结构假设本身。
几何刚度与弹性模量的拟合乘积约八点二五一七六亿帕，Maxwell时间约八万八千六百零二点五秒。取真实公转频率与零点零零五偏心率，这份拟合响应预测约一百点七四零吉瓦。它与无噪默认例子略有不同，因为拟合来自保存的带噪记录。
python3 analyze-response.py --input data --output rerun-output对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1
看清求解器最小化的对象
两个拟合坐标分别是无量纲有效刚度和Maxwell时间的对数，指数还原后始终为正。脚本先在固定粗网格中找初始点，再用中心差分计算标准化残差的导数，解一个二乘二线性系统得到局部改正；如果平方和没有下降，就缩短步长。
正规方程来自一个局部步骤：在当前参数附近线性展开残差向量，把长度平方展开，对改正量求导并令其为零。它是数值优化工具，不是新加的物理定律。即使最终残差很小，仍可能有多个内部，因为它们映到完全相同的有效坐标。
真实行星分析可能包含相关无线电数据、伴随参数和系统误差，那时需要完整协方差，不能继续使用这里的对角权重。小实验把噪声模型明确写出，才能在不一致时分别追查记录、符号、物理模型或优化步骤。
对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1
把替代内部逐行构造出来
打开shell_solid_family.csv，壳厚从十千米一直列到五百千米全固态。每个厚度先计算论文中的有限壳几何因子，再用拟合有效模量除以它，得到弹性模量；随后乘以拟合Maxwell时间得到黏度。外部预测在每个频率上保持不变。
带噪拟合的全固态行得到约零点二一七吉帕模量、十九点二四万亿帕秒黏度；五十千米行模量约三点五五吉帕。所有行的残差平方和相同。十千米替代体需要约十八点八吉帕，超出交互控件范围，也没有被声称为合理冰材料。独立材料信息的价值，正在这里出现。
再看conditional_thickness_profile.csv。严格固定三点五吉帕模量后，脚本在每个四分之一千米网格点对黏度最优化。最佳网格点为五十点五千米，局部剖面增量为一的区间为四十九点五至五十一点七五千米。读这个窄区间时，也要读它旁边的条件；放开模量后，恢复的是整个自由族，而不只是稍宽一点的误差棒。
对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 · Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025)
用失败记录检查符号与模型
第二份输入的虚部为正，并远超声明误差；采用正时间虚指数时，被动Maxwell天体无法产生它。脚本拒绝这份故意不相容的记录，并在summary.json写明原因。它不会裁掉响应、为了拟合翻转符号，也不会把它重新解释成炽热内部。
另一篇论文中的正虚数并不自动表示主动材料。作者可能使用相反时间指数，也可能把滞后大小报成正数。比较前应先翻译约定；2025年土卫六论文报告的是正耗散系数，导读会说明转换。这里的失败记录是在符号约定已经固定之后才不相容。
迁移练习可以先把两列误差都加倍再运行：自由族仍严格退化，有条件剖面则会放宽。随后只翻转虚部符号：被动模型应失败，而不是发现另一种海洋。这两次改动分别检验误差和相位究竟提供了什么信息。
对应阅读：Beuthe, Tides on Europa: the membrane paradigm, arXiv:1410.4735v1 · Petricca et al., Titan’s strong tidal dissipation precludes a subsurface ocean, Nature 648 (2025)