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.

Theory of multidimensional parametric band-gap simultons

Multidimensional spatiotemporal parametric simultons (simultaneous solitary waves) are possible in a nonlinear chi((2)) medium with a Bragg grating structure, where large effective dispersion occurs near two resonant band gaps for the carrier and second-harmonic field, respectively. The enhanced dispersion allows much reduced interaction lengths, as compared to bulk medium parametric simultons. The nonlinear parametric band-gap medium permits higher-dimensional stationary waves to form. In addition, solitons can occur with lower input powers than conventional nonlinear Schrodinger equation gap solitons. In this paper, the equations for electromagnetic propagation in a grating structure with a parametric nonlinearity are derived from Maxwell's equation using a coupled mode Hamiltonian analysis in one, two, and three spatial dimensions. Simultaneous solitary wave solutions are proved to exist by reducing the equations to the coupled equations describing a nonlinear parametric waveguide, using the effective-mass approximation (EMA). Exact one-dimensional numerical solutions in agreement with the EMA solutions are also given. Direct numerical simulations show that the solutions have similar types of stability properties to the bulk case, providing the carrier waves are tuned to the two Bragg resonances, and the pulses have a width in frequency space less than the band gap. In summary, these equations describe a physically accessible localized nonlinear wave that is stable in up to 3 + 1 dimensions. Possible applications include photonic logic and switching devices. [S1063-651X(98)06109-1].

Measurements on trapped laser-cooled ions using quantum computations

There are some interesting connections between the theory of quantum computation and quantum measurement. As an illustration, we present a scheme in which an ion trap quantum computer can be used to make arbitrarily accurate measurements of the quadrature phase variables for the collective vibrational motion of the ion. We also discuss some more general aspects of quantum computation and measurement in terms of the Feynman-Deutsch principle.

A relative gradient theory for layered materials

It is possible to remedy certain difficulties with the description of short wave length phenomena and interfacial slip in standard models of a laminated material by considering the bending stiffness of the layers. If the couple or moment stresses are assumed to be proportional to the relative deformation gradient, then the bending effect disappears for vanishing interface slip, and the model correctly reduces to an isotropic standard continuum. In earlier Cosserat-type models this was not the case. Laminated materials of the kind considered here occur naturally as layered rock, or at a different scale, in synthetic layered materials and composites. Similarities to the situation in regular dislocation structures with couple stresses, also make these ideas relevant to single slip in crystalline materials. Application of the theory to a one-dimensional model for layered beams demonstrates agreement with exact results at the extremes of zero and infinite interface stiffness. Moreover, comparison with finite element calculations confirm the accuracy of the prediction for intermediate interfacial stiffness.

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.

Predictors of local recurrence after treatment of ductal carcinoma in situ - A meta-analysis

lBACKGROUND. Management of patients with ductal carcinoma in situ (DCIS) is a dilemma, as mastectomy provides nearly a 100% cure rate but at the expense of physical and psychologic morbidity. It would be helpful if we could predict which patients with DCIS are at sufficiently high risk of local recurrence after conservative surgery (CS) alone to warrant postoperative radiotherapy (RT) and which patients are at sufficient risk of local recurrence after CS + RT to warrant mastectomy. The authors reviewed the published studies and identified the factors that may be predictive of local recurrence after management by mastectomy, CS alone, or CS + RT. METHODS. The authors examined patient, tumor, and treatment factors as potential predictors for local recurrence and estimated the risks of recurrence based on a review of published studies. They examined the effects of patient factors (age at diagnosis and family history), tumor factors (sub-type of DCIS, grade, tumor size, necrosis, and margins), and treatment (mastectomy, CS alone, and CS + RT). The 95% confidence intervals (CI) of the recurrence rates for each of the studies were calculated for subtype, grade, and necrosis, using the exact binomial; the summary recurrence rate and 95% CI for each treatment category were calculated by quantitative meta-analysis using the fixed and random effects models applied to proportions. RESULTS, Meta-analysis yielded a summary recurrence rate of 22.5% (95% CI = 16.9-28.2) for studies employing CS alone, 8.9% (95% CI = 6.8-11.0) for CS + RT, and 1.4% (95% CI = 0.7-2.1) for studies involving mastectomy alone. These summary figures indicate a clear and statistically significant separation, and therefore outcome, between the recurrence rates of each treatment category, despite the likelihood that the patients who underwent CS alone were likely to have had smaller, possibly low grade lesions with clear margins. The patients with risk factors of presence of necrosis, high grade cytologic features, or comedo subtype were found to derive the greatest improvement in local control with the addition of RT to CS. Local recurrence among patients treated by CS alone is approximately 20%, and one-half of the recurrences are invasive cancers. For most patients, RT reduces the risk of recurrence after CS alone by at least 50%. The differences in local recurrence between CS alone and CS + RT are most apparent for those patients with high grade tumors or DCIS with necrosis, or of the comedo subtype, or DCIS with close or positive surgical margins. CONCLUSIONS, The authors recommend that radiation be added to CS if patients with DCIS who also have the risk factors for local recurrence choose breast conservation over mastectomy. The patients who may be suitable for CS alone outside of a clinical trial may be those who have low grade lesions with little or no necrosis, and with clear surgical margins. Use of the summary statistics when discussing outcomes with patients may help the patient make treatment decisions. Cancer 1999;85:616-28. (C) 1999 American Cancer Society.

