Expokit: A software package for computing matrix exponentials

Expokit provides a set of routines aimed at computing matrix exponentials. More precisely, it computes either a small matrix exponential in full, the action of a large sparse matrix exponential on an operand vector, or the solution of a system of linear ODEs with constant inhomogeneity. The backbone of the sparse routines consists of matrix-free Krylov subspace projection methods (Arnoldi and Lanczos processes), and that is why the toolkit is capable of coping with sparse matrices of large dimension. The software handles real and complex matrices and provides specific routines for symmetric and Hermitian matrices. The computation of matrix exponentials is a numerical issue of critical importance in the area of Markov chains and furthermore, the computed solution is subject to probabilistic constraints. In addition to addressing general matrix exponentials, a distinct attention is assigned to the computation of transient states of Markov chains.

A six-year longitudinal study of the relationship of physical activity to bone mineral accrual in growing children: The University of Saskatchewan bone mineral accrual study

To investigate the influence of physical activity on bone mineral accrual during the adolescent years, we analyzed 6 years of data from 53 girls and 60 boys. Physical activity, dietary intakes, and anthropometry were measured every 6 months and dual-energy X-ray absorptiometry scans of the total body (TB), lumbar spine (LS), and proximal femur (Hologic 2000, array mode) were collected annually. Distance and velocity curves for height and bone mineral content (BMC) were fitted for each child at several skeletal sites using a cubic spline procedure, from which ages at peak height velocity (PHV) and peak BMC velocity (PBMCV) were identified. A mean age- and gender-specific standardized activity (Z) score was calculated for each subject based on multiple yearly activity assessments collected up until age of PHV. This score was used to identify active (top quartile), average (middle 2 quartiles), or inactive (bottom quartile) groups. Two-way analysis of covariance, with height and weight at PHV controlled for, demonstrated significant physical activity and gender main effects (but no interaction) for PBMCV, for BMC accrued for 2 years around peak velocity, and for BMC at 1 year post-PBMCV for the TB and femoral neck and for physical activity but not gender at the LS (all p < 0.05). Controlling for maturational and size differences between groups, we noted a 9% and 17% greater TB BMC for active boys and girls, respectively, over their inactive peers 1 year after the age of PBMCV. We also estimated that, on average, 26% of adult TB bone mineral was accrued during the 2 years around PBMCV.

A hybrid, inverse approach to the design of magnetic resonance imaging magnets

This paper describes a hybrid numerical method of an inverse approach to the design of compact magnetic resonance imaging magnets. The problem is formulated as a field synthesis and the desired current density on the surface of a cylinder is first calculated by solving a Fredholm equation of the first, kind. Nonlinear optimization methods are then invoked to fit practical magnet coils to the desired current density. The field calculations are performed using a semi-analytical method. The emphasis of this work is on the optimal design of short MRI magnets. Details of the hybrid numerical model are presented, and the model is used to investigate compact, symmetric MRI magnets as well as asymmetric magnets. The results highlight that the method can be used to obtain a compact MRI magnet structure and a very homogeneous magnetic field over the central imaging volume in clinical systems of approximately 1 m in length, significantly shorter than current designs. Viable asymmetric magnet designs, in which the edge of the homogeneous region is very close to one end of the magnet system are also presented. Unshielded designs are the focus of this work. This method is flexible and may be applied to magnets of other geometries. (C) 2000 American Association of Physicists in Medicine. [S0094-2405(00)00303-5].

Acyclic permutants of naturally occurring cyclic proteins - Characterization of cystine knot and beta-sheet formation in the macrocyclic polypeptide kalata B1

