A monolithic smoothing-gap algorithm for contact-impact based on the signed distance function

A new algorithm has been developed for smoothing the surfaces in finite element formulations of contact-impact. A key feature of this method is that the smoothing is done implicitly by constructing smooth signed distance functions for the bodies. These functions are then employed for the computation of the gap and other variables needed for implementation of contact-impact. The smoothed signed distance functions are constructed by a moving least-squares approximation with a polynomial basis. Results show that when nodes are placed on a surface, the surface can be reproduced with an error of about one per cent or less with either a quadratic or a linear basis. With a quadratic basis, the method exactly reproduces a circle or a sphere even for coarse meshes. Results are presented for contact problems involving the contact of circular bodies. Copyright (C) 2002 John Wiley Sons, Ltd.

Molecular phylogenetics of Lentibulariaceae inferred from plastid rps16 Intron and trnL-F DNA sequences: implications for character evolution and biogeography

Phylogenetic relationships among 75 species of Lentibulariaceae, representing the three recognized genera, were assessed by cladistic analysis of DNA sequences from the plastid rps16 intron and the trnL-F region. Sequence data from the two loci were analyzed both separately and in combination. Consensus trees from all analyses are congruent, and parsimony jackknife results demonstrate strong support for relationships both between and within each of the three demonstrably monophyletic genera. The genus Pinguicula is sister to a Genlisea-Utricularia clade, the phylogenetic structure within this clade closely follows Taylor's recent sectional delimitations based on morphology. Three principal clades are shown within Utricularia, with the basal sections Polypoinpholyx and Pleiochasia together forming the sister lineage of the remaining Utricularia species. Of the fundamental morphological specializations, the stoloniferous growth form apparently arose independently within Genlisea and Utricularia three times, and within Utricularia itself, perhaps more than once. The epiphytic habit has evolved independently at least three times, in Pinguicula, in Utricularia section Phyllaria, and within the two sections Orchidioides and Iperua (in the latter as bromeliad tank-epiphytes). The suspended aquatic habit may have evolved independently within sections Utricularia and Vesiculina. Biogeographic optimization on the phylogeny demonstrates patterns commonly associated with the boreotropics hypothesis and limits the spatial origin of Lentibulariaceae to temperate Eurasia or tropical America.

Universal quantum computation and simulation using any entangling Hamiltonian and local unitaries

What interactions are sufficient to simulate arbitrary quantum dynamics in a composite quantum system? We provide an efficient algorithm to simulate any desired two-body Hamiltonian evolution using any fixed two-body entangling n-qubit Hamiltonian and local unitary operations. It follows that universal quantum computation can be performed using any entangling interaction and local unitary operations.

Universal simulation of Hamiltonian dynamics for quantum systems with finite-dimensional state spaces

What interactions are sufficient to simulate arbitrary quantum dynamics in a composite quantum system? Dodd [Phys. Rev. A 65, 040301(R) (2002)] provided a partial solution to this problem in the form of an efficient algorithm to simulate any desired two-body Hamiltonian evolution using any fixed two-body entangling N-qubit Hamiltonian, and local unitaries. We extend this result to the case where the component systems are qudits, that is, have D dimensions. As a consequence we explain how universal quantum computation can be performed with any fixed two-body entangling N-qudit Hamiltonian, and local unitaries.

Evolution of the secondary structure of the rRNA internal transcribed spacer 2 (ITS2) in hard ticks (Ixodidae, Arthropoda)

ITS2 sequences are used extensively in molecular taxonomy and population genetics of arthropods and other animals yet little is known about the molecular evolution of ITS2. We studied the secondary structure of ITS2 in species from each of the six main lineages of hard ticks (family Ixodidae). The ITS2 of these ticks varied in length from 679 bp in Ixodes scapularis to 1547 bp in Aponomma concolor. Nucleotide content varied also: the ITS2 of ticks from the Prostriata lineage (Ixodes spp.) had 46-49% GC whereas ITS2 sequences of ticks from the Metastriata lineage (all other hard ticks) had 61-62% GC. Despite variation in nucleotide sequence, the secondary structure of the ITS2 of all of these ticks apparently has five domains. Stems 1, 3, 4 and 5 of this secondary structure were obvious in all of the species studied. However, stem 2 was not always obvious despite the fact that it is flanked by highly conserved sequence motifs in the adjacent stems, stems 1 and 3. The ITS2 of hard ticks has apparently evolved mostly by increases and decreases in length of the nucleotide sequences, which caused increases, and decreases in the length of stems of the secondary structure. This is most obvious when stems of the secondary structures of the Prostriata (Ixodes spp.) are compared to those of the Metastriata (all other hard ticks). Increases in the size of the ITS2 may have been caused by replication slippage which generated large repeats, like those seen in Haemaphysalis humerosa and species from the Rhipicepalinae lineage, and the small repeats found in species from the other lineages of ticks.

