We propose and demonstrate an excellent linearly frequency-swept laser diode (LD) for sensing system utilizing frequency-moduleted continuous-wave (FMCW) technique. In order to linearly sweep the optical frequency, we adopt a reference interferometer and an electric phase comparator. The interference beat signal of the reference interferometer is phase-compared with an external reference rectangular signal having a fixed frequency near the interference beat signal frequency by a lock-in amplifier. The error signal from the lock-in amplifier is fed back to the modulating signal of the injection current of the LD. Thus, a phase-locked loop composed of optical and electric circuits can be established, and the beat signal frequency is locked to the frequency of the reference signal. The optical frequency of the LD is, therefore, excellently linearly swept in time. In order to experimentally confirm the linearlity of the proposed method, we apply this light source to the FMCW reflectometry. Resultingly, the improvement of the linearity is estimated to be about 10 dB. And the theoretically limited spatial resolution of the FMCW reflectometry is achieved.

A new polarization splitter at optical frequencies is proposed. The basic structure of the device is a tapered velocity coupler which is composed of a straight and a dimensionally tapered slab waveguide on a LiNbO₃ substrate. The numerical results obtained with the finite difference method indicate that extinction ratios of polarization less than 2% for both TE and TM modes are possible of realization under moderate control voltages and that the splitting characteristics are stable over a wide range of frequencies.

A new simple formulation of absorbing boundary conditions with higher order approximation is proposed for the solution of Maxwell's equations with the finite-difference time-domain (FD-TD) method. Although this higher order approximation is based on the third order approximation of the one-way wave equations, we have succeeded in reducing it to an equation in a form quite similar to the second order appoximation. Numerical tests exhibit smaller reflection errors than the prevalent second order approximation.

A new type of variable beam splitter at optical frequencies is proposed. The basic structure of the device utilizes a tapered velocity coupler which is composed of a center slab waveguide of constant-thickness, constant-index type and two identical outer slab waveguides of constant-thickness, variable-index type. The coupler is assumed to be fabricated on a LiNbO₃ substrate, whith an external electric field applied in parallel with the optical axis. The numerical results obtained with the finite difference method show that a wide range of splitting ratios can be obtained with moderate drive voltages and that the splitting characteristics are stable over a wide range of frequencies.

This paper presents a new numerical procedure for solving the scattering wave by the moving surface. Recently, the author and her colleagues have proposed a novel numerical procedure of grid generation having a coordinate line coincident with an arbitrarily shaped moving boundary. The time dependent curvilinear coordinate system which coincides with a contour of moving boundary in a physical region is transformed into fixed rectangular coordinate system. The simple form for the transformation is made possible to choose the function between the physical region and the rectangular computational region. The FD-TD algorithm is exactly solved in this fixed rectangular coordinate system. In this paper, for the application of the FD-TD algorithm to the body fitted grid generation with moving boundary, the stability criterion of FD-TD algorithm for the body fitted grid generation with moving boundary in three dimensions is derived. The stability criterion is shown an upper bound for time step assuring stable numerical solutions. The study of scattering of electromagnetic and acoustic wave from a moving or a rotating body is very important for electromagnetic probing of moving body. The problem has been investigated in the past by numerous authors. One of them is solved by FD-TD method, where the linear interpolation is introduced for the movement. By using the presented transformation technique where time component is added to the grid generation, the time depending coordinate system can be transformed into fixed rectangular coordinate system. Then the problems are directly solved by FD-TD method in the transformed coordinate system. To verify this numemical technique, scattered field is evaluated in the case when a plane wave is normally incident on a moving perfectly conducting flat plate. The numerical results are compared with the exact ones and excellent agreement between them is obtained.

In this paper, the dispersion characteristics of magnetostatic surface waves (MSSW) in a nonlinear gyromagnetic medium are analytically investigated. Assuming the nonlinearity as the first order perturbation in permeability tensor of the gyromagnetic medium, the perturbation technique based on the multiple scale method is used to deduce the nonlinear dispersion relations for the MSSW. It is observed that for a given propagation constant of the MSSW the frequency decreases with microwave power. It is also observed that group velocity decreases, and as a consequence, delay time increases with power of the microwave. By evaluating the dependence of frequency on power and variation of group velocity on propagation constant within the accuracy of the perturbation theory, it is confirmed that the conditions for formation of bright soliton are not satisfied for the MSSW.

