忆阻微机电系统中准周期与混沌动力学的表征：参数诱导转变与热敏调制

《Journal of Computational Electronics》：Characterization of quasi-periodic and chaotic dynamics in a memristive MEMS: parameter-induced transitions and thermo-sensitive modulation

　　

： 本研究提出了集成惠普忆阻器（HP memristor）并与电学与机械谐振器耦合的微机电系统（MEMS）的数学建模与数值表征。基于线性漂移惠普忆阻器模型，通过李雅普诺夫指数（Lyapunov exponents）、分岔图（bifurcation diagrams）、庞加莱截面（Poincaré sections）、递归图（recurrence plots）、相图（phase portraits）和傅里叶谱（Fourier spectra）研究了非线性机电动力学。数值结果表明，忆阻器内在参数，即低阻态Ron、高阻态Roff、氧化物厚度D和离子迁移率μv，直接控制机电振荡的振幅、频率成分和动力学状态，诱导准周期与混沌行为之间的转变。忆阻参数随温度的变化可调制振荡幅度以及准周期和混沌状态的分布，突出了所提出架构的热敏特性。研究提供了与基于约瑟夫森结（Josephson junction）的MEMS的定性比较，强调惠普忆阻器展现出类似的非线性动力学特征，同时在室温下完全可运行。这些结果表明，忆阻热敏MEMS构成了自适应非线性振荡器、热敏传感设备和混沌辅助机电应用的有前景平台。

论文解读

一、研究背景与问题

自蔡少棠（L. Chua）于1971年预言忆阻器这一缺失的电路元件以来，忆阻器因其记忆特性和非线性特性受到了广泛的科学关注，并已在理论上、数值上和实验上得到验证。忆阻器件凭借其微纳尺寸和非线性阻变特性，被认为是微机电系统（MEMS）中极具潜力的候选元件。MEMS本身因其超小尺寸和对环境扰动的强敏感性，在高精度传感与驱动应用中具有独特优势。

研究人员的先前研究曾采用约瑟夫森结（Josephson junction）开发无源和有源热敏MEMS执行器与传感器。然而，基于约瑟夫森结的架构存在一个主要局限：其运行需要严格的超导条件，尤其是低温环境，这严重限制了其在环境条件下的实际应用。尽管如此，约瑟夫森结的显著非线性和动力学特性已被广泛研究和充分记录。

受上述局限性的驱动，本研究探索将惠普（HP）忆阻器集成到热敏MEMS架构中，作为在标准环境条件下运行的替代非线性记忆元件。选择HP忆阻器模型的原因在于：它提供了一个简单、具有物理依据且被广泛采用的框架来描述非线性记忆效应，同时可在标准室温条件下运行。与通常需要低温环境维持超导性的约瑟夫森结器件不同，HP忆阻器在室温下完全可运行，使其适用于热敏MEMS应用。此外，HP模型捕捉了忆阻器的基本特征，包括状态相关电阻和收缩滞回行为，同时保持数学复杂性可控。其紧凑的公式便于与热-电-机械方程耦合，并允许进行高效的分岔、稳定性和李雅普诺夫分析。

本研究所研究的机电系统由一个HP忆阻器与微梁谐振器和电气振荡器r₀L₀C₀耦合组成。在该配置中，忆阻器作为核心非线性元件，影响振荡幅度、共振频率和动力学转变。该工作的主要贡献在于开发了一种能够在耦合系统的机电响应中展现复杂非线性动力学的热敏MEMS架构。与仅用于混沌生成的简化忆阻器电路相比，本系统旨在建模一个物理耦合的热-电-机械平台，可能适用于自适应传感、非线性驱动和混沌辅助MEMS应用。该论文发表在《Journal of Computational Electronics》。

二、关键技术方法

研究人员在数值研究中采用了以下主要技术方法：基于基尔霍夫定律（Kirchhoff's laws）、电路网络关系和经典欧拉-伯努利梁理论（Euler–Bernoulli beam theory）建立系统的微分方程组，推导出无量纲化的六维动力学系统；采用Python中scipy.integrate.solve_ivp实现的龙格-库塔RK45（Dormand–Prince）方法进行数值求解；使用经典贝内廷（Benettin）重规范化算法计算最大李雅普诺夫指数；通过分岔图、庞加莱截面、递归图、傅里叶幅度谱、相图和时间演化分析表征系统动力学状态。研究中采用的物理参数包括：梁长L=200 μm，磁场强度B=1 T，初始条件为α₀=0.2，x₀=0.05，θ₀=0.05，y₀=0.35，γ₀=?1，z₀=0.25。

