第二节非线性光学极化率！.docVIP

第二节 非线性光学极化率 一 密度矩阵表述法 （一）刘维方程: 非线性光学极化率是介质的特征性质――与介质的电子和分子结构的细节有关――量子力学计算――密度矩阵表述法――最方便的方法，特别当必须处理激发的弛豫时. 令是在电磁场影响下物质系统的波函数. 密度矩阵算符： （2.1.1） 物理量P的系综平均由下式给出： （2.1.2） （2.1.3） 该方程称作刘维方程（Liouville’s equation）. 哈密顿算符是由三部分组成： （2.1.4） 1)是未受扰动的物质系统的哈密顿算符，其本征态是，而本征能量是，； 2)是描述光与物质相互作用的相互作用哈密顿算符； 3)而是描述系统周围的热库施于该系统随机的扰动的哈密顿算符. Hint 在电偶极矩近似下，相互作用哈密顿算符由下式给定： （2.1.5） 在这里将只考察电子对极化率的贡献. 对于离子的贡献，就必须用—代替,其中qi和分别是第i个离子的电荷和位置. H随机 哈密顿算符是造成物质激发的弛豫的原因，或者换言之，它是造成被扰动了的ρ弛豫回到热平衡的原因. 于是我们可以把式（2.1.3）表示成 （2.1.6） 其中 ρ的矩阵元的物理意义: 将本征态作为基矢，并把写成的线性组合： ，那么，ρ的矩阵元的物理意义就十分清楚了. 矩阵元表示系统在态中的布居， 而非对角矩阵元表明系统的态具有和的相干混合. 在和有混合的情况下，如果与的相对相位是随机的（或不相干的），那么，通过系综平均后就有。 寻找（）弛豫表达式. 布居的弛豫是系统与热库的相互作用引起的态之间的跃迁的结果.令Wn-n’是由热引起的丛态到态的跃迁的速率.于是，中的过剩布居的弛豫速率应是 弛豫= （2.1.8） 在热平衡时，就有 （2.1.9） 因此，也可以把式（2.1.8）写成 （2.1.10） 非对角元的弛豫更复杂. 然而，在一些简单的情况中，预期相位相干性指数的衰减到零.这样，对于nn’,我们有 （2.1.11） 这里是态与之间的特征弛豫时间.在磁共振中，布居的弛豫称作纵向弛豫，而非对角矩阵元的弛豫称作横向弛豫. 在某些情况下，态的纵向弛豫能用下式来近似： （2.1.12） （二）微扰法解刘维方程 在计算中采用微扰展开. 令 （2.1.13） 其中 （2.1.14） 式中是热平衡的系统的密度矩阵算符，而且我们假设在介质中没有固有极化，因而. 把的级数展开式代入式（2.1.6），再把视为一级微扰，相同级的相收集在一起，就得到 （2.1.15）  我们在这里感兴趣的是对能分解成傅立叶分量的场 ?i 的响应. 于是，由于 和 算符也能展开成傅立叶级数 当时，就能从式（2.1.15）具体的逐级解出.第一级解是 （2.1.16） 这里我们采用了记号. 可以很容易得到更高级的解，尽管这种推倒是冗长乏味的，每当在推导中出现对角元时，为了得到一个封闭的解，常常必须对式（2.1.8）中的作进一步的近似. 我们还需提及，只要式（2.1.16）中的表达式即使在n=n’时也是适用的，因为那时可在计算机中略去这一项. 二． 非线性极化率的微观表达式 非线性极化强度和非线性极化率的完全的微观表达式得到的. 在式（2.1.14）和（2.1.16）中，当Hint=e和时，很容易得到由电子贡献引起的一阶和二阶极化率. 用明显的笛卡儿张量标记，这些极化率就由下列各式给出： 一阶： χij(1)=pi1(1)(ω)/Ej(ω)= 注意：ij＝1，2，3 共有9个分量。 二阶： (2.2.) 在中有两项，而在中有8项. 注意： 有27个分量 三阶：（），它总共48项. 在文献（5）中给出了的完全表达式，这里就不在重述了. 的共振结构以后要在第十四章里讨论. 在非共振的情况下，可以忽略式（2.1.17）的