实现负频率_声子晶体- 2018 Nature 外尔晶体中声表面波的拓扑负折射-CSDN博客

本文链接：https://blog.csdn.net/weixin_39607474/article/details/112734530

本文导读：

1，模数哥按；

2，猫哥指点江山，膜拜一下国内声学诸侯；

3，COMSOL仿真三维声子晶体中的拓扑边界态；

模数哥按：我们的公众号创建半年多了，鼓励大家“原创复制品”并投稿也很久了。在模数哥恩威并施之下，今天终于迎来了第一份投稿，贡献了我们公众号第二篇声学类帖子（第一篇是Acoustics- 2018 PRL Comsol仿真单向“声聚焦”和“声透明”），作者是我的学长&师兄&老弟：小熊猫。小熊猫多年来一直关注于人工带隙材料的理论研究，尤其在光子晶体、声子晶体能带计算方面颇有建树，发表了多篇不错的论文。猫哥为人色而不淫，黄而不奸，一直是模数哥的role model。有意购买此原模型的朋友请直接联系我，或者加入我们的QQ群366619650，寻找小熊猫。模数哥非常鼓励和欢迎大家将自己拿手的仿真作品拿出来和大家分享或者交易，Comsol科研助手就是您的舞台：Comsol科研助手招聘科研助理。

今天猫哥和大家聊聊拓扑声学。要知道在国内，以第一单位发Nature、Science难如登天。当然最近几年情况好了很多，也侧面印证了中国科研的崛起嘛。不过今天要介绍的这篇工作，可以说是大陆高校第一篇声学Nature正刊，为什么强调大陆呢，因为通讯作者刘正猷老师之前在港科还发过一篇Science。刘老师手下Nature子刊、PR系列一篇接一篇，之前让很多国内的声学同行大跌眼镜：声学这个传统工程学科还能发Science，PRL？以猫哥的理解，这是因为搭上了光子晶体，光学超材料，非厄密系统，以及拓扑绝缘体的和谐号啊。这些热点领域，近十几年来以电磁波为载体的NS杂志文章层出不穷，声波作为另一个重要的counterpart当然也可以热闹一把。

国内做声学的院校其实不多，而像港科大，武大，南大，华南理工，同济这些学校能拿声学发大文章，套路总的来说就是凝聚态+声学（模数哥按：做科研怎么玩概念，参加往期例文非线性光学- 2015 OL 仿真Metasurface中的光学双稳）。如果非要在上述高校中排个第一，猫哥觉得首推港科大，那边真是人才济济，猫哥本人也是心向往之。多少人在那边进修一圈回来另起炉灶，玩的是风生水起。

今天猫哥要介绍的，就是今年八月新鲜出炉的Nature雄文一篇，Topological negative refraction of surface acoustic waves in a Weyl phononic crystal，点击文末阅读原文获得连接。光看这个标题，拓扑、负折射、surface wave Weyl crystal, 都是妥妥的时下热点，这几个热点凑一起，再加上有实验，基本就是一区文章的标配了。两江科技评论第一时间也做了相关报道：发现拓扑表面态负折射效应——《自然》刊登武汉大学刘正猷-邱春印团队最新研究成果。

在用Comsol克隆这篇大作之前，先简单介绍下这几个概念：

1，声子晶体：能带来源于周期势对电子的散射问题，那么材料参数或者结构参数的周期性变化对声波的散射，也可以产生能带，这就是声子晶体了。顺便多说一句为什么刘正猷老师两千年的声子晶体工作能发在science上呢，因为之前的散射都是布拉格散射，刘老师在这篇文章中提出了局域共振的新机制（模数哥按：最近这个概念被用到了metalens的设计中，后面会发帖介绍）。

2，负折射：顾名思义，折射是负的。这个负的折射，不代表有负的折射率，里面的水也是深似海，当然负折射率肯定会有负折射。实现负折射的方法也是层出不穷，总的来说，离不开调控色散关系，而这一点，超材料和声子晶体都很擅长！

3，上述两个概念是怎么和拓扑扯上关系的呢？这个就说来话长了。粗浅地讲，一维SSHmodel预言了界面态，二维拓扑会有边界态，三维的自然就有表面态了。装逼地讲，二维系统考虑单个Dirac点，打开它会产生非平庸拓扑。而三维系统，可以把一个Dirac 点拆成两个Weyl点，这两个Weyl点都是非平庸的，都带Charge，一个是N，一个就是负N。想象在动量空间有一个切面经过某个Weyl点，另一个切面经过另外一个weyl点，那么参考二维情形，两个切面之间就有手性的边界态了。这就是最近很火的Weyl point和Fermi arc。而这篇工作的重点在于，发现了这个Fermi arc 可以用来做负折射，而且是没有反射的负折射。

说了这么一大堆，赶紧让COMSOL飞一会。首先第一步，先算能带，找出Weyl Point。

这是个三维的系统，网格划分比较有学问，猫哥在这里只是粗浅地划分了下，结果和文献差不多，不过再精细地画，还是有一点点区别的，可见之前的精度还不够，但是考虑到时间成本，这个网格也够用了。值得一提的是，这里猫哥采用的unit cell和文章里展示的有一点点差别，我认为这种差别在comsol的计算里还是很有必要的（细节省去一万字），结果如何有图为证（原文Fig.2g）：

内行人一眼就看出图中有好几个Weyl point（是什么？在哪里？有什么用？），这是由系统的一些对称性决定的（废话么这不是）。具体细节可以参考文章的一篇参考文献。

我们在声子晶体系统里面找到Weyl point后，下面就要算Fermi Arc了。这篇文章里一共算了三个表面的Fermi Arc，所以我们首先构造好这三个supercell，这里以正文中的YZ界面为例（原文Fig.2j）：

从右下角场图可以看出，我们可以算出两个表面态。至于为什么我们算出两支，因为我们的supercell包含YZ1和YZ2两个界面，两个界面的ky是反向的，考虑这个因素，实际上YZ1界面只有其中的一支，和文献吻合的相当好了。

同样的道理，我们也可以构造一个supercell同时有XZ1和XZ2界面，这样一个模型就可以得到正文中的两个图了（原文Fig.2i&k）：

到了这里，故事才刚开始，走到这一步仍然不知道为什么这个表面态可以做负折射是不是？猫哥也看不出来啊，必须得算一个叫等频线（efc）的东西。参考正文中的图1，知道等频线的样子才能知道这个系统有没有负折射。在之前的计算中，我们输入k，得到频率。现在我们需要知道频率找k，猫哥不知道有没有直接的方法将这个系统变成输入频率求解K的本征值问题。但是，通过Comsol with Matlab可以实现！具体操作呢就是做一个k的二维循环，得到很多频率，用颜色来标定，找颜色一样的点连成线就可以。而这个文章中，Fig.2中的lmn三图中绿线的箭头就代表着频率由低到高的方向，这个就是重点。。画重点啊。这个箭头方向蕴含着这个系统有没有负折射：