三、研究结果

3.1 HP忆阻器关键参数的物理描述

研究人员详细描述了HP忆阻器的物理结构，其总电阻由两个不同材料区域的串联连接决定：导电掺杂的TiO_2?x区域（厚度为w(t)）和绝缘未掺杂的TiO₂区域（厚度为D?w(t)）。R_on和R_off的值仅由材料的固有电子电阻率、器件几何形状（厚度D、表面积S）、电子载流子密度和导电机制决定。离子迁移率μ_v描述了氧空位在电场下的漂移速度，仅参与动力学方程中边界位置w(t)的时间变化率，而不参与瞬时电子传导。因此，在理想的线性HP模型中，R_on和R_off不依赖于μ_v，而μ_v仅控制离子迁移的动力学。

3.2 无热效应忆阻MEMS的分岔分析与非线性动力学

研究人员首先分析了系统主电流随R_on和R_off的分布。结果表明，对于非常低的R_on值，忆阻器中流过最大电流，而流过电容C₀的电流i_b最小。随着R_on增加，流过梁和给电容C₀充电的电流i_b也增加。当R_off在10 kΩ至15 kΩ之间时，流过梁的电流i_b对R_off/R_on比值呈现近似线性依赖，这是HP忆阻器周期性振荡的首选工作范围。

3.2.1 忆阻MEMS中的混沌与准周期行为

3.2.1.1 正李雅普诺夫指数区域中交替的机电动力学状态

通过分岔图和最大李雅普诺夫指数图揭示了两个主要动力学区域：第一个是中心区域，以混沌和准周期振荡交替为特征，其宽度随R_on增加而逐渐减小；第二个区域关于中心轴近似对称，主要与准周期机电振荡相关。在正李雅普诺夫指数区域选取参数点后，忆阻器电流表现出不规则和非周期的时间振荡，庞加莱截面呈现高度分散和无结构的点分布，证实了混沌机电振荡的存在。时间演化显示忆阻器电流呈现相对有组织的振荡区间与不规则幅度调制交替出现的现象，庞加莱映射中局部集中的轨迹与高度分散的点分布交替，表明准周期与更复杂的混沌机电动力学之间存在交替转变。这种交替可能与有效忆阻电阻R(t)的时间演化有关，其动力学持续调制非线性机电耦合。

3.2.1.2 均匀李雅普诺夫指数域中的准周期动力学

在准周期区域选取参数点后，递归图在所有情况下都表现出主导主对角线，确认了确定性行为。随着R_on增加，振荡逐渐变得更加结构化：递归图显示出明确界定的次级对角线和更精细的线条模式，傅里叶幅度谱变得更加集中于主导频率附近，反映了向更相干和准单色振荡的转变。这些观察表明，增加R_on促进了更大的动力学规律性，减少了频谱离散并加强了忆阻器电流的准周期特征。

对于流过机械和电气谐振器的电流，相空间图在所有考虑的R_on值下都表现出双周期且高度规则的振荡轨迹。电荷γ和充电电流γ?～i_b均随R_on成比例增加。随着R_on增加，忆阻器变得更加电阻性，因此吸取更少的电流，导致更大比例的电气能量被转移到机电谐振器。傅里叶幅度谱揭示了单一主导共振频率，该频率基本不随R_on变化，确认系统在稳定谐振状态下运行。在较高R_on值下，频谱内容变得更干净，高频噪声成分减少。

四、讨论与结论

该研究通过定性比较分析表明，HP忆阻器展现出与约瑟夫森结类似的非线性动力学特征，同时保持室温可运行性，克服了约瑟夫森结架构需要低温环境的根本限制。忆阻器内在参数(R_on, R_off, D, μ_v)可直接控制机电振荡的幅度、频率内容和动力学状态，诱导准周期与混沌行为之间的转变。忆阻参数随温度的变化可调制振荡幅度以及准周期和混沌状态的分布，突出了所提出架构的热敏特性。这些结果表明，忆阻热敏MEMS构成了自适应非线性振荡器、热敏传感设备和混沌辅助机电应用的有前景平台，为在室温条件下实现复杂非线性动力学提供了新的可能性。