Analysis of the wind field and heat budget in an alpine lake basin during summertime fair weather conditions

Observational data collected in the Lake Tekapo hydro catchment of the Southern Alps in New Zealand are used to analyse the wind and temperature fields in the alpine lake basin during summertime fair weather conditions. Measurements from surface stations, pilot balloon and tethersonde soundings, Doppler sodar and an instrumented light aircraft provide evidence of multi-scale interacting wind systems, ranging from microscale slope winds to mesoscale coast-to-basin flows. Thermal forcing of the winds occurred due to differential heating as a consequence of orography and heterogeneous surface features, which is quantified by heat budget and pressure field analysis. The daytime vertical temperature structure was characterised by distinct layering. Features of particular interest are the formation of thermal internal boundary layers due to the lake-land discontinuity and the development of elevated mixed layers. The latter were generated by advective heating from the basin and valley sidewalls by slope winds and by a superimposed valley wind blowing from the basin over Lake Tekapo and up the tributary Godley Valley. Daytime heating in the basin and its tributary valleys caused the development of a strong horizontal temperature gradient between the basin atmosphere and that over the surrounding landscape, and hence the development of a mesoscale heat low over the basin. After noon, air from outside the basin started flowing over mountain saddles into the basin causing cooling in the lowest layers, whereas at ridge top height the horizontal air temperature gradient between inside and outside the basin continued to increase. In the early evening, a more massive intrusion of cold air caused rapid cooling and a transition to a rather uniform slightly stable stratification up to about 2000 m agl. The onset time of this rapid cooling varied about 1-2 h between observation sites and was probably triggered by the decay of up-slope winds inside the basin, which previously countered the intrusion of air over the surrounding ridges. The intrusion of air from outside the basin continued until about mid-night, when a northerly mountain wind from the Godley Valley became dominant. The results illustrate the extreme complexity that can be caused by the operation of thermal forcing processes at a wide range of spatial scales.

Polarised infrared and differential scanning calorimetry studies on oriented vinyl pipe materials

The molecular orientation in a conventionally extruded PVC pipe, a uniaxially oriented PVC pipe and a biaxially oriented PVC pipe has been studied via Infrared dichroism. The degree of order or crystallinity has also been studied by Differential Scanning Calorimetry and also via Infrared Spectroscopy. The fundamental structural difference between the conventional and oriented pipes was that polymer chains were preferentially aligning in the hoop direction for oriented pipes whereas they were fairly isotropic in the conventional pipe with a slight preferential alignment in the axial direction. Analysis of the C-Cl stretching mode indicated that the uniaxially oriented pipe had much higher alignment of the C-Cl bond in the axial direction than the biaxial pipe, which correlates with higher fracture toughness for circumferential cracking in the biaxial pipe. Both DSC and Infrared spectroscopy detected little change in the crystallinity or order in the oriented pipes compared to the conventionally extruded pipes. (C) 2002 Kluwer Academic Publishers.

Differential regulation of twitching motility and elastase production by Vfr in Pseudomonas aeruginosa

Vfr, a homolog of Escherichia coli cyclic AMP (cAMP) receptor protein, has been shown to regulate quorum sensing, exotoxin A production, and regA transcription in Pseudomonas aeruginosa. We identified a twitching motility-defective mutant that carries a transposon insertion in vfr and confirmed that vfr is required for twitching motility by construction of an independent allelic deletion-replacement mutant of vfr that exhibited the same phenotype, as well as by the restoration of normal twitching motility by complementation of these mutants with wild-type vfr. Vfr-null mutants exhibited severely reduced twitching motility with barely detectable levels of type IV pili, as well as loss of elastase production and altered pyocyanin production. We also identified reduced-twitching variants of quorum-sensing mutants (PAK lasl::Tc) with a spontaneous deletion in vfr (S. A. Beatson, C. B. Whitchurch, A. B. T. Semmler, and J. S. Mattick, J. Bacteriol., 184:3598-3604,2002), the net result of which was the loss of five residues (EQERS) from the putative cAMP-binding pocket or Vfr. This allele (VfrDeltaEQERS) was capable of restoring elastase and pyocyanin production to wild-type levels in vfr-null mutants but not their defects in twitching motility. Furthermore, structural analysis of Vfr and VfrDeltaEQERS in relation to E. coli CRP suggests that Vfr is capable of binding both cAMP and cyclic GMP whereas VfrDeltaEQERS is only capable of responding to cAMP. We suggest that Vfr controls twitching motility and quorum sensing via independent pathways in response to these different signals, bound by the same cyclic nucleotide monophosphate-binding pocket.

