B2 ErAg与ErAu金属间化合物的弹性各向异性及方向依赖超声衰减：热声效率、敏感性分析与GPR代理建模

《Results in Engineering》：Elastic Anisotropy and Directional Acoustic Response of B2 ErAg and ErAu Intermetallics: A Thermoacoustic Perspective

　　

原文未提供独立摘要段落；以下依据高亮要点进行专业翻译：研究人员分析了B2 ErAg与ErAu金属间化合物中方向依赖的超声衰减（ultrasonic attenuation）。沿?111?方向的ErAg表现出最低衰减和最高热声效率（thermo-acoustic efficiency）。敏感性分析表明，剪切波速度（shear-wave velocity，v_S）主导超声衰减。Monte Carlo（MC）分析证实该声学排序在所有方向上的稳健性。物理信息高斯过程回归（GPR）可准确插值温度依赖的超声衰减。

论文主体解读

稀土与贵金属金属间化合物可作为研究键合、弹性各向异性（elastic anisotropy，弹性响应随晶体方向变化）与声学响应关系的模型体系。已有研究讨论了金属间化合物延性、B2铝化物内聚与电子性质，以及部分稀土B2化合物在特定成分和变形条件下的韧性；Er基化合物、CeAg、稀土–金相形成及镧系–金电子结构提供了相关结构背景。ErAg与ErAu的直接第一性原理基准为Sekkal et al.采用密度泛函理论（DFT，电子结构第一性原理方法）和准谐Debye处理（quasi-harmonic Debye treatment，含温度效应的Debye近似）的研究；相关超声模型也用于稀土化合物等。但原文指出，这些文献不能构成ErAg或ErAu独立的衰减验证。因此，研究人员要解决的问题是：在温度依赖弹性张量基础上，从热声学视角分析B2 ErAg和ErAu的弹性各向异性、方向依赖声速、标量Debye热性质及超声衰减诊断，并评估敏感性、统计稳健性和代理插值。研究人员采用二阶/三阶弹性常数输入，经力学平均、Christoffel方程（Christoffel equation，由弹性张量和传播方向求声速本征值）、标量声学Debye重建和衰减诊断开展计算。结论是：两化合物在0–500 K满足立方弹性稳定不等式，300 K模量分别为ErAg B=52.937 GPa、G=31.041 GPa、E=77.896 GPa，ErAu B=46.730 GPa、G=27.033 GPa、E=67.988 GPa；Pugh比（Pugh ratio，B/G，经验韧性指标）接近常规阈值，但不建立宏观断裂或塑性；纵向速度最大沿[111]，而方向反立方平均速度v_m最大沿[100]；Debye温度Θ_D分别为201.815 K和165.850 K；衰减分支和A_Σ及条件声学比在[111]呈诊断最优，但受归一化和参数限制，不能验证低损耗取向。其意义在于给出方向依赖热声学诊断框架及其可复现性边界。

关键技术方法：以二阶和三阶弹性常数（SOECs/TOECs）为输入，采用Voigt–Reuss–Hill（VRH，多晶弹性模量上下界平均）平均计算力学模量；由Christoffel方程求方向依赖声速，并用全立体角48点Gauss–Legendre×96方位角积分的标量声学Debye模型重建热容与热能量；将频率归一化衰减分量A_th、A_L、A_S1、A_S2的分支和A_Σ作为诊断量；用条件归一化敏感性、20,000次Monte Carlo（MC，随机扰动统计模拟）扰动总结分析；以高斯过程回归（GPR，非参数代理插值模型）和随机森林重要性评估代理插值。无实验样本队列；采用已有报道的0 K、0 GPa固定参考密度（ErAg 10,784.51 kg m^?3，ErAu 13,104.65 kg m^?3）。

研究结果

3.1 弹性常数与稳定性：由所列二阶和三阶弹性常数，两化合物的三个SOECs在0–500 K降低；C₁₁为[100]波纵向刚度。立方弹性稳定不等式C₁₁?C₁₂>0、C₄₄>0和C₁₁+2C₁₂>0在温度网格成立，但原文强调这些条件不建立动力学或热力学相稳定。与DFT比较的0 K偏差为ErAg的C₁₁、C₁₂、C₄₄分别+25.30%、?2.09%、+1.07%，ErAu分别?22.36%、?39.99%、+44.42%。

3.2 机械性能：由二阶弹性常数经VRH平均得到力学性质。300 K时ErAg为B=52.937 GPa、G=31.041 GPa、E=77.896 GPa、B/G=1.705、ν=0.255；ErAu为B=46.730 GPa、G=27.033 GPa、E=67.988 GPa、B/G=1.729、ν=0.258。两者Pugh比略低于1.75，Cauchy压力（C_P=C₁₂?C₄₄）小且为负（?0.31和?0.06 GPa）。这些弹性指标接近常规脆–韧阈值，但不独立建立宏观断裂或塑性变形行为。

