非线性光学极化率的经典描述

给研究生的建议：
1.不像教学由一般到特殊的演绎法，做研究可以从特殊到一般，比如由分立频率到连续频率
2.不要总想着记已知的（应试思维）“尽信书，则不如无书”，而要思考怎么来的，并且学会创造。比如从简谐振子模型开始不断添加项，假设是线性关系…

一维振子的线性响应

一维自由振动简谐振子

质量为m的质点固定在弹性系数为k的弹簧的一端，弹簧另一端固定，忽略弹簧的质量，任何摩擦以及重力
由胡克定律，质点受力$F=-kx$
由牛顿第二定律，得运动方程$F=m\ddot{x}=-kx=-m{\omega }_0^{2}x$
微分方程通解$$x=A\cos ({\omega }_{0}t+{\varphi }_0)\left({\omega }_0=\sqrt{k/m}\right)$$

带阻尼的简谐振子（damped oscillator）

设阻力正比于速度$f=\alpha \dot{x}$
则阻尼振动运动方程$F=m\ddot{x}=-m{\omega }_0^{2}x-\alpha \dot{x}$
即$\ddot{x}+\frac{\alpha }{m}\dot{x}+{\omega }_0^{2}x=0$
定义阻尼系数$h=\frac{\alpha }{2m}$
简谐运动之阻尼振动 - 知乎 (zhihu.com)
[朗道力学笔记]（六）自由振动、强迫振动、共振、听一听“拍”的声音 - 知乎 (zhihu.com)

介质振子模型

设介质是一个含有固有振动频率为${\omega }_0$的振子的集合，振子模型是描述原子中电子运动的一种粗略模型。假设介质中每一个原子中的电子都受一个弹性恢复力的作用而保持在平衡位置上；阻尼对应能量耗散（光的辐射强度）；当原子受外加光场作用，原子中电子作受外力的强迫振动
$F=m\ddot{x}=-m{\omega }_0^{2}x-2mh\dot{x}+(-e)E$
运动方程（r是原子核与电子之间的距离）
$\frac{d^{2}r}{d{t}^2}+2h\frac{dr}{dt}+{\omega }_0^{2}r=-\frac{e}{m}E$
运动方程解法不唯一

这里采用傅里叶分析($\omega$是光场频率)
$\begin{array}{rl}& r(t)={\int }_{-\infty }^{\infty }r(\omega )e^{-i\omega t}d\omega \\ & \\ & E(t)={\int }_{-\infty }^{\infty }E(\omega )e^{-i\omega t}d\omega \end{array}$
对方程两边同时傅里叶变换
$-{\omega }^{2}r(\omega )-2ih\omega r(\omega )+{\omega }_0^{2}r(\omega )=-\frac{e}{m}E(\omega )$

即$r(\omega )=-\frac{e}{m}E(\omega )\frac{1}{{\omega }_0^2-{\omega }^2-2ih\omega }$

单位体积电偶极矩复振幅
$$P(\omega )=-ner(\omega )=\frac{n{e}^2}{m}E(\omega )\frac{1}{{\omega }_0^2-{\omega }^2-2\text{i}h\omega }\propto E\left(\omega \right)$$
频率域内线性极化强度复振幅$P^{(1)}(\omega )$与光电场复振幅$E(\omega )$的定义关系式$P^{(1)}(\omega )={\epsilon }_0{X}^{(1)}(\omega )\bullet E(\omega )$
$${\chi }^{(1)}(\omega )=\frac{P(\omega )}{{\epsilon }_{0}E(\omega )}=\frac{n{e}^2}{{\epsilon }_{0}m}\frac{1}{{\omega }_0^2-{\omega }^2-2ih\omega }$$

为什么要引入极化率？∵能反映物质本身的性质(强度量？) ，如$m=\rho V$，$F=Eq$

引入$F(\omega )=\frac{1}{{\omega }_0^2-{\omega }^2-2\mathrm{i}h\mathbf{\omega }}$
则${\chi }^{(1)}(\omega )=\frac{n{e}^2}{{\epsilon }_{0}m}F(\omega )={\chi }^{\prime }(\omega )+\mathrm{i}{\chi }^{\prime \prime }(\omega )$
其中
$$\begin{array}{rl}{\chi }^{\prime }\left(\omega \right)& =\frac{n{e}^2}{{\epsilon }_{0}m}\frac{{\omega }_0^2-{\omega }^2}{{\left({\omega }_0^2-{\omega }^2\right)}^2+4h^2{\omega }^2}\\ {\chi }^{\prime \prime }(\omega )& =\frac{n{e}^2}{{\epsilon }_{0}m}\frac{2h\omega }{{\left({\omega }_0^2-{\omega }^2\right)}^2+4h^2{\omega }^2}\end{array}$$
阻尼h引起虚部，虚部描述介质对光传输的吸收（或放大）
当频率$\omega$远离共振频率${\omega }_0$，频率$\omega$的光波在介质中无吸收地传输

一维振子的非线性响应

振子恢复力添加小的非简谐项（考虑到三次）
$F=-m{\omega }_0^{2}r+mA{r}^2+mB{r}^3$
$F=m\ddot{r}=-m{\omega }_0^{2}r+mA{r}^2+mB{r}^3-2mh\dot{r}-eE$
振子运动方程
$\frac{{\mathrm{d}}^{2}r}{\mathrm{d}{t}^2}+2h\frac{\mathrm{d}r}{\mathrm{d}t}+{\omega }_0^{2}r-A{r}^2-B{r}^3=-\frac{e}{m}E$

单个频率光场

用微扰方法求解一维振子非线性效应
假设频率为$\omega$的光电场表示为
$$E(t)=E(\omega )e^{-i\omega t}+E^{\ast }(\omega )e^{i\omega t}=E(\omega )e^{-i\omega t}+c.c.$$
线性响应有 $\vec{r}(t)\propto \vec{E}(t)$
$$r(t)=r\left(\omega \right)e^{-i\omega t}+r^{\ast }\left(\omega \right)e^{i\omega t}=r\left(\omega \right)e^{-i\omega t}+c.c.$$

包含多个频率分量光电场

非简谐振子模型反映线性光学与非线性光学的区别

线性光学	非线性光学
观测量正比于外加光场光强
不同频率的光波之间不交换能量	可能产生新的光波频率
电极化强度正比于外加光场