从COMSOLMultiphysics?软件5.6版本开始，半导体模块的薛定谔方程物理场接口新增了处理多分量波函数的功能。本文我们将以均匀磁场中的硅量子点模型为例，继续探索这项新功能。
量子点简介量子点是纳米技术中必不可少的组成部分，在太阳能电池、发光二极管（LEDs）、显示器，光电探测器和量子计算中都具有潜在应用前景。Jock等人最近发表了一篇与自旋轨道量子位的应用领域相关的论文（参考文献1）。他们在该文的补充说明1中，提供了描述均匀磁场中硅量子点的公式，并在补充图1中显示了数值解。今天，我们将通过仿真的方法来重现该数值解。
硅量子点的薛定谔方程在参考文献1的补充说明中，方程1给出了均匀磁场中硅量子点的单电子哈密顿量，不包括自旋轨道耦合：
(1)滑动查看完整公式
其中，和是分别在横向和垂直方向上的有效质量；是量子点的约束势能；是玻尔磁子；是Pauli矩阵的向量；根据该论文所述，假定旋磁比张量是值为2的标量；动量由下式给出：
(2)式中，是基本电荷，是给定的磁矢势，并且虚数单元前面没有减号，因为COMSOLMultiphysics中的所有物理场接口都采用工程符号而不是。约束势能项由论文中的等式9给出：()滑动查看完整公式
其中，前两项表示横向各向异性谐波捕获势，其中和分别表示在x和y方向上的捕获势的角频率，第三项描述了一个在z方向上的电场势，是粒子电荷（电子）；最后一项表示在z=0的氧化硅界面处存在一个高度为的潜在势垒。
参数和变量
建模参数如下表所示：
除磁场使用了一个推测值1[T]外，所有的参数在论文中均给出。
磁矢势使用变量指定，假设均匀磁场指向y方向：
名称表达式单元描述Axz*B/2Wb/m矢量势Ay0[Wb/m]Wb/m矢量势Az-x*B/2Wb/m矢量势物理场设置薛定谔方程物理场接口设置了两个波函数，psiu和psid，分别对应自旋向上和自旋向下分量：
默认的有效质量域条件用来指定横向和垂直有效质量：
在此模型中，尽管两个波函数分量共享相同的有效质量，但在一般情况下用户界面仍允许不同波动函数分量具有不同的值。
为横向谐波陷波使用三个电子势能域条件设置约束势能项（方程）：
对于电场：
对于氧化物阻挡层：
向量势（方程2中的第二项）对动量的贡献是使用洛伦兹力域条件来模拟的：
最后，利用零阶哈密顿域条件实现自旋向上和自旋向下分量的磁耦合（方程1中的最后一项）：
请注意，由于零阶哈密顿域条件会自动乘以因子，所以在输入表达式时需要除去该因子，以便得到的公式和方程1的最后一项保持一致：(4)滑动查看完整公式
上面的屏幕截图为零阶哈密顿域条件设置窗口，哈密顿输入表的第一行实现了方程4矩阵的（1,2）元素，第二行实现了（2,1）元素。我们还创建了一个简单的扫掠网格，并使用了一个特征值研究来查找前几个特征状态。
理解结果通过绘制前两个本征态的自旋波函数分量的实部和虚部，我们发现它们具有相似的大小和形状，除了一些符号翻转：
与自旋向上分量相比，自旋向下波函数分量在它们之间也具有这种趋势。
为了理解这一观察，我们评估了峰值密度点附近的波函数分量：
特征值
向上（1），点：（0，0，-2）
向下（1），点：（0，0，-2）
0.
1.+1.i
1.-1.i
0.
1.-1.i
1.+1.i
可以看出，对于基态（上面结果表的第一行），直到总缩放比例约为.6e11时，自旋分量的幅度为，并且向下旋转分量的幅度为。因此，由两个分量形成的向量与成比例，它被认为是y自旋算子的自旋向下本征态。这与直观的图像一致，即磁场中电子的较低能态的自旋磁矩与磁场平行，因此自旋与磁场反平行。
类似地，对于第一个激发态，直到总缩放比例约为.6e11，自旋向上分量的幅值为，并且自旋向下分量的幅值为。因此，由两个分量形成的向量与成比例，它被认为是y自旋算子的自旋本征态。这与直观的图片一致，即磁场中电子的高能态的自旋磁矩反平行于磁场，因此自旋平行于磁场。
通过将计算得到的本征态之间的能量差与在均匀磁场下自旋向上和自旋向下的电子之间的预期能量差（）进行比较来进一步确认这一观察结果，利用全局计算得到：
特征值
计算能量差（meV）
预期能量差（meV）
0.
0.
0.
我们看到前两个本征态之间计算出的能量差和预期的能量差非常吻合——两者的估计值都为0.meV。
概率密度和动量密度下图显示了xz平面上基态的概率密度（蓝灰色渐变）和动量密度（红色箭头）。它与论文中补充图1较好地吻合。
下图显示了模型缩略图的第八个激发态的概率密度（等值面）和动量密度（箭头）。使用COMSOLMultiphysics5.6版本中的透明度子节点可以绘制透明的概率密度的等值面。
动手尝试本文我们讨论了多分量波函数的薛定谔方程物理场接口的功能。您也可以自己动手尝试使用这些功能进行建模。
点击底部“阅读原文”尝试在均匀的磁场中对硅量子点建模，您也可以下载MPH文件。
参考文献R.M.Jocketal.,“Asiliconmetal-oxide-semiconductorelectronspin-orbitqubit,”NatureCommunications,9:,.预览时标签不可点收录于话题#个上一篇下一篇