基于移位算子时域有限差分的色散薄层节点修正算法

提出了一种时域有限差分（FDTD）计算中色散介质薄层问题处理的新算法.对于厚度小于一个元胞尺度的电小尺寸色散介质薄层问题,采用将元胞内电位移矢量和磁感应强度加权平均的方法,求得薄层所在元胞内修正点处的等效介质参数.然后根据常见色散介质模型,包括Debye模型、Lorenz模型、Drude模型等,介电常数和磁导率可以表示为jω分式多项式的特点,结合频域到时域的转换关系（即用/t代替jω）和移位算子方法得到了修正点处的时域本构关系,进而获得时域递推计算式.数值结果表明,该方法具有通用、节省计算时间、节省内存和计算精度良好等优点. 关键词： 色散介质薄层 节点修正 移位算子 时域有限差分     

磁各向异性色散介质散射的Pad啨时域有限差分方法分析

根据矩阵Pad啨逼近理论,把磁化色散介质的相对磁导率张量表示成以jω为自变量的矩阵函数形式,用/t代替jω后过渡到时域,再引入离散时域移位算子代替时间微分算子.进而导出磁化色散介质中的磁感应强度B和磁场强度H在离散时域的色散关系,并将其具体应用于旋磁介质,得到了这种磁化色散介质的Pad啨时域有限差分方法的递推表达式.作为验证,用这种方法计算了磁化铁氧体球的后向雷达散射截面,所得结果与文献结果一致.理论推导及算例表明,该方法是正确和有效的.     

磁各向异性色散介质散射的Padé时域有限差分方法分析

根据矩阵Padé逼近理论，把磁化色散介质的相对磁导率张量表示成以jω为自变量的矩阵函数形式，用/t代替jω后过渡到时域，再引入离散时域移位算子代替时间微分算子.进而导出磁化色散介质中的磁感应强度B和磁场强度H在离散时域的色散关系，并将其具体应用于旋磁介质，得到了这种磁化色散介质的Padé时域有限差分方法的递推表达式.作为验证，用这种方法计算了磁化铁氧体球的后向雷达散射截面，所得结果与文献结果一致.理论推导及算例表明，该方法是正确和有效的. 关键词： 各向异性介质 色散介质 矩阵Padé逼近 时域有限差分方法     

介电球体和无限长圆柱体沉浸在另一介质中的极化模和色散关系

本文导出了介电球体和无限长圆柱体沉浸在另一极性或非极性介质中的极化本征模矢量和色散关系。计算表明,每一介面振动对应两个振动频率,它们分别局限在两种介质的TO频率(ωTo)和LO频率(ωLO)之间。但是,当介质球沉浸在另一非晶介质中时,每一个介面振动只对应一个振动频率,位于内球介质的ωTo和ωLO之间。我们还将Cds微晶无规地分布在另一非晶介质中的实验结果同计算值进行了比较,两者是基本符合的。     

色散介质电磁特性时域有限差分分析的Newmark方法

根据Debye模型、Drude模型和Lorentz模型3种常见色散介质模型频域极化率的特点,利用频域到时域的转换关系jω→?/?t,将极化矢量P与电场强度E的频域关系转换成时域内关于P的二阶微分方程,其对3种色散介质模型皆适用,具有统一的形式.然后采用相比于中心差分具有更高精度的Newmark两步算法(Newmark-β-γ法)求解该方程,进而得到E→P的递推公式,再结合本构关系得到D→E的时域递推式.实现了色散介质电磁场量的时域有限差分迭代计算.数值计算结果表明该方法是适用于3种色散介质模型的通用算法,并且相比于移位算子时域有限差分方法等以中心差分为基础的离散方案具有更高的计算精度.     

一种处理色散介质问题的通用时域有限差分方法

色散介质的介电系数是频率的函数,使本构关系在时域成为卷积关系.这就给用时域有限差分方法计算色散介质中波的散射和传播带来了困难.现有算法往往要针对不同色散介质模型推导相应的递推公式,算法的通用性较差.本文完善和发展了移位算子-时域有限差分方法,使之成为一种处理色散介质电磁问题的通用方法.首先,证明了常见的三种色散介质模型(德拜模型、洛伦兹模型和德鲁模型)的介电系数均可以写成适于移位算子法计算的有理分式函数形式.然后,用/t代替jω,过渡到时域,再引入时域移位算子z_t代替时间微分算子来处理有理分式函数形式的介电系数,给出离散时域本构关系的表示式,进而导出时域有限差分方法当中电位移矢量和电场强度之间的关系.最后,计算了几种色散介质的电磁散射,数值结果表明了本文方法和程序的通用性和正确有效性. 关键词： 时域有限差分方法 色散介质 移位算子     

