Simulation of confined and enhanced optical near fields for a long narrow aperture in a pyramidal structure on a thick metallic screen

A new type of long narrow aperture in a pyramidal structure on a thick metallic screen is proposed, and optical wave scattering by this structure is simulated. This aperture structure provides high emission intensity and small spot size simultaneously through excitation of the surface plasmon polaritons on the sidewalls of the pyramidal structure. Scattering of optical waves by this structure in the thick metallic screen is solved numerically with a volume integral equation by generalized conjugate residual iteration and fast Fourier transformation. The basic characteristics of the near-field intensities of the aperture are investigated in detail.     

Single-scattering theory of light diffraction by a circular subwavelength aperture in a finitely conducting screen

A perturbation theory based on a single-scattering approximation is developed from the rigorous differential theory of diffraction in cylindrical coordinates. It results in analytical formulas in the inverse space for the field amplitudes providing results that are in quantitative agreement with the results of the rigorous method, in both the near- and far-field regions, when a proper correction to the incident field inside the aperture is made by using the renormalized Born approximation. When working in reflection by a screen having permittivity high in modulus, the method proposes an equivalence with the simple model consisting of the emission by a single magnetic dipole excited inside the pierced layer, emission that is then transferred back into the cladding following the Fresnel's coefficients of transmission from the layer into the cladding. The theory predicts a directivity of the radiation pattern that increases for smaller values of modulus of permittivity, both for dielectrics and metals, thus independently of the possibility of plasmon surface wave excitation along the interface. The theory can take into account such surface wave resonances, as well as the waveguide supported by a dielectric slab, but cannot implicitly recognize the modes carried out by the cylindrical waveguide corresponding to the aperture. This fact limits its domain of validity when used in transmission, although the far- and near-field maps can be reconstructed sufficiently well within a multiplicative factor corresponding to the enhanced transmission due to the excitation of these modes.     

An integral equation method for the diffraction of oblique waves by an infinite cylinder

A hybrid integral equation method is formulated to study the diffraction of oblique waves by an infinite cylinder. The water depth and the geometry of the floating cylinder are assumed to be uniform in the y-direction (one of the horizontal axes). Numerical discretization and integrations are performed in the vertical plane. Analytical solutions are used in far fields such that radiation boundary conditions are satisfied. Numerical results are obtained for the case of wave scattering by a floating rectangular cylinder in a constant water depth. The accuracy and efficiency of present method are compared with those obtained by other numerical techniques.     

Phase retrieval from diffraction intensities by use of a scanning slit aperture

A noniterative method of retrieving the phase of a wave field from intensities measured during scanning of a slit aperture is proposed. In the measurements, one records the diffraction intensities of wave fields transmitted through a slit that is scanned along two directions in the Fresnel-zone plane of an object's field. From these intensities, the phase in the Fresnel-zone plane can be retrieved by a method in which a novel phase calculation technique that uses Fourier transforms is included. Because the method does not require lens systems, it provides a potentially useful means for coherent imaging by use of x rays, electrons, or nuclear particles.     

Analysis of the polarization-dependent diffraction from a metallic grating by use of a three-dimensional combined vectorial method

A three-dimensional combined vectorial method, which is based on the finite-difference time-domain algorithm and vectorial diffraction formulation, is introduced to analyze the interaction of a realistic focused beam with a metallic grating in an optical storage system. The diffracted field patterns and the detected signals are calculated for gratings with different geometries, and the polarization-dependent characteristics are studied. The combined method can give accurately the detected signals and the field pattern at any desired position of the optical storage system.     

Fresnel diffraction by a circular aperture with off-axis illumination and its use in deconvolution of microscope images

The Fresnel approximation for off-axis illumination of a circular aperture is reexamined. The point-spread function for an aberration-free system can be expressed in terms of redefined optical coordinates. An improved expression is given for contours of constant intensity in the focal plane. The variation in axial width of the focal spot with angle of offset is discussed. The predictions are compared with exact calculations of the Rayleigh-Sommerfeld diffraction integral. Limitations for application in deconvolution of microscope images formed with objectives of finite tube length are discussed.     

Characterization of optical diffraction gratings by use of a neural method.

Optical scatterometry by use of a neural network is now recognized as an efficient method for retrieving dimensions of gratings in semiconductors or glasses. For an on-line control, a small number of measurements and a rapid data treatment are needed. We demonstrate that these requirements can be met by combining data preprocessing and a proper neural learning method. A good accuracy is attainable with the measurement of only a few orders, even in the presence of experimental errors, with a reduction in learning and computing time.     

Analysis and observation of finite beam Bragg diffraction by a thick planar phase grating

