Wei E. I. Sha

Yongpin P. Chen, Wei E. I. Sha, Li Jun Jiang et al.

A novel unified Hamiltonian approach is proposed to solve Maxwell-Schrodinger equation for modeling the interaction between classical electromagnetic (EM) fields and particles. Based on the Hamiltonian of electromagnetics and quantum mechanics, a unified Maxwell-Schrodinger system is derived by the variational principle. The coupled system is well-posed and symplectic, which ensures energy conserving property during the time evolution. However, due to the disparity of wavelengths of EM waves and that of electron waves, a numerical implementation of the finite-difference time-domain (FDTD) method to the multiscale coupled system is extremely challenging. To overcome this difficulty, a reduced eigenmode expansion technique is first applied to represent the wave function of the particle. Then, a set of ordinary differential equations (ODEs) governing the time evolution of the slowly-varying expansion coefficients are derived to replace the original Schrodinger equation. Finally, Maxwell's equations represented by the vector potential with a Coulomb gauge, together with the ODEs, are solved self-consistently. For numerical examples, the interaction between EM fields and a particle is investigated for both the closed, open and inhomogeneous electromagnetic systems. The proposed approach not only captures the Rabi oscillation phenomenon in the closed cavity but also captures the effects of radiative decay and shift in the open free space. After comparing with the existing theoretical approximate models, it is found that the approximate models break down in certain cases where a rigorous self-consistent approach is needed. This work is helpful for the EM simulation of emerging nanodevices or next-generation quantum electrodynamic systems.

Ming Fang, Zhixiang Huang, Wei E. I. Sha et al.

The interaction between the electromagnetic field and plasmonic nanostructures leads to both the strong linear response and inherent nonlinear behavior. In this paper, a time-domain hydrodynamic model for describing the motion of electrons in plasmonic nanostructures is presented, in which both surface and bulk contributions of nonlinearity are considered. A coupled Maxwell-hydrodynamic system capturing full-wave physics and free electron dynamics is numerically solved with the parallel finite-difference time-domain (FDTD) method. The validation of the proposed method is presented to simulate linear and nonlinear responses from a plasmonic metasurface. The linear response is compared with the Drude dispersion model and the nonlinear terahertz emission from a difference-frequency generation process is validated with theoretical analyses. The proposed scheme is fundamentally important to design nonlinear plasmonic nanodevices, especially for efficient and broadband THz emitters.

Jing Shen, Wei E. I. Sha, Xiaojing Kuang et al.

A symplectic pseudospectral time-domain (SPSTD) scheme is developed to solve Schrodinger equation. Instead of spatial finite differences in conventional finite-difference time-domain (FDTD) method, the fast Fourier transform is used to calculate the spatial derivatives. In time domain, the scheme adopts high-order symplectic integrators to simulate time evolution of Schrodinger equation. A detailed numerical study on the eigenvalue problems of 1D quantum well and 3D harmonic oscillator is carried out. The simulation results strongly confirm the advantages of the SPSTD scheme over the traditional PSTD method and FDTD approach. Furthermore, by comparing to the traditional PSTD method and the non-symplectic Runge-Kutta (RK) method, the explicit SPSTD scheme which is an infinite order of accuracy in space domain and energy-conserving in time domain, is well suited for a long-term simulation.

Zi He, Ji Hong Gu, Wei E. I. Sha et al.

Metallic nanoparticles (NPs) support localized surface plasmon resonances (LSPRs), which enable to concentrate sunlight at the active layer of solar cells. However, full-wave modeling of the plasmonic solar cells faces great challenges in terms of huge computational workload and bad matrix condition. It is tremendously difficult to accurately and efficiently simulate near-field multiple scattering effects from plasmonic NPs embedded into solar cells. In this work, a preconditioned volume integral equation (VIE) is proposed to model plasmonic organic solar cells (OSCs). The diagonal block preconditioner is applied to different material domains of the device structure. As a result, better convergence and higher computing efficiency are achieved. Moreover, the calculation is further accelerated by two-dimensional periodic Green's functions. Using the proposed method, the dependences of optical absorption on the wavelengths and incident angles are investigated. Angular responses of the plasmonic OSCs show the super-Lambertian absorption on the plasmon resonance but near-Lambertian absorption off the plasmon resonance. The volumetric method of moments and explored physical understanding are of great help to investigate the optical responses of OSCs.

Qi I. Dai, Jun Wei Wu, Ling Ling Meng et al.

Characteristic mode (CM) analysis poses challenges in computational electromagnetics (CEM) as it calls for efficient solutions of dense generalized eigenvalue problems (GEP). Multilevel fast multipole algorithm (MLFMA) can greatly reduce the computational complexity and memory cost for matrix-vector product operations, which is powerful in iteratively solving large scattering problems. In this article, we demonstrate that MLFMA can be easily incorporated into the implicit restarted Arnoldi (IRA) method for the calculation of CMs, where MLFMA with the sparse approximate inverse (SAI) preconditioning technique is employed to accelerate the construction of Arnoldi vectors. This work paves the way of CM analysis for large-scale and complicated three-dimensional ($3$-D) objects with limited computational resources.

Xiao Z. Wang, Wei E. I. Sha

This report mainly focused on the basic concepts and the recovery methods for the random sampling. The recovery methods involve the orthogonal matching pursuit algorithm and the gradient-based total variation strategy. In particular, a fast and efficient observation matrix filling technique was implemented by the classic Shannon interpolation and Poisson summation formulae. The numerical results for the trigonometric signal, the Gaussian-modulated sinusoidal pulse, and the square wave were demonstrated and discussed. The work may give some help for future work in theoretical study and practical implementation of the random sampling.