Kalata B1 is a prototypic member of the unique cyclotide family of macrocyclic polypeptides in which the major structural features are a circular peptide backbone, a triple stranded beta-sheet, and a cystine knot arrangement of three disulfide bonds. The cyclotides are the only naturally occurring family of circular proteins and have prompted us to explore the concept of acyclic permutation, i.e. opening the backbone of a cross-linked circular protein in topologically permuted ways. We have synthesized the complete suite of acyclic permutants of kalata B1 and examined the effect of acyclic permutation on structure and activity. Only two of six topologically distinct backbone loops are critical for folding into the native conformation, and these involve disruption of the embedded ring in the cystine knot. Surprisingly, it is possible to disrupt regions of the p-sheet and still allow folding into native-like structure, provided the cystine knot is intact. Kalata B1 has mild hemolytic activity, but despite the overall structure of the native peptide being retained in all but two cases, none of the acyclic permutants displayed hemolytic activity. This loss of activity is not localized to one particular region and suggests that cyclization is critical for hemolytic activity.

Identifying rate-limiting nodes in large-scale cortical networks for visuospatial processing: an illustration using fMRI

With the advent of functional neuroimaging techniques, in particular functional magnetic resonance imaging (fMRI), we have gained greater insight into the neural correlates of visuospatial function. However, it may not always be easy to identify the cerebral regions most specifically associated with performance on a given task. One approach is to examine the quantitative relationships between regional activation and behavioral performance measures. In the present study, we investigated the functional neuroanatomy of two different visuospatial processing tasks, judgement of line orientation and mental rotation. Twenty-four normal participants were scanned with fMRI using blocked periodic designs for experimental task presentation. Accuracy and reaction time (RT) to each trial of both activation and baseline conditions in each experiment was recorded. Both experiments activated dorsal and ventral visual cortical areas as well as dorsolateral prefrontal cortex. More regionally specific associations with task performance were identified by estimating the association between (sinusoidal) power of functional response and mean RT to the activation condition; a permutation test based on spatial statistics was used for inference. There was significant behavioral-physiological association in right ventral extrastriate cortex for the line orientation task and in bilateral (predominantly right) superior parietal lobule for the mental rotation task. Comparable associations were not found between power of response and RT to the baseline conditions of the tasks. These data suggest that one region in a neurocognitive network may be most strongly associated with behavioral performance and this may be regarded as the computationally least efficient or rate-limiting node of the network.

A family of perfect factorisations of complete bipartite graphs

A 1-factorisation of a graph is perfect if the union of any two of its 1-factors is a Hamiltonian cycle. Let n = p(2) for an odd prime p. We construct a family of (p-1)/2 non-isomorphic perfect 1-factorisations of K-n,K-n. Equivalently, we construct pan-Hamiltonian Latin squares of order n. A Latin square is pan-Hamiltoilian if the permutation defined by any row relative to any other row is a single Cycle. (C) 2002 Elsevier Science (USA).

Coaxial flattening at deep levels of orogenic belts: evidence from blueschists and eclogites on Syros and Sifnos (Cyclades, Greece)

This work presents new Structural data from a high-pressure/low-temperature (HP/LT) metamorphic terrane exposed on the islands of Syros and Sifnos (Cyclades, Greece). The structure and the metamorphism of a relatively coherent HP/LT rock section were studied in order to elucidate how strain was accommodated at deep crustal levels during the formation and exhumation of HP/LT rocks. At least three deformation phases associated with eclogite- and blueschist-facies conditions (P = 8-15 kbar; T = 400-550 degreesC) were recognised. The earliest deformation fabric (S1), preserved as inclusion trails within garnet porphyroblasts, is aligned to define a sub-vertical schistosity (at present orientation), which is frequently orthogonal to the flat matrix schistosity (S2), and may indicate that deep crustal thickening involved upright folding. The currently dominant fabric in the HP rock section, S2, is Usually moderately dipping and locally contains NW-trending glaucophane lineations, symmetric pressure-shadows and eclogitic boudins. The symmetric structures associated with this fabric seem to indicate coaxial vertical thinning, although the existence of non-coaxial structures out of the study area cannot be excluded. Glaucophane-bearing shear bands (S3), with top-to-NW sense of shearing, locally crosscut the earlier structures. The latest recognised fabric (D4) is scarce and often absent within the HP rocks. It is associated with top-to-NE kinematic criteria that formed at greenschist-facies conditions (P = 4-7 kbar; T = 400-450 degreesC). Based on these observations, it is suggested that partitioning of strain occurred at different crustal levels and at different times. Deep crustal deformation was governed by thickening via upright folding followed by coaxial vertical thinning, whereas non-coaxial shearing occurred when the rocks were already exhumed to relatively shallow crustal levels. The earliest fabrics (D1 to D3) pertain to Alpine orogenesis and possibly to syn-orogenic extension, whereas the latest correspond to whole-crust back-are extension. (C) 2002 Elsevier Science Ltd. All rights reserved.