A numerical-impulse-response technique for studying the propagation and diffraction of finite-width beams in planar phase gratings is described. It can account for both symmetric and asymmetric diffractions. The grating-length-to-beam-width ratio is shown to govern the extent of beam-profile distortion and selectivity sidelobe suppression. Trade-offs between diffraction efficiency and beam profile have also been demonstrated. Theoretical results have been verified by experimental observations in a planar waveguide geometry of diffracted beams that change from a single diffraction peak to multiple peaks as the grating-length-to-beam-width ratio increases.     

Elastic analysis of a half-plane containing an inclusion and a void using a mixed volume and boundary integral equation method

A mixed volume and boundary integral equation method is used to calculate the plane elastostatic field in an isotropic elastic half-plane containing an isotropic or anisotropic inclusion and a void subject to remote loading parallel to a traction-free boundary. A detailed analysis of the stress field is carried out for three different geometries of the problem. It is demonstrated that the method is very accurate and effective for investigating local stresses in an isotropic elastic half-plane containing multiple isotropic or anisotropic inclusions and multiple voids.     

A hybrid integral equation method for diffraction and radiation of water waves by three-dimensional bodies

A hybrid integral equation method is formulated for the analysis of water wave diffraction and radiation by arbitrarily-shaped three-dimensional bodies. It is based upon the direct application of Green's second identity and uses the simple fundamental solution (i. e. a simple source) rather than the special Green's function. The boundary element idealisation is used only in an inner fluid region close to the body and local depth irregularities, while an analytical solution is employed in the outer region of constant depth extending to infinity. The two representations are matched on a fictitious vertical cylindrical surface. Special treatment is discussed to reduce computational efforts when the body has one or two planes of symmetry. Numerical results are presented for a variety of geometries to illustrate the applicability and potentialities of the method.This is a revised and expanded version of the paper originally presented at the 7th International Conference on Boundary Element Methods in Engineering, Villa Olmo, Lake Como, Italy, September 1985     

Analysis of an optical depth converter used in a three-dimensional integral imaging system

An optical depth converter that uses a lens array pair is analyzed theoretically and experimentally. We present a theory of depth conversion and explain the effects of the system parameters in the optical depth converter by using wave-optical analysis. Ray-optical analysis is applied to the investigation of the tendencies of the system parameter effects. We also show that the optical depth converter can be used for the three-dimensional screen in projection-type integral imaging systems.     

Modeling microlenses by use of vectorial field rays and diffraction integrals

A nonparaxial vector-field method is used to describe the behavior of low-f-number microlenses by use of ray propagation, Fresnel coefficients and the solution of Maxwell equations to determine the field propagating through the lens boundaries, followed by use of the Rayleigh-Sommerfeld method to find the diffracted field behind the lenses. This approach enables the phase, the amplitude, and the polarization of the diffracted fields to be determined. Numerical simulations for a convex-plano lens illustrate the effects of the radii of curvature, the lens apertures, the index of refraction, and the wavelength on the variations of the focal length, the focal plane field distribution, and the cross polarization of the field in the focal plane.     

Improving the accuracy of the magnetic field integral equation with the use of a new impedance matrix element formulation

The electric field integral equation (EFIE) and the magnetic field integral equation (MFIE) are widely used in conjunction with method of moments for electromagnetic scattering analysis of three-dimensional conducting objects with closed surfaces. However, the MFIE suffers from an accuracy problem compared with the EFIE with the use of the Rao?Wilton?Glisson (RWG) basis function. This accuracy problem is more serious for objects with sharp edges or corners. To solve this problem, a new technique to compute the impedance matrix elements (IME) of the MFIE using an RWG basis function is presented here. Details to compute the IME and the advantage of this new formulation are displayed. In addition, the relationship between this new IME formulation and the formulation using the low-order curl-conforming basis function for the MFIE is given. Through the computation of the RCS of several relatively small sharp-edged conducting objects, it is shown that the accuracy of the MFIE can be greatly improved by the use of the new IME formulation.     

Transmission of an inhomogeneous plane wave through an electrically small aperture in a perfectly conducting plane screen

Most solutions for electromagnetic diffraction by a circular aperture in a perfectly conducting plane screen are for an incident homogeneous (propagating) plane wave. When the aperture is electrically small (dimensions small compared to the wavelength), the well-known transmission coefficient behaves as the fourth power of the diameter/wavelength. We consider the case in which the incident field is an inhomogeneous (evanescent) plane wave. Numerical calculations for the electrically small circular aperture show that the transmission coefficient for an inhomogeneous plane wave can be substantially greater than for a homogeneous plane wave at the same frequency. This observation may be helpful in explaining the increased transmission recently reported for electrically small apertures in plane screens with modifications. The numerical calculations for the electrically small aperture are in agreement with results from approximate analytical expressions that are based on the equivalent electric and magnetic dipole moments for the electrically small complementary disk.     

