基于描述符引导的反应–扩散声子介质中带选择性弹性波频率对比筛选

《Materials Today Communications》：Descriptor-Guided Screening of Band-Selective Elastic-Wave Frequency Contrast in Reaction–Diffusion Phononic Media

　　

基于描述符引导的反应–扩散（reaction–diffusion, RD）声子介质中带选择性弹性波频率对比筛选。反应–扩散（RD）形态能够实现面内弹性波筛选。弹性波各向异性局域于狭窄的分支–半径热点区域。长度尺度和界面密度呈现最清晰的描述符趋势。有限元（finite element method, FEM）计算再现代表性的面内能带趋势。

Hyeon Sik Won等围绕反应–扩散（reaction–diffusion, RD；由反应与扩散耦合生成空间图案的动力学模型）形态与声子超胞（supercell；周期性重复的异质单元）中带选择性（band-selective）弹性波频率对比之间的关系开展了探索性研究。研究背景方面，声子晶体（phononic crystals）和弹性超材料（elastic metamaterials）可通过材料对比、几何和结构层级调控机械波传播；该领域已从完全带隙扩展到等频轮廓（equal-frequency contours）、群速度场（group-velocity fields）、方向输运（directional transport）和频率选择性各向异性（frequency-selective anisotropy）等更广泛的色散工程。与此同时，非晶格和无序波介质受到关注，既有研究涉及无序声子色散、超均匀无序声子结构中的带隙与波导、软超均匀弹性超材料中的振动集中、光子系统中的短程有序与结构关联，以及RD生成的Turing图案（Turing patterns）与超均匀性之间的联系。问题在于，这类非晶格介质不能由少量晶体学参数自然描述，哪些形态特征实际支配弹性波各向异性仍难以识别。RD系统特别是Gray–Scott模型（Gray–Scott model；由两个控制参数生成丰富可重复形态的双变量反应–扩散模型）可由低维参数变化产生斑点、条纹、迷宫和混合过渡态，因而可作为结构化波介质的生成设计空间。研究人员因此将二维Gray–Scott图案族映射到共同Bloch波（Bloch wave；周期介质中的本征波）框架下的异质弹性介质，比较形态描述符（descriptor）与面内弹性能带计算得到的频率对比度量。研究结论是：主波长（dominant wavelength）和界面密度（interface density）在全49例数值数据集中具有最大的绝对相关系数，但原始场振幅审计、均匀各向同性参考和互易基研究说明这些关联必须保持探索性；报告的最大值只是规定扫描上的频率排序统计，而非已证明的宽带性质或物理各向异性最优。其意义在于提出描述符引导的探索性计算筛选流程，并明确后续形态引导波设计需要接受的原始输入、一致材料映射、数值解析且物理可解释的响应以及单独输运验证，而非器件级验证。

方法概括（不超过250字）：研究人员采用二维Gray–Scott反应–扩散模型在7×7的(F,k)参数网格上生成128×128场，归一化并阈值化为两相弹性介质，赋予密度1200与2600 kg·m^?3及两组Lamé参数(λ, μ)=(1.5, 0.7) GPa和(5.0, 2.8) GPa；计算主波长、界面密度、谱熵、环滤波谱各向异性、向列序、连通性评分、骨架长度密度、迷宫指数和成熟度评分等描述符；采用二维面内弹性平面波展开（plane-wave expansion, PWE）进行Bloch波计算，扫描前六个频率排序分支、六个归一化波数半径和36个角度；保留代表性COMSOL Multiphysics有限元法（FEM）能带沿Γ–X–M–Γ比较；另做固定探针PWE基组细化至NG=20、原始振幅审计和均匀各向同性参考。本研究无样本队列来源。

