A retrieval method for nonlinear metamaterials

Metamaterials are engineered periodic structures for which it is possible to assign effective homogenized constitutive properties. In recent years, metamaterials in which the constituent elements are integrated with inherently nonlinear materials or electronic components have been considered for their potential impact on nonlinear wave propagation. As is the case with their linear counterparts, nonlinear metamaterials can also be assigned homogenized effective properties. The effective constitutive parameters of a metamaterial can be determined by a retrieval method applied to full-wave numerical simulations of a single layer of the structure. In this work, we present a transfer matrix approach that extends the retrieval of metamaterial properties to include the effective nonlinear susceptibilities. Comparisons with time-domain finite element simulations of continuous nonlinear slabs confirm the validity of this approach. The proposed approach is also applied to determine the nonlinear susceptibility of a simple nonlinear metamaterial.     

具有电与磁非线性效应超常介质中光束的聚焦特性

基于Drude模型,采用解析法与数值模拟法研究了具有电与磁非线性效应超常介质中光束的聚焦特性。结果表明超常材料中光束聚焦现象不但可以发生在自聚焦非线性介质中(正折射区)而且可以发生在自散焦非线性介质中(负折射区),且光束中心频率离介质的等离子频率越近,其聚焦效应越明显。特别是当电与磁非线性极化率为异号时,超常介质的非线性符号可以通过调控电与磁等离子频率的相对大小控制,这为主动操控光束的传输提供一种新的自由度。     

超常介质中超短电磁脉冲的传输特性研究

将超常介质的色散磁导率合并到非线性极化项中,借鉴常规介质中超短脉冲传输方程的推导方法,得到了非线性超常介质中超短脉冲的传输方程。在德鲁德(Drude)色散模型下,根据脉冲中心频率的不同在传输方程中出现了可正、可负、可为零的自陡峭系数,以及高阶非线性色散项。此外,利用矩方法对传输方程进行分析,得到了超常介质中超短脉冲传输方程的能量守恒定律表达式,揭示了色散磁导率导致的超短脉冲传输的新特性,发现二阶非线性色散使超短脉冲的能量、脉冲频移、脉冲宽度、中心位置和啁啾都随传输距离呈现振荡式变化。     

Direct observation of high-order optical susceptibilities via angularly-resolved multiwave mixing

A new method of monitoring specific high-order optical susceptibilities of nonlinear media via angularly-resolved emission in a phase-matched multiwave mixing configuration is presented. This technique, which can be viewed as the creation by two optical beams of a strongly anharmonic absorption or index grating, and the simultaneous exploration of higher Bragg-diffraction orders by a third wave, is demonstrated here in the case of gas media near optical resonance.     

饱和非线性正折射异向介质中的调制不稳定性

直接从异向介质中包含饱和非线性、自陡峭和二阶非线性色散效应的非线性扩展传输方程出发,采用线性稳定性分析法,导出了调制不稳定性的色散关系、不稳定条件、无量纲的临界扰动频率和增益谱。计算和讨论了异向介质正折射区无量纲的增益谱随归一化角频率和入射功率密度的变化关系。结果表明,在异向介质正折射区,随归一化角频率和入射功率密度的不同,增益谱将会出现扰动频率大于零、扰动频率大于某个非零临界值、扰动频率大于零而小于某个非零临界值3种形式。在归一化角频率较小时,调制不稳定性可能出现阈值入射功率密度;且在较小的入射功率密度时,调制不稳定性只能出现在大于某个临界频率时。     

超常介质中部分相干光束传输的调制不稳定性研究

超常介质与常规光学介质的一个最重要的区别是前者具有色散磁导率,其折射率呈现出可正可负,这样调制不稳定性出现新特点。借助于Wigner变换,从一般的非线性薛定谔方程出发,得到了部分相干光满足Wigner-Moyal的方程,获得了色散关系和增益谱,揭示了部分相干度不仅可抑制调制不稳定性,而且还可以缩小发生不稳定性的频率范围。     

Effective material parameter retrieval for thin sheets: Theory and application to graphene,thin silver films,and single-layer metamaterials

An important tool in the field of metamaterials is the extraction of effective material parameters from simulated or measured scattering parameters of a sample. Here we discuss a retrieval method for thin-film structures that can be approximated by a two-dimensional scattering sheet. We determine the effective sheet conductivity from the scattering parameters and we point out the importance of the magnetic sheet current to avoid an overdetermined inversion problem. Subsequently, we present two applications of the sheet retrieval method. First, we determine the effective sheet conductivity of thin silver films and we compare the resulting conductivities with the sheet conductivity of graphene. Second, we apply the method to a cut-wire metamaterial with an electric dipole resonance. The method is valid for thin-film structures such as two-dimensional metamaterials and frequency-selective surfaces and can be easily generalized for anisotropic or chiral media.     

