An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 1. Methodology for residual generation

For dynamic simulations to be credible, verification of the computer code must be an integral part of the modelling process. This two-part paper describes a novel approach to verification through program testing and debugging. In Part 1, a methodology is presented for detecting and isolating coding errors using back-to-back testing. Residuals are generated by comparing the output of two independent implementations, in response to identical inputs. The key feature of the methodology is that a specially modified observer is created using one of the implementations, so as to impose an error-dependent structure on these residuals. Each error can be associated with a fixed and known subspace, permitting errors to be isolated to specific equations in the code. It is shown that the geometric properties extend to multiple errors in either one of the two implementations. Copyright (C) 2003 John Wiley Sons, Ltd.

An approach to verifying and debugging simulation models governed by ordinary differential equations: Part 2. Residuals analysis and a case study

In Part 1 of this paper a methodology for back-to-back testing of simulation software was described. Residuals with error-dependent geometric properties were generated. A set of potential coding errors was enumerated, along with a corresponding set of feature matrices, which describe the geometric properties imposed on the residuals by each of the errors. In this part of the paper, an algorithm is developed to isolate the coding errors present by analysing the residuals. A set of errors is isolated when the subspace spanned by their combined feature matrices corresponds to that of the residuals. Individual feature matrices are compared to the residuals and classified as 'definite', 'possible' or 'impossible'. The status of 'possible' errors is resolved using a dynamic subset testing algorithm. To demonstrate and validate the testing methodology presented in Part 1 and the isolation algorithm presented in Part 2, a case study is presented using a model for biological wastewater treatment. Both single and simultaneous errors that are deliberately introduced into the simulation code are correctly detected and isolated. Copyright (C) 2003 John Wiley Sons, Ltd.

Differential Priors for Elastic Nets

The elastic net and related algorithms, such as generative topographic mapping, are key methods for discretized dimension-reduction problems. At their heart are priors that specify the expected topological and geometric properties of the maps. However, up to now, only a very small subset of possible priors has been considered. Here we study a much more general family originating from discrete, high-order derivative operators. We show theoretically that the form of the discrete approximation to the derivative used has a crucial influence on the resulting map. Using a new and more powerful iterative elastic net algorithm, we confirm these results empirically, and illustrate how different priors affect the form of simulated ocular dominance columns.

Gauge fields, geometric phases, and quantum adiabatic pumps

Quantum adiabatic pumping of charge and spin between two reservoirs (leads) has recently been demonstrated in nanoscale electronic devices. Pumping occurs when system parameters are varied in a cyclic manner and sufficiently slowly that the quantum system always remains in its ground state. We show that quantum pumping has a natural geometric representation in terms of gauge fields (both Abelian and non-Abelian) defined on the space of system parameters. Tunneling from a scanning tunneling microscope tip through a magnetic atom could be used to demonstrate the non-Abelian character of the gauge field.

Geometric and electronic information from the spectroscopy of six-coordinate copper(II) compounds

The ground and excited state geometry of the six-coordinate copper(II) ion is examined in detail using the CuF64- and Cu(H2O)(6)(2+) complexes as examples. A variety of spectroscopic techniques are used to illustrate the relations between the geometric and electronic properties of these complexes through the characterization of their potential energy surfaces.

Semi-analytical solution for nonlinear vibration of laminated FGM plates with geometric imperfections

This paper investigates the nonlinear vibration of imperfect shear deformable laminated rectangular plates comprising a homogeneous substrate and two layers of functionally graded materials (FGMs). A theoretical formulation based on Reddy's higher-order shear deformation plate theory is presented in terms of deflection, mid-plane rotations, and the stress function. A semi-analytical method, which makes use of the one-dimensional differential quadrature method, the Galerkin technique, and an iteration process, is used to obtain the vibration frequencies for plates with various boundary conditions. Material properties are assumed to be temperature-dependent. Special attention is given to the effects of sine type imperfection, localized imperfection, and global imperfection on linear and nonlinear vibration behavior. Numerical results are presented in both dimensionless tabular and graphical forms for laminated plates with graded silicon nitride/stainless steel layers. It is shown that the vibration frequencies are very much dependent on the vibration amplitude and the imperfection mode and its magnitude. While most of the imperfect laminated plates show the well-known hard-spring vibration, those with free edges can display soft-spring vibration behavior at certain imperfection levels. The influences of material composition, temperature-dependence of material properties and side-to-thickness ratio are also discussed. (C) 2004 Elsevier Ltd. All rights reserved.

