comsol如何定义狄利克雷边界_COMSOL模型分享-2020 Science声子晶体中的Klein隧穿效应...

说明：模数哥免费提供本文模型，点击此处，查看获取方式。

本文导读：

1， 狄拉克线性色散——无质量粒子；

2， COMSOL仿真声子晶体中的Klein隧穿效应；

---

昨天，Science上线了UCBerkeley翔叔组的最新文章：Direct observation of Klein tunneling in phononic crystals，具体内容可见两江科技评论今天的报道：声子晶体中的Klein隧穿效应的直接实验观测。文章的一作江雪，是南大声学所10级校友，是程建春老师、梁彬老师的得意门生，目前供职于复旦大学。模数哥今天看了一下这篇文章，顺手重复了一下仿真模型，放到了我们的QQ群文件里面，感兴趣的童鞋请自取，拿来玩玩。

---

一、狄拉克线性色散——无质量粒子

关于Klein隧穿，两江科技评论的报道里面有专业的说明：有别于非相对论效应的传统量子隧穿现象，具备相对论效应的狄拉克粒子在隧穿效应中会呈现不受势垒宽度和能量高度影响的全透射现象。模数哥缺乏凝聚态物理的功底，就不班门弄斧，歪曲解读了。而实验所需的狄拉克粒子，倒是能带理论中的常客，比如这篇文章也引用了赖耘老师早年的大作，我们也做过仿真说明：光子晶体- 2011 Nat. Mat. 光子晶体随机简并的狄拉克锥和零折射。需要注意的是赖老师这里的狄拉克粒子是在布里渊区中心设计出来的，后来由模数哥的师兄猫哥拓展到了声子晶体领域：声学拓扑-2016 Nat. Phys. 基于能带反转的声拓扑绝缘体(下)，进一步由何程老师通过打破对称性(打开狄拉克锥)来实现了量子自旋霍尔效应：声学拓扑-2016 Nat. Phys. 基于能带反转的声拓扑绝缘体(上)。而江雪他们这里，使用的是三角晶格在布里渊区边界处的狄拉克锥，类似于刘正猷和陆久阳老师他们2016年的谷自旋工作：声拓扑-2017 Nat.Phys.声子晶体中的拓扑谷传输(二)。

二、COMSOL仿真声子晶体中的Klein隧穿效应

如上图所示，我们首先在COMSOL中画一个三角晶格的元胞，这里我们跟原文类似，画成菱形。然后通过分段函数来定义布里渊区Γ-M-K-Γ的路径，使用Eigenfrequency求解器，就可以仿真出其能带(原图Fig. 1c)：

如何具体定义Floquet周期边界，压力声学场(Pressure Acoustics, Frequency Domain, acpr)，我们就不细说了。通过改变中间刚性柱的尺寸，我们可以对K点的狄拉克锥进行自由设计(原图Fig. s1)：

模数哥此时需要停下来做点点评：

1，需要指出的是，实验中的刚性柱比波导稍微矮一点，为了与实验结果相匹配，原文采用的是三维仿真，模数哥这里偷懒采用了2D仿真，所以数值上和原文有些差别。但是道理和方法都是一样的。

2，我们这里的K点看上去不像线性的狄拉克锥，这是由于我们对布里渊区的定义又偷懒了，ΓK的长度应该是MK的两倍。

当我们采用两种不同刚性柱半径R0的声子晶体构造一个异质结时，只要保证声波频率属于两种声子晶体狄拉克锥的线性色散区域内，Klein隧穿就产生了(原图Fig. 2b)：

这里的仿真仍然是使用压力声学acpr，求解步则换成了Frequency Domain，由于要对散射的声波进行吸收，我们在x方向使用了完美匹配层(Perfectly Matched Layer，PML)，采用背景场进行入射激发就非常方便了。图所示声波从左向右传播经过该异质结时，能量上非常的均匀，并没有发生反射现象。并且该现象在很宽的频段内都能保持(原图Fig. 3c)：

文章大概的内容就重复出来了。至于具体深层次的物理原因，欢迎大家细看文章。