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.展开更多
To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed...To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Edge structures are ubiquitous in the processing and fabrication of various optoelectronic devices.Novel physical properties and enhanced light–matter interactions are anticipated to occur at crystal edges due to the...Edge structures are ubiquitous in the processing and fabrication of various optoelectronic devices.Novel physical properties and enhanced light–matter interactions are anticipated to occur at crystal edges due to the broken spatial translational symmetry.However,the intensity of first-order Raman scattering at crystal edges has been rarely explored,although the mechanical stress and edge characteristics have been thoroughly studied by the Raman peak shift and the spectral features of the edge-related Raman modes.Here,by taking Ga As crystal with a well-defined edge as an example,we reveal the intensity enhancement of Raman-active modes and the emergence of Raman-forbidden modes under specific polarization configurations at the edge.This is attributed to the presence of a hot spot at the edge due to the redistributed electromagnetic fields and electromagnetic wave propagations of incident laser and Raman signal near the edge,which are confirmed by the finite-difference time-domain simulations.Spatially-resolved Raman intensities of both Raman-active and Raman-forbidden modes near the edge are calculated based on the redistributed electromagnetic fields,which quantitatively reproduce the corresponding experimental results.These findings offer new insights into the intensity enhancement of Raman scattering at crystal edges and present a new avenue to manipulate light–matter interactions of crystal by manufacturing various types of edges and to characterize the edge structures in photonic and optoelectronic devices.展开更多
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.展开更多
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.展开更多
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.展开更多
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.展开更多
Metal-semiconductor Janus nanostructures with asymmetry and directionality have recently aroused significant interest,both in fundamental light-matter interactions mechanism and in technological applications.Here we r...Metal-semiconductor Janus nanostructures with asymmetry and directionality have recently aroused significant interest,both in fundamental light-matter interactions mechanism and in technological applications.Here we report the synthesis of different Au-ZnO Janus nanostructures via a facile one-pot colloid method.The growth mechanism is revealed by a series of designed synthesis experiments.The light absorption properties are determined by both the decrease of dipole oscillations of the free electrons and the plasmon-induced hot-electron transfer.Moreover,the finite-difference time-domain(FDTD)simulation method is used to elucidate the electric field distributions of these Janus nanostructures.展开更多
The influence ot oxidation aperture on the output characteristics ot the circularly symmetric vertical-cavity-surtaceemitting laser(VCSEL) structure is investigated.To do so,VCSELs with different oxide aperture sizes ...The influence ot oxidation aperture on the output characteristics ot the circularly symmetric vertical-cavity-surtaceemitting laser(VCSEL) structure is investigated.To do so,VCSELs with different oxide aperture sizes are simulated by the finite-difference time-domain(FDTD) method.The relationships among the field distribution of mode superposition,mode wavelength,output spectra,and far-field divergence with different oxide apertures are obtained.Further,VCSELs respectively with oxide aperture sizes of 2.7 μm,4.4 μm,5.9 μm,7 μm,8 μm,9 μm,and 18.7 μm are fabricated and characterized.The maximum output power increases from 2.4 mW to 5.7 mW with oxide aperture increasing from 5.9 μm to 9 μm.Meanwhile,the wavelength tuning rate decreases from 0.93 nm/mA to 0.375 nm/mA when the oxide aperture increases from 2.7 μm to 9 μm.The thermal resistance decreases from 2.815℃/mW to 1.015℃/mW when the oxide aperture increases from 4.4 μm to 18.7μm.It is demonstrated theoretically and experimentally that the wavelength spacing between adjacent modes increases with the augment of the injection current and the spacing becomes smaller with the oxide aperture increasing.Thus it can be reported that the aperture size can effectively reduce the mode overlaying but at the cost of the power decreasing and the wavelength tuning rate and thermal resistance increasing.展开更多
Theoretical studies show that a Hertzian-conical crack can be considered to be composed of double cone faces for simplicity. In the present study, the three-dimensional finite-difference time-domain method is employed...Theoretical studies show that a Hertzian-conical crack can be considered to be composed of double cone faces for simplicity. In the present study, the three-dimensional finite-difference time-domain method is employed to quantify the electric-field distribution within the subsurface in the presence of such a defect under normal incidence irradiation. Both impurities （inside the crack） and the chemical etching have been investigated. The results show that the maximum electric field amplitude │E│max is 9.57374 V/m when the relative dielectric constant of transparent impurity equals 8.5. And the near-field modulation will be improved if the crack is filled with the remainder polishing powders or water vapor/drops. Meanwhile, the laser-induced initial damage moves to the glass-air surface. In the etched section, the magnitude of intensification is strongly dependent on the inclination angle θ. There will be a highest modulation when θ is around π/6, and the maximum value of IEImax is 18.57314 V/m. When θ ranges from π/θ to π/4, the light intensity enhancement factor can easily be larger than 100, and the modulation follows a decreasing trend. On the other hand, the modulation curves become smooth when θ〉π/4 or θ 〈 π/8.展开更多
基金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 under (Grant No.107 72040,50709005 and 50921001)the Major National Science and Technology Projects of China under (Grant No.2008ZX05026-02)the Open Fund of State Key Laboratory of Ocean Engineering
文摘To study wave-current actions on 3-D bodies a time-domain numerical model was established using a higher-order boundary element method(HOBEM).By assuming small flow velocities,the velocity potential could be expressed for linear and higher order components by perturbation expansion.A 4th-order Runge-Kutta method was applied for time marching.An artificial damping layer was adopted at the outer zone of the free surface mesh to dissipate scattering waves.Validation of the numerical method was carried out on run-up,wave exciting forces,and mean drift forces for wave-currents acting on a bottom-mounted vertical cylinder.The results were in close agreement with the results of a frequency-domain method and a published time-domain method.The model was then applied to compute wave-current forces and run-up on a Seastar mini tension-leg platform.
基金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.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.
基金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 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 Key Research and Development Program of China(Grant No.2023YFA1407000)the Strategic Priority Research Program of Chinese Academy of Sciences(Grant No.XDB0460000)+4 种基金the National Natural Science Foundation of China(Grant Nos.12322401,12127807,and 12393832)CAS Key Research Program of Frontier Sciences(Grant No.ZDBS-LY-SLH004)Beijing Nova Program(Grant No.20230484301)Youth Innovation Promotion Association,Chinese Academy of Sciences(Grant No.2023125)CAS Project for Young Scientists in Basic Research(Grant No.YSBR-026)。
文摘Edge structures are ubiquitous in the processing and fabrication of various optoelectronic devices.Novel physical properties and enhanced light–matter interactions are anticipated to occur at crystal edges due to the broken spatial translational symmetry.However,the intensity of first-order Raman scattering at crystal edges has been rarely explored,although the mechanical stress and edge characteristics have been thoroughly studied by the Raman peak shift and the spectral features of the edge-related Raman modes.Here,by taking Ga As crystal with a well-defined edge as an example,we reveal the intensity enhancement of Raman-active modes and the emergence of Raman-forbidden modes under specific polarization configurations at the edge.This is attributed to the presence of a hot spot at the edge due to the redistributed electromagnetic fields and electromagnetic wave propagations of incident laser and Raman signal near the edge,which are confirmed by the finite-difference time-domain simulations.Spatially-resolved Raman intensities of both Raman-active and Raman-forbidden modes near the edge are calculated based on the redistributed electromagnetic fields,which quantitatively reproduce the corresponding experimental results.These findings offer new insights into the intensity enhancement of Raman scattering at crystal edges and present a new avenue to manipulate light–matter interactions of crystal by manufacturing various types of edges and to characterize the edge structures in photonic and optoelectronic devices.
文摘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 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 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.
基金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.
基金Project supported by the National Natural Science Foundation of China(Grant Nos.11274370 and 51471185)the National Key Research and Development Program of China(Grant Nos.2016YFJC020013 and 2018FYA0305800)
文摘Metal-semiconductor Janus nanostructures with asymmetry and directionality have recently aroused significant interest,both in fundamental light-matter interactions mechanism and in technological applications.Here we report the synthesis of different Au-ZnO Janus nanostructures via a facile one-pot colloid method.The growth mechanism is revealed by a series of designed synthesis experiments.The light absorption properties are determined by both the decrease of dipole oscillations of the free electrons and the plasmon-induced hot-electron transfer.Moreover,the finite-difference time-domain(FDTD)simulation method is used to elucidate the electric field distributions of these Janus nanostructures.
文摘The influence ot oxidation aperture on the output characteristics ot the circularly symmetric vertical-cavity-surtaceemitting laser(VCSEL) structure is investigated.To do so,VCSELs with different oxide aperture sizes are simulated by the finite-difference time-domain(FDTD) method.The relationships among the field distribution of mode superposition,mode wavelength,output spectra,and far-field divergence with different oxide apertures are obtained.Further,VCSELs respectively with oxide aperture sizes of 2.7 μm,4.4 μm,5.9 μm,7 μm,8 μm,9 μm,and 18.7 μm are fabricated and characterized.The maximum output power increases from 2.4 mW to 5.7 mW with oxide aperture increasing from 5.9 μm to 9 μm.Meanwhile,the wavelength tuning rate decreases from 0.93 nm/mA to 0.375 nm/mA when the oxide aperture increases from 2.7 μm to 9 μm.The thermal resistance decreases from 2.815℃/mW to 1.015℃/mW when the oxide aperture increases from 4.4 μm to 18.7μm.It is demonstrated theoretically and experimentally that the wavelength spacing between adjacent modes increases with the augment of the injection current and the spacing becomes smaller with the oxide aperture increasing.Thus it can be reported that the aperture size can effectively reduce the mode overlaying but at the cost of the power decreasing and the wavelength tuning rate and thermal resistance increasing.
基金Project supported by the Young Scientists Fund of the National Natural Science Foundation of China (Grant No.10904008)the Joint Fund of the National Natural Science Foundation of China and the China Academy of Engineering Physics (GrantNo. 11076008)the Scientific Research Foundation for the Central Universities of China (Grant Nos.ZYGX2009X007,ZYGX2010J045, and ZYGX2011J043)
文摘Theoretical studies show that a Hertzian-conical crack can be considered to be composed of double cone faces for simplicity. In the present study, the three-dimensional finite-difference time-domain method is employed to quantify the electric-field distribution within the subsurface in the presence of such a defect under normal incidence irradiation. Both impurities （inside the crack） and the chemical etching have been investigated. The results show that the maximum electric field amplitude │E│max is 9.57374 V/m when the relative dielectric constant of transparent impurity equals 8.5. And the near-field modulation will be improved if the crack is filled with the remainder polishing powders or water vapor/drops. Meanwhile, the laser-induced initial damage moves to the glass-air surface. In the etched section, the magnitude of intensification is strongly dependent on the inclination angle θ. There will be a highest modulation when θ is around π/6, and the maximum value of IEImax is 18.57314 V/m. When θ ranges from π/θ to π/4, the light intensity enhancement factor can easily be larger than 100, and the modulation follows a decreasing trend. On the other hand, the modulation curves become smooth when θ〉π/4 or θ 〈 π/8.