Predictor-corrector methods of Runge-Kutta type for stochastic differential equations

In this paper we construct predictor-corrector (PC) methods based on the trivial predictor and stochastic implicit Runge-Kutta (RK) correctors for solving stochastic differential equations. Using the colored rooted tree theory and stochastic B-series, the order condition theorem is derived for constructing stochastic RK methods based on PC implementations. We also present detailed order conditions of the PC methods using stochastic implicit RK correctors with strong global order 1.0 and 1.5. A two-stage implicit RK method with strong global order 1.0 and a four-stage implicit RK method with strong global order 1.5 used as the correctors are constructed in this paper. The mean-square stability properties and numerical results of the PC methods based on these two implicit RK correctors are reported.

A variable stepsize implementation for stochastic differential equations

Stochastic differential equations (SDEs) arise from physical systems where the parameters describing the system can only be estimated or are subject to noise. Much work has been done recently on developing higher order Runge-Kutta methods for solving SDEs numerically. Fixed stepsize implementations of numerical methods have limitations when, for example, the SDE being solved is stiff as this forces the stepsize to be very small. This paper presents a completely general variable stepsize implementation of an embedded Runge Kutta pair for solving SDEs numerically; in this implementation, there is no restriction on the value used for the stepsize, and it is demonstrated that the integration remains on the correct Brownian path.

The differential role of alcohol expectancies, drinking refusal self-efficacy and coping resources in predicting alcohol consumption in community and clinical samples

This study explored the relationship between coping, alcohol expectancies and drinking refusal self-efficacy in predicting drinking behaviour in both community and clinical samples. These variables were found to have differential effects in their association with frequency and volume of alcohol consumption across the two samples. Generally, drinking refusal self-efficacy was a more salient factor in relation to frequency and volume of community drinking, while coping and expectancies were more strongly associated with frequency of drinking sessions by problem drinkers. The interaction between expectancies and drinking refusal self-efficacy was related to volume of consumption in both groups, while coping and expectancies interacted in their association with frequency in the clinical group. The findings are discussed with regard to the different patterns of cognitive variables governing the decision to drink and the amount consumed in each drinking session, which may differentiate community and problem drinkers.

Differential sorting and fate of endocytosed GPI-anchored proteins

In this paper, we studied the fate of endocytosed glycosylphosphatidyl inositol anchored proteins (GPI-APs) in mammalian cells, using aerolysin, a bacterial toxin that binds to the GPI anchor, as a probe. We find that GPI-APs are transported down the endocytic pathway to reducing late endosomes in BHK cells, using biochemical, morphological and functional approaches. We also find that this transport correlates with the association to raft-like membranes and thus that lipid rafts are present in late endosomes (in addition to the Golgi and the plasma membrane). In marked contrast, endocytosed GPI-APs reach the recycling endosome in CHO cells and this transport correlates with a decreased raft association. GPI-APs are, however, diverted from the recycling endosome and routed to late endosomes in CHO cells, when their raft association is increased by clustering seven or less GPI-APs with an aerolysin mutant. We conclude that the different endocytic routes followed by GPI-APs in different cell types depend on the residence time of GPI-APs in lipid rafts, and hence that raft partitioning regulates GPI-APs sorting in the endocytic pathway.

Robust model-order reduction of complex biological processes

This paper addresses robust model-order reduction of a high dimensional nonlinear partial differential equation (PDE) model of a complex biological process. Based on a nonlinear, distributed parameter model of the same process which was validated against experimental data of an existing, pilot-scale BNR activated sludge plant, we developed a state-space model with 154 state variables in this work. A general algorithm for robustly reducing the nonlinear PDE model is presented and based on an investigation of five state-of-the-art model-order reduction techniques, we are able to reduce the original model to a model with only 30 states without incurring pronounced modelling errors. The Singular perturbation approximation balanced truncating technique is found to give the lowest modelling errors in low frequency ranges and hence is deemed most suitable for controller design and other real-time applications. (C) 2002 Elsevier Science Ltd. All rights reserved.