《Bioengineering》：AI-Driven Robotic PCI: A Perception-to-Action Framework for Coronary Guidewire Control
编辑推荐：
经皮电刺激（TES）受到皮肤不适的限制，这是由于在向深部靶点输送电流时不可避免地激活了浅表感觉神经。虽然心理物理学研究已凭经验确定了能减少皮肤感觉的波形，但现有的计算模型通常采用准静态近似，忽略了生物组织显著的介电色散，导致其生物物理机制尚不明确。研究人员构建
经皮电刺激（TES）受到皮肤不适的限制，这是由于在向深部靶点输送电流时不可避免地激活了浅表感觉神经。虽然心理物理学研究已凭经验确定了能减少皮肤感觉的波形，但现有的计算模型通常采用准静态近似，忽略了生物组织显著的介电色散，导致其生物物理机制尚不明确。研究人员构建了一个人类前臂的有限元模型，该模型包含了皮肤、脂肪、肌肉、骨骼和神经组织（直流至1MHz）的频率依赖介电特性，并结合了基于系统传递函数H(f)的线性时不变信号处理框架。该模型定量复现了Hsu等人报告的波形依赖感觉趋势。研究人员引入了一个穿透比来量化深部与浅表神经的激活程度，并发现在生物组织中所有时变波形均表现出比直流（DC）更低的穿透性，揭示了电流在生物组织中的类似趋肤效应行为。此外，深部和浅表神经的激活在波形参数变化下共变，表明仅靠波形优化无法改善深度选择性。研究表明，电极空间配置提供了一种补充策略：将主动电极靠近目标肌肉神经，同时避开浅表皮肤神经，并结合足够的横向间距，可增强深度选择性。这项工作将心理物理学观察与组织电生理学联系起来，为TES中的波形和电极优化提供了一种计算工具。
**论文解读：基于频率依赖介电特性的经皮电刺激波形与电极优化计算框架**
**研究背景与问题**
经皮电刺激（TES）是一种非侵入性神经调控技术，通过表面电极输送受控电流以调节外周神经、脊髓回路或皮层区域的活动，广泛应用于慢性疼痛管理、中风康复和脊髓损伤康复等领域。然而，TES面临一个根本性挑战：刺激电流在到达深部靶点神经前必须穿过皮肤，在此过程中不可避免地激活浅表感觉神经，引起刺痛、灼烧、瘙痒等不适，直接限制了可耐受的刺激强度并降低治疗效果。如何有效激活深部靶点神经同时减少浅表感觉，是TES领域长期关注的核心问题。
Hsu等人通过心理物理学实验系统比较了单相方波、单相正弦波和调幅（AM）正弦波在前臂上的感知强度，发现恒定直流（DC）产生最弱感觉；增加占空比或频率可降低感觉，但从未低于DC水平；对于双相AM波，更高载波频率和更低振幅产生更弱感觉。这些结果为波形选择提供了重要经验基础，但背后的生物物理机制尚不清楚。现有计算模型通常采用准静态近似或单频率组织电导率值，忽略了生物组织在直流至1MHz频段内显著的介电色散——电导率和介电常数可变化数个数量级，这从根本上改变了不同频率成分的电流分布。当比较包含不同频率成分的刺激波形时，忽略这种色散可能导致预测偏差。
**研究内容与结论**
为弥补这一空白，研究人员构建了一个包含皮肤（角质层和真皮）、脂肪、肌肉、骨骼、血管和主要神经的人类前臂有限元模型，为所有组织赋予频率依赖介电特性（DC至1MHz）。在频域分析中，计算单位电流激励下每条神经的激活函数（activating function，即沿神经轴的电势二阶空间导数）；定义系统传递函数H(f)为激活函数与输入电流之比，该函数完整表征了中间组织的频率滤波效应。在线性时不变（LTI）假设下，通过FFT将任意输入波形变换到频域，乘以H(f)后再通过逆FFT重建时域响应y(t)，取其均方根（RMS）值作为刺激强度的定量指标。该框架通过复现Hsu等人的实验观察得到验证，并通过7×6电极阵列的25种配置检验了电极空间排列对深度选择性的影响。此外，引入“穿透比”来量化不同波形和电极配置下深部与浅表神经的相对激活程度。
本研究的主要结论包括：第一，在相同平均电流下，单相波形产生比DC更高的刺激强度；在相同峰值电流下，DC产生比所有双相波形更高的刺激强度。第二，所有时变波形的穿透比均低于DC，反映出电流在生物组织中的类似趋肤效应行为——时变电流优先集中在浅表组织层，该行为归因于脂肪-肌肉界面的介电弛豫。第三，任何减少浅表激活的波形参数调整都会同时降低深部肌肉神经的刺激强度，表明仅靠波形优化无法解决深度选择性问题。电极空间配置提供了一种补充策略：将主动电极靠近目标肌肉神经并远离皮肤神经，同时采用足够大的横向间距以降低浅表神经上的电流密度，可增强深度选择性。
**主要关键技术方法**
研究人员采用以下关键方法：1）构建基于UGNX 2206的多层圆柱体前臂几何模型，包含皮肤（角质层、真皮）、脂肪、肌肉、皮质骨、骨髓、血管及七条代表性神经（三条肌神经、四条皮肤神经），神经为三层同轴圆柱体（内膜、束膜、外膜）；2）在COMSOL Multiphysics 6.1中利用AC/DC模块的“电流（ec）”物理场接口进行频域有限元仿真（0Hz至1MHz），所有组织赋予频率依赖介电特性（基于Cole-Cole弛豫模型或Gabriel等的数据），肌肉和皮质骨设为各向异性，使用扫掠网格和自由四面体网格混合剖分；3）基于线性时不变系统理论，通过Python程序计算系统传递函数H(f)（定义为激活函数与输入电流之比），对任意时域波形经FFT后乘以H(f)再逆FFT得到时域神经响应，以RMS值量化刺激强度，并定义为穿透比ρ = y
桡神经/y
第二皮肤神经；4）通过7×6电极阵列的25种空间配置（改变主动电极位置及横向、纵向间距）评估深度选择性。
**研究结果**
**3.1 神经的频率响应特性**
电流密度分布显示，电流主要沿高电导率路径（肌肉>真皮>脂肪）流动。所有神经的传递函数幅值响应随频率增加呈近似指数衰减，表明虽然组织电导率随频率略有增加，但神经感知的刺激强度主要由介电常数的急剧下降主导。相位响应从低频接近0°渐近趋向-90°（高频极限），表明所有路径在高频下由位移电流（电容耦合）主导。不同神经的相位曲线呈现三种模式：单调下降（第二皮肤神经）、非单调带一个极小值（桡神经、尺神经、正中神经、第三皮肤神经）、非单调带两个极值（第一、第四皮肤神经），分别对应电流路径是否跨越脂肪-肌肉界面以及是否存在平行通路。实部响应在部分神经上出现暂态负值，归因于脂肪-肌肉界面显著的介电弛豫效应使感应电场方向暂时反向。虚部响应从接近零的负值下降至波谷再回升至零，波谷对应特定神经通路的最大特征介电弛豫频率，应避免在此频率附近刺激以减少能量耗散。
**3.2 实验验证**
通过复现Hsu等人的五组实验验证模型：实验1（5kHz单相方波，占空比10%/20%/30%，基线电流0.5/1/2mA，净电流3mA）——模拟结果精确复现了刺激强度随占空比或基线电流增加而单调下降但始终高于DC水平的趋势，时域响应呈近似三角形状（低通滤波效应），穿透比随占空比或基线电流增加而提高但始终低于DC。实验2（20%占空比单相方波，频率5/10/20kHz）——刺激强度随频率增加而下降，穿透比提高但低于DC。实验3（单相正弦波，1-50kHz，3mA振幅+3mA DC偏移）——刺激强度随频率指数衰减并趋近DC水平，穿透比在约8kHz处出现明显波谷，与桡神经和第二皮肤神经的特征弛豫频率算术平均值（约7.55kHz）吻合。实验4（双相AM波，5Hz调制频率，100Hz-400kHz载波频率，0.5-2mA峰峰值振幅）——刺激强度随载波频率增加而单调下降，随振幅增加而增加，但穿透比单调下降（与单相正弦波的非单调行为不同，因AM波不含DC分量）。实验5（AC波形比较，方波/正弦波/AM波，载波频率500Hz和1kHz，最大振幅3mA）——模拟结果正确再现了方波刺激强度最高，但发现正弦波始终高于AM波，与Hsu等实验结果中的排序存在差异（可能原因待讨论）。所有AC波形在相同最大振幅下刺激强度均低于DC。补充参数研究（固定1mA振幅）显示，方波刺激强度随占空比单调增加，随频率单调下降，穿透比在约8kHz处再次出现波谷，证实其由介电弛豫特性决定而非波形形状。
**3.3 主动电极配置的影响**
电极-神经距离效应：移动主动电极远离皮肤神经可降低其激活函数幅值，但若肌神经恰好位于皮肤神经下方（如正中神经在第一皮肤神经下），则存在解剖矛盾。增大主动电极与返回电极之间的横向间距，可有效降低覆盖在深部靶点上方的皮肤神经的激活，而对深部神经的影响较小，从而解决此矛盾。纵向间距减小仅轻微增加激活函数幅值，过小间距可能增加短路风险。多电极阵列测试（恒定总电流）：横向扩展电极数量（垂直于前臂方向）对皮肤神经和肌神经的影响呈神经依赖性，虽对单个肌神经效果有限，但可增强所有肌神经的累积激活；纵向扩展电极数量（平行于前臂方向）则一致地降低所有七条神经的刺激强度（因肌肉纵向电导率高，增加了电流路径横截面积）。通过对25种电极配置的综合评估（累加所有皮肤神经和所有肌神经的激活函数幅值），确定最佳配置标准为：主动电极靠近肌神经聚集区并远离皮肤神经，同时增大横向间距以扩展覆盖范围。
**总结与讨论**
本研究提出了一个用于优化TES的计算框架，解决了先前模型忽略频率依赖组织介电特性的关键局限。通过将频率依赖介电特性纳入前臂有限元模型，并结合基于传递函数的信号处理管道，建立了能够预测任意刺激波形下神经响应的方法。定量复现Hsu等人的心理物理学发现验证了该方法，并确认频率依赖介电色散是波形依赖神经激活的关键决定因素。三重主要发现：波形优化不能单独解决深度选择性问题；电极空间配置提供补充策略；最佳配置为主动电极靠近靶肌神经、远离皮肤神经并增大横向间距。局限性包括：简化圆柱体模型未捕捉个体解剖变异；激活函数作为线性代理未模拟动作电位起源的非线性动力学，但由于其线性依赖性及固定几何下传递函数的一致性，相对激活强度排序得以保留，足以满足本研究比较性目标。总之，本研究证明了中间组织的频率依赖介电特性是TES中波形依赖皮肤感觉的主要生物物理机制，基于传递函数的框架为筛选刺激波形和电极配置提供了计算高效的工具，并可适用于其他解剖靶点。
**结论部分翻译**：本研究提出了一个用于优化经皮电刺激（TES）的计算框架，解决了先前模型的一个关键局限：忽视频率依赖的组织介电特性。通过将频率依赖介电特性纳入前臂有限元模型，并结合基于传递函数的信号处理管道，建立了一种能够预测任意刺激波形下神经响应的方法。定量复现Hsu等人的心理物理学发现验证了该方法，并确认频率依赖介电色散是波形依赖神经激活的关键决定因素。波形分析得出三个主要发现：第一，在相同平均电流下，单相波形产生比DC更高的刺激强度；在相同峰值电流下，DC产生比所有双相波形更高的刺激强度。第二，所有时变波形均表现出比DC更低的穿透比，反映了电流在生物组织中的类似趋肤效应行为，即浅表组织层优先聚集。第三，减少浅表感觉激活的波形参数调整不可避免地衰减深部神经激活，表明仅靠波形优化无法解决深度选择性问题。电极空间配置提供了一种补充策略：最佳配置将主动电极靠近靶肌肉神经，同时保持与皮肤神经的距离，并采用足够大的横向跨度以降低浅表神经上的电流密度。