受迫薛罗格双匣化学反应模型的随机共振

为了在理论上揭示高斯白噪声激励的薛罗格双匣化学反应模型对弱周期扰动的线性与非线性响应 ,分四态近似和两态近似两种情形 ,基于绝热近似与速率方程方法 ,解析导出线性的和非线性的敏感性以及信噪比的表达式 ,并与数值模拟结果进行比较 ,在一次谐波的意义上得到了解析结果与数值模拟结果的定量一致性 .理论上讲 ,该模型只能表现出奇次谐波的随机共振 ,但数值模拟结果也出现了二次谐波的随机共振 ,其原因可能归结为在数值模拟中有限频率的截断引入了误差 ,也可能归结为信号的高次谐波与背景噪声难以区分所致 .     

基于电磁超材料的两种等效参数提取算法的比较分析

基于传统反演算法和Kramers-Kronig关系改进算法, 分别提取弱耦合和强耦合超材料渔网结构模型的等效参数, 并对两种算法的有效性及普适性进行了探讨.理论分析及计算结果表明, 传统的反演算法可以准确地反演弱、强耦合情形下电磁超材料结构的等效参数, 但计算复杂度较高;而基于Kramers-Kronig关系的改进算法巧妙地降低了计算的复杂度, 能简单准确地提取弱耦合情形下电磁超材料结构的等效参数, 但对于强耦合情形则不适用, 原因在于强耦合情形破坏了Kramers-Kronig关系的解析且连续性要求. 研究结果拓展了等效媒质理论并可为新的电磁超材料的设计提供理论参考. 关键词： 超材料 反演算法 Kramers-Kronig关系 等效参数     

同时实现介电常数和磁导率为负的H型结构单元左手材料

基于用一种结构同时实现材料介电常数和磁导率为负的思想，提出了H型结构单元左手材料模型.采用矩形波导法测试表明，H型结构单元在微波频率范围出现左手透射峰，并且可以由参数t对左手特性区域进行调控.同时利用相位法、棱镜折射法和散射参量法从实验和理论证明了在左手透射峰区域材料的折射率为负值，介电常数和磁导率亦同时为负.相对由金属开口谐振环与金属线两种结构组合实现的左手材料，H型结构集磁谐振与电谐振于一体，结构简单、制备方便，对微波器件左手材料表现出更多的优越性.     

Numerical demonstration of three-dimensional terahertz metamaterials based on the causality principle

A parameter retrieval algorithm based on the causality principle and Kramers-Kronig （KK） relations is employed to calculate the effective parameters of three-dimensional （3D） metamaterials. Using KK relations, the branch selecting problem, which is the challenge of effective parameter retrieval method, can be removed. To reveal the validity of the proposed algorithm, the constitutive refractive index of a homogeneous polymide cube is extracted. The result is in excellent agreement with the intrinsic refractive index of the polymide. Finally, the two terahertz metamaterials with 3D structures are designed and their effective parameters are then retrieved using the proposed algorithm. Numerical simulations are performed using the fuiI-wave electromagnetic solver, CST Microwave Studio.     

Wide-band tuneability,nonlinear transmission,and dynamic multistability in SQUID metamaterials

Superconducting metamaterials comprising rf Superconducting QUantum Interference Devices (SQUIDs) have been recently realized and investigated with respect to their tuneability, permeability, and dynamic multistability properties. These properties are a consequence of intrinsic nonlinearities due to the sensitivity of the superconducting state to external stimuli. SQUIDs, made of a superconducting ring interrupted by a Josephson junction, possess yet another source of nonlinearity, which makes them widely tuneable with an applied dc dlux. A model SQUID metamaterial, based on electric equivalent circuits, is used in the weak coupling approximation to demonstrate the dc flux tuneability, dynamic multistability, and nonlinear transmission in SQUID metamaterials comprising non-hysteretic SQUIDs. The model equations reproduce the experimentally observed tuneability patterns and predict tuneability with the power of an applied ac magnetic field. Moreover, the results indicate the opening of nonlinear frequency bands for energy transmission through SQUID metamaterials, for sufficiently strong ac fields.     

White noise nonlinear system analysis in nuclear magnetic resonance spectroscopy

