Rational Function Interpolation of Electromagnetic Transfer Functions of High-Speed Interconnect Systems from Discrete Time-Domain and Frequency-Domain Data

We present new methodologies to generate rational function approximations of broadband electromagnetic responses of linear and passive networks of high-speed interconnects, and to construct SPICE-compatible, equivalent circuit representations of the generated rational functions. These new methodologies are driven by the desire to improve the computational efficiency of the rational function fitting process, and to ensure enhanced accuracy of the generated rational function interpolation and its equivalent circuit representation. Toward this goal, we propose two new methodologies for rational function approximation of high-speed interconnect network responses. The first one relies on the use of both time-domain and frequency-domain data, obtained either through measurement or numerical simulation, to generate a rational function representation that extrapolates the input, early-time transient response data to late-time response while at the same time providing a means to both interpolate and extrapolate the used frequency-domain data. The aforementioned hybrid methodology can be considered as a generalization of the frequency-domain rational function fitting utilizing frequency-domain response data only, and the time-domain rational function fitting utilizing transient response data only. In this context, a guideline is proposed for estimating the order of the rational function approximation from transient data. The availability of such an estimate expedites the time-domain rational function fitting process. The second approach relies on the extraction of the delay associated with causal electromagnetic responses of interconnect systems to provide for a more stable rational function process utilizing a lower-order rational function interpolation. A distinctive feature of the proposed methodology is its utilization of scattering parameters. For both methodologies, the approach of fitting the electromagnetic network matrix one element at a time is applied. It is shown that, with regard to the computational cost of the rational function fitting process, such an element-by-element rational function fitting is more advantageous than full matrix fitting for systems with a large number of ports. Despite the disadvantage that different sets of poles are used in the rational function of different elements in the network matrix, such an approach provides for improved accuracy in the fitting of network matrices of systems characterized by both strongly coupled and weakly coupled ports. Finally, in order to provide a means for enforcing passivity in the adopted element-by-element rational function fitting approach, the methodology for passivity enforcement via quadratic programming is modified appropriately for this purpose and demonstrated in the context of element-by-element rational function fitting of the admittance matrix of an electromagnetic multiport.

Design, construction, and applications of a high-resolution terahertz time-domain spectrometer

This thesis reports on the design, construction, and initial applications of a high-resolution terahertz time-domain ASOPS spectrometer. The instrument employs asynchronous optical sampling (ASOPS) between two Ti:sapphire ultrafast lasers operating at a repetition rate of approximately 80 MHz, and we thus demonstrate a THz frequency resolution approaching the limit of that repetition rate. This is an order of magnitude improvement in resolution over typical THz time-domain spectrometers. The improved resolution is important for our primary effort of collecting THz spectra for far-infrared astronomy. We report on various spectroscopic applications including the THz rotational spectrum of water, where we achieve a mean frequency error, relative to established line centers, of 27.0 MHz. We also demonstrate application of the THz system to the long-duration observation of a coherent magnon mode in a anti-ferromagnetic yttrium iron oxide (YFeO3) crystal. Furthermore, we apply the all-optical virtual delay line of ASOPS to a transient thermoreflectance experiment for quickly measuring the thermal conductivity of semiconductors.

Exploration of Ground Penetrating Radar and Time Domain Reflectometry Methods for the Determination of Pavement Dielectric Constant

Traditionally, densities of newly built roadways are checked by direct sampling (cores) or by nuclear density gauge measurements. For roadway engineers, density of asphalt pavement surfaces is essential to determine pavement quality. Unfortunately, field measurements of density by direct sampling or by nuclear measurement are slow processes. Therefore, I have explored the use of rapidly-deployed ground penetrating radar (GPR) as an alternative means of determining pavement quality. The dielectric constant of pavement surface may be a substructure parameter that correlates with pavement density, and can be used as a proxy when density of asphalt is not known from nuclear or destructive methods. The dielectric constant of the asphalt can be determined using ground penetrating radar (GPR). In order to use GPR for evaluation of road surface quality, the relationship between dielectric constants of asphalt and their densities must be established. Field measurements of GPR were taken at four highway sites in Houghton and Keweenaw Counties, Michigan, where density values were also obtained using nuclear methods in the field. Laboratory studies involved asphalt samples taken from the field sites and samples created in the laboratory. These were tested in various ways, including, density, thickness, and time domain reflectometry (TDR). In the field, GPR data was acquired using a 1000 MHz air-launched unit and a ground-coupled unit at 200 and 500 MHz. The equipment used was owned and operated by the Michigan Department of Transportation (MDOT) and available for this study for a total of four days during summer 2005 and spring 2006. The analysis of the reflected waveforms included “routine” processing for velocity using commercial software and direct evaluation of reflection coefficients to determine a dielectric constant. The dielectric constants computed from velocities do not agree well with those obtained from reflection coefficients. Perhaps due to the limited range of asphalt types studied, no correlation between density and dielectric constant was evident. Laboratory measurements were taken with samples removed from the field and samples created for this study. Samples from the field were studied using TDR, in order to obtain dielectric constant directly, and these correlated well with the estimates made from reflection coefficients. Samples created in the laboratory were measured using 1000 MHz air-launched GPR, and 400 MHz ground-coupled GPR, each under both wet and dry conditions. On the basis of these observations, I conclude that dielectric constant of asphalt can be reliably measured from waveform amplitude analysis of GJPR data, based on the consistent agreement with that obtained in the laboratory using TDR. Because of the uniformity of asphalts studied here, any correlation between dielectric constant and density is not yet apparent.