Measuring geometric phases of scattering states in nanoscale electronic devices

We show how a quantum property, a geometric phase, associated with scattering states can be exhibited in nanoscale electronic devices. We propose an experiment to use interference to directly measure the effect of this geometric phase. The setup involves a double-path interferometer, adapted from that used to measure the phase evolution of electrons as they traverse a quantum dot (QD). Gate voltages on the QD could be varied cyclically and adiabatically, in a manner similar to that used to observe quantum adiabatic charge pumping. The interference due to the geometric phase results in oscillations in the current collected in the drain when a small bias across the device is applied. We illustrate the effect with examples of geometric phases resulting from both Abelian and non-Abelian gauge potentials.

Geometric phases of scattering states in a ring geometry: adiabatic pumping in mesoscopic devices

Geometric phases of scattering states in a ring geometry are studied on the basis of a variant of the adiabatic theorem. Three timescales, i.e., the adiabatic period, the system time and the dwell time, associated with adiabatic scattering in a ring geometry play a crucial role in determining geometric phases, in contrast to only two timescales, i.e., the adiabatic period and the dwell time, in an open system. We derive a formula connecting the gauge invariant geometric phases acquired by time-reversed scattering states and the circulating (pumping) current. A numerical calculation shows that the effect of the geometric phases is observable in a nanoscale electronic device.

The importance of wet-powder dynamic mechanical properties in understanding granulation

Granule impact deformation has long been recognised as important in determining whether or not two colliding granules will coalesce. Work in the last 10 years has highlighted the fact that viscous effects are significant in granulation. The relative strengths of different formulations can vary with strain rate. Therefore, traditional strength measurements made at pseudo-static conditions give no indication, even qualitatively, of how materials will behave at high strain rates, and hence are actually misleading when used to model granule coalescence. This means that new standard methods need to be developed for determining the strain rates encountered by granules inside industrial equipment and also for measuring the mechanical properties of granules at these strain rates. The constitutive equations used in theoretical models of granule coalescence also need to be extended to include strain-rate dependent components.

Spectral properties of the tandem Jackson network, seen as a quasi-birth-and-death process

Quasi-birth-and-death (QBD) processes with infinite “phase spaces” can exhibit unusual and interesting behavior. One of the simplest examples of such a process is the two-node tandem Jackson network, with the “phase” giving the state of the first queue and the “level” giving the state of the second queue. In this paper, we undertake an extensive analysis of the properties of this QBD. In particular, we investigate the spectral properties of Neuts’s R-matrix and show that the decay rate of the stationary distribution of the “level” process is not always equal to the convergence norm of R. In fact, we show that we can obtain any decay rate from a certain range by controlling only the transition structure at level zero, which is independent of R. We also consider the sequence of tandem queues that is constructed by restricting the waiting room of the first queue to some finite capacity, and then allowing this capacity to increase to infinity. We show that the decay rates for the finite truncations converge to a value, which is not necessarily the decay rate in the infinite waiting room case. Finally, we show that the probability that the process hits level n before level 0 given that it starts in level 1 decays at a rate which is not necessarily the same as the decay rate for the stationary distribution.

A protein with protective properties against the cellular defense reactions in insects

The molecular mechanism of how insects recognize intruding microorganisms and parasites and distinguish them from own body structures is not well known. We explored evolutionary adaptations in an insect parasitoid host interaction to identify components that interfere with the recognition of foreign objects and cellular encapsulation. Because some parasitoids provide protection for the developing wasp in the absence of an overt suppression of the insect host defense, we analyzed the surface of eggs and symbiotic viruses for protective properties. Here we report on the molecular cloning of a 32-kDa protein (Crp32) that is one of the major protective components. It is produced in the calyx cells of the female wasp ovaries and attached to the surface of the egg and other particles including polydnaviruses. The recombinant protein confers protection to coated objects in a cellular encapsulation assay suggesting that a layer of Crp32 may prevent cellular encapsulation reactions by a local inactivation of the host defense system.

