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.

Measurements of small ice crystals in arctic cirrus: observations from the Indirect and Semi-Direct Aerosol Campaign (ISDAC)

In-situ observations on the size and shape of particles in arctic cirrus are less common than those in mid-latitude and tropical cirrus with considerable uncertainty about the contributions of small ice crystals (maximum dimension D<50 µm) to the mass and radiative properties that impact radiative forcing. In situ measurements of small ice crystals in arctic cirrus were made during the Indirect and Semi-Direct Aerosol Campaign (ISDAC) in April 2008 during transits of the National Research Council of Canada Convair-580 between Fairbanks and Barrow, Alaska and during Mixed Phase Arctic Cloud Experiment (MPACE) in October 2004 with the University of North Dakota (UND) Citation over Barrow, Alaska. Concentrations of small ice crystals with D < 50 μm from a Cloud and Aerosol Spectrometer (CAS), a Cloud Droplet Probe (CDP), a Forward Scattering Spectrometer Probe (FSSP), and a two-dimensional stereo probe (2DS) were compared as functions of the concentrations of crystals with D > 100 μm measured by a Cloud Imaging Probe (CIP) and two-dimensional stereo probe (2DS) in order to assess whether the shattering of large ice crystals on protruding components of different probes artificially amplified measurements of small ice crystal concentrations. The dependence of the probe comparison on other variables as CIP N>100 (number concentrations greater than diameter D>100 μm),temperature, relative humidity respect to ice (RHice), dominant habit from the Cloud Particle Imager (CPI), aircraft roll, pitch, true air speed and angle of attack was examined to understand potential causes of discrepancies between probe concentrations. Data collected by these probes were also compared against the data collected by a CAS, CDP and CIP during the Tropical Warm Pool-International Cloud Experiment (TWP-ICE) and by a CAS and 2DS during the Tropical Composition, Cloud and Climate Coupling (TC4) missions. During ISDAC, the CAS and FSSP both overestimated measurements of small ice crystals compared to both the CDP and 2DS by 1-2 orders of magnitude. Further, the amount of overestimation increased with the concentrations from the CIP2 (N>100 > 0.1 L-1). There was an unexplained discrepancy in concentrations of small crystals between the CDP and 2DS during ISDAC. In addition, there was a strong dependence on RHice of the average ratios of the N3-50, CAS/N3-50,CDP, N3-50, FSSP096/N3-50,CDP, N3-50, CAS/N3-50,FSSP096, N10-50, CDP/N3-50,2DS, N10-50, FSSP096/N10-50,2DS. Continued studies are needed to understand the discrepancy of these probes.