Analysis of fractures in 3D piezoelectric media by a weakly singular integral equation method

A weakly singular, symmetric Galerkin boundary element method (SGBEM) is established to compute stress and electric intensity factors for isolated cracks in three-dimensional, generally anisotropic, piezoelectric media. The method is based upon a weak-form integral equation, for the surface traction and the surface electric charge, which is established by means of a systematic regularization procedure; the integral equation is in a symmetric form and is completely regularized in the sense that its integrand contains only weakly singular kernels of (hence allowing continuous interpolations to be employed in the numerical approximation). The weakly singular kernels which appear in the weak-form integral equation are expressed explicitly, for general anisotropy, in terms of a line integral over a unit circle. In the numerical implementation, a special crack-tip element is adopted to discretize the region near the crack front while the remainder of the crack surface is discretized by standard continuous elements. The special crack-tip element allows the relative crack-face displacement and electric potential in the vicinity of the crack front to be captured to high accuracy (even with relatively large elements), and it has the important feature that the mixed-mode intensity factors can be directly and independently extracted from the crack front nodal data. To enhance the computational efficiency of the method, special integration quadratures are adopted to treat both singular and nearly singular integrals, and an interpolation strategy is developed to approximate the weakly singular kernels. As demonstrated by various numerical examples for both planar and non-planar fractures, the method gives rise to highly accurate intensity factors with only a weak dependence on mesh refinement.     

Near-field anomalous spectral behavior in diffraction of a Gaussian pulsed beam from an annular aperture

Based on the Fresnel diffraction integral and by introducing a hard-aperture function into a finite sum of complex Gaussian functions, the approximate analytical expression for the near-field spectral intensity distribution of a space-time-dependent Gaussian pulsed beam passing through an annular aperture is derived, which permits us to study the on- and off-axis spectral anomalies that are near phase singularities of the diffracted Gaussian pulsed beam in the near-field. The expressions for a circular black screen and a circular aperture are given as special cases of the general results. The relative spectral shift of a space-time-dependent Gaussian pulsed beam versus the different values of the truncation parameters and the position parameters of observation points are also studied and illustrated with numerical calculations. It is shown that the spectral switch appears near phase singularities in the near-field, and the near-field spectral behavior depends on the truncation parameters, the pulse duration tau, and the position parameter. The results of this work have potential applications in free-space information encoding and transmission.     

Spectral changes and 1 X N spectral switches in the diffraction of partially coherent light by an aperture

The off-axis and on-axis spectra in the far zone of an aperture for the case in which a particular class of partially coherent light with a broad spectrum is diffracted by an aperture are studied. It is shown that the spectrum in the far zone is generally different from that at the aperture; i.e., the spectrum is split into two or more peaks. Moreover, the spectrum varies with the diffractive angle. For a fixed diffractive angle, the spectral shift, defined as the difference between the frequencies at which the observed spectrum and the spectrum at the aperture take their maximum, shows a gradual change with the change in the coherence at the aperture. However, as the coherence reaches some critical values, the spectral shift exhibits a rapid transition; i.e., spectral switch occurs. The coherence that causes the spectral switch to take place is different for different diffractive angles. Therefore we propose a new kind of 1 x N spectral switch, where N detectors (output ports) are placed at different diffractive angles in the far zone, and the spectral shifts at different detectors are measured. By adjusting the coherence of the aperture (input port) to the desired values, we obtain a rapid transition of the spectral shift in the desired output ports.     

Demonstration of an optical isolator by use of a nonreciprocal phase shift

The experimental study of an optical isolator by use of a nonreciprocal phase shift is demonstrated. The isolator has an optical interferometer composed of tapered couplers, nonreciprocal phase shifters, and a reciprocal phase shifter. The isolator, designed for a 1.55-mum wavelength, was fabricated to investigate the characteristics of each component. The branching and coupling characteristics of the tapered coupler were measured. The nonreciprocal and reciprocal phase shifts were also evaluated. By applying an external magnetic field to the interferometer, we confirmed the nonreciprocal phase shift in the interferometric isolator.     

Simultaneous characterization of the shape and refractive index of transparent living cells by an optical aperture

We report a method of using optical aperture diffraction to simultaneously detect the image and the refractive index of a living cell. By simulation, it is found that, when a suitable gap is introduced between the cell and the aperture, the image of the cell on the detection plane can be amplified hundreds to thousands of times, and a limit of detection of 3e-4 to 9e-5 can be reached for the refractive index of the cell. Experiments show that this method is feasible to realize, but the achievement of such a detection system is yet to be proved.