3.3 弹性表面：由柔度张量S=C^?1计算方向杨氏模量和剪切模量；剪切模量还依赖满足n·m=0的单位横向方向，并在每个加载方向上对全部允许横向偏振取最小和最大剪切模量。结果表明方向模量面及各主平面截面被给出；立方对称下体模量为标量并另行报告。该部分由方向模量计算展示弹性各向异性。

3.4 Christoffel速度：由Christoffel本征问题在固定参考密度下得到100–500 K声速。300 K时ErAg沿[100]为(v_L,v_S1,v_S2)=(2.786,1.874,1.874) km s^?1，沿[110]为(3.022,1.874,1.462) km s^?1，沿[111]为(3.097,1.611,1.611) km s^?1；ErAu相应纵向速度为2.350、2.578和2.650 km s^?1。两化合物v_L均沿[111]最大，而方向反立方平均速度v_m在300 K沿[100]最大。结论是最大纵向速度方向与最大方向平均速度方向不同。

3.5 标量Debye温度与热学量：对三个声学分支进行全立体角平均，300 K标量Debye速度v_D为1875.545 m s^?1（ErAg）和1585.970 m s^?1（ErAu），对应Θ_D为201.815 K和165.850 K。谐波晶格热容C_V,m为48.776和49.133 J mol^?1 K^?1，排除零点能的热能U_m为11,527.404和12,091.252 J mol^?1。在有效各向同性动力学近似κ_lat=C_V,Vv_D²τ/3下报告τ/κ_lat；该关系不独立建立分支分辨声子弛豫时间或绝对Akhiezer衰减（Akhiezer attenuation，声波与热声子相互作用引起的损耗）。

3.6 报告的衰减分量：300 K衰减分量及其和A_Σ=A_th+A_L+A_S1+A_S2被列出，分量和与所列分量一致。A_Σ是不同偏振贡献的分支和诊断量，不是单一传播模式的衰减。绝对物理解释仍受未解决的热归一化和弛豫时间设定限制。

3.7 条件声学比：条件声学比API在300 K于两化合物中均以[111]最大：ErAg为0.10638×10²¹ m² s^?3，ErAu为0.05204×10²¹ m² s^?3。ErAg/ErAu比值沿[100]、[110]、[111]约为1.915、1.785和2.044。原文强调这些代数比较保留A_Σ局限，不建立经验证的低损耗取向；因热导归一化未解决，未评估热声比。

3.8 条件敏感性与Monte Carlo总结：由归一化分量缩放表达式得到条件敏感性。六种化合物–方向情形平均后，最大绝对系数对应v_S1（1.713），E₀和τ均为0.982，v_S2为0.909，v_L为0.413，κ为0.018。这些值刻画所述偏导数，不是独立物理输入的因果排序。Monte Carlo总结基于20,000个报告扰动样本，描述所述扰动计算，不是实验误差条，也不是对条件声学比的不确定度。

3.9 代理插值及其局限：30对模型生成目标与留一预测给出R²=0.98711、MAE=2.57227、RMSE=5.67049；同系列分段线性留一比较器给出R²=0.96940、MAE=3.58378、RMSE=8.73663。在所列模型生成目标中，A_Σ从300到500 K沿[100]和[110]降低、沿[111]升高。这些趋势描述所报告的模型输出；热学和衰减规格不完整，不能赋予独立建立的损耗机制。训练矩阵和拟合模型规格不可得，不能验证GPR验证协议或外推到独立材料测量。

讨论与结论总结：原文将衰减分量、Monte Carlo总结和机器学习目标–预测对解释为模型诊断；在所列目标中A_Σ随温度的方向趋势不能建立独立验证的损耗机制。热导、弛豫时间和衰减归一化的不一致，以及原子势和机器学习过程规格不完整，限制其物理阐释和可重复性。

研究结论翻译：B2 ErAg和ErAu的温度依赖弹性张量在300 K给出ErAg体、剪切和杨氏模量为52.937、31.041和77.896 GPa，ErAu为46.730、27.033和67.988 GPa。两个张量在所研究温度网格内满足立方弹性稳定不等式。其Pugh比位于常规阈值附近，但不建立宏观断裂或塑性变形行为。计算采用已有报道的参考密度。Christoffel计算显示纵向速度沿[111]最大，而方向反立方平均速度在300 K沿[100]最大。标量Debye重建给出ErAg的Θ_D=201.815 K、ErAu的Θ_D=165.850 K；这些结果使用固定参考密度和谐波双原子模型。补充信息中的衰减分量、Monte Carlo总结和机器学习目标–预测对作为模型诊断解释。在所列目标中，衰减分支和从300到500 K沿[100]和[110]降低、沿[111]升高；这些趋势不建立独立验证的损耗机制。热导、弛豫时间和衰减归一化的不一致，以及原子间势和机器学习过程规格不完整，限制其物理阐释与可重复性。