A three-dimensional analysis of Whispering-Gallery modes (W. G. modes) in a coaxial dielectric resonator is proposed and presented. The coaxial dielectric resonator is constructed from a lossy dielectric disk and ring which have diameters of several tens times as large as wavelength. Eigenvalue equations of the W. G. modes are derived rigorously from field expressions and boundary conditions. The resonant frequencies, unloaded Q values and field distributions are calculated numerically from the eigenvalue equations. These calculated results are in good agreement with experimental ones for an X band model. As a result, it is shown that a considerable quantity of modal energy can be confined in a loss-less gap between the disk and ring, and then the unloaded Q value is higher than that of a conventional dielectric disk and ring resonator.

In this paper, a novel coupler-type wavelength demultiplexer composed of the K⁺-ion diffused waveguides, which has an adjustment function for optimizing the diffusion depth, is proposed to achieve reliably the high extinction ratio. The optimization in the diffusion depth is made by repeating the K⁺-ion diffusion and extinction-ratio measurement alternatively, and the high extinction ratios more than 20 dB are measured reliably at both operation wavelengths of 0.6328 and 0.83 µm. Experimental results on the polarization dependence in the extinction-ratio adjustment are also reported.

We describe the characteristics of scattering and diffraction of plane E-waves by a lossy dielectric elliptic cylinder. The computational programs for calculating the analytic solutions for the diffraction of a lossy dielectric elliptic cylinder can be achieved. From the calculated results of the backscattering cross section (BSCS) (usually the radar cross section: RCS) and the forward-scattering cross section (FSCS) due to the cross-sectional shape and complex dielectric constant of the elliptic cylinder, the features of the BSCS and FSCS can be clarified as follows. (1) There is a cross-sectional shape of the cylinder which results in a minimum BSCS with a complex dielectric constant of the cylinder. (2) The BSCS and FSCS of the lossy dielectric scatterer approach zero as the scatterer approaches a strip. This result means that no material composing such a strip exists, and the features are very different from those in a perfectly conducting strip. (3) The influence of conductivity, σ, of the cylinder on a scattered wave is small for the relative dielectric constant of ε_r6. (4) The total scattering cross section of the lossy dielectric elliptic cylinder which causes the minimum BSCS is not small. Hence, it may be considered that the minimum BSCS is determined mainly by interference based on the cross-sectional shape and complex dielectric constant of cylinder, and is not caused by incident wave absorption due to the lossy dielectric.

In this paper, the influences of the cross-sectional deformation on the coplanar waveguide (CPW) characteristics for the use of Ti: LiNbO₃optical modulator are presented based on quasi-static analysis. In particular, the influences of the changes in the thickness of Ti: LiNbO₃ substrate and the cross-sectional shape of electrodes are studied in detail by using the finite element method proposed previously. As a result, it is found that the propagation characteristics of the dominant mode change significantly with the thickness of LiNbO₃ substrate when it is less than 100 µm. It is also shown that an inverted trapezoidal deformation of the electrode cross section is promising because a wider electrode gap and thinner electrodes are available in the design of optical modulators.

A coupled-mode analysis of a symmetric planar nonlinear directional coupler (NLDC) is presented by using a singular perturbation scheme. The effects of linear coupling and nonlinear modification of refractive index are treated to be small perturbations, and the modal fields of isolated linear waveguides are employed as the basis of propagation model. The self-consistent first-order coupled-mode equations governing the transfer of optical power along the NLDC are obtained in analytically closed form. It is shown that tha critical power for optical switching derived from the coupled-mode equations is in close agreement with that obtained by the numerical analysis using the finite difference beam propagation mathod.

Coupling optics for a monolithic array LD module has been developed. High efficient and small sized confocal optics with aspheric lenses based on hyperbolic surfaces has been designed to achieve the uniformity of coupling loss. A small sized (7.2 cc) 4 channel array LD module with this optics was fabricated. This module has low (4.2 dB) and uniform (0.5 dB) coupling loss.

A method is presented for reconstructing the surface profile of a perfectly conducting rough surface boundary from the measurements of the scattered far-field. The proposed inversion algorithm is based on the use of the Kirchhoff approximation and in order to determine the surface profile, the Fletcher-Powell optimization procedure is applied. A number of numerical results illustrating the method are presented.