Parity measurement of one- and two-electron double well systems

We outline a scheme to accomplish measurements of a solid state double well system (DWS) with both one and two electrons in nonlocalized bases. We show that, for a single particle, measuring the local charge distribution at the midpoint of a DWS using a SET as a sensitive electrometer amounts to performing a projective measurement in the parity (symmetric/antisymmetric) eigenbasis. For two-electrons in a DWS, a similar configuration of SET results in close-to-projective measurement in the singlet/triplet basis. We analyze the sensitivity of the scheme to asymmetry in the SET position for some experimentally relevant parameter, and show that it is experimentally realizable.

The use of electrospray ionization tandem mass spectrometry on the structural characterization of novel asymmetric metallo-organic supermolecules, based on pentafluorophenylporphyrins and ruthenium complexes

The novel asymmetric metallo-organic triads cis- and trans-[B(4-py)BPFPH(2){Ru(3)O(Ac)(6)(py)(2)}(Ru(bpy)(2)Cl}](PF(6))(2) (5a,b) for which cis- and trans-B(4-py)BPFPH(2)=5,10-bis(pentafluorophenyl)-15,20-bis(4-pyridyl)porphyrin and 5,15-bis(pentafluorophenyl)-10,20-bis(4-pyridyl)porphyrin, respectively; Ac = acetate; py = pyridine and bpy = 2,2`-bipyridine, as well as their corresponding monosubstituted dyads cis- and trans-[B(4-py)BPFPH(2){Ru(3)O(Ac)(6)(py)(2)}]PF(6) (4a,b) have been structurally characterized via electrospray ionization mass spectrometry (ESI-MS and ESI-MS/MS). The ESI-MS of dyads 4a,b display two characteristic Ru-multicomponent clusters of isotopologue ions corresponding to singly charged ions 4a,b(+) of m/z 1629 and doubly charged ions [4a,b+H](2+) of m/z 815 and the triads 5a,b are detected by ESI-MS as the intact doubly charged cluster of isotopologue ions of m/z 1039 [5a,b](2+). The ESI-MS/MS of 4a,b(+), [4a,b+H](2+) and [5a,b](2+) reveal characteristic dissociation pathways, which confirm the structural assignments providing additional information on the intrinsic binding strengths of the gaseous ions. Although the gas-phase behavior of each pair of isomers was rather similar, the less symmetric dyads 4a,b are distinguished via the (1)H NMR spectral profile of the pyrrolic signals. Exploratory photophysical assays have shown that both modifying motifs alter the porphyrinic core emission profile, opening the possibility to use these asymmetric systems as photophysical devices. (C) 2008 Elsevier Ltd. All rights reserved.

Drinfield twists and alegebraic Bethe ansatz of the supersymmetric t-J model

We construct the Drinfeld twists (factorizing F-matrices) for the supersymmetric t-J model. Working in the basis provided by the F-matrix (i.e. the so-called F-basis), we obtain completely symmetric representations of the monodromy matrix and the pseudo-particle creation operators of the model. These enable us to resolve the hierarchy of the nested Bethe vectors for the gl(2\1) invariant t-J model.

Binary games with many players

We examine a problem with n players each facing the same binary choice. One choice is superior to the other. The simple assumption of competition - that an individual's payoff falls with a rise in the number of players making the same choice, guarantees the existence of a unique symmetric equilibrium (involving mixed strategies). As n increases, there are two opposing effects. First, events in the middle of the distribution - where a player finds itself having made the same choice as many others - become more likely, but the payoffs in these events fall. In opposition, events in the tails of the distribution - where a player finds itself having made the same choice as few others - become less likely, but the payoffs in these events remain high. We provide a sufficient condition (strong competition) under which an increase in the number of players leads to a reduction in the equilibrium probability that the superior choice is made.

Corruption and auctions

We investigate how corruption affects the outcome of a first-price auction (bidding behavior, efficiency and the seller's expected revenue). The auctioneer approaches the winner to offer the possibility of a reduction in his bid in exchange for a bribe. The bribe can be a percentage of the difference between the winning and the second-highest bid or a fixed amount. We show that there exists a symmetric bidding strategy equilibrium that is monotone, i.e., higher valuation buyers bid higher. Corruption does not affect efficiency but both the auctioneer's expected bribe and the seller's expected revenue depend on the format of the bribe payments. We also find the optimal bribe scheme.

Periodic electric field enhanced transport through membranes

We examine the mean flux across a homogeneous membrane of a charged tracer subject to an alternating, symmetric voltage waveform. The analysis is based on the Nernst-Planck flux equation, with electric field subject to time dependence only. For low frequency electric fields the quasi steady-state flux can be approximated using the Goldman model, which has exact analytical solutions for tracer concentration and flux. No such closed form solutions can be found for arbitrary frequencies, however we find approximations for high frequency. An approximation formula for the average flux at all frequencies is also obtained from the two limiting approximations. Numerical integration of the governing equation is accomplished by use of the numerical method of lines and is performed for four different voltage waveforms. For the different voltage profiles, comparisons are made with the approximate analytical solutions which demonstrates their applicability. (c) 2005 Elsevier B.V. All rights reserved.

Reconstruction of large cranial defects in nonimmunosuppressed experimental design with human dental pulp stem cells

The main aim of this study is to evaluate the capacity of human dental pulp stem cells (hDPSC), isolated from deciduous teeth, to reconstruct large-sized cranial bone defects in nonimmunosuppressed (NIS) rats. To our knowledge, these cells were not used before in similar experiments. We performed two symmetric full-thickness cranial defects (5 x 8 mm) on each parietal region of eight NIS rats. In six of them, the left side was supplied with collagen membrane only and the right side (RS) with collagen membrane and hDPSC. In two rats, the RS had collagen membrane only and nothing was added at the left side (controls). Cells were used after in vitro characterization as mesenchymal cells. Animals were euthanized at 7, 20, 30, 60, and 120 days postoperatively and cranial tissue samples were taken from the defects for histologic analysis. Analysis of the presence of human cells in the new bone was confirmed by molecular analysis. The hDPSC lineage was positive for the four mesenchymal cell markers tested and showed osteogenic, adipogenic, and myogenic in vitro differentiation. We observed bone formation 1 month after surgery in both sides, but a more mature bone was present in the RS. Human DNA was polymerase chain reaction-amplified only at the RS, indicating that this new bone had human cells. The us e of hDPSC in NIS rats did not cause any graft. rejection. Our findings suggest that hDPSC is an additional cell resource for correcting large cranial defects in rats and constitutes a promising model for reconstruction of human large cranial defects in craniofacial surgery.

High-breakdown linear discriminant analysis

The classification rules of linear discriminant analysis are defined by the true mean vectors and the common covariance matrix of the populations from which the data come. Because these true parameters are generally unknown, they are commonly estimated by the sample mean vector and covariance matrix of the data in a training sample randomly drawn from each population. However, these sample statistics are notoriously susceptible to contamination by outliers, a problem compounded by the fact that the outliers may be invisible to conventional diagnostics. High-breakdown estimation is a procedure designed to remove this cause for concern by producing estimates that are immune to serious distortion by a minority of outliers, regardless of their severity. In this article we motivate and develop a high-breakdown criterion for linear discriminant analysis and give an algorithm for its implementation. The procedure is intended to supplement rather than replace the usual sample-moment methodology of discriminant analysis either by providing indications that the dataset is not seriously affected by outliers (supporting the usual analysis) or by identifying apparently aberrant points and giving resistant estimators that are not affected by them.