Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open...Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain（ECT-CFDTD） method is used to simulate the wave propagation inside a perfect electric conductor（PEC） waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.展开更多
This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the...This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.展开更多
The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the ...The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.展开更多
An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D tra...An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.展开更多
A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spheric...A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.展开更多
The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order sy...The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.展开更多
To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE met...To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.展开更多
Composite electromagnetic scattering from a two-dimensional （2D） ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain （FDTD） method. A uniaxial perfectly mat...Composite electromagnetic scattering from a two-dimensional （2D） ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain （FDTD） method. A uniaxial perfectly matched layer is adopted for truncation of FDTD lattices.The FDTD updated equations can be used for the total computation domain by choosing the uniaxial parameters properly. To validate the proposed numerical technique,a 2D infinitely long cylinder over the sea surface is taken into account first.The variation of angular distribution of the scattering changing with incident angle is calculated. The results show good agreement with the conventional moment method. Finally,the influence of the incident angle,the polarization,and the size of the ship-like target on the composite scattering coefficient is discussed in detail.展开更多
In this paper, a modified sub-gridding scheme that hybridizes the conventional finite-difference time-domain（FDTD）method and the unconditionally stable locally one-dimensional（LOD） FDTD is developed for analyzing ...In this paper, a modified sub-gridding scheme that hybridizes the conventional finite-difference time-domain（FDTD）method and the unconditionally stable locally one-dimensional（LOD） FDTD is developed for analyzing the periodic metallic nanoparticle arrays. The dispersion of the metal, caused by the evanescent wave propagating along the metal-dielectric interface, is expressed by the Drude model and solved with a generalized auxiliary differential equation（ADE） technique.In the sub-gridding scheme, the ADE–FDTD is applied to the global coarse grids while the ADE–LOD–FDTD is applied to the local fine grids. The time step sizes in the fine-grid region and coarse-grid region can be synchronized, and thus obviating the temporal interpolation of the fields in the time-marching process. Numerical examples about extraordinary optical transmission through the periodic metallic nanoparticle array are provided to show the accuracy and efficiency of the proposed method.展开更多
The frequency dependent permittivity for dusty plasmas is provided by introducing the charging response factor and charge relaxation rate of airborne particles. The field equations that describe the characteristics of...The frequency dependent permittivity for dusty plasmas is provided by introducing the charging response factor and charge relaxation rate of airborne particles. The field equations that describe the characteristics of Terahertz（THz） waves propagation in a dusty plasma sheath are derived and discretized on the basis of the auxiliary differential equation（ADE） in the finite difference time domain（FDTD） method. Compared with numerical solutions in reference, the accuracy for the ADE FDTD method is validated. The reflection property of the metal Aluminum interlayer of the sheath at THz frequencies is discussed. The effects of the thickness, effective collision frequency, airborne particle density, and charge relaxation rate of airborne particles on the electromagnetic properties of Terahertz waves through a dusty plasma slab are investigated.Finally, some potential applications for Terahertz waves in information and communication are analyzed.展开更多
The piecewise linear recursive convolution （PLRC） finite-different time-domain （FDTD） method improves accuracy over the original recursive convolution （RC） FDTD approach and current density convolution （JEC） b...The piecewise linear recursive convolution （PLRC） finite-different time-domain （FDTD） method improves accuracy over the original recursive convolution （RC） FDTD approach and current density convolution （JEC） but retains their advantages in speed and efficiency. This paper describes a revised piecewise linear recursive convolution PLRC-FDTD formulation for magnetized plasma which incorporates both anisotropy and frequency dispersion at the same time, enabling the transient analysis of magnetized plasma media. The technique is illustrated by numerical simulations of the reflection and transmission coefficients through a magnetized plasma layer. The results show that the revised PLRC-FDTD method has improved the accuracy over the original RC FDTD method and JEC FDTD method.展开更多
In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is i...In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint.展开更多
The surface plasmon resonance( SPR) of gold( Au) nanorod can be tuned in a large visible-near infrared( Vis-NIR) region by changing the aspect ratio of nanorod. Compared with the SPR of isolated Au nanorod,assembly of...The surface plasmon resonance( SPR) of gold( Au) nanorod can be tuned in a large visible-near infrared( Vis-NIR) region by changing the aspect ratio of nanorod. Compared with the SPR of isolated Au nanorod,assembly of Au nanorods exhibits strong coupling effect in the nanogap and wealthy changes in the optical spectra. The SPR coupling effects and localized electronic fields for end-to-end( E-E) and side-by-side( S-S) assembled Au nanorods dimers are studied through finite-difference time-domain( FDTD) simulation. With decreasing the gap spacing,the longitudinal SPR( SPRL) red-shifts for the E-E orients Au nanorods dimer and blue-shifts for S-S orientes Au nanorods dimer. The transverse SPR( SPRT) has slight red-shifting for S-S assembly and no shifting for the E-E assembly. Moreover,a new coupling SPR appears for the E-E assembly in a long wavelength in the NIR region,blue-shifting and enhancing with decreasing the gap spacing. Based on the spring oscillator model and the polarization of nanoparticles with incident electric field,the SPR shifting and the appearance of new coupling SPR of assembled Au nanorods are proposed.展开更多
The analysis of electromagnetic propagation in a dispersive medium is complicated in the time-domain because its dielectric constant is frequency-dependent. In this paper, the dielectric constant of the dispersive med...The analysis of electromagnetic propagation in a dispersive medium is complicated in the time-domain because its dielectric constant is frequency-dependent. In this paper, the dielectric constant of the dispersive medium is written as a rational polynomial function, and the relationship between D and E is derived in the time-domain. It is referred to as the shift operator finite-different time-domain （SO-FDTD） method. Compared to an analytical solution and a piecewise linear current density recursive convolution （PLJERC） solution, the high accuracy and efl%iency of this method is verified by calculating the reflectance of the electromagnetic wave through a cold plasma slab. As the electron density in plasma is distributed as the Epstein formula, the effect of distribution grads and electron collision frequency on the reflectance is calculated by using the SO-FDTD method. The result shows that the increase in the distribution grads coefficient affects the reflectance sharply. When it comes to a smaller distribution grads coelBcient, the increase of the collision frequency showed a significant effect on the reflectance, but on the contrary, there is actually less and less effect till it disappears.展开更多
Based on FDTD difference expressions and eigenfunctions of Maxwell functions in cylindrical coordinates, mesh wave impedances (MWIs) in 2D and 3D cylindrical coordinates were introduced. Combined with the concept of p...Based on FDTD difference expressions and eigenfunctions of Maxwell functions in cylindrical coordinates, mesh wave impedances (MWIs) in 2D and 3D cylindrical coordinates were introduced. Combined with the concept of perfectly matched layer (PML), MWI PML absorbing boundary condition (ABC) algorithm was deduced in 2D cylindrical coordinates. Numerical experiments were done to investigate the validity of MWI and its application in cylindrical coordinates FDTD algorithm. The results showed that MWI in cylindrical coordinates can be used to accurately calculate the numerical reflection error caused by different mesh increments in non uniform FDTD. MWI can also provide theoretical criterion to define the permitted variable range of mesh dimension. MWI PML ABC is easy to be applied and reduces low numerical reflection, which only causes a little higher reflection error compared with Teixeira's PML.展开更多
A novel method based on the relevance vector machine（RVM） for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.T...A novel method based on the relevance vector machine（RVM） for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.The nonlinearity is embodied in the relation between the scattered field and the target property,which can be obtained through the RVM training process.Besides,rather than utilizing regularization,the ill-posed nature of the inversion is naturally accounted for because the RVM can produce a probabilistic output.Simulation results reveal that the proposed RVM-based approach can provide comparative performances in terms of accuracy,convergence,robustness,generalization,and improved performance in terms of sparse property in comparison with the support vector machine（SVM） based approach.展开更多
The propagation characteristics of oblique incidence terahertz(THz) waves through non-uniform plasma are investigated by the shift-operator finite-difference time-domain(SO-FDTD) method combined with the phase matchin...The propagation characteristics of oblique incidence terahertz(THz) waves through non-uniform plasma are investigated by the shift-operator finite-difference time-domain(SO-FDTD) method combined with the phase matching condition.The electron density distribution of the non-uniform plasma is assumed to be in a Gaussian profile. Validation of the present method is performed by comparing the results with those obtained by an analytical method for a homogeneous plasma slab.Then the effects of parameters of THz wave and plasma layer on the propagation properties are analyzed. It is found that the transmission coefficients greatly depend on the incident angle as well as on the thickness of the plasma, while the polarization of the incident wave has little influence on the propagation process in the range of frequency considered in this paper. The results confirm that the THz wave can pass through the plasma sheath effectively under certain conditions,which makes it a potential candidate to overcome the ionization blackout problem.展开更多
In this study, the propagation of electromagnetic waves in one-dimensional plasma photonic crystals （PPCs）, namely, superlattice structures consisting alternately of a homogeneous unmagnetized plasma and dielectric ...In this study, the propagation of electromagnetic waves in one-dimensional plasma photonic crystals （PPCs）, namely, superlattice structures consisting alternately of a homogeneous unmagnetized plasma and dielectric material, is simulated numerically using the finite-difference time-domain （FDTD） algorithm. A perfectly matched layer （PML） absorbing technique is used in this simulation. The reflection and transmission coefficients of electromagnetic （EM） waves through PPCs are calculated. The characteristics of the photonic band gap （PBG） are discussed in terms of plasma density, dielectric constant ratios, number of periods, and introduced layer defect. These may provide some useful information for designing plasma photonic crystal devices.展开更多
In this paper, a subwavelength metal-grating assisted sensor of Kretschmann style that is capable of detecting the sample with a refractive index higher than that of the substrate is proposed. The sensor configuration...In this paper, a subwavelength metal-grating assisted sensor of Kretschmann style that is capable of detecting the sample with a refractive index higher than that of the substrate is proposed. The sensor configuration is similar to the traditional Kretschmann structure, but the metal film is pattered into a grating. As a TM-polarized laser beam impinges from the substrate, a resonant dip point in reflectance curve is produced at a certain incident angle. Our studies indicate that the sensing sensitivity and resolution are affected by the grating＇s gap and period, and after these parameters havelaeen optimized, a sensing sensitivity of 51.484°/RIU is obtained with a slightly changing resolution.展开更多
Three-dimensional photonic crystal （PC） heterostructures with high quality are fabricated by using a pressure controlled isothermal heating vertical deposition technique. The formed heterostructures have higher qual...Three-dimensional photonic crystal （PC） heterostructures with high quality are fabricated by using a pressure controlled isothermal heating vertical deposition technique. The formed heterostructures have higher quality, such as deeper band gaps and sharper band edges, than the heterostructures reported so far. Such a significant improvement in quality is due to the introduction of a thin TiO2 buffer layer between the two constitutional PCs. It is revealed that the disorder caused by lattice mismatch is successfully removed if the buffer layer is used once. As a result, the formed heterostructures possess the main features in the band gap of constitutional PCs. The crucial role of the thin buffer layer is also verified by numerical simulations based on the finite-difference time-domain technique.展开更多
基金supported by the National Natural Science Foundation of China(Grant No.61231003)
文摘Based on conformal construction of physical model in a three-dimensional Cartesian grid,an integral-based conformal convolutional perfectly matched layer（CPML） is given for solving the truncation problem of the open port when the enlarged cell technique conformal finite-difference time-domain（ECT-CFDTD） method is used to simulate the wave propagation inside a perfect electric conductor（PEC） waveguide.The algorithm has the same numerical stability as the ECT-CFDTD method.For the long-time propagation problems of an evanescent wave in a waveguide,several numerical simulations are performed to analyze the reflection error by sweeping the constitutive parameters of the integral-based conformal CPML.Our numerical results show that the integral-based conformal CPML can be used to efficiently truncate the open port of the waveguide.
基金Project supported by Tianjin Research Program Application Foundation and Advanced Technology,China(Grant No.15JCQNJC01100)
文摘This paper investigates the phenomenon of three-pulse photon echo in thick rare-earth ions doped crystal whose thickness is far larger than 0.002 cm which is adopted in previous works.The influence of thickness on the three-pulse photon echo＇s amplitude and efficiency is analyzed with the Maxwell-Bloch equations solved by finite-difference timedomain method.We demonstrate that the amplitude of three-pulse echo will increase with the increasing of thickness and the optimum thickness to generate three-pulse photon echo is 0.3 cm for Tm^（3＋）：YAG when the attenuation of the input pulse is taken into account.Meanwhile,we find the expression 0.09 exp（α＇L）,which is previously employed to describe the relationship between echo＇s efficiency and thickness,should be modified as 1.3 · 0.09 exp（2.4 ·α＇L） with the propagation of echo considered.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11304074,61475042,and 11274088)the Natural Science Foundation of Hebei Province,China(Grant Nos.A2015202320 and GCC2014048)the Key Subject Construction Project of Hebei Province University,China
文摘The finite-difference time-domain method is used to simulate the optical characteristics of an in-plane switching blue phase liquid crystal display.Compared with the matrix optic methods and the refractive method,the finite-difference timedomain method,which is used to directly solve Maxwell＇s equations,can consider the lateral variation of the refractive index and obtain an accurate convergence effect.The simulation results show that e-rays and o-rays bend in different directions when the in-plane switching blue phase liquid crystal display is driven by the operating voltage.The finitedifference time-domain method should be used when the distribution of the liquid crystal in the liquid crystal display has a large lateral change.
基金supported by the National Natural Science Foundation of China(Grant Nos.61331007 and 61471105)
文摘An efficient conformal locally one-dimensional finite-difference time-domain（LOD-CFDTD） method is presented for solving two-dimensional（2D） electromagnetic（EM） scattering problems. The formulation for the 2D transverse-electric（TE） case is presented and its stability property and numerical dispersion relationship are theoretically investigated. It is shown that the introduction of irregular grids will not damage the numerical stability. Instead of the staircasing approximation, the conformal scheme is only employed to model the curve boundaries, whereas the standard Yee grids are used for the remaining regions. As the irregular grids account for a very small percentage of the total space grids, the conformal scheme has little effect on the numerical dispersion. Moreover, the proposed method, which requires fewer arithmetic operations than the alternating-direction-implicit（ADI） CFDTD method, leads to a further reduction of the CPU time. With the total-field/scattered-field（TF/SF） boundary and the perfectly matched layer（PML）, the radar cross section（RCS） of two2 D structures is calculated. The numerical examples verify the accuracy and efficiency of the proposed method.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61301063 and 41305017)
文摘A higher-order finite-difference time-domain（HO-FDTD） in the spherical coordinate is presented in this paper. The stability and dispersion properties of the proposed scheme are investigated and an air-filled spherical resonator is modeled in order to demonstrate the advantage of this scheme over the finite-difference time-domain（FDTD） and the multiresolution time-domain（MRTD） schemes with respect to memory requirements and CPU time. Moreover, the Berenger＇s perfectly matched layer（PML） is derived for the spherical HO-FDTD grids, and the numerical results validate the efficiency of the PML.
基金supported by the National Natural Science Foundation of China(Grant Nos.60931002 and 61101064)the Universities Natural Science Foundation of Anhui Province,China(Grant Nos.KJ2011A002 and 1108085J01)
文摘The method of splitting a plane-wave finite-difference time-domain （SP-FDTD） algorithm is presented for the initiation of plane-wave source in the total-field / scattered-field （TF/SF） formulation of high-order symplectic finite- difference time-domain （SFDTD） scheme for the first time. By splitting the fields on one-dimensional grid and using the nature of numerical plane-wave in finite-difference time-domain （FDTD）, the identical dispersion relation can be obtained and proved between the one-dimensional and three-dimensional grids. An efficient plane-wave source is simulated on one-dimensional grid and a perfect match can be achieved for a plane-wave propagating at any angle forming an integer grid cell ratio. Numerical simulations show that the method is valid for SFDTD and the residual field in SF region is shrinked down to -300 dB.
基金Project supported by the National Natural Science Foundation of China(Grant No.51025622)
文摘To deal with the staircase approximation problem in the standard finite-difference time-domain（FDTD） simulation,the two-dimensional boundary condition equations（BCE） method is proposed in this paper.In the BCE method,the standard FDTD algorithm can be used as usual,and the curved surface is treated by adding the boundary condition equations.Thus,while maintaining the simplicity and computational efficiency of the standard FDTD algorithm,the BCE method can solve the staircase approximation problem.The BCE method is validated by analyzing near field and far field scattering properties of the PEC and dielectric cylinders.The results show that the BCE method can maintain a second-order accuracy by eliminating the staircase approximation errors.Moreover,the results of the BCE method show good accuracy for cylinder scattering cases with different permittivities.
基金Project supported by the National Natural Science Foundation of China (Grant No 60571058)the Specialized Research Fund for the Doctoral Program of Higher Education,China (Grant No 20070701010)
文摘Composite electromagnetic scattering from a two-dimensional （2D） ship-like target on a one-dimensional sea surface is investigated by using the finite-difference time-domain （FDTD） method. A uniaxial perfectly matched layer is adopted for truncation of FDTD lattices.The FDTD updated equations can be used for the total computation domain by choosing the uniaxial parameters properly. To validate the proposed numerical technique,a 2D infinitely long cylinder over the sea surface is taken into account first.The variation of angular distribution of the scattering changing with incident angle is calculated. The results show good agreement with the conventional moment method. Finally,the influence of the incident angle,the polarization,and the size of the ship-like target on the composite scattering coefficient is discussed in detail.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61471105 and 61331007)
文摘In this paper, a modified sub-gridding scheme that hybridizes the conventional finite-difference time-domain（FDTD）method and the unconditionally stable locally one-dimensional（LOD） FDTD is developed for analyzing the periodic metallic nanoparticle arrays. The dispersion of the metal, caused by the evanescent wave propagating along the metal-dielectric interface, is expressed by the Drude model and solved with a generalized auxiliary differential equation（ADE） technique.In the sub-gridding scheme, the ADE–FDTD is applied to the global coarse grids while the ADE–LOD–FDTD is applied to the local fine grids. The time step sizes in the fine-grid region and coarse-grid region can be synchronized, and thus obviating the temporal interpolation of the fields in the time-marching process. Numerical examples about extraordinary optical transmission through the periodic metallic nanoparticle array are provided to show the accuracy and efficiency of the proposed method.
基金supported by National Natural Science Foundation of China(Nos.41104097,11504252,61201007,41304119)the Fundamental Research Funds for the Central Universities(Nos.ZYGX2015J039,ZYGX2015J041)the Specialized Research Fund for the Doctoral Program of Higher Education of China(No.20120185120012)
文摘The frequency dependent permittivity for dusty plasmas is provided by introducing the charging response factor and charge relaxation rate of airborne particles. The field equations that describe the characteristics of Terahertz（THz） waves propagation in a dusty plasma sheath are derived and discretized on the basis of the auxiliary differential equation（ADE） in the finite difference time domain（FDTD） method. Compared with numerical solutions in reference, the accuracy for the ADE FDTD method is validated. The reflection property of the metal Aluminum interlayer of the sheath at THz frequencies is discussed. The effects of the thickness, effective collision frequency, airborne particle density, and charge relaxation rate of airborne particles on the electromagnetic properties of Terahertz waves through a dusty plasma slab are investigated.Finally, some potential applications for Terahertz waves in information and communication are analyzed.
基金National Natural Science Foundation of China (No. 60471002) and the Natural Science Foundation ofJiangxi Province (No. 0412014)
文摘The piecewise linear recursive convolution （PLRC） finite-different time-domain （FDTD） method improves accuracy over the original recursive convolution （RC） FDTD approach and current density convolution （JEC） but retains their advantages in speed and efficiency. This paper describes a revised piecewise linear recursive convolution PLRC-FDTD formulation for magnetized plasma which incorporates both anisotropy and frequency dispersion at the same time, enabling the transient analysis of magnetized plasma media. The technique is illustrated by numerical simulations of the reflection and transmission coefficients through a magnetized plasma layer. The results show that the revised PLRC-FDTD method has improved the accuracy over the original RC FDTD method and JEC FDTD method.
基金Supported by the Specialized Research Fund for the Doctoral Program of Higher Education(No.20010614003)
文摘In order to eliminate Courant-Friedrich-Levy(CFL) condition restraint and improvecomputational efficiency,a new finite-difference time-domain(FDTD)method based on the alternating-direction implicit(ADI) technique is introduced recently.In this paper,a theoretical proof of the stabilityof the three-dimensional(3-D)ADI-FDTD method is presented.It is shown that the 3-D ADI-FDTDmethod is unconditionally stable and free from the CFL condition restraint.
基金Supported by the National Natural Science Foundation of China(11274173,11374159)the Fundamental Research Funds for the Central Universities(NZ2013304,NJ20140005)Qing Lan Project of Jiangsu Province
文摘The surface plasmon resonance( SPR) of gold( Au) nanorod can be tuned in a large visible-near infrared( Vis-NIR) region by changing the aspect ratio of nanorod. Compared with the SPR of isolated Au nanorod,assembly of Au nanorods exhibits strong coupling effect in the nanogap and wealthy changes in the optical spectra. The SPR coupling effects and localized electronic fields for end-to-end( E-E) and side-by-side( S-S) assembled Au nanorods dimers are studied through finite-difference time-domain( FDTD) simulation. With decreasing the gap spacing,the longitudinal SPR( SPRL) red-shifts for the E-E orients Au nanorods dimer and blue-shifts for S-S orientes Au nanorods dimer. The transverse SPR( SPRT) has slight red-shifting for S-S assembly and no shifting for the E-E assembly. Moreover,a new coupling SPR appears for the E-E assembly in a long wavelength in the NIR region,blue-shifting and enhancing with decreasing the gap spacing. Based on the spring oscillator model and the polarization of nanoparticles with incident electric field,the SPR shifting and the appearance of new coupling SPR of assembled Au nanorods are proposed.
基金supported by National Natural Science Foundation of China(Nos.60271005 and 60431010)the National Science Fund for Distinguished Young Scholars of China(No.60325103)
文摘The analysis of electromagnetic propagation in a dispersive medium is complicated in the time-domain because its dielectric constant is frequency-dependent. In this paper, the dielectric constant of the dispersive medium is written as a rational polynomial function, and the relationship between D and E is derived in the time-domain. It is referred to as the shift operator finite-different time-domain （SO-FDTD） method. Compared to an analytical solution and a piecewise linear current density recursive convolution （PLJERC） solution, the high accuracy and efl%iency of this method is verified by calculating the reflectance of the electromagnetic wave through a cold plasma slab. As the electron density in plasma is distributed as the Epstein formula, the effect of distribution grads and electron collision frequency on the reflectance is calculated by using the SO-FDTD method. The result shows that the increase in the distribution grads coefficient affects the reflectance sharply. When it comes to a smaller distribution grads coelBcient, the increase of the collision frequency showed a significant effect on the reflectance, but on the contrary, there is actually less and less effect till it disappears.
文摘Based on FDTD difference expressions and eigenfunctions of Maxwell functions in cylindrical coordinates, mesh wave impedances (MWIs) in 2D and 3D cylindrical coordinates were introduced. Combined with the concept of perfectly matched layer (PML), MWI PML absorbing boundary condition (ABC) algorithm was deduced in 2D cylindrical coordinates. Numerical experiments were done to investigate the validity of MWI and its application in cylindrical coordinates FDTD algorithm. The results showed that MWI in cylindrical coordinates can be used to accurately calculate the numerical reflection error caused by different mesh increments in non uniform FDTD. MWI can also provide theoretical criterion to define the permitted variable range of mesh dimension. MWI PML ABC is easy to be applied and reduces low numerical reflection, which only causes a little higher reflection error compared with Teixeira's PML.
基金Project supported by the National Natural Science Foundation of China (Grant No. 61071022)the Graduate Student Research and Innovation Program of Jiangsu Province,China (Grant No. CXZZ11-0381)
文摘A novel method based on the relevance vector machine（RVM） for the inverse scattering problem is presented in this paper.The nonlinearity and the ill-posedness inherent in this problem are simultaneously considered.The nonlinearity is embodied in the relation between the scattered field and the target property,which can be obtained through the RVM training process.Besides,rather than utilizing regularization,the ill-posed nature of the inversion is naturally accounted for because the RVM can produce a probabilistic output.Simulation results reveal that the proposed RVM-based approach can provide comparative performances in terms of accuracy,convergence,robustness,generalization,and improved performance in terms of sparse property in comparison with the support vector machine（SVM） based approach.
基金Project supported by the National Basic Research Program of China(Grant No.2014CB340203)the National Natural Science Foundation of China(Grant Nos.61431010 and 61501350)the Natural Science Foundation of Shaanxi Province,China(Grant Nos.2018JM6016 and 2016JM1001)
文摘The propagation characteristics of oblique incidence terahertz(THz) waves through non-uniform plasma are investigated by the shift-operator finite-difference time-domain(SO-FDTD) method combined with the phase matching condition.The electron density distribution of the non-uniform plasma is assumed to be in a Gaussian profile. Validation of the present method is performed by comparing the results with those obtained by an analytical method for a homogeneous plasma slab.Then the effects of parameters of THz wave and plasma layer on the propagation properties are analyzed. It is found that the transmission coefficients greatly depend on the incident angle as well as on the thickness of the plasma, while the polarization of the incident wave has little influence on the propagation process in the range of frequency considered in this paper. The results confirm that the THz wave can pass through the plasma sheath effectively under certain conditions,which makes it a potential candidate to overcome the ionization blackout problem.
基金supported by the Program for New Century Excellent Talents in University(No.NCET-05-0575)the Education Science Foundation of Jiangxi Province(No.Z-03510)
文摘In this study, the propagation of electromagnetic waves in one-dimensional plasma photonic crystals （PPCs）, namely, superlattice structures consisting alternately of a homogeneous unmagnetized plasma and dielectric material, is simulated numerically using the finite-difference time-domain （FDTD） algorithm. A perfectly matched layer （PML） absorbing technique is used in this simulation. The reflection and transmission coefficients of electromagnetic （EM） waves through PPCs are calculated. The characteristics of the photonic band gap （PBG） are discussed in terms of plasma density, dielectric constant ratios, number of periods, and introduced layer defect. These may provide some useful information for designing plasma photonic crystal devices.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.61137005 and 61178067)the Science Foundation of Shanxi Province,China(Grant No.2013021004-3/2014021021-1)+1 种基金the Pre-studied Project on Weapon Equipment,China(Grant No.201262401090404)the Specialized Research Foundation for Doctor of School,China(Grant No.20122027)
文摘In this paper, a subwavelength metal-grating assisted sensor of Kretschmann style that is capable of detecting the sample with a refractive index higher than that of the substrate is proposed. The sensor configuration is similar to the traditional Kretschmann structure, but the metal film is pattered into a grating. As a TM-polarized laser beam impinges from the substrate, a resonant dip point in reflectance curve is produced at a certain incident angle. Our studies indicate that the sensing sensitivity and resolution are affected by the grating＇s gap and period, and after these parameters havelaeen optimized, a sensing sensitivity of 51.484°/RIU is obtained with a slightly changing resolution.
基金Project supported by the National Natural Science Foundation of China (Grant No 10674051)the Program for Innovative Research Team of the Higher Education of Guangdong Province,China (Grant No 06CXTD005)the Key Program of Extracurricular Research in South China Normal University (SCNU),China (Grant No 08GDKC02)
文摘Three-dimensional photonic crystal （PC） heterostructures with high quality are fabricated by using a pressure controlled isothermal heating vertical deposition technique. The formed heterostructures have higher quality, such as deeper band gaps and sharper band edges, than the heterostructures reported so far. Such a significant improvement in quality is due to the introduction of a thin TiO2 buffer layer between the two constitutional PCs. It is revealed that the disorder caused by lattice mismatch is successfully removed if the buffer layer is used once. As a result, the formed heterostructures possess the main features in the band gap of constitutional PCs. The crucial role of the thin buffer layer is also verified by numerical simulations based on the finite-difference time-domain technique.