The problem of electromagnetic scattering by inductive discontinuities located in rectangular waveguides, in particular when dealing with discontinuous conductors of finite thickness, is analyzed using the modified residue-calculus method, and form of the equation suitable for a numerical calculation is derived. The incident wave is taken to be the dominant mode, and reflection and transmission properties of an asymmetric inductive iris are discussed. After the modal representation of the filed, the modal matching is apply to satisfy the boundary conditions at the discontinuity. And using the modified residue-calculus method, simultaneous infinite equations, which are concerned with the scattered mode coefficients, are derived. Then they are approximated at the thick diaphragm. The solutions obtained take on the form of an infinite product, and a numerical solution based on the method of successive approximations is presented as a technique for concretely determining the reflection coefficients. As confirmation, experiments are also carried out in the X-band and close agreement is shown between the calculated and experimental values.

A new approach using radiation mode expansions is presented for calculating radiated fields from arbitrary distribution of electromagnetic sources in the half space region partitioned by a dielectric layer with a ground conductor. This method is applied to the calculation of radiation from microstrip-type antennas with a dielectric substrate of theoretically infinite extent. To be able to use this method, it is necessary to obtain first the field distribution around antenna patches, which is accomplished rather easily by using the FD-TD method. Radiation pattern calculations are presented for a rectangular patch antenna to verify the feasibility of this approach.

A particle trajectory estimation method from far electromagnetic fields are discussed in this paper. Authors have already presented a trajectory estimation method for single particle system and good agreements between a source particle trajectory and an estimated one have been obtained. For this, this paper discusses twin particles system as an examples of multi-particles systems for simplicity. First of all, it is pointed out that far electromagnetic fields from the twin particles system show quite different aspect from the single particle system using an example, radiation patterns produced by two particles which carry out circular motion. This result tells us that any trajectory estimations for general multi-particles system are almost impossible. However, it is shown that when the distance between the particles is small, the estimation method for the single particle system can be applied to the twin particles system, and that twin particles effects appear as disturbance of estimated trajectory.

In an image fiber containing a large number of cores, a certain class of crosstalk has been found to decrease with the distance along the fiber axis. This crosstalk is absolutely distinguished from the usual crosstalk that increases with the distance. A theoretical model is presented based on the power transfer between three groups of modes supported by each core. The process of power transfer is described by coupled power equations. Values of the coupling coefficients can be determined from the measurement of the crosstalk. The equations are solved numerically for the transmission of a point image. The results are in good agreement with measurement results.

In this paper, we developed the analytical method for the radiation field from a line current source placed in a stratified slab waveguide. This method is applicable to the analysis of excitation problem of inhomogeneous slab waveguide by increasing the number of layers. The numerical results are given for the cases of five layers, such as W and M type waveguides, and the inhomogeneous slab waveguide. The influence of guided and leaky modes on the radiation field are studied.

The propagation characteristics of dielectric waveguides with slanted grating structure are analyzed by using the combination of the improved Fourier series expansion method and the approximated multilayer method. The slanted grating region is appoximated by a structure with stratified thin modulated index layers. This method is effective to the guiding problems of the planar slanted grating, because the electromagnetic fields in each layer can be expressed by shifting the phase of the solution in the first layer. In this paper, numerical results are given for the grating with the rectangular and the sinusoidal profile for arbitrary slant angle. The radiation efficiencies for the grating with negative and positive slant angle are also discussed.

Numerical simulations for the (3+1)-dimensional optical-field dynamics of nonstationary pulsed beams that propagate down Kerr-like nonlinear channel waveguides are carried out for what is to our knowledge the first time. Time-resolved intrapulse switching due to spontaneous symmetry breaking of optical fields from a quasilinear symmetric field to a nonlinear asymmetric field is analyzed. A novel instability phenomenon triggered by the symmetry breakdown is found.

A simple numerical method for calculating paths of creeping rays around an arbitrary convex object is presented. The adventage of this method is that the path of creeping ray is iteratively determined from initial values of incident point and incident direction of the creeping ray without solving differential equation of geodesic path. As the numerical examples, the path of creeping ray on the prolate spheroid and the resonance path of natural modes are shown.

We present a novel algorithm to reconstruct the refractive-index profile of a circularly symmetric object from measurements of the electromagnetic field scattered when the object is illuminated by a plane wave. The reconstruction algorithm is besed on an iterative procedure of matching the scattered field calculated from a certain refractive-index distribution with the measured scattered field on the boundary of the object. In order to estimate the convergence of the reconstruction, the mean square error between the calculated and measured scattered fields is introduced. It is shown through reconstructing several examples of lossy dielectric cylinders that the algorithm is quite stable and is applicable to high-contrasty models in situations where the Born approximation is not valid.