有机材料在各种波长的折射率计算

使用Vogel方法和介质的色散关系,计算有机化合物在各波长的折射率n(λ).用Vogel方法把化合物分子分解成碎片,加和得到化合物在某几个波长的折射率.假设化合物在所计算的波长范围内为正常色散关系,使用Voyel方法得到的两个折射率,通过Cauchy公式得到色散曲线,从而得到各波长的折射率n(λ).通过对有机化合物芪盐和高分子材料聚苯乙烯的各波长折射率进行计算,并与材料在某一特定波长的实际折射率对比,证明了计算结果比较准确.并对双折射液晶材料的折射率计算进行了探索.     

一种适用于任意阶空间差分时域有限差分方法的色散介质通用吸收边界条件算法

基于伸展坐标完全匹配层方程和辅助微分方程 方法, 给出了一种在时域有限差分(FDTD) 计算中适用于常见色散介质的通用边界条件算法. 该算法适用于任意阶的FDTD空间差分, 并且由于所采用的D-H方程独立于计算区域, 所以可以被用于截断任意电介质. 数值试验结果表明, 与卷积完全匹配层 算法相比较, 所提出的吸收边界条件算法不仅通用性强、 计算复杂度低、 计算时间短, 并且吸收效果有明显的提高.     

宽带三倍频混频过程的群速匹配关系

基于非线性耦合波理论分析了三倍频过程中群速度之间的关系，推导出三波混频条件下的群速匹配关系式.对超短脉冲在Ⅱ类KDP晶体中的混频过程进行了模拟分析，结果表明：当三波(基波ω，谐波2ω和3ω)的群速度满足该匹配关系式时，三倍频(3ω)的带宽和转换效率均可达最大值；而随着三波群速度对该关系式的偏离逐渐增大，带宽变窄、转换效率下降.导出的三波群速匹配关系式对寻找合适的色散晶体和选择有效的匹配措施、实现宽带三次谐波转换具有指导意义.     

一种可用于分析任意介电常数张量的各向异性波导导模的二维时域有限差分模型

为分析具有任意介电常数张量的各向异性波导的导模,本文通过把简化的二维时域有限差分(2-DFDTD)法扩展至任意各向异性介质,提出了一种以D、E和H场为基础的统一的简化2-DFDTD模型.利用该模型,研究了简化的复数2-DFDTD方法与实变数2-DFDTD方法之间的关系.文中还讨论了复变数方法和实变数方法的激励技术.     

Theoretical analysis and simulation of conjugate heights for dual-conjugate AO system in lidar

A multi-conjugate adaptive optics （MCAO） can offer a possibility of widening field of view （FOV） characterized by the isoplanatic angle, and the choose of conjugate height becomes a basic problem for MCAO, which influences the size of iosplanatic angle. Considering the application of lidar, the isoplanatic angle＇s expressions of two deformable mirrors （DMs） MCAO for uplink and downlink are deduced. The effects of conjugate heights for dual-conjugate AO are thoughtfully discussed, and the isoplanatic angles are further analyzed. The results show that the isopanatic angle varies with the conjugate height and reaches the maximum as the conjugate height is at the optimal altitude. Moreover, the optimal conjugate height changes with the propagation distance.     

Quick single-photon detector with many avalanche photo diodes working on the time division

Due to the limit of response speed of the present single-photon detector, the code rate is still too low to come into practical use for the present quantum key distribution (QKD) system.A new idea is put up to design a quick single-photon detector.This quick single-photon detector is composed of a multi-port optic-fiber splitter and many avalanche photo diodes (APDs).Au of the ports with APDs work on the time division and cooperate with a logic discriminating and deciding unit driven by the clock signal.The operation frequency lies on the number N of ports, and can reach N times of the conventional single-photon detector.The single-photon prompt detection can come true for high repetition-rate pulses.The applying of this detector will largely raise the code rate of the QKD, and boost the commercial use.     

Characterization of PolyA and PolyC mismatches by Raman spectroscopy