Nonlinear noise excitation in nuclear magnetic resonance is a form of nonlinear spectroscopy which exploits the nonlinear susceptibilities in a very direct way. The nonlinear susceptibilities are defined by perturbation theory in the frequency domain. In nonlinear system analysis, on the other hand, the system response is described by a Volterra series in the time domain. The kernels of the Volterra functionals carry the information about the system and are to be determined by experiment.The series expansion of a molecular, atomic or nuclear system response is derived in quantum mechanics by time dependent perturbation theory, leading to a Volterra series with time ordered, triangular kernels. The kernels are multi-dimensional products of decaying exponentials, which describe coherence decays of particular density matrix elements. The Fourier transforms of the triangular Volterra kernels are the susceptibilies, which are formally identical in NMR spectroscopy and nonlinear optical spectroscopy. The nonlinear susceptibilities are multi-dimensional spectra, which in NMR spectroscopy reveal the spin communication pathways. These are established by various forms of single quantum coherence connectivities, such as indirect coupling, chemical exchange, cross-relaxation, dipolar and quadrupolar coupling.If the functionals of the Volterra series are orthogonalized with respect to Gaussian white noise excitation, the Wiener series results. The Wiener kernels can be derived by multi-dimensional cross-correlation of the system response with different powers of the Gaussian white noise excitation.Cross-correlation of the transverse magnetization response to noise excitation in NMR leads to multi-dimensional time functions, the Fourier transforms of which closely resemble the nonlinear susceptibilities. The cross-correlation spectra differ from the susceptibilities in the governing Liouvillean and the dynamic density matrix, which are affected by saturation for continuous excitation. Cross-correlation spectra and susceptibilities converge for vanishing excitation power. Therefore the cross-correlation spectra are referred to as stochastic susceptibilities.In stochastic NMR spectroscopy only odd order susceptibilities exist for transverse magnetization. The first nonlinear order is the third, and the nonlinear spectral information is derived from the third order susceptibility. Higher order susceptibilities are not feasible to derive from experimental data. An important share of the nonlinear information is found on the six subdiagonal 2D cross-sections through the third order susceptibility. These cross-sections arise in three pairs, which carry distinct information, separated according to longitudinal magnetization and population effects, zero quantum coherences, and double quantum coherences.In practice a nonlinear 3D spectrum is computed from experimental data by an algorithm in the frequency domain, which yields access to selected regions in the 3D spectrum. This spectrum is the symmetrized stochastic third order susceptibility. All its sub-diagonal 2D cross-sections are equivalent. They are the average of the six different sub-diagonal 2D cross-sections through the asymmetric third order susceptibility.The stochastic excitation technique in NMR is characterized by several unique attributes. (1) There is no minimum time for a data acquisition cycle, so that, at the expense of signal-to-noise ratio, strong samples can be investigated faster with stochastic NMR than with pulsed FT NMR. (2) Stochastic excitation tests the sample extensively, and measures a maximum amount of information in a single experiment. This feature is of particular interest for investigation of short-lived samples and of samples with little a priori information. (3) An experiment with stochastic excitation is simple to perform, but the data processing is more complex than in FT spectroscopy. (4) The nonlinear information about spin communication pathways is derived for individual frequency regions only, which are identified in the stochastic ID spectrum. This information is located primarily on the sub-diagonal 2D cross-sections through the third order susceptibility. (5) Stochastic NMR spectra derived from random noise excitation are contaminated by systematic noise. In the sub-diagonal 2D cross-sections the noise is reduced by filtering and symmetrization during data processing. (6) Sub-diagonal 2D cross-sections are sensitive to experimental phase distortions in one direction only. They are readily adjusted in phase with the same parameters as the ID spectrum. (7) Stochastic multi-dimensional spectra can be computed at variable resolution from one and the same set of raw data.So far stochastic NMR spectroscopy is not applied routinely in analytical spectroscopy. More practical experience is needed to evaluate its merits in comparison with Fourier transform NMR.Stochastic excitation is distinguished from continuous wave and sparsely pulsed excitation by low input power in connection with large bandwidth. This important property cannot be exploited in high resolution NMR in liquids, because excitation power is not a restricting factor in this case. The situation is different in NMR imaging, where large field gradients require large bandwidths and the excitation power becomes a point of concern. For this reason stochastic RF excitation is being investigated in NMR imaging.The multi-dimensional cross-correlation functions obtained from random noise excitation generally are contaminated by systematic noise. The occurrence of systematic noise can be avoided if pseudo-random excitation is used in combination with a transformation of the system response to obtain the kernels. This technique is used successfully in Hadamard spectroscopy, where the linear Volterra kernel is the Hadamard transform of the linear response functional. Nonlinear transformations^(220,221) for retrieval of nonlinear kernels have not yet been realized in NMR spectroscopy.The cross-correlation technique underlying the data evaluation in stochastic nonlinear system analysis is equivalent to interferometry in optical spectroscopy. The Michelson interferometer is the most prominent optical correlator. The time resolution of the kernels derived by cross-correlation is determined by the inverse bandwidth of the excitation. With the Michelson interferometer a time resolution of 10⁻¹⁴ s is achieved in IR spectroscopy. Since the IR correlogramm is Fourier transformed for spectral analysis, the time resolution cannot be exploited otherwise. For analysis of fast time dependent processes a two-dimensional interferometer should be constructed, which performs a 2D cross-correlation of the system response to two in general different noise inputs. One input pumps the time dependent process, the other is used to investigate the time dependence spectroscopically. This technique is introduced by the name of ‘two-dimensional interferometry’. It uses low excitation power, but provides high time resolution at large response energy. Related work is pursued in nonlinear optical spectroscopy with incoherent excitation. In this area the use of broad band lasers is investigated for generation of echoes and for correlation based measurements of relaxation times.     

