Focusing of stochastic electromagnetic Gaussian Schell-model beams through a high numerical aperture objective

Based on vectorial Debye diffraction theory, the focusing properties of stochastic electromagnetic Gaussian Schell-model beams in the focal region of high numerical aperture objective are investigated. Expressions for the intensity distribution and the degree of polarization are derived near the focus. Numerical calculations are performed to analyze the influences of varying corresponding parameters on the changes in the intensity distribution and in the degree of polarization in the focal region.     

静态电磁场边值问题计算方法

介绍常用的电磁场分析方法——解析法、有限差分法、有限元法.以静态电磁场边值问题为例,分别用自己编写的MATLAB有限差分程序及大型有限元分析软件ANSYS求出了数值解,并与解析法得到的精确解进行比较,得出了用数值法求解电磁场问题基本满足工程需要的结论.采用数值计算时,对如何减小计算机内存,取消冗余数据,也做了进一步的讨论.     

Enhanced transmission of electromagnetic waves through metamaterials

This paper reviews our recent experimental and simulation results regarding the electromagnetic wave transmission through three configurations of sandwiching structures of metamaterials: a metallic mesh sandwiched between two identical layers composed of split rings, metallic fractals, and fractal slits, respectively. We observed the enhanced transmission of the waves through these three types of sandwiching composites with respect to the opaque metallic mesh. The locations of the transmission peaks in the spectrum are associated closely with the band characteristics of the sandwiching layer by appearing either on the left- or the right-hand side of its band. PACS 41.20.Jb; 42.70.Qs; 78.20.-e; 42.25.Bs; 47.53.+n     

Invariants and numerical methods for ODEs

The preservation of some invariants is important during the numerical integration of ODEs. In some cases, failure to maintain certain invariants leads to physically impossible solutions, in other cases to instability. Some authors report that the long term macroscopic characteristics of a solution are better represented if certain invariants are maintained. This paper considers two types of invariants, point invariants which are functions along a particular trajectory, and differential invariants which are invariant relations between neighboring trajectories, such as symplectic invariants. Conventional numerical methods usually introduce (h^p+1) errors in an invariant at each numerical step and these will accumulate over the integration interval. Two approaches to preventing this accumulation are (1) to find methods that maintain the invariants within round-off error, or (2) to find methods that maintain another invariant which is no more than (h^q) different from the invariant satisfied by the ODE. These approaches will be illustrated with some examples. The major objective of this paper is to consider the potential of general methods (that is, methods that do not depend on the specific differential equation) to maintain the invariants.     

Studies of scattering theory using numerical methods

Abstract

Numerical simulations, using both exact and approximate methods, are used to study rough surface scattering in both the smd and large roughness regimes. This study is limited lo scattcring lrom rough one-dimensional surfaces that obey the Dirichlet boundary condition and have a Gaussian roughness spectrum. For surfdces with small roughness (kh?1, where k is the radiation wavenumber and h is the root-mean-square (RMS) Surface height), perturbation theory is known to be valid. However, it is shown numerically that when kh?1 and kl?6 (where I is the surface correlation length) the Kirchhoffapprorimation is valid except at low grazing angles, and one must sum the first three orders of perturbation theory obtain the correct result. For kh?1 and kl?1, first-order perturbation theory is accurate. In this region, the accuracy of the first two terms of the iterative series solution of the exact integral equation is examined; the first term a1 this series is the Kirchhoff approximation, It is shown numerically that lor very small kh these first two terms reduce to first-order perturbation theory. However, lor this reduction to occur, kh must be made smaller than necessdry lor first-order perturbation theory to be accurate. In the regime of large roughness (kh?1) backscattering enhancement occurs when the RMS slope is on the order of unity. Several investigators have recently shown that the second term of the iterative series solution (the double-scattering term) replicates the properties of backscattering enhancement reasonably well. However, the double-scattering term has a lundamental flaw: predictions lor the scattering cross section per unit length based on the double-scattering term increase as the surfdce length is increased. This is shown here with numerical simulations and with an approximate analytical result based on the high frequency limit. The physical significance of this finding is also discussed. The final topic is the use of the double-scattering approximation to study the mechanism for backscattering enhancement with the Dirichlet boundary condition. This mechanism is usually assumed to be interference between reciprocal scattering paths. When the interlerence between reciprocal scattering paths is removed, the enhancement is eliminated. This shows that interference between reciprocal paths is almost certainly the dominant mechanism for backscattering enhancement in the scattering regime studied.     

Modulation in the statistical properties of stochastic electromagnetic beams through an electromagnetic induced transparency atomic vapor

We report analytical expressions for the elements of the 2 × 2 cross-spectral density matrix of a stochastic electromagnetic beam passing through an electromagnetic induced transparency (EIT) atomic vapor. By use of the derived formulas the changes in the spectral density, the spectral degree of coherence, and the spectral degree of polarization of such a beam on propagation can be studied in detail. Numerical examples show that the statistical properties of the stochastic electromagnetic beam can be modulated by the Rabi frequency of the control light when the beam propagates through the EIT atomic vapor.     

Fast quantum algorithm for numerical gradient estimation

Given a black box for f, a smooth real scalar function of d real variables, one wants to estimate [symbol: see text]f at a given point with n bits of precision. On a classical computer this requires a minimum of d + 1 black box queries, whereas on a quantum computer it requires only one query regardless of d. The number of bits of precision to which f must be evaluated matches the classical requirement in the limit of large n.     