A modified dual-band microstrip monopole antenna

A modified microstrip-fed planar monopole antenna with open circuited coupled line is presented in this paper. The operational bandwidth of the proposed antenna covers the 2.4 GHz ISM band (2.42-2.48 GHz) and the 5 GHz WLAN band (5 GHz to 6 GHz). The radiating elements occupy a small area of 23×8 mm2. The Finite Difference Time Domain method is used to predict the input impedance of the antenna. The calculated return loss shows very good agreement with measured data. Reasonable antenna gain is observed across the operating band. The measured radiation patterns are similar to those of a simple monopole antenna.

Analysis of ECG by multipulse excited linear prediction model

He propose a new time domain method for efficient representation of the KCG and delineation of its component waves. The method is based on the multipulse Linear prediction (LP) coding which is being widely used in speech processing. The excitation to the LP synthesis filter consists of a few pulses defined by their locations and amplitudes. Based on the amplitudes and their distribution, the pulses are suitably combined to delineate the component waves. Beat to beat correlation in the ECG signal is used in QRS periodicity prediction. The method entails a data compression of 1 in 6. The method reconstructs the signal with an NMSE of less than 5%.

Propagation of an arbitrary elliptically polarized few-cycle ultrashort laser pulse in resonant two-level quantum systems

We investigate the propagation of an arbitrary elliptically polarized few-cycle ultrashort laser pulse in resonant two-level quantum systems using an iterative predictor-corrector finite-difference time-domain method. It is shown that when the initial effective area is equal to 2 pi, the effective area will remain invariant during the course of propagation, and a complete Rabi oscillation can be achieved. However, for an elliptically polarized few-cycle ultrashort laser pulse, polarization conversion can occur. Eventually, the laser pulse will evolve into two separate circularly polarized laser pulses with opposite helicities.

Propagation of a few-cycle laser pulse in a V-type three-level system

Propagation of a few-cycle laser pulse in a V-type three-level system (fine structure levels of rubidium) is investigated numerically. The full three-level Maxwell-Bloch equations without the rotating wave approximation and the standing slowly varying envelope approximation are solved by using a finite-difference time-domain method. It is shown that, when the usual unequal oscillator strengths are considered, self-induced transparency cannot be recovered and higher spectral components can be produced even for small-area pulses. (c) 2005 Pleiades Publishing, Inc.

Polarization evolution of a few-cycle ultrashort laser pulse propagating in a degenerate three-level medium

The propagation of an arbitrary polarized few-cycle ultrashort laser pulse in a degenerate three-level medium is investigated by using an iterative predictor-corrector finite-difference time-domain method. It is found that the polarization evolution of the ultrashort laser pulse is dependent not only on the initial atomic coherence of the medium but also on the polarization condition of the incident laser pulse. When the initial effective area is equal to 2 pi, complete linear-to-circular and circular-to-linear polarization conversion of few-cycle ultrashort laser pulses can be achieved due to the quantum interference effects between the two different transition paths.

Polarization-dependent Talbot effect

The term "polarization-dependent Talbot effect" means that the Talbot self-imaging intensity of a high-density grating is different for TE and TM polarization modes. Numerical simulations with the finite-difference time-domain method show that the polarization dependence of the Talbot images is obvious for gratings with period d between 2 lambda and 3 lambda. Such a polarization-dependent difference for TE and TM polarization of, a high-density grating of 630 lines/mm (corresponding to d/lambda = 2.5) is verified through experiments with the scanning near-field optical microscopy technique, in which a He-Ne laser is used as its polarization is changed from the TE mode to the TM mode. The polarization-dependent Talbot effect should help us to understand more clearly the diffraction behavior of a high-density grating in nano-optics and contribute to wide application of the Talbot effect. (c) 2006 Optical Society of America.

