A simulated annealing algorithm for finding consensus sequences

Motivation: A consensus sequence for a family of related sequences is, as the name suggests, a sequence that captures the features common to most members of the family. Consensus sequences are important in various DNA sequencing applications and are a convenient way to characterize a family of molecules. Results: This paper describes a new algorithm for finding a consensus sequence, using the popular optimization method known as simulated annealing. Unlike the conventional approach of finding a consensus sequence by first forming a multiple sequence alignment, this algorithm searches for a sequence that minimises the sum of pairwise distances to each of the input sequences. The resulting consensus sequence can then be used to induce a multiple sequence alignment. The time required by the algorithm scales linearly with the number of input sequences and quadratically with the length of the consensus sequence. We present results demonstrating the high quality of the consensus sequences and alignments produced by the new algorithm. For comparison, we also present similar results obtained using ClustalW. The new algorithm outperforms ClustalW in many cases.

Reasoning about real-time repetitions: terminating and nonterminating

It is common for a real-time system to contain a nonterminating process monitoring an input and controlling an output. Hence, a real-time program development method needs to support nonterminating repetitions. In this paper we develop a general proof rule for reasoning about possibly nonterminating repetitions. The rule makes use of a Floyd-Hoare-style loop invariant that is maintained by each iteration of the repetition, a Jones-style relation between the pre- and post-states on each iteration, and a deadline specifying an upper bound on the starting time of each iteration. The general rule is proved correct with respect to a predicative semantics. In the case of a terminating repetition the rule reduces to the standard rule extended to handle real time. Other special cases include repetitions whose bodies are guaranteed to terminate, nonterminating repetitions with the constant true as a guard, and repetitions whose termination is guaranteed by the inclusion of a fixed deadline. (C) 2002 Elsevier Science B.V. All rights reserved.

The fluorolysis assay, a highly sensitive method for measuring the cytolytic activity of T cells at very low numbers

We have developed a highly sensitive cytolysis test, the fluorolysis assay, as a simple nonradioactive and inexpensive alternative to the standard Cr-51-release assay. P815 cells were stably transfected with a plasmid expressing the enhanced green fluorescent protein (EGFP) gene. These target cells were coated with or without cognate peptide or anti-CD3 Ab and then incubated with CD8(+) T cells to allow antigen-specific or nonspecific lysis. The degree of target cell lysis was measured using flow cytometry to count the percentage of viable propidium iodide(-) EGFP(+) cells, whose numbers were standardized to a reference number of fluorochrome-linked beads. By using small numbers of target cells (200-800 per reaction) and extended incubation times (up to 2 days), the antigen-specific cytolytic activity of one to two activated CD8(+) T cells of a CTL line could be detected. The redirected fluorolysis assay also measured the activity of very few ( greater than or equal to6) primary CD8(+) T cells following polyclonal activation. Importantly, antigen-specific lysis by small numbers ( greater than or equal to 25) of primary CD8(+) T cells could be directly measured ex vivo. This exquisite sensitivity of the fluorolysis assay, which was at least 8-33-folds higher than an optimized 51 Cr-release assay, allows in vitro and ex vivo studies of immune responses that would otherwise not be possible due to low CTL numbers or frequencies. (C) 2002 Elsevier Science B.V. All rights reserved.

Marx as a Critical Discourse Analyst: The genesis of a critical method and its relevance to the critique of global capital

In this paper we identify elements in Marx´s economic and political writings that are relevant to contemporary critical discourse analysis (CDA). We argue that Marx can be seen to be e n gaging in a form of discourse analysis. We identify the elements in Marx´s historical materialist method that support such a perspective, and exemplify these in a longitudinal comparison of Marx´s texts.

Analysis of an array of monopoles with the use of a radial waveguide approach

A radial guide field matching method (RGFMM) is used to analyze a circular array antenna consisting of one active monopole surrounded by a concentric array of passive monopoles terminated in arbarary loads. An equivalent admittance matrix for this antenna system is determined to study the input admittance of the active monopole when the peripheral elements are terminated in open or short circuits. RGFMM results are compared with free-space method of moments (FS-MoM) results for a small switched-beam array a seven monopoles. Good agreement is noted. (C) 2002 Wiley Periodicals, Inc.

Analysis of a circular patch antenna radiating in a parallel-plate radial guide

A field matching method is described to analyze a recessed circular cavity radiating into a radial waveguide. Using the wall impedance approach, the analysis is divided into two separate problems of the cavity and its external environment. Based on this analysis, a computer algorithm is developed for determining wall admittances as seen at the edge of the patch in the cavity, the radial admittance matrix for the two-probe feed arrangement, and the input impedance as observed from the coaxial line feeding the cavity. This algorithm is tested against the general-purpose Hewlett-Packard finite-element High Frequency Structure Simulator as well as against measured results. Good agreement in all considered cases is noted.

Violation of multiparticle Bell inequalities for low- and high-flux parametric amplification using both vacuum and entangled input states