Theoretical and numerical analyses of convective instability in porous media with upward throughflow

Exact analytical solutions have been obtained for a hydrothermal system consisting of a horizontal porous layer with upward throughflow. The boundary conditions considered are constant temperature, constant pressure at the top, and constant vertical temperature gradient, constant Darcy velocity at the bottom of the layer. After deriving the exact analytical solutions, we examine the stability of the solutions using linear stability theory and the Galerkin method. It has been found that the exact solutions for such a hydrothermal system become unstable when the Rayleigh number of the system is equal to or greater than the corresponding critical Rayleigh number. For small and moderate Peclet numbers (Pe less than or equal to 6), an increase in upward throughflow destabilizes the convective flow in the horizontal layer. To confirm these findings, the finite element method with the progressive asymptotic approach procedure is used to compute the convective cells in such a hydrothermal system. Copyright (C) 1999 John Wiley & Sons, Ltd.

Groups with exponent six

Burnside asked questions about periodic groups in his influential paper of 1902. The study of groups with exponent six is a special case of the study of the Burnside questions on which there has been significant progress. It has contributed a number of worthwhile aspects to the theory of groups and in particular to computation related to groups. Finitely generated groups with exponent six are finite. We investigate the nature of relations required to provide proofs of finiteness for some groups with exponent six. We give upper and lower bounds for the number of sixth powers needed to define the largest 2-generator group with exponent six. We solve related questions about other groups with exponent sis using substantial computations which we explain.

Effect of disorder on quantum phase transitions in anisotropic XY spin chains in a transverse field

We present some exact results for the effect of disorder on the critical properties of an anisotropic XY spin chain in a transverse held. The continuum limit of the corresponding fermion model is taken and in various cases results in a Dirac equation with a random mass. Exact analytic techniques can then be used to evaluate the density of states and the localization length. In the presence of disorder the ferromagnetic-paramagnetic or Ising transition of the model is in the same universality class as the random transverse field Ising model solved by Fisher using a real-space renormalization-group decimation technique (RSRGDT). If there is only randomness in the anisotropy of the magnetic exchange then the anisotropy transition (from a ferromagnet in the x direction to a ferromagnet in the y direction) is also in this universality class. However, if there is randomness in the isotropic part of the exchange or in the transverse held then in a nonzero transverse field the anisotropy transition is destroyed by the disorder. We show that in the Griffiths' phase near the Ising transition that the ground-state energy has an essential singularity. The results obtained for the dynamical critical exponent, typical correlation length, and for the temperature dependence of the specific heat near the Ising transition agree with the results of the RSRODT and numerical work. [S0163-1829(99)07125-8].

Transport properties of strongly correlated metals: A dynamical mean-field approach