Linear phase retrieval for wave-front sensing

Phase retrieval from one or more intensity measurements is a potentially powerful and appealing technique for real-time adaptive-optics wave-front sensors. Under the assumption of small wave-front phase excursions, one is able to derive an exact solution to the inverse problem given three or more intensity measurements with known phase offsets. Applications include a high-order wave-front sensor to correct for residual aberrations in an adaptive-optics system in tandem with a low-resolution Hartmann-Shack wave-front sensor. The formula can also furnish mathematical insights into the full nonlinear phase-retrieval task.     

非线性左手材料中的二次谐波

基于电磁场理论，推导了无耗非线性左手材料中二次谐波的曼利－诺关系，及相位匹配条件下正向基波与逆向二次谐波的能量转换过程及其空间分布.验证了可将无耗非线性左手材料的入射面作为反射镜，把能量以二次谐波的形式反射.同时分别给出有限厚度介质板中基波和二次谐波的场强分布数值结果，验证了结论的正确性.最后从相位失配角度说明了相位匹配是分析非线性左手材料二次谐波的重要条件.这为研究左手材料的非线性理论奠定了基础.     

Nonlinear phenomena of left-handed nonlinear split-ring resonators

We theoretically analyze a two-dimensional periodic structure consisting of period nonlinear split-ring resonators (SRRs) with varactor diode. The diode is loaded into the slit of the SRRs. Then, we demonstrate nonlinear phenomena of left-handed nonlinear SRRs. This paper introduces nonlinear SRRs based on left-handed (LH) media and simulates self-tuning mechanisms, bistable effects in the microwave frequency range. In addition, with the increase of input power, the nonlinear response of the SRRs becomes multi-valued, paving a way for creating bistable tunable metamaterials.     

人工表面等离激元超材料

以微带为代表的传统微波传输线无法精细操控电磁模式,因此传统电子信息系统在空间耦合、动态响应和性能鲁棒性等方面存在瓶颈。人工表面等离激元(SSPP)超材料可打破上述瓶颈,是光学与信息领域的研究热点之一。人工表面等离激元超材料是一类模拟光频段表面等离激元特性的新型超材料,可在微波和太赫兹频段精细操控表面波,具有与平面电路相似的构型特性,可用于制备下一代集成电路的基础传输线。人工表面等离激元分为传输型和局域型两类。传输型人工表面等离激元超材料始于三维立体结构,后发展成超薄梳状金属条带构型。学者们构建了以其为基础的微波电路新体系,研制了人工表面等离激元滤波器、天线、放大器和倍频器等典型的无源和有源器件,并将其集成为可实现亚波长间距多通道信号非视距传输的无线通信系统。人工局域表面等离激元(SLSP)超材料也经历了从三维立体构型到超薄构型的发展历程,并通过螺旋构型、链式构型、高阶模式和杂化模式等为电磁波的亚波长尺度调控提供了更多自由度。系统讨论了人工表面等离激元超材料在微波电路中的相关理论和应用,包括人工表面等离激元超材料的基本概念、构型发展、无源/有源器件以及无线通信系统。     

Subpicosecond pulse compression in nonlinear photonic crystal waveguides based on the formation of high-order optical solitons

We investigate by numerical simulation the compression of subpicosecond pulses in two-dimensional nonlinear photonic crystal (PC) waveguides. The compression originates from the generation of high-order optical solitons through the interplay of the huge group-velocity dispersion and the enhanced self-phase modulation in nonlinear PC waveguides.Both the formation of Bragg grating solitons and gap solitons can lead to efficient pulse compression. The compression factors under different excitation power densities and the optimum length for subpicosecond pulse compression have been determined. As a compressor, the total length of the nonlinear PC waveguide is only ten micrometres and therefore can be easily incorporated into PC integrated circuits.