基于Kramers-Kronig关系建立金属太赫兹色散模型

提出了一种基于测量反射率谱、使用Kramers-Kronig(KK)关系建立金属太赫兹色散模型的方法.结合合金铝和合金铜4—40 THz的测量反射率谱,通过反射系数振幅和相位的KK关系,采用高频端指数外推,低频端常数外推的方法,反演金属复折射率.以KK反演的复折射率作为实验值,以拟合复折射率和实验值误差最小为准则,使用遗传优化算法,拟合了合金铝和合金铜的Drude色散参数(等离子频率和碰撞频率).基于优化的Drude模型计算了0.1—40 THz材料的复折射率,与椭偏仪的实测结果符合,验证了模型的准确性.该方法理论与实验相互验证,以测量的复折射率作为实验定标,将远红外频段的色散信息拓展到太赫兹频域,确定了太赫兹频段金属的微观物理参数,提供了太赫兹频段色散和散射机理的研究依据.

Establishment of THz dispersion model of metals based on Kramers-Kronig relation

The extraction of terahertz dispersion parameters is confined in a limited region due to the limitation of the existing THz techniques. A method of studying the dispersion model of metals from the measurements of reflection spectrum and analysis of Kramers-Kronig (KK) relation is proposed. The reflection spectrum is measured by Vertex 80V Fourier transform spectrometer. In order to eliminate the signal noise of measured reflection spectrum, the measured spectrum is smoothed by Drude estimation. Using the smoothed reflection spectra of copper (Cu) alloy and aluminum (Al) alloy in a range of 4–40 THz, the complex refractivities are inversed based on the KK relation of amplitude and phase of reflective coefficient. The constant extrapolations at lower frequencies and the exponential extrapolation at higher frequencies are adopted in the KK integration. The exponential extrapolation index is adjusted according to the calibrating complex refractivity measured from far-infrared ellipsometer. According to the inversed complex refractivity, the plasma frequency and damping frequency in Drude model are optimized using the genetic algorithm. The objective function is defined as the error between the fitted complex refractivity and KK inversion. Since the optimal plasma frequency and damping frequency are different for different fitting frequencies, the obtained Drude parameters are averaged in order to reduce the influences of errors from KK inversion, measured reflection spectrum and calibrations. The complex refractivity indexes in a range from 15 THz to 40 THz, calculated by the established Drude model, are in good agreement with the measured calibrations from ellipsometer, which demonstrates the accuracy of the established Drude dispersion model. The reflection spectra below 4 THz are greatly distorted due to the signal noise, and the calibrating refractivity is located in the far infrared region, thus the complex refractivity is inversed in a region of 4–40 THz by KK algorithm. The complex refractivity indexes in a range of 0.1–20 THz, obtained by the proposed scheme, are for the vacancy, which will provide great support for the dispersion analysis in the whole terahertz gap. The procedures are helpful for extrapolating the dispersion information to terahertz band from the far infrared region. The scheme takes the advantage of the spectrometer and ellipsometer, and it requires high experimental precisions of reflection spectrum and calibrating refractivity. In addition, the scheme is adaptive to both metals and nonmetals by applying proper dispersion model which depends on the property of the reflection spectrum. The established model determines the microscopic dispersion parameters of material, which provides great support for the investigation of terahertz dispersion analysis, scattering mechanisms and imaging processes.