We show how polarization measurements on the output fields generated by parametric down conversion will reveal a violation of multiparticle Bell inequalities, in the regime of both low- and high-output intensity. In this case, each spatially separated system, upon which a measurement is performed, is comprised of more than one particle. In view of the formal analogy with spin systems, the proposal provides an opportunity to test the predictions of quantum mechanics for spatially separated higher spin states. Here the quantum behavior possible even where measurements are performed on systems of large quantum (particle) number may be demonstrated. Our proposal applies to both vacuum-state signal and idler inputs, and also to the quantum-injected parametric amplifier as studied by De Martini The effect of detector inefficiencies is included, and weaker Bell-Clauser-Horne inequalities are derived to enable realistic tests of local hidden variables with auxiliary assumptions for the multiparticle situation.

R-matrices and the tensor product graph method

A systematic method for constructing trigonometric R-matrices corresponding to the (multiplicity-free) tensor product of any two affinizable representations of a quantum algebra or superalgebra has been developed by the Brisbane group and its collaborators. This method has been referred to as the Tensor Product Graph Method. Here we describe applications of this method to untwisted and twisted quantum affine superalgebras.

Modelling of Lamb waves in composite laminated plates excited by interdigital transducers

The technique of permanently attaching interdigital transducers (IDT) to either flat or curved structural surfaces to excite single Lamb wave mode has demonstrated great potential for quantitative non-destructive evaluation and smart materials design, In this paper, the acoustic wave field in a composite laminated plate excited by an IDT is investigated. On the basis of discrete layer theory and a multiple integral transform method, an analytical-numerical approach is developed to evaluate the surface velocity response of the plate due to the IDTs excitation. In this approach, the frequency spectrum and wave number spectrum of the output of IDT are obtained directly. The corresponding time domain results are calculated by applying a standard inverse fast Fourier transformation technique. Numerical examples are presented to validate the developed method and show the ability of mode selection and isolation. A new effective way of transfer function estimation and interpretation is presented by considering the input wave number spectrum in addition to the commonly used input frequency spectrum. The new approach enables the simple physical evaluation of the influences of IDT geometrical features such as electrode finger widths and overall dimension and excitation signal properties on the input-output characteristics of IDT. Finally, considering the convenience of Mindlin plate wave theory in numerical computations as well as theoretical analysis, the validity is examined of using this approximate theory to design IDT for the excitation of the first and second anti-symmetric Lamb modes. (C) 2002 Elsevier Science Ltd. All rights reserved.

Control of chaotic instabilities in a spinning spacecraft with dissipation using Lyapunov's method

Control of chaotic instability in a simplified model of a spinning spacecraft with dissipation is achieved using an algorithm derived using Lyapunov's second method. The control method is implemented on a realistic spacecraft parameter configuration which has been found to exhibit chaotic instability for a range of forcing amplitudes and frequencies when a sinusoidally varying torque is applied to the spacecraft. Such a torque, may arise in practice from an unbalanced rotor or from vibrations in appendages. Numerical simulations are performed and the results are studied by means of time history, phase space, Poincare map, Lyapunov characteristic exponents and bifurcation diagrams. (C) 2002 Elsevier Science Ltd. All rights reserved.

A four-stage index 2 Diagonally Implicit Runge-Kutta method

High index Differential Algebraic Equations (DAEs) force standard numerical methods to lower order. Implicit Runge-Kutta methods such as RADAU5 handle high index problems but their fully implicit structure creates significant overhead costs for large problems. Singly Diagonally Implicit Runge-Kutta (SDIRK) methods offer lower costs for integration. This paper derives a four-stage, index 2 Explicit Singly Diagonally Implicit Runge-Kutta (ESDIRK) method. By introducing an explicit first stage, the method achieves second order stage calculations. After deriving and solving appropriate order conditions., numerical examples are used to test the proposed method using fixed and variable step size implementations. (C) 2001 IMACS. Published by Elsevier Science B.V. All rights reserved.

Single mode Lamb waves in composite laminated plates generated by piezoelectric transducers

The technique of permanently attaching piezoelectric transducers to structural surfaces has demonstrated great potential for quantitative non-destructive evaluation and smart materials design. For thin structural members such as composite laminated plates, it has been well recognized that guided Lamb wave techniques can provide a very sensitive and effective means for large area interrogation. However, since in these applications multiple wave modes are generally generated and the individual modes are usually dispersive, the received signals are very complex and difficult to interpret. An attractive way to deal with this problem has recently been introduced by applying piezoceramic transducer arrays or interdigital transducer (IDT) technologies. In this paper, the acoustic wave field in composite laminated plates excited by piezoceramic transducer arrays or IDT is investigated. Based on dynamic piezoelectricity theory, a discrete layer theory and a multiple integral transform method, an analytical-numerical approach is developed to evaluate the input impedance characteristics of the transducer and the surface velocity response of the plate. The method enables the quantitative evaluation of the influence of the electrical characteristics of the excitation circuit, the geometric and piezoelectric properties of the transducer array, and the mechanical and geometrical features of the laminate. Numerical results are presented to validate the developed method and show the ability of single wave mode selection and isolation. The results show that the interaction between individual elements of the piezoelectric array has a significant influence on the performance of the IDT, and these effects can not be neglected even in the case of low frequency excitation. It is also demonstrated that adding backing materials to the transducer elements can be used to improve the excitability of specific wave modes. (C) 2002 Elsevier Science Ltd. All rights reserved.