-
轨道选择Mott相变(orbital-selective Mott transition, OSMT)是过渡金属材料中的一种新奇物理现象: 多轨道体系中部分轨道上的电子进入局域化状态, 而其他轨道的电子仍然具有巡游特性. 实验上在Ca
$ _{2-x} $ SrxRuO4[1], 过渡金属硫化物[2]以及铁基超导材料[3]中均观测到了OSMT, 其对了解强关联系统的本质很有帮助. 双轨道Hubbard模型是研究OSMT最简单的模型[4-7]. 动力学平均场理论(dynamical mean-field theory, DMFT)在处理单轨道Hubbard 模型的相关问题中成效显著[8-11], 并且它拓展到双轨道Hubbard 模型中也十分有效[5,12,13].DMFT中处理Anderson杂质模型所用的杂质求解器对于最终的计算结果有很大的影响, 多种多样的理论以及数值方法都能用作动力学平均场的杂质求解器, 如量子蒙特卡罗方法(quantum Monte Carlo simulations, QMC)[14-17], 重整化群理论(renormalization-group theory)[18-20]以及隶算符表示方法 (slave-variable representations)[21-23]等. 当DMFT用于多轨道体系的求解时, 不同的杂质求解器会面临不同的困难: 在处理掺杂的费米体系时, QMC方法会出现负几率问题; 重整化群理论能很好地处理单带问题, 但却难以拓展到多带体系; 隶算符表示方法无法精确处理多轨道体系中的电子关联相互作用. 除了上述杂质求解器外, 还有以Lanczos方法为代表的精确对角化方法, 它能精确地处理多轨道体系中包括轨道内相互作用U、轨道间相互作用U' 以及洪特耦合
$ J_{\rm{H}} $ 在内的电子关联作用, 但是对于含有非局域轨道间跃迁的系统, 轨道间跃迁所引入的非对角格林函数将会使自洽参数求解过程的复杂程度显著增加, 并需要耗费巨大的计算资源.轨道间跃迁会诱发轨道杂化, 其在具有多轨道关联效应的过渡金属中至关重要[24-26], 但同时它所带来的非对角格林函数为处理多轨道Hubbard模型带来了巨大挑战. 本文引入正则变换将双轨道Hubbard 模型中的紧束缚部分对角化[27,28], 得到了没有轨道间跃迁的有效杂化轨道. 同时为了更精确地计算轨道间相互作用对体系的影响, 选取Lanczos精确对角化方法作为杂质求解器, 比起早先的工作[27,28], 该方法能更加精确地描述库仑相互作用以及洪特耦合效应对体系量子相变的影响, 因此所得相图的相变点也更加准确. 除了DMFT理论中常用的Bethe晶格外, 本文还基于正方晶格能带结构计算了双轨道体系的相图, 并发现DMFT中选用的晶格能带结构不同时, 虽然所得相变点存在差异, 但相图展现出来的轨道选择Mott相变的物理图像基本一致. 最后应用文中的方法研究了最近发现的新型铜氧化物超导材料Ba2CuO
$ _{4-\delta} $ [29]. Ba2CuO$ _{4-\delta} $ 是高空穴掺杂浓度下的铜氧化物超导材料, 其超导转变温度约为73 K. 由于其具有极其扁平的CuO6八面体结构, 铜原子的$ {\rm{d}}_{x^2-y^2} $ 轨道和$ {\rm{d}}_{3 z^2-r^2} $ 轨道均对费米面附近的电子态有贡献[29-32], 这种特殊的双轨道特性是其区别于传统铜氧化物超导材料的地方, 同时密度泛函理论(density functional theory, DFT) 的计算结果也对双带特性进行了印证[33-35]. 有理论研究表明这种特殊的多轨道关联特征极有可能是Ba2CuO$ _{4-\delta} $ 的超导转变温度高于大部分高温超导材料的原因[34,36]. 对于半满的Ba2CuO$ _{4-\delta} $ 体系, 我们考虑非局域轨道间跃迁的动量空间各向异性, 通过Ba2CuO$ _{4-\delta} $ 的完整能带, 得到了材料在半满条件下的相图.本文主要结构如下: 第2部分对含轨道间跃迁的双带Hubbard模型中的正则变换进行介绍; 第3部分对研究所用的理论方法—以Lanczos为求解器的DMFT进行详细的说明; 第4部分介绍基于上述方法得到的双轨道Hubbard模型中轨道选择Mott相变的物理图像以及轨道间跃迁对其的影响, 并对比用正方晶格能带代入DMFT中计算与Bethe晶格的差异; 第5部分将方法应用于双轨道超导材料Ba2CuO3.5中, 并介绍了所得的结果; 第6部分为文章的主要结论.
-
双轨道体系的哈密顿量可以分为两部分: 紧束缚部分
$ H_{\rm{t}} $ 和相互作用部分$ H_{\rm{I}} $ , 紧束缚部分$ H_{\rm{t}} $ 的形式为$ \begin{split} H_{\rm{t}} =\;& -\sum\limits_{\langle ij \rangle}\sum\limits_{l \sigma} t_ld^{\dagger}_{il\sigma}d_{jl\sigma} -\sum\limits_{\langle ij \rangle}\sum\limits_{l\neq l', \sigma} t_{ll'}d^{\dagger}_{il\sigma}d_{jl'\sigma}\\ &-\mu\sum\limits_{il\sigma} d^{\dagger}_{il\sigma}d_{il\sigma}, \\[-18pt]\end{split} $ $ \begin{split} H_{\rm{I}} =\;& \frac{U}{2}\sum\limits_{il\sigma}n_{il\sigma}n_{il\bar{\sigma}} + \sum\limits_{i,l < l',\sigma\sigma'} (U' - \delta_{\sigma\sigma'}J_{\rm{H}})n_{il\sigma}n_{il'\sigma'} \\ &+\frac{J_{\rm{H}}}{2}\sum\limits_{i,l\neq l',\sigma} d^{\dagger}_{il\sigma}d^{\dagger}_{il\bar{\sigma}}d_{il'\bar{\sigma}}d_{il'\sigma}\\ &+\frac{J_{\rm{H}}}{2}\sum\limits_{i,l\neq l',\sigma\sigma'} d^{\dagger}_{il\sigma}d^{\dagger}_{il'\sigma'}d_{il\sigma'}d_{il'\sigma}, \\[-18pt]\end{split} $ 其中
$ d^{\dagger}_{il\sigma} $ ($ d_{il\sigma} $ )是i 格点的l轨道上自旋为σ电子的产生(湮灭)算符,$ \langle ij \rangle $ 表示遍布全空间最近邻格点的求和.$ t_l $ 为最近邻同轨道的跃迁系数,$ t_{ll'} $ 为最近邻不同轨道间的跃迁系数. U 和U' 分别表示轨道内电子相互作用和轨道间电子相互作用,$ J_{\rm{H}} $ 为洪特耦合强度. 本文研究的系统具有自旋旋转对称性, 相互作用部分$ H_{\rm{I}} $ 的系数满足$ U= U'+2 J_{\rm{H }}$ 的关系.通过正则变换消除轨道间跃迁, 进而引入两个退耦合的有效轨道α和β [27,28], 正则变换的形式如下:
$ \begin{split} &d_{i1\sigma } = u\alpha _{i\sigma } + v\beta _{i\sigma }, \\ &d_{i2\sigma } = -v\alpha _{i\sigma } + u\beta _{i\sigma }, \end{split} $ 其中
$ \begin{split} &{u} = -\frac{t_{12}}{\sqrt{2}\left|t_{12}\right|}\left[ {1 + \sqrt {\frac{{{{\left[ {{t_{1}} - {t_{2}}} \right)}^2}}}{{{(t_{12})^2} + {{\left( {{t_{1}} - {t_{2}}} \right)}^2}}}} }\; \right]^{1/2},\\ &{v} = \frac{1}{\sqrt 2}\left[ {1 - \sqrt {\frac{{{{\left( {{t_{1}} - {t_{2}}} \right)}^2}}}{{{(t_{12})^2} + {{\left( {{t_{1}} - {t_{2}}} \right)}^2}}}} } \;\right]^{1/2}.\ \end{split} $ (3)式中,
$ \alpha_{i\sigma} $ 和$ \beta_{i\sigma} $ 分别表示i格点上的有效轨道α和β中自旋为σ的电子湮灭算符, 由(4)式决定的u值和v值代入哈密顿量(1)中将会使轨道间跃迁项消失. 通过正则变换, (1)式和(2)式描述的原双轨道的哈密顿量会转换为有效轨道的双轨道哈密顿量$ H^{{\rm{eff}}} $ , 其同样包含紧束缚部分$ H_{\rm{t}}^{{\rm{eff}}} $ 和相互作用部分$ H_{\rm{I}}^{{\rm{eff}}} $ , 紧束缚部分$ H_{\rm{t}}^{{\rm{eff}}} $ 为$ \begin{split} H_{\rm{t}}^{{\rm{eff}}} = \;& - \sum\limits_{\langle {i,j} \rangle \sigma } {\left( {{t_\alpha }\alpha _{i\sigma }^+ {\alpha _{j\sigma }} + {t_\beta }\beta _{i\sigma }^ + {\beta _{j\sigma }}} \right)}\\ &-\mu \sum\limits_{i\sigma}{\left({\alpha _{i\sigma }^ + {\alpha _{i\sigma }} + \beta _{i\sigma }^ + {\beta _{i\sigma }}} \right)}, \end{split} $ 相互作用部分
$H_{\rm{I}}^{{\rm{eff}}}$ 为$ \begin{split}\;& H_{\rm{I}}^{{\rm{eff}}} = \frac{U}{2} \sum\limits_{i \sigma}\left(n_{i \alpha \sigma} n_{i \alpha \bar{\sigma}}+n_{i \beta \sigma} n_{i \beta \bar{\sigma}}\right)\\ &~~+\sum\limits_{i \sigma \sigma^{\prime}}\left(U^{\prime}-\delta_{\sigma \sigma^{\prime}} J_{\rm{H}}\right) n_{i \alpha \sigma} n_{i \beta \sigma^{\prime}}\\ &~~+\frac{J_{{\rm{H}}}}{2} \sum\limits_{i, \sigma}\left(\alpha_{i \sigma}^{\dagger} \alpha_{i \bar{\sigma}}^{\dagger} \beta_{i \bar{\sigma}} \beta_{i \sigma}+\beta_{i \sigma}^{\dagger} \beta_{i \bar{\sigma}}^{\dagger} \alpha_{i \bar{\sigma}} \alpha_{i \sigma}\right)\\ &~~+\frac{J_{{\rm{H}}}}{2} \sum\limits_{i, \sigma \sigma^{\prime}} \left(\alpha_{i \sigma}^{\dagger} \beta_{i \sigma^{\prime}}^{\dagger} \alpha_{i \sigma^{\prime}}\beta_{i \sigma} + \beta_{i \sigma}^{\dagger} \alpha_{i \sigma^{\prime}}^{\dagger} \beta_{i \sigma^{\prime}} \alpha_{i \sigma}\right) . \end{split} $ 根据正则变换(3)可知紧束缚部分
$H_{\rm{I}}^{{\rm{eff}}}$ 的系数为$ \begin{split} &t_{\alpha}=t_1 u^{2}+t_2 v^{2}-t_{12}uv,\\ &t_{\beta}=t_1 v^{2}+t_2 u^{2}+t_{12}uv, \end{split} $ 结合(4)式可以看出, 正则变换后的跃迁积分
$ t_{\alpha} $ 和$ t_{\beta} $ 与$ t_{12} $ 的取值正负无关. 因为本文研究的系统具有自旋旋转对称性$ U=U'+2 J_{\rm{H}} $ , 相互作用部分$ H_{\rm{I}}^{{\rm{eff}}} $ 的系数在正则变换后保持不变. -
正则变换之后去除轨道杂化的双轨道模型, 可以通过DMFT有效地处理.
在DMFT的基本框架中[8], 我们将多体问题中格点数量巨大的格点哈密顿量(5)和(6) 映射成格点数目有限的Anderson杂质哈密顿量
$H_{{\rm{imp}}}$ ,$ \begin{split} H_{{\rm{imp}}} =\;& \sum\limits_{m\sigma}( \varepsilon _{m\sigma}^{\alpha} c^{\dagger}_{\alpha m\sigma} c_{\alpha m\sigma} + \varepsilon _{m\sigma}^{\beta} c^{\dagger}_{\beta m\sigma}c_{\beta m\sigma}) \\ &+\sum\limits_{m\sigma}V_{m\sigma}^{\alpha}(c^{\dagger}_{\alpha m\sigma}\alpha_{\sigma} +\alpha^{\dagger}_{\sigma}c_{\alpha m\sigma}) \\ &+\sum\limits_{m\sigma}V_{m\sigma}^{\beta}(c^{\dagger}_{\beta m\sigma}\beta_{\sigma} +\beta^{\dagger}_{\sigma}c_{\beta m\sigma}) -\mu\sum\limits_{\sigma}\alpha^{\dagger}_{\sigma}\alpha_{\sigma} \\ &-\mu\sum\limits_{\sigma}\beta^{\dagger}_{\sigma}\beta_{\sigma} +H^{{\rm{eff}}}_{{\rm{I}}}(\alpha, \beta), \\[-18pt]\end{split} $ 其中
$ c^{\dagger}_{\gamma m\sigma} $ ($ c_{\gamma m\sigma} $ ) 表示外场格点(浴格点) m上γ轨道($ \gamma=\alpha, \beta $ )中自旋为σ 的电子产生(湮灭)算符,$ \varepsilon _{m\sigma}^{\gamma} $ 为第m个外场格点上轨道γ的能量,$ V_{m\sigma}^{\gamma} $ 表示外场格点m的轨道γ和本文研究的杂质格点轨道γ之间的耦合系数. 本文计算选取的外场格点数目$ n_{\rm{b}}=3 $ , 因为DMFT的相关研究证明当外场格点数$n_{\rm{b}}\geqslant3$ 时, 以Lanczos作为杂质求解器的DMFT得到的OSMT相变点基本一致[37].正则变换之后的哈密顿量去除了轨道间的杂化, 因此格林函数和自能都不含关于轨道的非对角项, 各轨道的格林函数可以独立地计算. 杂质模型的Weiss场格林函数(无相互作用格林函数)可以由杂质模型的系数直接计算:
$ {{{\rm{G}}}}_{0\gamma \sigma}^{-1}\left({\rm{i}}\omega_n\right) ={\rm{ i}}\omega_n +\mu -\varepsilon _{\gamma}-\sum\limits_{m}\frac{\left(V_{m\sigma}^{\gamma}\right)^2} {{\rm{i}}\omega_n- \varepsilon _{m\sigma}^{\gamma}}. $ 通过Lanczos杂质求解器可以得到杂质模型格林函数
$G_{{\rm{i m p}}}^{(\gamma)}$ [38-40], 在Lanczos方法中$ \begin{array}{l} {\rm{G}}_{{\rm{i m p}}}^{(\gamma)}\left( {\rm{i}}\omega_n\right)={\rm{G}}_{\gamma}^{(+)}\left( {\rm{i}}\omega_n\right)+{\rm{G}}_{\gamma}^{(-)}\left( {\rm{i}}\omega_n\right), \end{array} $ 其中
$ {\rm{G}}_{\gamma}^{(+)}\left( {\rm{i}}\omega_n\right)=\dfrac{\langle\phi_{0}\left|\gamma \gamma^{\dagger}\right| \phi_{0}\rangle}{ {\rm{i}}\omega_n-a_{0}^{(+)}-\dfrac{b_{1}^{(+) 2}}{ {\rm{i}}\omega_n-a_{1}^{(+)}-\dfrac{b_{2}^{(+) 2}}{ {\rm{i}}\omega_n-a_{2}^{(+)}-\cdots}}}, $ $ {\rm{G}}_{\gamma}^{(-)}\left( {\rm{i}}\omega_n\right)=\dfrac{\langle\phi_{0}\left|\gamma^{\dagger} \gamma\right| \phi_{0}\rangle}{ {\rm{i}}\omega_n+a_{0}^{(-)}-\dfrac{b_{1}^{(-) 2}}{ {\rm{i}}\omega_n+a_{1}^{(-)}-\dfrac{b_{2}^{(-) 2}}{ {\rm{i}}\omega_n+a_{2}^{(-)}-\cdots}}}. $ DMFT通过自洽方程来使Weiss场格林函数和无相互作用的格点格林函数相等, 进而可以用杂质模型来表征格点模型的物理性质. 如果考虑具有半圆形态密度的Bethe晶格, DMFT中的无相互作用的格点格林函数可以化简为非常简洁且便于计算的形式[8]:
$ \begin{array}{l} {\rm{g}}_{0}^{(\gamma)}\left({\rm{i}}\omega_n\right)^{-1}={\rm{i}}\omega_n+\mu-t_{\gamma}^{2} {\rm{G}}_{{\rm{i m p}}}^{(\gamma)}\left( {\rm{i}}\omega_n\right), \end{array} $ 其中
${\rm{g}}_{0}( {\rm{i}}\omega_n)$ 为无相互作用的格点模型格林函数. (13)式中${\rm{g}}_{0}( {\rm{i}}\omega_n)$ 以及由(9)式计算的Weiss场格林函数${\cal {{\rm{G}}}}_{0}( {\rm{i}}\omega_n)$ 做差即为DMFT的自洽方程:$ \begin{equation} d=\frac{1}{n_{\max }+1} \sum\limits_{\gamma=1}^{2}\sum\limits_{n=0}^{n_{\max }}\left|{\rm{g}}_{0}^{(\gamma)}\left({\rm{i}} \omega_{n}\right)^{-1}-{\cal {{\rm{G}}}}_{0\gamma}^{-1}\left({\rm{i}} \omega_{n}\right)\right|^{2}, \end{equation} $ 其中
$ n_{{\rm{max}}} $ 为截断系数, 计算中取nmax = 2048 . 利用共轭梯度算法调整自洽参数$ V_{m}^{\gamma} $ 和$ \varepsilon _{m}^{\gamma} $ 使(14)式中的d 趋于0, 我们工作中体系自洽完成的判断标准为d值小于$ 10^{-5} $ .如果计算中选用其他晶格的能带结构, 则在计算局域的格点格林函数
${\rm{G}}_{{\rm{lat}}}\left({\rm{i}} \omega_{n}\right)$ 的过程中需要在动量空间对动量$ {\boldsymbol{k}} $ 求和, 以正方晶格为例,$ \begin{equation} {\rm{G}}_{{\rm{lat}}}\left({\rm{i}} \omega_{n}\right)=\frac{1}{N} \sum\limits_{{\boldsymbol{k}}} \frac{1}{{\rm{i}} \omega_{n}+\mu-\varepsilon_{{\boldsymbol{k}}}-\varSigma\left({\rm{i}} \omega_{n},{\boldsymbol{k}}\right)}, \end{equation} $ 其中
$ \varepsilon_{{\boldsymbol{k}}}=2 t(\cos k_x+\cos k_y) $ 代表正方晶格的色散关系,$\varSigma\left({\rm{i}} \omega_{n}, {\boldsymbol{k}}\right)$ 为体系自能(因为各轨道独立, 所以这里略去了轨道指标γ). 在DMFT中, 当体系的维度d趋于无穷大时空间的涨落将会被冻结, 此时系统自能将不再依赖动量$ {\boldsymbol{k}} $ 而局域化[8]:$ \begin{equation} \lim\limits _{d \rightarrow \infty}\varSigma\left({\rm{i}} \omega_{n},{\boldsymbol{k}}\right)=\varSigma_{{\rm{loc}}}\left({\rm{i}} \omega_{n}\right). \end{equation} $ 通过Dyson方程可以计算杂质模型的自能:
$ \begin{equation} \varSigma_{{\rm{loc}}}\left({\rm{i}} \omega_{n}\right)={\cal {{\rm{G}}}}_{0}^{{\rm{imp}}}\left({\rm{i}} \omega_{n}\right)^{-1}-{\rm{G}}_{{\rm{imp}}}\left({\rm{i}} \omega_{n}\right)^{-1}, \end{equation} $ 其中,
${\cal {{\rm{G}}}}^{{\rm{imp}}}_{0}{\rm{(i}} \omega_{n})$ 和${\rm{G}}_{{\rm{imp}}}({\rm{i}} \omega_{n})$ 分别为利用(9)式和 (10)式计算的杂质模型的Weiss场格林函数以及考虑相互作用的格林函数. DMFT中认为杂质格林函数的自能与格点格林函数的局域自能相等, 因此可以通过(15)式—(17)式计算格点格林函数${\rm{G}}_{{\rm{lat}}}({\rm{i }}\omega_{n})$ . 对于正方晶格, 本文选取的自洽标准为${\rm{G}}_{{\rm{imp}}}\left({\rm{i}} \omega_{n}\right)={\rm{G}}_{{\rm{lat}}}\left({\rm{i}} \omega_{n}\right)$ , 同样地, 在自洽过程通过共轭梯度算法调整自洽参数V和ε来使${\rm{G}}_{{\rm{imp}}}\left({\rm{i}} \omega_{n}\right)$ 和$ G_{{\rm{lat}}}\left({\rm{i}} \omega_{n}\right) $ 相等.利用体系的格林函数虚部来计算各轨道的态密度:
$ \begin{equation} A_{\gamma}(\omega)=-\frac{1}{\pi}\rm{Im}\it{{\rm{G}}^{\gamma}_{ii}(\omega+{\rm{i}}\eta)}, \end{equation} $ 其中
$ \eta $ 为能量展宽因子, γ为轨道指标. 各轨道的准粒子权重可以通过自能计算[41]:$ \begin{equation} Z_{\gamma}=\left[1-\frac{\partial}{\partial\omega}{\rm{Re}} \varSigma_{\gamma}(\omega+{\rm{i}}\eta)|_{\omega=0}\right]^{-1}. \end{equation} $ -
首先研究相互作用U和轨道间跃迁
$ t_{12} $ 对Bethe晶格双带Hubbard模型相变的影响. 研究中选取轨道间跃迁$ t_{12} > 0 $ , 因为正则变换后的跃迁积分$ t_{\alpha} $ 和$ t_{\beta} $ 与$ t_{12} $ 的取值正负无关, 所以之后的研究结果能够自然地拓展到$ t_{12} < 0 $ 的情况中. 在本文研究的半满体系中, 化学势保持$ \mu=U/2+ U^{\prime}- J_{\rm{H}}/2 $ 以满足粒子空穴对称性. 图1比较了当两轨道带宽比$ t_2/t_1=0.4 $ , 相互作用U以及轨道间跃迁$ t_{12} $ 不同时体系准粒子权重的变化. 当$ t_{12}=0 $ 时, 从准粒子权重曲线可以看出, 宽有效轨道α金属绝缘体相变 (metal-insulator transition, MIT)的相变点$ U_{c\alpha}=3.7 $ , 而窄有效轨道β在相互作用强度$ U_{c\beta}=2.3 $ 时发生金属绝缘体相变, 因此体系在相互作用强度$ 2.3\leqslant U < 3.7 $ 时处于轨道选择Mott相(orbital-selective Mott phase, OSMP). 当轨道间跃迁系数$ t_{12}=0.4 $ 时, 轨道选择Mott相对应的相互作用区间扩大为$ 2.1\leqslant U < 4.1 $ , 轨道间跃迁的引入使窄有效轨道β 发生MIT所需的相互作用强度减弱, 而使宽有效轨道α的临界相互作用$ U_{c\alpha} $ 增强. 由(7)式可知, 轨道间相互作用通过调整两个有效轨道的跃迁系数来促进OSMT的发生, 如当$ t_{12}=0.4 $ 时, 两个有效轨道的跃迁系数比值为$ t_{\beta}/t_{\alpha}=0.32 $ , 因此体系更容易发生轨道选择Mott相变.图 1 不同轨道间跃迁系数下两个有效轨道的准粒子权重
$Z$ 随相互作用U 的变化曲线. 轨道α 和轨道β分别用黑色圆圈和红色三角标识, 实心标志和空心标志分别代表$t_{12}=0$ 和$t_{12}=0.4$ 的情况. 体系中的洪特耦合作用$ J_{\rm{H}}=0.25 U$ , 两轨道带宽比$t_2/t_1=0.4$ , 所选用的能量单位为宽带跃迁系数$t_1$ Figure 1. The interaction dependencies of the quasiparticle weight Z of the effective orbital α (black circles) and β (red triangles) with interorbital hopping
$t_{12}=0.4$ (hollow symbols) and without interorbital hopping$t_{12}=0$ (solid symbols) when$J_{\rm{H}}=0.25 U$ and$t_2/t_1=0.4$ . The energy is in unit$t_1$ 图2给出了当轨道间跃迁
$ t_{12}=0.4 $ 、带宽比$ t_2/t_1=0.4 $ 以及洪特耦合$ J_{\rm{H}}=0.25 U $ 时双轨道体系的态密度随相互作用强度U变化的细节. 当相互作用强度$ U=1.5 $ 时, 有效轨道α (图2(a))和β(图2(b)) 在费米面附近都存在明显的准粒子峰, 这意味着两个轨道都处于金属态, 因此整个体系表现出金属相的特征. 当相互作用强度增加至$ U=3.0 $ 时, 轨道α在费米面附近仍存在明显的准粒子峰(图2(c)), 但此时轨道β 的费米面附近打开了Mott能隙(图2(d)), 这是OSMP的典型特征. 继续增加相互作用强度至$ U=4.5 $ , 图2(e)和图2(f)中两个有效轨道均在费米面附近打开了Mott能隙, 这也意味着体系进入了Mott绝缘相. 图2给出的物理图像与图1中$ t_{12}=0.4 $ 的曲线相符合, 这二者也是本文判定体系物态以及绘制相图的依据.图 2 两个有效轨道的态密度随相互作用U 增加的演化图像, 体系中的轨道间跃迁
$t_{12}=0.4$ , 带宽比$t_2/t_1=0.4$ 以及相互作用强度$J_{\rm{H}}=0.25 U.$ 左边的图像为α轨道的结果, 右边图像为β轨道的结果, 计算中选取的能量展宽因子为$\eta=0.05$ Figure 2. Evolution of the orbital-resolved spectral density
$A(\omega)$ with the increasing Coulomb interaction U when$t_{12}=0.4$ ,$t_2/t_1=0.4$ and$J_{\rm{H}}=0.25 U.$ Left panels show the results of effective orbital α and right panels are for effective orbital β. The energy broadening factor in our calculation takes$\eta=0.05$ 图3(a)所示为轨道间跃迁系数
$ t_{12} $ 不同时, 两轨道带宽比$ t_2/t_1 $ 和相互作用强度U 影响下的体系相图. 在没有轨道间跃迁的双轨道体系中, 除了$ t_2/t_1=1 $ 的情况外体系中均存在OSMP, 这主要是由于洪特耦合的作用[42], 相图中OSMP区域的范围随着带宽比$ t_2/t_1 $ 的增加而收缩, 说明两轨道差异越大体系越容易发生轨道选择Mott相变, 这与早先对OSMP的研究结果相符合[43-45]. 在双轨道体系中引入轨道间跃迁之后, 两个有效轨道之间的差异会增大, 这将导致相图中轨道选择Mott相的区域扩大, 图3(b) 直观地反映了当$ t_{12}=0.4 $ 时相图中OSMP区域的扩大, 和图3(a)对比可以看出, 所有的$ t_2/t_1 $ 情况下轨道选择Mott相的区域都扩大了, 并且两轨道差异越小轨道间跃迁对OSMP的促进作用越明显.图 3 轨道间跃迁系数不同时, 相互作用U和带宽比
$t_2/t_1$ 影响下的双轨道体系相图 (a) 轨道间跃迁$t_{12}=0$ ; (b) 轨道间跃迁$t_{12}=0.4$ . 体系的洪特耦合$J_{\rm{H}}=0.25 U$ . 黑色圆圈和红色三角分别代表轨道α和β 的MIT 相变点Figure 3. Phase diagrams in the plane of interaction U and hopping integral
$t_2/t_1$ of the effective two-orbital Hubbard model with different interorbital hopping, when$J_{\rm{H}}=0.25 U$ : (a)$t_{12}=0$ ; (b)$t_{12}=0.4$ . The black circles (red triangles) denote the critical points of MIT for effective orbital α(β)根据(15)式及正方晶格的能带结构, 研究了正方晶格双轨道Hubbard模型中的OSMT. 图4给出了正方晶格
$ t_2/t_1=0.4 $ 时轨道间跃迁$ t_{12} $ 及相互作用U影响的准粒子权重变化, 以及$ t_{12}=0.4 $ 时相互作用U和带宽比$ t_2/t_1 $ 决定的相图, 可与Bethe晶格的结果(图1和图3(b))做对比. 如图4(a)所示, 当$ t_{12}=0 $ 时体系在相互作用强度$ 3.9\leqslant U < 6.6 $ 时处于轨道选择Mott相; 引入轨道间跃迁$ t_{12}= 0.4 $ , OSMP存在范围扩大至$ 3.5\leqslant U < 6.9 $ , 尽管相变点的具体数值与Bethe晶格的结果存在差异, 但是$ t_{12}=0.4 $ 对OSMP影响的物理图像是一致的. 图4(b)进一步印证了这一点,$ t_{12}=0.4 $ 时体系相图中的OSMT的相变点与图3(b)中的存在差异, 但二者的OSMP的变化趋势基本一致, 都反映出两轨道带宽差异越小轨道间相互作用对OSMP的促进作用越明显. 综上, Bethe晶格是一种理想的晶格, 其态密度作为一种近似手段来定性研究实际晶格中关联效应驱动的相变是可行的, 这也正是Bethe晶格被广泛应用于DMFT的原因.图 4 DMFT结合正方晶格能带结构所得到的结果 (a)轨道间跃迁系数不同时, 准粒子权重Z 随相互作用U 的变化曲线, 体系的洪特耦合
$J_{\rm{H}}$ = 0.25U, 带宽比$t_2$ /$t_1$ = 0.4; (b)$t_{12}=0.4$ 以及$J_{\rm{H}}=0.25 U$ 时, 相互作用U 和带宽比$t_2/t_1$ 影响下的有效轨道模型相图Figure 4. DMFT results of the square lattice: (a) The interaction dependencies of the quasiparticle weight Z with different interorbital hoppings when
$J_{\rm{H}}$ = 0.25U and$t_2$ /$t_1$ = 0.4; (b) phase diagrams in the plane of interaction U and hopping integral$t_2/t_1$ of the effective two-orbital Hubbard model when$t_{12}=0.4$ and$J_{\rm{H}}=0.25 U$ -
把结合了正则变换的扩展DMFT运用到最近发现的铜氧化物超导材料Ba2CuO
$ _{4-\delta} $ 中[29], 由于Ba2CuO$ _{4-\delta} $ 晶格中具有极其扁平的CuO6八面体结构, 所以可以用双轨道Hubbard模型描述其电子结构与性质. 根据DFT计算半满情况下Ba2CuO3.5所得到的电子能带结构可知, 费米面附近的电子态主要由Cu的$ {\rm{d}}_{x^2-y^2} $ 轨道和$ {\rm{d}}_{3 z^2-r^2} $ 轨道贡献[33-35]. 将DFT所得的能带结果近似处理, 只考虑最近邻跃迁而忽略较小的次近邻跃迁, 并忽略非局域轨道间跃迁的动量空间各向异性. 之后采用适用于动量空间的各向异性的扩展DMFT研究了Ba2CuO$ _{4-\delta} $ 半满情况下的相图.选取如下的跃迁模型系数代入(1)式:
$ t_{1} = 0.504 {\rm{eV}}$ ,$ t_{2} =0.196 {\rm{eV}}$ , 及$ t_{12}=-0.302 {\rm{eV}} $ ,$ {\rm{d}}_{x^2-y^2} $ 轨道和$ {\rm{d}}_{3 z^2-r^2} $ 轨道的能量分别为$ \mu_1= -0.222 {\rm{eV}}$ 和$ \mu_2=0.661 {\rm{eV}} $ [33]. 第2部分中采用正则变换后所得的跃迁积分$ t_{\alpha} $ 和$ t_{\beta} $ 与$ t_{12} $ 的取值正负无关, 因此可以将第4部分的结果拓展至$ t_{12} < 0 $ 的情况. 结合本文的Ba2CuO$ _{4-\delta} $ 紧束缚模型, 体系中轨道间跃迁$ t_{12} $ 将促进体系的有效轨道的OSMP, 同时由于Ba2CuO$ _{4-\delta} $ 材料中具有晶体场劈裂$\varepsilon _{\rm{d}}= \mu_1-\mu_2$ , 不同轨道化学势的差异会导致有效轨道的格林函数矩阵形式非对角元中出现常数项, 其中常数项的处理采用与之前工作中同样的方法[28].图5(a)给出了
$J_{\rm{H}}=0.25 U$ 时各轨道的准粒子权重随相互作用U的变化曲线. 根据准粒子权重的变化曲线可知, 当相互作用$ U < 2.1 $ eV时, 双轨道体系为金属态; 当相互作用强度$ U\geqslant 2.3 $ eV时, 体系转变为绝缘态; 在这两个态之间存在一个狭窄的OSMP区域, 即作用强度$ 2.1\; {\rm{eV}}$ $ \leqslant U < 2.3\;{\rm{eV}} $ 时, 宽有效轨道α处于金属态而窄有效轨道β呈现出绝缘态的特征. 图5(b)—(d)进一步展现了$ J_{\rm{H}}= 0.25 U $ 时相互作用强度U对两个有效轨道α和β电子态密度的影响. 当相互作用强度U较弱, 如图5(b)所示,$ U=1.0 $ eV时两个轨道在费米能级附近的电子态密度都不为零, 说明体系为金属态; 随着相互作用强度U增加到2.5 eV, 如图5(d)所示, 两个轨道的态密度图像中均有Mott能隙打开, 此时体系为绝缘态; 而当相互作用强度$ U=2.2 $ eV时, 如图5(c)所示, 宽轨道费米能级附近的态密度仍为有限值的同时窄轨道打开了Mott能隙, 说明此时体系处于轨道选择Mott相[44-46].图 5 (a)
$J_{\rm{H}}=0.25 U$ 时有效轨道的准粒子权重随相互作用强度的变化曲线. 不同相互作用强度U下体系各有效轨道的态密度图像 (b)$U=1.0$ eV; (c)$U=2.2$ eV; (d)$U=2.5$ eV. 当相互作用2.1 eV$\leqslant U < $ 2.3 eV时体系处于OSMPFigure 5. (a) Quasiparticle weight Z as a function of interaction U when
$J_{\rm{H}}=0.25 U$ . The orbital-resolved spectral density$A(\omega)$ with different intraorbital interaction: (b)$U=1.0$ eV; (c)$U=2.2$ eV; (d)$U=2.5$ eV for the two effective orbitals. An OSMP occurs in a narrow interaction region with 2.1 eV$\leqslant U < $ 2.3 eV图6给出了在相互作用U 和洪特耦合
$ J_{\rm{H}} $ 作用下的相图, 在强电子关联作用导致的Mott绝缘相和相互作用强度较弱时体系的金属相之间存在着一定范围的OSMP区域, 随着洪特耦合$ J_{\rm{H}} $ 的减小以及相互作用U的增加OSMP区域逐渐缩小, 并在$ J_{\rm{H}}=0.34 $ eV以及$ U=2.7 $ eV时在相图上消失. 根据之前对电子关联效应的研究, 洪特耦合作用有利于OSMP的形成[12,42,47], 这也解释了相图中当$ J_{\rm{H}} > 0.34 $ eV时OSMP区域随洪特耦合增加而扩大的原因. 而当洪特耦合$ J_{\rm{H}} < 0.34 $ eV 时, 晶体场劈裂对体系的影响压制住了库仑相互作用以及洪特耦合的影响, 因此OSMP从相图上消失[12,42,47].图 6 相互作用U和洪特耦合
$J_{\rm{H}}$ 影响下体系的相图. 随着相互作用U的增加以及洪特耦合$J_{\rm{H}}$ 的减小, OSMP区域逐渐缩小, 大约在$J_{\rm{H}}=0.34$ eV以及$U=2.7$ eV时体系中的OSMT消失Figure 6. The phase diagram of the effective two-orbital Hubbard model with interaction U and Hund's rule coupling
$J_{\rm{H}}$ . The region of OSMP becomes narrower with the decreasing$J_{\rm{H}}$ and increasing U, and OSMT vanishes around$J_{\rm{H}}=0.34$ eV and$U=2.7$ eV图7(a)给出了本文研究的紧束缚模型的能带图像, 并将其与根据DFT模拟Ba2CuO
$ _{4-\delta} $ 的双带模型能带图像(图7(b))进行了对比.图 7 由不同模型得到的半满Ba2CuO
$_{4-\delta}$ 体系的能带图 (a)$t_{1}$ =0.504 eV,$t_{2}$ =0.196 eV,$t_{12}=-0.302$ eV,$\mu_1=-0.222$ eV 和$\mu_2=0.661$ eV的紧束缚模型所得到的能带图; (b) 根据DFT模拟的双带模型能带图像[33]Figure 7. Bandstructure for half-filling Ba2CuO3.5 with different model: (a) Tight binding model with
$t_{1}$ =0.504 eV,$t_{2}$ =0.196 eV,$t_{12}=-0.302$ eV,$\mu_1=-0.222$ eV and$\mu_2=0.661$ eV; (b) two-orbital model according to DFT results[33].模拟DFT计算结果所得的半满情况下Ba2CuO
$ _{4-\delta} $ 的能带结构如下[33]:$H_{0}(k)=\sum\limits_{k \sigma} \sum\limits_{l l^{\prime}}\left\{\xi_{l l^{\prime}}(k)+\left( \varepsilon _{l}-\mu\right) \delta_{l l^{\prime}}\right\} d_{l \sigma}^{\dagger}(k) d_{l^{\prime} \sigma}(k), $ 其中
$ \begin{split} \xi_{11 / 22}(k) = \;&2 t_{11 / 22}\left(\cos k_{x}+\cos k_{y}\right)\\ &+4 t_{11 / 22}^{\prime} \cos k_{x} \cos k_{y}\\ &+2 t_{11 / 22}^{\prime \prime}\left[\cos (2 k_{x})+\cos (2 k_{y})\right], \end{split} $ 以及
$ \begin{split} &\xi_{12}(k) = \xi_{21}(k)\\ = \;&2 t_{12}\left(\cos k_{x}-\cos k_{y}\right) \\ & +2 t_{12}^{\prime \prime}\left[\cos (2 k_{x})-\cos (2 k_{y})\right], \end{split} $ 式中
$ t_{ab}(a, b=1, 2) $ 表示由轨道$ a $ 跃迁至$ b $ 的跃迁系数, 下标1和2分别表示$ {\rm{d}}_{x^2-y^2} $ 轨道和$ {\rm{d}}_{3 z^2-r^2} $ 轨道,$ t $ ,$ t^{\prime} $ 和$ t^{\prime \prime} $ 分别表示最近邻跃迁、次近邻跃迁以及第三近邻跃迁系数, 具体的数值见表1. 由于本文选取的紧束缚模型忽略了非局域轨道间跃迁的动量空间各向异性, 因此其能带结构和根据DFT模拟的双带模型能带结构存在明显差异, 各向同性的紧束缚模型中两个能带存在分离的特征, 因此需要采用和材料能带相符合的双带模型开展研究. 本文基于图7(b)的能带结构, 研究了电子关联效应对Ba2CuO$ _{4-\delta} $ 中量子相变的影响, 并在附录中对如何用DMFT处理图7(b)的能带结构进行了详细的介绍.on-site energy (ε) 1st hopping (t) 2nd hopping ($t^\prime$) 3rd hopping ($t^{\prime\prime}$) Orbital ${\rm{d}}_{x^2-y^2}$ –0.222 0.504 –0.067 0.130 Orbital ${\rm{d}}_{3 z^2-r^2}$ 0.661 0.196 0.026 0.029 inter-orbital 0 –0.302 0 –0.051 考虑非局域轨道间跃迁项的动量空间的各向异性, 将DFT模拟所得完整能带结构代入DMFT中, 计算半满Ba2CuO
$ _{4-\delta} $ 体系[36]的相图. 考虑完整能带结构的相图如图8所示, 相图中的相变点与相图6中存在差异, 主要是因为本文选用的紧束缚模型的能带结构与Ba2CuO$ _{4-\delta} $ 的能带结构差异比较大, 图8给出的相图能更加精确地描述半满Ba2CuO$ _{4-\delta} $ 体系的量子相变. 同时能带结构及相图的对比也说明了利用正则变换处理轨道间跃迁的局限性. -
运用以Lanczos方法为求解器的DMFT对相互作用U及非局域轨道间跃迁
$ t_{12} $ 在双轨道Hubbard模型中轨道选择Mott相变的影响进行了研究, 结果表明轨道间跃迁通过调整两有效轨道的带宽比来促进双轨道体系中轨道选择Mott相的发生. 对比了动力学平均场中使用不同晶格能带结构得到的相图, 尽管相变点存在差异但其描述的轨道选择Mott相变的物理图像具有一致性, 进一步说明了Bethe晶格作为近似的理论工具在研究关联效应驱动的量子相变中的泛用性. 最后将拓展的DMFT应用于双轨道铜氧化物超导材料Ba2CuO$ _{4-\delta} $ , 并采用动量空间各向异性的非局域轨道间跃迁的双带模型绘制了体系的相图, 发现半满的Ba2CuO$ _{4-\delta} $ 为轨道选择Mott材料.感谢中国科学院固体物理研究所邹良剑教授在本工作中的讨论.
-
在DMFT中考虑(20)式的紧束缚哈密顿量
$, $ 则映射之后的杂质模型为$ \begin{split} {H_{{\rm{imp}}}} =\;& \sum\limits_{ml\sigma} {{\varepsilon_{ml\sigma}}} c_{ml\sigma}^{\dagger} {c_{ml\sigma}}+ \mathop \sum \limits_{l\sigma} ( \varepsilon _l-\mu)d_{l\sigma}^{\dagger} {d_{l\sigma}} \\ &+ \sum\limits_{ll'm\sigma} {{V_{ll'm\sigma}}\big( {d_{l\sigma}^{\dagger} {c_{ml'\sigma}} + c_{ml'\sigma}^{\dagger} {d_{l\sigma}}} \big)} + {H_{\rm{I}}^{{\rm{imp}}}}, \end{split} \tag{A1}$ 其中l (l = 1, 2)为轨道指标, 与(8)式不同的地方在于杂质格点与外场格点的耦合系数
$ V_{ll'm\sigma} $ 含不同轨道间耦合的部分,$ V_{ll'm\sigma} $ 表示外场格点m的l轨道与杂质格点l' 轨道间的耦合. (A1)式对应的格林函数具有$ 2\times 2 $ 的矩阵形式:$ \begin{array}{l} {{\boldsymbol{{\rm{G}}}}_{{\rm{imp}}}}\left( {{\rm{i}}{\omega _n}} \right) = \left( {\begin{array}{*{20}{c}} {{G_{11}}\left( {{\rm{i}}{\omega _n}} \right)}&{{G_{12}}\left( {{\rm{i}}{\omega _n}} \right)}\\ {{G_{21}}\left( {{\rm{i}}{\omega _n}} \right)}&{{G_{22}}\left( {{\rm{i}}{\omega _n}} \right)} \end{array}} \right).\ \end{array} \tag{A2}$ 通过Lanczos精确对角化方法可以直接求得这个矩阵的对角元部分
$ G_{ll} $ [37-40], 求解方式与(10)式一致. 非对角元$ G_{l, l'} $ 可由对角元的线性组合求出:$ G_{12}^{}\left( {{\rm{i}}{\omega _n}} \right) = G_{21}^{}\left( {{\rm{i}}{\omega _n}} \right) =\frac{1}{2}\left[ {{\rm{G}}_{1+2,1+2}^{}\left( {{\rm{i}}{\omega _n}} \right) - \sum\limits_l {G_{ll}^{}\left( {{\rm{i}}{\omega _n}} \right)} } \right], \tag{A3} $ 其中
${\rm{G}}_{1+2, 1+2}$ 为两轨道结合的格林函数, 在Lanczos 算法中将构建三对角矩阵的初始矢量分别替换为$ \left( {d_{1}^{\dagger} + d_{2}^{\dagger} } \right)\left| {{\varPsi _0}} \right\rangle $ 和$ \left( {d_{1} + d_{2}} \right)\left| {{\varPsi _0}} \right\rangle $ 即可得到.杂质模型的Weiss场格林函数
${{{\cal {{\rm{G}}}}}}_0^{{\rm{imp}}}\left( {{\rm{i}}{\omega _n}} \right)$ 也为$ 2\times 2 $ 的矩阵, 可以通过杂质模型无相互作用部分的系数直接求出,$ \begin{array}{l} {\cal {{\rm{G}}}}_0^{{\rm{imp}}}{\left( {{\rm{i}}{\omega _n}} \right)^{ - 1}} = {\rm{i}}{\omega _n} + \mu - {\boldsymbol{\varDelta}} ({\rm{i}}{\omega _n}), \end{array} \tag{A4}$ 其中杂化函数
${\boldsymbol{\varDelta}} \left( {{\rm{i}}{\omega _n}} \right)$ 也为$ 2\times 2 $ 的矩阵:$ \begin{array}{l} {\boldsymbol{\varDelta}} \left( {{\rm{i}}{\omega _n}} \right) = \left( {\begin{array}{*{20}{c}} {{\varDelta _{11}}\left( {{\rm{i}}{\omega _n}} \right)}&{{\varDelta _{12}}\left( {{\rm{i}}{\omega _n}} \right)} \\ {{\varDelta _{21}}\left( {{\rm{i}}{\omega _n}} \right)}&{{\varDelta _{22}}\left( {{\rm{i}}{\omega _n}} \right)} \end{array}} \right), \end{array}\tag{A5} $ $ {\varDelta _{{l_1}{l_2}}}\left( {{\rm{i}}\omega _n} \right) $ 的表达形式为$ {\varDelta _{{l_1}{l_2}}}\left( {{\rm{i}}\omega _n} \right) \equiv \sum\limits_{ml} {\frac{{{V_{{l_1}lm}}{V_{{l_2}lm}}}}{{{\rm{i}}\omega - {\varepsilon_{ml}}}}}.\tag{A6}$ 将Lanczos方法得到的
${{\boldsymbol{{\rm{G}}}}_{{\rm{imp}}}}\left( {{\rm{i}}{\omega _n}} \right)$ 与Weiss场格林函数${\cal {{\rm{G}}}}_0^{{\rm{imp}}}\left( {{\rm{i}}{\omega _n}} \right)$ 相结合, 并利用Dyson方程即可求出杂质模型自能的矩阵形式:$ {{\boldsymbol{\varSigma}} _{{\rm{imp}}}}({\rm{i}}\omega _n) = {\cal {{\rm{G}}}}_0^{{\rm{imp}}}{\left( {{\rm{i}}{\omega _n}} \right)^{ - 1}}-{{\boldsymbol{{\rm{G}}}}_{{\rm{imp}}}}{\left( {{\rm{i}}{\omega _n}} \right)^{ - 1}}.\tag{A7} $ DMFT中杂质模型自能与格点模型的局域自能相等, 由此可计算格点格林函数
${{\boldsymbol{{\rm{G}}}}_{{\rm{lat}}}}\left( {{\rm{i}}\omega _n} \right)$ $ \begin{split} {{\boldsymbol{{\rm{G}}}}_{{\rm{lat}}}}\left( {{\rm{i}}\omega _n} \right) =\;& \frac{1}{N}\sum\limits_k {{\boldsymbol{G}}\left( {{\rm{i}}\omega _n,k} \right)} \\ = \;&\frac{1}{N}\sum\limits_k {\frac{1}{{{\rm{i}}\omega _n - {{\boldsymbol{H}}_0}\left( k \right) - {{\boldsymbol{\varSigma}} _{{\rm{imp}}}}({\rm{i}}\omega _n)}}}, \end{split}\tag{A8} $ 其中
$ {{\boldsymbol{H}}_0} $ 表示体系能带结构, 对半满情况下的Ba2CuO$ _{4-\delta} $ , 代入(20)式中的能带表达式即可计算. 最后通过调节自洽参数$ \varepsilon_{ml} $ 和$ V_{ll'm} $ 使得自洽条件${\bf G}_{{\rm{imp}}}({\rm{i}} \omega _n)= {\bf G}_{{\rm{lat}}}({\rm{i}}\omega _n)$ 成立后, 即可根据格林函数来了解Ba2CuO$ _{4-\delta} $ 体系的电子性质.
-
基于以Lanczos方法为杂质求解器的动力学平均场理论, 研究了非局域轨道间跃迁对于双轨道强关联体系中轨道选择Mott相变的影响. 计算了轨道间跃迁系数不同的双轨道Hubbard模型的准粒子权重和态密度, 并构建了它们在相互作用强度U 和轨道带宽比t2/t1 影响下的相图. 通过正则变换引入两个有效的退耦和轨道, 在一定条件下轨道间跃迁会有利于轨道选择Mott相变的发生. 还比较了Bethe晶格和正方晶格的相图, 虽然基于两种不同晶格能带结构得到的轨道选择Mott相变的相变点存在一定的差异, 但其中关于轨道选择Mott相变的基本物理图像具有一致性. 并将方法拓展到半满的Ba2CuO4–δ材料的研究中, 与根据密度泛函理论得到的能带对比, 我们发现各向同性的轨道间跃迁对能带结构影响较大, 进一步采用动量空间各向异性的非局域轨道间跃迁项, 得到了材料的相图, 在半满条件下Ba2CuO4–δ应为轨道选择Mott材料.
-
关键词:
- 动力学平均场理论 /
- Mott相变 /
- 双轨道Hubbard模型
The effect of interorbital hopping on orbital selective Mottness in a two-band correlation system is investigated by using the dynamical mean-field theory with the Lanczos method as impurity solver. The phase diagrams of the two-orbital Hubbard model with non-local interorbital hopping (t12) , where the orbital selective Mott phases (OSMP) region is expanded by the increasing of the interorbital hopping. We compare the results obtained by self-consistent relations of Bethe lattice and squate lattice based on DMFT procedure, and the general OSMP physics of Bethe lattice is consistent with that of the square lattice while the critical points of two kinds of lattices are different. We extend the method to the study of half-filled Ba2CuO4–δ materials. By comparing with the band structure obtained from the density functional theory (DFT), it is found that the isotropic inter-orbital hopping has a great influence on the bandstructure. The DFT bandstructure in DMFT is considered, and the phase diagram of the material is obtained. The half-filled Ba2CuO4–δ should be orbital-selective Mott compound according to our results.[1] Anisimov V I, Nekrasov I A, Kondakov D E, Rice T M, Sigrist M 2002 Eur. Phys. J. B 25 191
[2] Zhong Z, Hansmann P 2017 Phys. Rev. X 7 011023
[3] Huang J, Yu R, Xu Z, Zhu J X, Jiang Q, Wang M, Wu H, Chen T, Denlinger J D, Mo S K, Hashimoto M, Gu G, Dai P, Chu J H, Lu D, Si Q, Birgeneau R J, Yi M 2021 arXiv: 2010.13913
[4] Oleś A M 1983 Phys. Rev. B 28 327
Google Scholar
[5] Koga A, Imai Y, Kawakami N 2002 Phys. Rev. B 66 165107
Google Scholar
[6] Koga A, Kawakami N, Rice T M, Sigrist M 2005 Phys. Rev. B 72 045128
Google Scholar
[7] Kubo K 2007 Phys. Rev. B 75 224509
Google Scholar
[8] Georges A, Kotliar G, Krauth W, Rozenberg M J 1996 Rev. Mod. Phys. 68 13
Google Scholar
[9] Held K, Keller G, Eyert V, Vollhardt D, Anisimov V I 2001 Phys. Rev. Lett. 86 5345
Google Scholar
[10] Kotliar G, Savrasov S Y, Haule K, Oudovenko V S, Parcollet O, Marianetti C A 2006 Rev. Mod. Phys. 78 865
Google Scholar
[11] Zhao J Z, Zhuang J N, Deng X Y, Bi Y, Cai L C, Fang Z, Dai X 2012 Chin. Phys. B 21 057106
Google Scholar
[12] Liebsch A 2005 Phys. Rev. Lett. 95 116402
Google Scholar
[13] Peters R, Kawakami N, Pruschke T 2011 Phys. Rev. B 83 125110
Google Scholar
[14] Hirsch J E 1983 Phys. Rev. Lett. 51 1900
Google Scholar
[15] Hirsch J E 1985 Phys. Rev. B 31 4403
Google Scholar
[16] Chang C C, Zhang S 2008 Phys. Rev. B 78 165101
Google Scholar
[17] Tocchio L F, Becca F, Sorella S 2016 Phys. Rev. B 94 195126
Google Scholar
[18] Hanke W, Hirsch J E, 1982 Phys. Rev. B 25 6748
Google Scholar
[19] Shankar R 1994 Rev. Mod. Phys. 66 129
Google Scholar
[20] Hille C, Kugler F B, Eckhardt C J, He Y Y, Kauch A, Honerkamp C, Toschi A, Andergassen S 2020 Phys. Rev. Research 2 033372
Google Scholar
[21] Kotliar G, Ruckenstein A E 1986 Phys. Rev. Lett. 57 1362
Google Scholar
[22] Lechermann F, A Georges, Kotliar G, Parcollet O 2007 Phys. Rev. B 76 155102
Google Scholar
[23] Yu R, Si Q 2012 Phys. Rev. B 86 085104
Google Scholar
[24] Imada M, Fujimori A, Tokura Y 1998 Rev. Mod. Phys. 70 1039
Google Scholar
[25] Lee P A, Nagaosa N, Wen X G 2006 Rev. Mod. Phys. 78 17
Google Scholar
[26] Rohringer G, Hafermann H, Toschi A, Katanin A A, Antipov A E, Katsnelson M I, Lichtenstein A I, Rubtsov A N, Held K 2018 Rev. Mod. Phys. 90 025003
Google Scholar
[27] Song Y, Zou L J 2005 Phys. Rev. B 72 085114
Google Scholar
[28] Song Y, Zou L J 2009 Eur. Phys. J. B 72 59
Google Scholar
[29] Li W M, Zhao J F, Cao L P, Hu Z, Huang Q Z, Wang X C, Y. Liu, Zhao G Q, Zhang J, Liu Q Q, Yu R Z, Long Y W, Wu H, Lin H J, Chen C T, Li Z, Gong Z Z, Guguchia Z, Kim J S, Stewart G R, Uemura Y J, Uchida S, Jin C Q 2019 PNAS 116 12156
Google Scholar
[30] Liu K, Lu Z, Xiang T 2019 Phys. Rev. Mater. 3 044802
Google Scholar
[31] Li Y, Du S, Weng Z, Liu Z 2020 Phys. Rev. Mater. 4 044801
Google Scholar
[32] Wang Z, Zhou S, Chen W, Zhang F 2020 Phys. Rev. B 101 180509
Google Scholar
[33] Maier T, Berlijn T, Scalapino D J 2019 Phys. Rev. B 99 224515
Google Scholar
[34] Jiang K, Le C, Li Y, Qin S, Wang Z, Zhang F, Hu J 2021 Phys. Rev. B 103 045108
Google Scholar
[35] Jin H, Pickett W E, Lee K 2021 Phys. Rev. B 104 054516
Google Scholar
[36] Ni Y, Quan Y M, Liu J, Song Y, Zou L J 2021 Phys. Rev. B 103 214510
Google Scholar
[37] Niu Y, Sun J, Ni Y, Liu J, Song Y, Feng S 2019 Phys. Rev. B 100 075158
Google Scholar
[38] Dagotto E 1994 Rev. Mod. Phys. 66 763
Google Scholar
[39] Caffarel M, Krauth W 1994 Phys. Rev. Lett. 72 1545
Google Scholar
[40] Capone M, de’Medici L, Georges A 2007 Phys. Rev. B 76 245116
Google Scholar
[41] Perroni C A, Ishida H, Liebsch A 2007 Phys. Rev. B 75 045125
Google Scholar
[42] Georges A, de’Medici L, Mravlje J 2013 Annu. Rev. Condens. Matter Phys. 4 137
Google Scholar
[43] Koga A, Kawakami N, Rice T M, Sigrist M 2004 Phys. Rev. Lett. 92 216402
Google Scholar
[44] Costi T A, Liebsch A 2007 Phys. Rev. Lett. 99 236404
Google Scholar
[45] de'Medici L, Mravlje J, Georges A 2011 Phys. Rev. Lett. 107 256401
Google Scholar
[46] 孙健, 刘洋, 宋筠 2015 物理学报 64 247101
Google Scholar
Sun J, Liu Y, Song Y 2015 Acta Phys. Sin. 64 247101
Google Scholar
[47] de’Medici L 2011 Phys. Rev. B 83 205112
Google Scholar
-
图 1 不同轨道间跃迁系数下两个有效轨道的准粒子权重
$Z$ 随相互作用U 的变化曲线. 轨道α 和轨道β分别用黑色圆圈和红色三角标识, 实心标志和空心标志分别代表$t_{12}=0$ 和$t_{12}=0.4$ 的情况. 体系中的洪特耦合作用$ J_{\rm{H}}=0.25 U$ , 两轨道带宽比$t_2/t_1=0.4$ , 所选用的能量单位为宽带跃迁系数$t_1$ Fig. 1. The interaction dependencies of the quasiparticle weight Z of the effective orbital α (black circles) and β (red triangles) with interorbital hopping
$t_{12}=0.4$ (hollow symbols) and without interorbital hopping$t_{12}=0$ (solid symbols) when$J_{\rm{H}}=0.25 U$ and$t_2/t_1=0.4$ . The energy is in unit$t_1$ 图 2 两个有效轨道的态密度随相互作用U 增加的演化图像, 体系中的轨道间跃迁
$t_{12}=0.4$ , 带宽比$t_2/t_1=0.4$ 以及相互作用强度$J_{\rm{H}}=0.25 U.$ 左边的图像为α轨道的结果, 右边图像为β轨道的结果, 计算中选取的能量展宽因子为$\eta=0.05$ Fig. 2. Evolution of the orbital-resolved spectral density
$A(\omega)$ with the increasing Coulomb interaction U when$t_{12}=0.4$ ,$t_2/t_1=0.4$ and$J_{\rm{H}}=0.25 U.$ Left panels show the results of effective orbital α and right panels are for effective orbital β. The energy broadening factor in our calculation takes$\eta=0.05$ 图 3 轨道间跃迁系数不同时, 相互作用U和带宽比
$t_2/t_1$ 影响下的双轨道体系相图 (a) 轨道间跃迁$t_{12}=0$ ; (b) 轨道间跃迁$t_{12}=0.4$ . 体系的洪特耦合$J_{\rm{H}}=0.25 U$ . 黑色圆圈和红色三角分别代表轨道α和β 的MIT 相变点Fig. 3. Phase diagrams in the plane of interaction U and hopping integral
$t_2/t_1$ of the effective two-orbital Hubbard model with different interorbital hopping, when$J_{\rm{H}}=0.25 U$ : (a)$t_{12}=0$ ; (b)$t_{12}=0.4$ . The black circles (red triangles) denote the critical points of MIT for effective orbital α(β)图 4 DMFT结合正方晶格能带结构所得到的结果 (a)轨道间跃迁系数不同时, 准粒子权重Z 随相互作用U 的变化曲线, 体系的洪特耦合
$J_{\rm{H}}$ = 0.25U, 带宽比$t_2$ /$t_1$ = 0.4; (b)$t_{12}=0.4$ 以及$J_{\rm{H}}=0.25 U$ 时, 相互作用U 和带宽比$t_2/t_1$ 影响下的有效轨道模型相图Fig. 4. DMFT results of the square lattice: (a) The interaction dependencies of the quasiparticle weight Z with different interorbital hoppings when
$J_{\rm{H}}$ = 0.25U and$t_2$ /$t_1$ = 0.4; (b) phase diagrams in the plane of interaction U and hopping integral$t_2/t_1$ of the effective two-orbital Hubbard model when$t_{12}=0.4$ and$J_{\rm{H}}=0.25 U$ 图 5 (a)
$J_{\rm{H}}=0.25 U$ 时有效轨道的准粒子权重随相互作用强度的变化曲线. 不同相互作用强度U下体系各有效轨道的态密度图像 (b)$U=1.0$ eV; (c)$U=2.2$ eV; (d)$U=2.5$ eV. 当相互作用2.1 eV$\leqslant U < $ 2.3 eV时体系处于OSMPFig. 5. (a) Quasiparticle weight Z as a function of interaction U when
$J_{\rm{H}}=0.25 U$ . The orbital-resolved spectral density$A(\omega)$ with different intraorbital interaction: (b)$U=1.0$ eV; (c)$U=2.2$ eV; (d)$U=2.5$ eV for the two effective orbitals. An OSMP occurs in a narrow interaction region with 2.1 eV$\leqslant U < $ 2.3 eV图 6 相互作用U和洪特耦合
$J_{\rm{H}}$ 影响下体系的相图. 随着相互作用U的增加以及洪特耦合$J_{\rm{H}}$ 的减小, OSMP区域逐渐缩小, 大约在$J_{\rm{H}}=0.34$ eV以及$U=2.7$ eV时体系中的OSMT消失Fig. 6. The phase diagram of the effective two-orbital Hubbard model with interaction U and Hund's rule coupling
$J_{\rm{H}}$ . The region of OSMP becomes narrower with the decreasing$J_{\rm{H}}$ and increasing U, and OSMT vanishes around$J_{\rm{H}}=0.34$ eV and$U=2.7$ eV图 7 由不同模型得到的半满Ba2CuO
$_{4-\delta}$ 体系的能带图 (a)$t_{1}$ =0.504 eV,$t_{2}$ =0.196 eV,$t_{12}=-0.302$ eV,$\mu_1=-0.222$ eV 和$\mu_2=0.661$ eV的紧束缚模型所得到的能带图; (b) 根据DFT模拟的双带模型能带图像[33]Fig. 7. Bandstructure for half-filling Ba2CuO3.5 with different model: (a) Tight binding model with
$t_{1}$ =0.504 eV,$t_{2}$ =0.196 eV,$t_{12}=-0.302$ eV,$\mu_1=-0.222$ eV and$\mu_2=0.661$ eV; (b) two-orbital model according to DFT results[33].表 1 半满Ba2CuO4–δ紧束缚哈密顿量中的参数(单位为eV)[33]
Table 1. Model parameters of the TB Hamiltonian of Ba2CuO4–δ at half-filling (unit: eV)[33]
on-site energy (ε) 1st hopping (t) 2nd hopping ($t^\prime$) 3rd hopping ($t^{\prime\prime}$) Orbital ${\rm{d}}_{x^2-y^2}$ –0.222 0.504 –0.067 0.130 Orbital ${\rm{d}}_{3 z^2-r^2}$ 0.661 0.196 0.026 0.029 inter-orbital 0 –0.302 0 –0.051 -
[1] Anisimov V I, Nekrasov I A, Kondakov D E, Rice T M, Sigrist M 2002 Eur. Phys. J. B 25 191
[2] Zhong Z, Hansmann P 2017 Phys. Rev. X 7 011023
[3] Huang J, Yu R, Xu Z, Zhu J X, Jiang Q, Wang M, Wu H, Chen T, Denlinger J D, Mo S K, Hashimoto M, Gu G, Dai P, Chu J H, Lu D, Si Q, Birgeneau R J, Yi M 2021 arXiv: 2010.13913
[4] Oleś A M 1983 Phys. Rev. B 28 327
Google Scholar
[5] Koga A, Imai Y, Kawakami N 2002 Phys. Rev. B 66 165107
Google Scholar
[6] Koga A, Kawakami N, Rice T M, Sigrist M 2005 Phys. Rev. B 72 045128
Google Scholar
[7] Kubo K 2007 Phys. Rev. B 75 224509
Google Scholar
[8] Georges A, Kotliar G, Krauth W, Rozenberg M J 1996 Rev. Mod. Phys. 68 13
Google Scholar
[9] Held K, Keller G, Eyert V, Vollhardt D, Anisimov V I 2001 Phys. Rev. Lett. 86 5345
Google Scholar
[10] Kotliar G, Savrasov S Y, Haule K, Oudovenko V S, Parcollet O, Marianetti C A 2006 Rev. Mod. Phys. 78 865
Google Scholar
[11] Zhao J Z, Zhuang J N, Deng X Y, Bi Y, Cai L C, Fang Z, Dai X 2012 Chin. Phys. B 21 057106
Google Scholar
[12] Liebsch A 2005 Phys. Rev. Lett. 95 116402
Google Scholar
[13] Peters R, Kawakami N, Pruschke T 2011 Phys. Rev. B 83 125110
Google Scholar
[14] Hirsch J E 1983 Phys. Rev. Lett. 51 1900
Google Scholar
[15] Hirsch J E 1985 Phys. Rev. B 31 4403
Google Scholar
[16] Chang C C, Zhang S 2008 Phys. Rev. B 78 165101
Google Scholar
[17] Tocchio L F, Becca F, Sorella S 2016 Phys. Rev. B 94 195126
Google Scholar
[18] Hanke W, Hirsch J E, 1982 Phys. Rev. B 25 6748
Google Scholar
[19] Shankar R 1994 Rev. Mod. Phys. 66 129
Google Scholar
[20] Hille C, Kugler F B, Eckhardt C J, He Y Y, Kauch A, Honerkamp C, Toschi A, Andergassen S 2020 Phys. Rev. Research 2 033372
Google Scholar
[21] Kotliar G, Ruckenstein A E 1986 Phys. Rev. Lett. 57 1362
Google Scholar
[22] Lechermann F, A Georges, Kotliar G, Parcollet O 2007 Phys. Rev. B 76 155102
Google Scholar
[23] Yu R, Si Q 2012 Phys. Rev. B 86 085104
Google Scholar
[24] Imada M, Fujimori A, Tokura Y 1998 Rev. Mod. Phys. 70 1039
Google Scholar
[25] Lee P A, Nagaosa N, Wen X G 2006 Rev. Mod. Phys. 78 17
Google Scholar
[26] Rohringer G, Hafermann H, Toschi A, Katanin A A, Antipov A E, Katsnelson M I, Lichtenstein A I, Rubtsov A N, Held K 2018 Rev. Mod. Phys. 90 025003
Google Scholar
[27] Song Y, Zou L J 2005 Phys. Rev. B 72 085114
Google Scholar
[28] Song Y, Zou L J 2009 Eur. Phys. J. B 72 59
Google Scholar
[29] Li W M, Zhao J F, Cao L P, Hu Z, Huang Q Z, Wang X C, Y. Liu, Zhao G Q, Zhang J, Liu Q Q, Yu R Z, Long Y W, Wu H, Lin H J, Chen C T, Li Z, Gong Z Z, Guguchia Z, Kim J S, Stewart G R, Uemura Y J, Uchida S, Jin C Q 2019 PNAS 116 12156
Google Scholar
[30] Liu K, Lu Z, Xiang T 2019 Phys. Rev. Mater. 3 044802
Google Scholar
[31] Li Y, Du S, Weng Z, Liu Z 2020 Phys. Rev. Mater. 4 044801
Google Scholar
[32] Wang Z, Zhou S, Chen W, Zhang F 2020 Phys. Rev. B 101 180509
Google Scholar
[33] Maier T, Berlijn T, Scalapino D J 2019 Phys. Rev. B 99 224515
Google Scholar
[34] Jiang K, Le C, Li Y, Qin S, Wang Z, Zhang F, Hu J 2021 Phys. Rev. B 103 045108
Google Scholar
[35] Jin H, Pickett W E, Lee K 2021 Phys. Rev. B 104 054516
Google Scholar
[36] Ni Y, Quan Y M, Liu J, Song Y, Zou L J 2021 Phys. Rev. B 103 214510
Google Scholar
[37] Niu Y, Sun J, Ni Y, Liu J, Song Y, Feng S 2019 Phys. Rev. B 100 075158
Google Scholar
[38] Dagotto E 1994 Rev. Mod. Phys. 66 763
Google Scholar
[39] Caffarel M, Krauth W 1994 Phys. Rev. Lett. 72 1545
Google Scholar
[40] Capone M, de’Medici L, Georges A 2007 Phys. Rev. B 76 245116
Google Scholar
[41] Perroni C A, Ishida H, Liebsch A 2007 Phys. Rev. B 75 045125
Google Scholar
[42] Georges A, de’Medici L, Mravlje J 2013 Annu. Rev. Condens. Matter Phys. 4 137
Google Scholar
[43] Koga A, Kawakami N, Rice T M, Sigrist M 2004 Phys. Rev. Lett. 92 216402
Google Scholar
[44] Costi T A, Liebsch A 2007 Phys. Rev. Lett. 99 236404
Google Scholar
[45] de'Medici L, Mravlje J, Georges A 2011 Phys. Rev. Lett. 107 256401
Google Scholar
[46] 孙健, 刘洋, 宋筠 2015 物理学报 64 247101
Google Scholar
Sun J, Liu Y, Song Y 2015 Acta Phys. Sin. 64 247101
Google Scholar
[47] de’Medici L 2011 Phys. Rev. B 83 205112
Google Scholar
计量
- 文章访问数: 1522
- PDF下载量: 93
- 被引次数: 0