研究结果如下。3.1 代表性形态区域（Representative morphology regimes）：研究人员以三个历史选择案例(F,k)=(0.038,0.050)、(0.038,0.060)、(0.042,0.060)说明从形态到波响应的探索流程，历史标签分别为“类各向同性成熟”“最强波各向异性”和“迷宫主导”；原始选择先采用成熟度分数（maturity score）阈值M≥0.35，再取保存的SH频率对比最小值和最大值及最大迷宫指数，并非面内最优搜索。第一例原始范围仅约1.67×10^?15，数值退化，不能作为成熟各向同性对照；另两例原始范围较大，仅作为固定数值例子，不声称是面内材料各向异性收敛最大值。3.2 描述符图（Descriptor maps）：在7×7参数空间上，环滤波谱各向异性（ring-filtered spectral anisotropy）、主波长/a、界面密度、向列序（nematic order）、连通性（connectivity）和迷宫指数（labyrinth index）呈现多描述符图像而非单参数转变；23个原始输入被数值拒绝，2个弱振幅，24个为暂定输入；有限描述符值不能证明图案已解析或成熟，且取向/拓扑描述符热点不完全重合。3.3 面内角度频率对比图（In-plane angular frequency-contrast maps）：通过共同分支–半径搜索，49个材料图的最大采样对比局域于(F,k)网格的受限区域而非均匀分布，半径图记录最大值出现位置；筛选输出依赖分支和半径，不能据此建立工作带宽或全色散面响应，仍需原始场有效性、基灵敏度和超胞带折叠检查。3.4 描述符与面内频率对比相关性（Descriptor correlations with in-plane frequency contrast）：最大绝对Pearson相关系数出现在主波长（r=?0.57）和界面密度（r=0.54）；环滤波谱各向异性（r=0.02）、连通性（r=?0.03）关联很小，向列序（r=0.22）和迷宫指数（r=0.28）较小；这些是49个数值输出（含23个被拒绝输入）的描述性关联，移除23个后约为主波长?0.70和界面密度0.65，但仍有2个弱振幅输入；该分析未控制相分数、未分离相关描述符，也未验证响应统计量就是材料各向异性，因此不是预测层级或唯一散射机制证据。3.5 面内PWE–FEM能带结构比较（In-plane PWE–FEM band-structure comparison）：FEM采用相同组分常数和案例坐标并沿Γ–X–M–Γ施加Floquet周期边界，但FEM将128×128原始浓度场双线性缩放到40×40、重归一化并再次阈值化，所得材料图与保存的PWE阵列不同；例如(0.038,0.060)的相2分数从0.500061变为0.550625，(0.042,0.060)从0.606262变为0.651250。保留的代表性PWE输出采用NG=3，FEM胞边长为0.4 m且归档比较将FEM频率乘0.4以在数值a=1参考长度下表达，但这不校正重映射材料场；因此该比较不能建立匹配输入的PWE–FEM一致，也不能证明族趋势独立于离散化。3.6 代表性面内热点极化特征（Polarization characteristics of representative in-plane hotspots）：保存的NG3代表性扫描角度调制相近，选定对比约0.2590、0.2767、0.2752；第一例因原始场数值退化不能作成熟各向同性对照。最大化频率排序分支为n^?=4、η^?=0.75，平均纵向傅里叶位移分数（longitudinal Fourier-displacement fraction）较低，仅支持“横向主导傅里叶位移”的有限描述；该分数未按异质密度加权，不是动能分配，Γ处刚性零模极化未定义，近简并时本征向量依赖所选本征空间基，故不推断纯物理模式身份或收敛极化剖面。3.7 形态引导弹性波设计启示（Implications for morphology-guided elastic-wave design）：RD场提供参数化内部纹理空间，周期重复使其成为声子超胞，可在同一框架中评估几何描述符与弹性本征频率；主波长和界面密度具有最大绝对相关，其他描述符描述互补特征。但当前数据集不能据此优先选择物理各向异性候选；分支和波数依赖有助于组织后续计算并区分选定数值区域与宽带性质，仍需区分超胞带折叠（band folding）与材料效应并识别相关物理模式。该工作流是探索性计算筛选，不是器件级验证；7×7扫描是结构化输入集，保留的归一化图像包括未通过原始场检查的案例；推广需要接受的图案及独立实现、阈值和胞尺寸研究，并需等频轮廓、模式与能流信息、有限域频率响应或传输计算连接带级统计与可测输运。文中未报告制造或实验测量；组分局部各向同性且互易，极性敏感或Willis型描述需额外均匀化或富集连续介质建模；与机电本构模型耦合的扩展需材料耦合系数、电热边界条件、相关偏置平衡和新波计算，功能梯度本身不意味着主动控制。3.8 数值检查和局限（Numerical checks and limitations）：原始振幅审计拒绝跨度不大于100 eps(max(1,max|H|))的场，其中22例仅在舍入尺度变化，另有(0.018,0.060)最大原始值约6.1×10^?244，共23例失败；另有约2.48×10^?11和1.41×10^?9两个弱振幅跨度，剩余24个暂定输入高于示例性10^?8截断，但该截断和成熟度分数均不证明物理图案形成。频率对比统计还与均匀各向同性参考比较：在分支4、η=0.75时，两相解析折叠谱在相同36个角度上给出A_f=0.2770666，说明大对比可由超胞折叠和频率排序产生而非材料各向异性；这不意味着所有异质角度变化都是折叠伪影。固定探针PWE基组序列采用NG=2,4,6,8,10,12,16,20、前六个频率排序频率和12个波矢（η=0.15,0.35,0.75；θ=0°,30°,45°,90°），相邻基变化取绝对频率差与较低NG频率点wise比值最大值；最终NG12→16与NG16→20变化分别为1.9441%/0.3321%、1.5911%/0.5283%、0.4747%/0.5137%，仅诊断案例(0.022,0.047)满足两次最终变化均低于1%；早期小变化后出现更大变化，说明粗基一致不足。诊断案例原始跨度较大但不满足原成熟度分数过滤，是固定比较而非认证替代代表；频率比较仅在规定探针上建立一致性，不证明角度极值、跟踪模式身份或族排名收敛，且未开展新RD演化、全族角度重扫、COMSOL求解或实验。