A fast directional algorithm for high-frequency electromagnetic scattering

This paper is concerned with the fast solution of high-frequency electromagnetic scattering problems using the boundary integral formulation. We extend the O(N log N) directional multilevel algorithm previously proposed for the acoustic scattering case to the vector electromagnetic case. We also detail how to incorporate the curl operator of the magnetic field integral equation into the algorithm. When combined with a standard iterative method, this results in an almost linear complexity solver for the combined field integral equations. In addition, the butterfly algorithm is utilized to compute the far field pattern and radar cross section with O(N log N) complexity.     

New hybrid FDTD algorithm for electromagnetic problem analysis

Since the time step of the traditional finite-difference time-domain(FDTD) method is limited by the small grid size, it is inefficient when dealing with the electromagnetic problems of multi-scale structures.Therefore, the explicit and unconditionally stable FDTD(US-FDTD) approach has been developed to break through the limitation of Courant–Friedrich–Levy(CFL) condition.However, the eigenvalues and eigenvectors of the system matrix must be calculated before the time iteration in the explicit US-FDTD.Moreover, the eigenvalue decomposition is also time consuming, especially for complex electromagnetic problems in practical application.In addition, compared with the traditional FDTD method, the explicit US-FDTD method is more difficult to introduce the absorbing boundary and plane wave.To solve the drawbacks of the traditional FDTD and the explicit US-FDTD, a new hybrid FDTD algorithm is proposed in this paper.This combines the explicit US-FDTD with the traditional FDTD, which not only overcomes the limitation of CFL condition but also reduces the system matrix dimension, and introduces the plane wave and the perfectly matched layer(PML) absorption boundary conveniently.With the hybrid algorithm, the calculation of the eigenvalues is only required in the fine mesh region and adjacent coarse mesh region.Therefore, the calculation efficiency is greatly enhanced.Furthermore, the plane wave and the absorption boundary introduction of the traditional FDTD method can be directly utilized.Numerical results demonstrate the effectiveness, accuracy, stability, and convenience of this hybrid algorithm.     

Propagation of electromagnetic dipole waves through dielectric interfaces

The problem of divergent electromagnetic dipole waves propagating through parallel dielectric interfaces is solved. The solution is obtained in an analytic form that can be readily evaluated numerically. The result is obtained as a solution to a boundary-value problem. Applications of the solution are described.     

海洋可控源三维电磁响应显式灵敏度矩阵的快速算法

借助电场耦合势三维有限体积法与直接求解技术,研究建立了一套海洋可控源三维电磁响应显式灵敏度矩阵(或称为Fréchet导数)高效算法.首先,利用Yee氏交错网格和有限体积法对电场混合势Helmholtz方程进行离散处理,建立与移动源电磁场正演模拟相对应的大型代数方程组,再应用直接法得到的逆矩阵和三维线性插值技术事先确定插值算子和投影算子,并利用投影算子与各个发射源离散向量的乘积计算多发射源电磁响应,极大地提高了多发射源电磁场正演模拟效率.在此基础上,根据块状模型和像素模型中异常体电导率分片常数分布特征,将电导率摄动产生的一次散射电流场表示成Yee氏剖分网格上散射电流元的叠加,由投影算子与散射电流元的离散向量的乘积快速计算出电场强度与磁场强度的显式灵敏度矩阵.最后,通过数值计算检验算法的有效性,并通过块状模型与像素模型分别研究海洋可控源电磁响应特征.     

Transmission of electromagnetic waves through thin Cu films

Optical transmission is used to study the electronic structure of nanosized Cu phases that are obtained on the surface of single-crystalline NaCl (100) using vacuum deposition. It is demonstrated that Cu nanofilms with a thickness of 1–1.5 monolayers exhibit properties typical of narrow-band-gap semiconductors.     

All-optical switching of the transmission of electromagnetic radiation through subwavelength apertures

Unprecedented optical control of the surface plasmon polariton assisted transmission of terahertz radiation through subwavelength apertures is rendered possible by carrier-induced changes to the dielectric properties of a semiconductor grating. Although the study presented is static, the extension of our approach to dynamic switching and tuning is deemed straightforward, opening the way for the realization of ultrafast surface plasmon based devices.     

Propagation of stochastic electromagnetic vortex beams through the turbulent biological tissues

The general analytical expression of the stochastic electromagnetic vortex beams through turbulent biological tissues is derived based on the fractal model. The statistical properties, including the spectral density, the spectral degree of coherence and the spectral degree of polarization are investigated in detail. It can be found that the normalized spectral density of the stochastic electromagnetic vortex beams with higher topological charge is less influenced by turbulence than that with lower topological charge. In addition, the change of the degree of polarization versus propagation distance of the anisotropic vortex beams in biological tissues differs from that of the isotropic vortex beams. The findings might be useful in the investigation of the structures of biological tissues and operation of communication and sensing systems involving biological tissues turbulence channels.     

Unwinding algorithm for numerical generation and writing of fullerenes

An efficient algorithm is proposed for the numerical generation of fullerenes of an arbitrary structure. The algorithm is a combination of the method for unwinding fullerenes into a triangular network and the topological-invariants method for describing quasi-two-dimensional closed clusters. Graphs of possible structures of C_n fullerenes are found satisfying the isolated-pentagon rule for numbers of atoms in the range 100<n≤150.