A numerical study of large sparse matrix exponentials arising in Markov chains

Krylov subspace techniques have been shown to yield robust methods for the numerical computation of large sparse matrix exponentials and especially the transient solutions of Markov Chains. The attractiveness of these methods results from the fact that they allow us to compute the action of a matrix exponential operator on an operand vector without having to compute, explicitly, the matrix exponential in isolation. In this paper we compare a Krylov-based method with some of the current approaches used for computing transient solutions of Markov chains. After a brief synthesis of the features of the methods used, wide-ranging numerical comparisons are performed on a power challenge array supercomputer on three different models. (C) 1999 Elsevier Science B.V. All rights reserved.AMS Classification: 65F99; 65L05; 65U05.

Pore structure tailoring of pillared clays with cation doping techniques

Techniques and mechanism of doping controlled amounts of various cations into pillared clays without causing precipitation or damages to the pillared layered structures are reviewed and discussed. Transition metals of great interest in catalysis can be doped in the micropores of pillared clay in ionic forms by a two-step process. The micropore structures and surface nature of pillared clays are altered by the introduced cations, and this results in a significant improvement in adsorption properties of the clays. Adsorption of water, air components and organic vapors on cation-doped pillared clays were studied. The effects of the amount and species of cations on the pore structure and adsorption behavior are discussed. It is demonstrated that the presence of doped Ca2+ ions can effectively aides the control of modification of the pillared clays of large pore openings. Controlled cation doping is a simple and powerful tool for improving the adsorption properties of pillared clay.

Modelling the large deformations in stratified media - the Cosserat continuum approach

Methods employing continuum approximation in describing the deformation of layered materials possess a clear advantage over explicit models, However, the conventional implicit models based on the theory of anisotropic continua suffers from certain difficulties associated with interface slip and internal instabilities. These difficulties can be remedied by considering the bending stiffness of the layers. This implies the introduction of moment (couple) stresses and internal rotations, which leads to a Cosserat-type theory. In the present model, the behaviour of the layered material is assumed to be linearly elastic; the interfaces are assumed to be elastic perfectly plastic. Conditions of slip or no slip at the interfaces are detected by a Coulomb criterion with tension cut off at zero normal stress. The theory is valid for large deformation analysis. The model is incorporated into the finite element program AFENA and validated against analytical solutions of elementary buckling problems in layered medium. A problem associated with buckling of the roof and the floor of a rectangular excavation in jointed rock mass under high horizontal in situ stresses is considered as the main application of the theory. Copyright (C) 1999 John Wiley & Sons, Ltd.

Chemical synthesis and folding pathways of large cyclic polypeptide: Studies of the cystine knot polypeptide kalata B1

Kalata B1 is a member of a new family of polypeptides, isolated from. plants, which have a cystine knot structure embedded within an amide-cyclized backbone. This family of molecules are the largest known cyclic peptides, and thus, the mechanism of synthesis and folding is of great interest. To provide information about both these phenomena, we have synthesized kalata B1 using two distinct strategies. In the first, oxidation of the cysteine residues of a linear precursor peptide to form the correct disulfide bonds results in folding of the three-dimensional structure and preorganization of the termini in close proximity for subsequent cyclization. The second approach involved cyclization prior to oxidation. In the first method, the correctly folded peptide was produced only in the presence of partially hydrophobic solvent conditions. These conditions are presumably required to stabilize the surface-exposed hydrophobic residues. However,; in the synthesis,involving cyclization prior to oxidation, the cyclic reduced peptide folded to a significant degree in the absence of hydrophobic solvents and even more efficiently in the presence of hydrophobic solvents. Cyclization clearly has a major effect on the folding pathway and facilitates formation of the correctly disulfide-bonded form in aqueous solution; In addition to facilitating folding to a compact stable structure cyclization has an important effect on biological activity as assessed by hemolytic activity.

Designs for an asymmetric gradient set and a compact superconducting magnet for neural magnetic resonance imaging

Imaging of the head and neck is the most commonly performed clinical magnetic resonance imaging (MRI) examination [R. G. Evans and J. R. G. Evans, AJR 157, 603 (1991)]. This is usually undertaken in a generalist MRI instrument containing superconducting magnet system capable of imaging all organs. These generalist instruments are large, typically having a bore of 0.9-1.0 m and a length of 1.7-2.5 m and therefore are expensive to site, somewhat claustrophobic to the patient, and offer little access by attending physicians. In this article, we present the design of a compact, superconducting MRI magnet for head and neck imaging that is less than 0.8 m in length and discuss in detail the design of an asymmetric gradient coil set, tailored to the magnet profile. In particular, the introduction of a radio-frequency FM modulation scheme in concert with a gradient sequence allows the epoch of the linear region of the gradient set to be much closer to the end of the gradient structure than was previously possible. Images from a prototype gradient set demonstrate the effectiveness of the designs. (C) 1999 American Institute of Physics. [S0034-6748(99)04910-2].

Two-dimensional nonlinear dynamics of evanescent-wave guided atoms in a hollow fiber

We describe the classical and quantum two-dimensional nonlinear dynamics of large blue-detuned evanescent-wave guiding cold atoms in hollow fiber. We show that chaotic dynamics exists for classic dynamics, when the intensity of the beam is periodically modulated. The two-dimensional distributions of atoms in (x,y) plane are simulated. We show that the atoms will accumulate on several annular regions when the system enters a regime of global chaos. Our simulation shows that, when the atomic flux is very small, a similar distribution will be obtained if we detect the atomic distribution once each the modulation period and integrate the signals. For quantum dynamics, quantum collapses, and revivals appear. For periodically modulated optical potential, the variance of atomic position will be suppressed compared to the no modulation case. The atomic angular momentum will influence the evolution of wave function in two-dimensional quantum system of hollow fiber.