Ester prodrugs of a potent analgesic, morphine-6-sulfate: syntheses, spectroscopic and physicochemical properties

The aim of this work is to develop 3-acyl prodrugs of the potent analgesic morphine-6-sulfate (M6S). These are expected to have higher potency and/or exhibit longer duration of analgesic action than the parent compound. M6S and the prodrugs were synthesized, then purified either by recrystallization or by semi-preparative HPLC and the structures confirmed by mass spectrometry, IR spectrophotometry and by detailed 1- and 2-D NMR studies. The lipophilicities of the compounds were assessed by a combination of shake-flask, group contribution and HPLC retention methods. The octanol-buffer partition coefficient could only be obtained directly for 3-heptanoylmorphine-6-sulfate, using the shake-flask method. The partition coefficients (P) for the remaining prodrugs were estimated from known methylene group contributions. A good linear relationship between log P and the HPLC log capacity factors was demonstrated. Hydrolysis of the 3-acetyl prodrug, as a representative of the group, was found to occur relatively slowly in buffers (pH range 6.15-8.01), with a small buffer catalysis contribution. The rates of enzymatic hydrolysis of the 3-acyl group in 10% rat blood and in 10% rat brain homogenate were investigated. The prodrugs followed apparent first order hydrolysis kinetics, with a significantly faster hydrolysis rate found in 10% rat brain homogenate than in 10% rat blood for all compounds. (C) 1998 Elsevier Science B.V. All rights reserved.

Effects of acidic treatments on the pore and surface properties of Ni catalyst supported on activated carbon

Activated carbon as catalyst support was treated with HCl, HNO3, and HF and the effects of acid treatments on the properties of the activated carbon support were studied by N-2 adsorption, mass titration, temperature-programmed desorption (TPD), and X-ray photoelectron spectrometry (XPS). Ni catalysts supported on untreated and treated activated carbons were prepared, characterized and tested for the reforming reaction of methane with carbon dioxide. It is found that acid treatment significantly changed the surface chemical properties and pore structure of the activated carbon. The surface area and pore volume of the carbon supports are generally enhanced upon acid treatment due to the removal of impurities present in the carbon. The adsorption capacity of Ni2+ on the carbon supports is also increased, and the increase can be closely correlated with the surface acidity. The impregnation of nickel salts decreases the surface area and pore volume of carbon supports both in micropores and mesopores. Acid treatment results in a more homogeneous distribution of the nickel salt in carbon. When the impregnated carbons are heated in inert atmosphere, there exists a redox reaction between nickel oxide and the carbon. Catalytic activity tests for methane reforming with carbon dioxide show that the activity of nickel catalysts based on the acid-treated carbon supports is closely related with the surface characteristics of catalysts. (C) 1998 Elsevier Science Ltd. All rights reserved.

Optical properties of metallic films for vertical-cavity optoelectronic devices

We present models for the optical functions of 11 metals used as mirrors and contacts in optoelectronic and optical devices: noble metals (Ag, Au, Cu), aluminum, beryllium, and transition metals (Cr, Ni, Pd, Pt, Ti, W). We used two simple phenomenological models, the Lorentz-Drude (LD) and the Brendel-Bormann (BB), to interpret both the free-electron and the interband parts of the dielectric response of metals in a wide spectral range from 0.1 to 6 eV. Our results show that the BE model was needed to describe appropriately the interband absorption in noble metals, while for Al, Be, and the transition metals both models exhibit good agreement with the experimental data. A comparison with measurements on surface normal structures confirmed that the reflectance and the phase change on reflection from semiconductor-metal interfaces (including the case of metallic multilayers) can be accurately described by use of the proposed models for the optical functions of metallic films and the matrix method for multilayer calculations. (C) 1998 Optical Society of America.