The temperature dependence of the transport properties of the metallic phase of a frustrated Hubbard model on the hypercubic lattice at half-filling is calculated. Dynamical mean-held theory, which maps the Hubbard model onto a single impurity,Anderson model that is solved self-consistently, and becomes exact in the limit of large dimensionality, is used. As the temperature increases there is a smooth crossover from coherent Fermi liquid excitations at low temperatures to incoherent excitations at high temperatures. This crossover leads to a nonmonotonic temperature dependence for the resistance, thermopower, and Hall coefficient, unlike in conventional metals. The resistance smoothly increases from a quadratic temperature dependence at low temperatures to large values which can exceed the Mott-Ioffe-Regel value ha/e(2) (where a is a lattice constant) associated with mean free paths less than a lattice constant. Further signatures of the thermal destruction of quasiparticle excitations are a peak in the thermopower and the absence of a Drude peak in the optical conductivity. The results presented here are relevant to a wide range of strongly correlated metals, including transition metal oxides, strontium ruthenates, and organic metals.

Finite element analysis of heat transfer and mineralization in layered hydrothermal systems with upward throughflow

We use the finite element method to solve the coupled problem between convective pore-fluid flow, heat transfer and mineralization in layered hydrothermal systems with upward throughflow. In particular, we present the improved rock alteration index (IRAI) concept for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in the systems. To validate the numerical method used in the computation, analytical solutions to a benchmark problem have been derived. After the numerical method is validated, it is used to investigate the pattern of pore-fluid Aom, the distribution of temperature and the mineralization pattern of gold minerals in a layered hydrothermal system with upward throughflow. The related numerical results have demonstrated that the present concept of IRAI is useful and applicable for predicting the most probable precipitation and dissolution regions of gold (Au) minerals in hydrothermal systems. (C) 2000 Elsevier Science S.A. All rights reserved.

A model of specification-based testing of interactive systems

In this paper we present a model of specification-based testing of interactive systems. This model provides the basis for a framework to guide such testing. Interactive systems are traditionally decomposed into a functionality component and a user interface component; this distinction is termed dialogue separation and is the underlying basis for conceptual and architectural models of such systems. Correctness involves both proper behaviour of the user interface and proper computation by the underlying functionality. Specification-based testing is one method used to increase confidence in correctness, but it has had limited application to interactive system development to date.

The pit organs of elasmobranchs: a review

Elasmobranchs have hundreds of tiny sensory organs, called pit organs, scattered over the skin surface. The pit organs were noted in many early studies of the lateral line, but their exact nature has long remained a mystery. Although pit organs were known to be innervated by the lateral line nerves, and light micrographs suggested that they were free neuromasts, speculation that they may be external taste buds or chemoreceptors has persisted until recently Electron micrographs have now revealed that the pit organs are indeed free neuromasts. Their functional and behavioural role(s), however, are yet to be investigated.

Transmission of ocular media in labrid fishes

Wrasses (Labridae) are the second largest family of fishes on the: Great Barrier Reef (after the Gobiidae) and, in terms of morphology and lifestyle, one of the most diverse. They occupy all zones of the reef from the very shadow reef flats to deep slopes, feeding on a variety of fauna. Many wrasses also have elaborately patterned bodies and reflect a range of colours from ultraviolet (UV) to far red. As a first step to investigating the visual system of these fishes we measured the transmission properties of the ocular media of 36 species from the Great Barrier Reef, Australia, and Hawaii, California and the Florida Keys, USA. Transmission measurements were made of whole eyes with a window cut into the back, and also of isolated lenses and corneas. Based on the transmission properties of the corneas the species could be split into two distinct groups within which the exact wavelength of the cut-off was variable. One group had visibly yellow corneas, while the corneas of the other group appeared clear to human observers. Five species had ocular media that transmitted wavelengths below 400 nm, making a perception of UV wavelengths for those species possible. Possible functional roles for the different filler types are discussed.

Teleportation using coupled oscillator states

We analyze the fidelity of teleportation protocols, as a function of resource entanglement, for three kinds of two-mode oscillator states: states with fixed total photon number, number states entangled at a beam splitter, and the two-mode squeezed vacuum state. We define corresponding teleportation protocols for each case including phase noise to model degraded entanglement of each resource.