光弹性常数

折射率(介电常数)相对于内部或施加的应变的变化构成了晶体的所谓光弹性或弹性光学布里渊散射，即光与布里渊区中心附近的声分支声子在互反空间中的相互作用，是测定弹性和弹性光学波克尔斯常数的有效实验技术。特别地，从布里渊散射的现象学理论出发，用布里渊分量的强度给出了弹性光学常数。对晶体这种特殊的、基本的光学特性的微观建模和理解与光学和微电子学领域有关:光弹性确实降低了α-石英光纤布拉格光栅的效率，并降低了用于光刻设计晶体管的纯硅光学透镜的分辨率。晶体化合物通常的各向异性增加了许多物理性质的信息量:介电常数成为二级张量ζ，并可能出现双折射现象，这对应于光沿着不同的晶体学方向以不同的速度传播弹性光学常数成为四阶波克尔斯张量P的分量，它将介电张量ζ−1的逆与二阶弹性应变张量η线性联系起来。一般来说，Pockels张量有21个独立的分量，对于立方晶体，系统的特定对称性可以进一步减少到至少3个。

从基本的观点来看，在20世纪60年代和70年代，大量的注意力被投入到非常简单的晶体的光弹性的实验表征上，如面心立方碱卤化物8 - 10，如MgO, NaCl, NaF, LiF, KCl, KBr和KI，立方体类金刚石半导体8,11,12，如Si, SiC, BN, AlN, GaN和六方绿柱石，α-CdS， α-ZnS和ZnO。13-15在更复杂的晶体中，如氧化钛TiO2金红石、钛酸钡BaTiO3钙钛矿具有临界极化波动效应、磷酸二氢钾KH2PO4具有异常居里-魏斯温度依赖性、硅酸盐玻璃、20、21 α-石英SiO2、22、23等，发现了技术应用的有趣特征和固体光弹性行为的异常。由于光弹性的大多数特殊方面表现在相对复杂的晶体结构中，因此预计需要相当复杂的理论计算来准确地确认甚至预测这些特征。

本文的目的是提出一种精确的、全自动的从头计算方案，用于模拟和预测晶体的光弹性常数，该方案是在周期边界条件的形式化框架内发展起来的，单电子哈密顿量[har树- fock (HF)或Kohn-Sham (KS)]，并使用原子中心高斯型函数(GTF)基集(全电子基集可以很容易地使用)。本实现的最有趣的特征如下:(i)在对称空间群上的一般性;(ii)充分利用晶体的点对称性和平移对称性;(iii)分离电子和核“束缚离子”贡献的可能性;(iv)通过解析耦合摄动HF或KS方案(分别为CPHF或CPKS)评估静态介电张量的电子贡献;(v)波克尔斯张量的电场频率依赖性。

迄今为止，文献中很少报道晶体光弹性的从头算研究:Levine和同事在伪势平面波方案中进行了局部密度近似(LDA)计算，并进行了剪刀校正，24,25 Detraux和Gonze报道了LDA中α-石英的密度-函式微扰理论研究，并进行了剪刀式修正;26 Donadio等人将相同的形式应用于晶体和非晶态二氧化硅的研究;27,28和Hounsome等人使用剪刀式修正的LDA计算了硅和金刚石的弹性光学常数与电场的依赖关系