The effects of change in lead stimulus modality on the modulation of acoustic blink startle

In two experiments we investigated the effect of generalized orienting induced by changing the modality of the lead stimulus on the modulation of blink reflexes elicited by acoustic stimuli. In Experiment 1 (n = 32), participants were presented with acoustic or visual change stimuli after habituation training with tactile lead stimuli. In Experiment 2 (n = 64), modality of the lead stimulus (acoustic vs. visual) was crossed with experimental condition (change vs. no change). Lead stimulus change resulted in increased electrodermal orienting in both experiments. Blink latency shortening and blink magnitude facilitation increased from habituation to change trials regardless of whether the change stimulus was presented in the same or in a different modality as the reflex-eliciting stimulus. These results are not consistent with modality-specific accounts of attentional startle modulation.

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.

H-reflex modulation during passive lengthening and shortening of the human triceps surae

1. The present study investigated the effects of lengthening and shortening actions on IT-reflex amplitude. H-reflexes were evoked in the soleus (SOL) and medial gastroenemius (MG) of human subject, during passive isometric, lengthening and shortening actions performed at angular velocities of 0, +/-2, +/-5 and +/- 15 deg s(-1). 2. H-reflex amplitude, in froth SOL and MG were significantly depressed during passive lengthening actions and facilitated during passive shortening actions, when compared with the isometric R-reflex amplitude. 3. Four experiments were performed in which the latencies front the onset of movement to delivery of the stimulus were altered. Passive H-reflex modulation during lengthening actions was found tee begin at latencies of less than 60 ms suggesting that this inhibition was due to peripheral and/or spinal mechanisms. 4. It is postulated that, the H-reflex modulation seen in the present study is related to the tunic discharge of muscle spindle afferents and the consequent effects of transmission within the la pathway. Inhibition of the H-reflex at less than 60 ms after the onset of muscle lengthening may he attributed to several mechanisms, which cannot be distinguished using the current protocol. These may include the inability to evoke volleys in la fibres that are refractory following muscle spindle discharge during; rapid muscle lengthening, a reduced probability of transmitter release front the presynaptic terminal (homosynaptic post.-activation depression) and presynaptic inhibition of la afferents from plantar flexor agonists. Short latency facilitation of the H-reflex may be attributed to temporal summation of excitatory postsynaptic potentials arising from muscle spindle afferents during rapid muscle lengthening. At longer latencies, presynaptic inhibition of Ia afferents cannot be excluded as a potential inhibitory mechanism.

Field-swept pulsed electron paramagnetic resonance of Cr3+-doped ZBLAN fluoride glass

Field-swept pulsed electron paramagnetic resonance (EPR) spectra of a ZBLAN fluoride glass doped with a low concentration of Cr3+ are obtained using echo-detected EPR and hole-burning free induction decay detection. We review the utility of the pulsed EPR technique in generating field-swept EPR spectra, as well as some of the distorting effects that are peculiar to the pulsed detection method. The application of this technique to Cr3+-doped ZBLAN reveals that much of the broad resonance extending from g(eff) = 5.1 to g(eff) = 1.97, characteristic of X-band continuous wave EPR of Cr3+ in glasses, is absent. We attribute this largely to the variation in nutation frequencies across the spectrum that result from sites possessing large fine structure interactions. The description of the spin dynamics of such sites is complicated and we discuss some possible approaches to the simulation of the pulsed EPR spectra.

Large mitochondrial repeats multiplied during the polymerase chain reaction

The number of repeats in repetitive DNA like micro- and minisatellites is often determined by polymerase chain reaction (PCR). When we counted repeats in an array of mitochondrial repeats in the cattle tick (Boophilus microplus) we found that the number of repeats increased during PCR. Multiplication of the repeats was independent of the primers used to amplify the region, the PCR annealing temperature and the length of the PCR product. The use of PCR to determine the number of repeats in arrays needs to be reassessed. For long repeats, a subset of samples should always be analysed by Southern blot hybridization to confirm the PCR results.

The VITATOPS (vitamins to prevent stroke) Trial: Rationale and design of an international, large, simple, randomised trial of homocysteine-lowering multivitamin therapy in patients with recent transient ischaemic attack or stroke

Background: Epidemiological studies suggest that raised plasma concentrations of total homocysteine (tHcy) may be a common, causal and treatable risk factor for atherothromboembolic ischaemic stroke. Although tHcy can be lowered effectively with small doses of folic acid, vitamin B-12 and vitamin B-6, it is not known whether lowering tHcy, by means of multivitamin therapy, can prevent stroke and other major atherothromboembolic vascular events. Purpose: To determine whether vitamin supplements (folic acid 2 mg, B-6 25 Mg, B-12 500 mug) reduce the risk of stroke, and other serious vascular events, in patients with recent stroke or transient ischaemic attacks of the brain or eye (TIA). Methods: An international, multi-centre, randomised, double-blind, placebo-controlled clinical trial. Results: As of November 2001, more than 1,400 patients have been randomised from 10 countries in four continents. Conclusion: VITATOPS aims to recruit and follow up 8,000 patients between 2000 and 2004, and provide a reliable estimate of the safety and effectiveness of dietary supplementation with folic acid, vitamin B-12, and vitamin B-6 in reducing recurrent serious vascular events among a wide range of patients with TIA and stroke. Copyright (C) 2002 S. Karger AG, Basel.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.