Talbot effect of a grating with different kinds of flaws

The Talbot effect of a grating with different kinds of flaws is analyzed with the finite-difference time-domain (FDTD) method. The FDTD method can show the exact near-field distribution of different flaws in a high-density grating, which is impossible to obtain with the conventional Fourier transform method. The numerical results indicate that if a grating is perfect, its Talbot imaging should also be perfect; if the grating is distorted, its Talbot imaging will also be distorted. Furthermore, we evaluate high-density gratings by detecting the near-field distribution with the scanning near-field optical microscopy technique. Experimental results are also given. (c) 2005 Optical Society of America.

Phase interpretation for polarization-dependent near-field images of high-density gratings

It has been described that the near-field images of a high-density grating at the half self-imaging distance could be different for TE and TM polarization states. We propose that the phases of the diffraction orders play an important role in such polarization dependence. The view is verified through the coincidence of the numerical result of finite-difference time-domain method and the reconstructed results from the rigorous coupled-wave analysis. Field distributions of TE and TM polarizations are given numerically for a grating with period d = 2.3 lambda, which are verified through experiments with the scanning near-field optical microscopy technique. The concept of phase interpretation not only explains the polarization dependence at the half self-imaging distance of gratings with a physical view, but also, it could be widely used to describe the near-field diffraction of a variety of periodic diffractive optical elements whose feature size comparable to the wavelength. (C) 2008 Elsevier B.V. All rights reserved.

Optical near-field simulation of Sb thin film thermal lens and its application in optical recording

Using the finite-difference-time-domain method, the near-field optical distribution and properties of Sb thin film thermal lens are calculated and simulated. The results show as follows. Within the near-field distance to the output plane of thermal lens, the spot size is approximately 100 nm, and its intensity is greatly enhanced, which is higher than that of incident light. The spot shape gradually changes from ellipse to round at the distance of more than 12 nm to the output plane. The above-simulated results are further demonstrated by the static optical recording experiment. (C) 2005 American Institute of Physics.

数值模拟探针诱导表面等离子体共振耦合纳米光刻

采用有损耗介质和色散介质的二维时域有限差分方法,数值模拟了以光波长514.5 nm的p偏振基模高斯光束为入射光源,激发Kretschmann型表面等离子体共振,并通过探针的局域场增强效应实现纳米光刻的新方法——探针诱导表面等离子体共振耦合纳米光刻.分别就探针与记录层的间距以及探针针尖大小,模拟分析了不同情况下探针的局域场增强效应和记录层表面的相对电场强度振幅分布.结果表明,探针工作在接触模式时,探针的局域场增强效应最明显,记录层表面的相对电场强度振幅的对比度最大;当探针针尖距记录层5 nm时,针尖下方记录层表面的相对电场强度振幅大于光刻临界值的分布宽度与针尖尺寸相近.

Design of a compact polarizing beam splitter based on a photonic crystal ring resonator with a triangular lattice

We propose and simulate a new kind of compact polarizing beam splitter (PBS) based on a photonic crystal ring resonator (PCRR) with complete photonic bandgaps. The two polarized states are separated far enough by resonant and nonresonant coupling between the waveguide modes and the microring modes. Some defect holes are utilized to control the beam propagation. The simulated results obtained by the finite-difference time-domain method show that high transmission (over 95%) is obtained and the polarization separation is realized with a length as short as 3.1 mu m. The design of the proposed PBS can be flexible, thanks to the advantages of PCRRs.

Proposal of an Ultracompact Triplexer Using Photonic Crystal Waveguide with an Air Holes Array

We propose an ultracompact triplexer based on a shift of the cutoff frequency of the fundamental mode in a planar photonic crystal waveguide (PCW) with a triangular lattice of air holes. The shift is realized by modifying the radii of the border holes adjacent to the PCW core. Some defect holes are introduced to control the beam propagation. The numerical results obtained by the finite-difference time-domain method show that the presented triplexer can separate three specific wavelengths, i.e. 1310, 1490 and 1550 nm with the extinction ratios higher than - 18 dB. The designed device with a size as compact as 12 mu m x 6.5 mu m is feasible for the practical application, and can be utilized in the system of fiber to the home.