讨论部分总结：研究人员强调RD衍生场提供可参数化的内部纹理生成空间，周期重复后可作为声子超胞；描述符图与圆形PWE扫描可按几何特征和选定分支–半径响应组织保存输入族。主波长和界面密度虽显示最大绝对相关，但原始场失败、弱振幅、均匀各向同性参考的带折叠以及基组敏感性使这些相关只能作为探索性数值诊断，不能形成物理预测层级或材料各向异性排序。分支–半径依赖的筛选输出有助于组织后续计算，但必须与超胞折叠、物理模式识别和直接波计算区分；后续设计要求接受的原始输入、一致材料映射、数值解析且物理可解释的响应以及单独输运验证。

研究结论翻译：本研究提出了一种探索性描述符引导工作流，用于Gray–Scott反应–扩散超胞中的带选择性弹性波频率对比。描述符图和圆形PWE扫描按几何特征及选定分支–半径响应组织保存的输入族。所报告的最大值是规定扫描上的频率排序统计，而非已证明的宽带性质或物理各向异性最优。在整个49例数值数据集中，主波长（r=?0.57）和界面密度（r=0.54）具有最大的绝对相关；其余报告的系数较小。附加检查说明了这些关联为何必须保持探索性：23个原始场未通过数值振幅守卫，另有两个场振幅较弱，且一个均匀各向同性超胞通过带折叠给出A_f=0.2771。剩余24个场为暂定候选，当前相关结果不能推出物理描述符层级。新增固定探针PWE研究量化了三个保存材料图至NG=20的基灵敏度。最终NG16→20最大变化为0.3321–0.5283%，但仅诊断案例满足规定的两次连续细化准则；角度极值和物理排名仍未验证。保留的PWE–FEM比较使用不同材料图，不能提供匹配验证。综合而言，原始探索输出和附加检查界定了后续形态引导波设计的要求：接受的原始输入、一致的材料映射、数值解析且物理可解释的响应，以及单独的输运验证。