A.C mismatches are studied by Raman spectral characterization of PolyA, PolyC, and their equimolar complex in solution of 0.14 mol/L Na ,pH7.0.Experimental results show that A·C mismatches occur to be A/B (mainly A) conformers, and unlike Watson-Crick base pairing, this kind of mismatches is stabilized by only one hydrogen bond involving cytosine N4H2 and adenine N7.The formation of A·C complex makes the base stacking interactions much stronger, and conformation of the backbone more ordered, which leads to obvious Raman hypochromic effect with some shifts in corresponding bands.     

Improvement of four-wave mixing-based wavelength conversion efficiency in dispersion shifted fiber by 40-GHz clock pumping

40-GHz clock modulated signal as a pump to improve the efficiency of four-wave mixing (FWM)-based wavelength conversion in a 26.5-km dispersion shifted fiber (DSF) is investigated. The experimental results demonstrate that the conjugated FWM component has higher intensity with the clock pumping than that with the continuous-wave (CW) light pumping. The improvement of FWM-based wavelength conversion efficiency is negligible when the pump power is less than Brillouin threshold. But when the pump power is greater than Brillouin threshold, the improvement becomes significant and increases with the increment of pump power. The improvement can increase up to 9 dB if pump power reaches 17 dBm.     

Multi-pulse operation of a Kerr-lens mode-locked femtosecond laser

Our experimental results show that the presence of a proper amount of negative group velocity dispersion is essential to multi-pulse operation of a Kerr-lens mode-locked femtosecond laser. We demonstrate that the pulse separations and the number of pulses contained within a cavity round trip are strongly dependent on the initial perturbations. The results allow us to get a better understanding on the influences of the convoluted self-phase modulation and intra-cavity dispersions on the stable multi-pulse oscillation in a Kerr-lens mode-locked femtosecond laser.     

Propagation of Helmholtz-Gauss beams in weak turbulent atmosphere

Based on the Rytov approximation of light propagation in weak turbulent atmosphere,the closed-form expressions of field and average irradiance of each one of the four fundamental families of Helmholtz-Gauss (HzG)beams:cosine-Gauss beams,stationary Mathieu-Gauss beams,stationary parabolic-Gauss beams,and Bessel-Gauss beams,which are propagating in weak turbulent atmosphere,are obtained.The results show that the field and average irradiance can be written as the product of four factors:complex amplitude depending on the z-coordinate only,a Ganssian beam.a factor of complex phase perturbation induced by atmospheric turbulence,and a complex scaled version of the transverse shape of the non-diffracting beam.The effect of weak atmospheric turbulence on irradiance distribution of the HzG beam can be ignored.     

An effective method for reducing speckle noise in digital holography

An effective method for reducing the speckle noise in digital holography is proposed in this paper.Different from the methods based on classical filtering technique,it utilizes the multiple holograms which are generated by rotating the illuminating light continuously.The intensity images reconstructed by a series of holograms generated by rotating the illuminating light possess different speckle patterns.Hence by properly averaging the reconstructed intensity fields,the speckle noises can be reduced greatly.Experimental results show that the proposed method is simple and effective to reduce speckle noise in digital holography.     

Influence of purity of HfO2 on reflectance of ultraviolet multilayer

The impurities in two kinds of HfO2 materials and in their corresponding single layer thin films were determined through glow discharge mass spectrum technology and secondary ion mass spectrometry(SIMS) equipment respectively.It was found that ZrO2 was the main impurity in the two kinds of HfO2 either in the original HfO2 materials or in the electron beam deposited films.In addition,the difference of Zr content in the two kinds of HfO2 single laver films was much laxger than that of the other impurities such as Ti and Fe.which showed that it was just ZrO2 that made the difference between the optical performance of the film products including the two kinds of HfO2.With these two kinds of HfO2 and the same kind of SiO2.we deposited HfO2/SiO2 multilayer reflective coatings at the wavelength of 266 nm.Experimental results showed that the reflectances of these two mirrors were about 99.85% and 99.15% respectively,which agreed well with the designed results what were based on the optical constants obtained from the corresponding single layer thin films.     

Speech enhancement by array crosstalk resistant ANC and spectrum subtraction

Microphone array-based speech enhancement has great importance for speech communications and speech recognition. To reduce the aperture of the microphone array and to increase the effect of the speech enhancement will greatly broaden the application areas of the microphone array. An array crosstalk resistant adaptive noise cancellation method is therefore presented. And then an improved spectral subtraction algorithm is further cascaded to obtain better enhancement results. Theoretic analysis and experiments indicate that the proposed scheme needs only a very small microphone array while it simultaneously achieves a higher SNR improvement. Besides, the proposed scheme can be used in many noisy environments and is easy for real-time implementation.