Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Dependence of transient dynamics in a class-C laser upon variation of inversion with time

The transient statistics of a gain-switched coherently pumped class-C laser displays a linear correlation between the first passage time and subsequent peak intensity. Measurements are reported showing a positive or negative sign of this linear correlation, controlled through the switching time and the laser detuning. Further measurements of the small-signal laser gain combined with calculations involving a three-level laser model indicate that this sign fundamentally depends upon the way the laser inversion varies during the gain switching, despite the added dynamics of the laser polarization in the class-C laser. [S1050-2947(97)07112-6].

Anomalous resonance fluorescence and dressed-state inversion by squeezed light in a Fabry-Perot microcavity

It has been suggested that phased atomic decay in a squeezed vacuum could be detected in the fluorescence spectrum emitted from a driven two-level atom in a cavity. Recently, the existence of other very distinctive features in the fluorescence spectra arising from the nonclassical features of the squeezed vacuum has been reported. In this paper, we investigate the possibility of experimental observation of these spectra. The main obstacle to the experimentalist is ensuring an effective squeezed-vacuum-atom coupling. To overcome this problem we propose the use of a Fabry-Perot microcavity. The analysis involves a consideration of the three-dimensional nature of the electromagnetic held, and the possibility of a mismatch between the squeezed and cavity modes. The problem of squeezing bandwidths is also addressed. We show that under experimentally realistic circumstances many of the spectral anomalies predicted in free space also occur in this environment. In addition, we report large population inversions in the dressed states of the two-level atom. [S1050-2947(98)02301-4].

Matrix elements of U(2n) generators in a multishell spin-orbit basis. I. The del-operator MEs in a two-shell composite Gelfand-Paldus basis

This is the first in a series of three articles which aimed to derive the matrix elements of the U(2n) generators in a multishell spin-orbit basis. This is a basis appropriate to many-electron systems which have a natural partitioning of the orbital space and where also spin-dependent terms are included in the Hamiltonian. The method is based on a new spin-dependent unitary group approach to the many-electron correlation problem due to Gould and Paldus [M. D. Gould and J. Paldus, J. Chem. Phys. 92, 7394, (1990)]. In this approach, the matrix elements of the U(2n) generators in the U(n) x U(2)-adapted electronic Gelfand basis are determined by the matrix elements of a single Ll(n) adjoint tensor operator called the del-operator, denoted by Delta(j)(i) (1 less than or equal to i, j less than or equal to n). Delta or del is a polynomial of degree two in the U(n) matrix E = [E-j(i)]. The approach of Gould and Paldus is based on the transformation properties of the U(2n) generators as an adjoint tensor operator of U(n) x U(2) and application of the Wigner-Eckart theorem. Hence, to generalize this approach, we need to obtain formulas for the complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The nonzero shift coefficients are uniquely determined and may he evaluated by the methods of Gould et al. [see the above reference]. In this article, we define zero-shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis which are appropriate to the many-electron problem. By definition, these are proportional to the corresponding two-shell del-operator matrix elements, and it is shown that the Racah factorization lemma applies. Formulas for these coefficients are then obtained by application of the Racah factorization lemma. The zero-shift adjoint reduced Wigner coefficients required for this procedure are evaluated first. All these coefficients are needed later for the multishell case, which leads directly to the two-shell del-operator matrix elements. Finally, we discuss an application to charge and spin densities in a two-shell molecular system. (C) 1998 John Wiley & Sons.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. II. The two-shell nonzero shift ACCs and U(2n) generator MEs

This is the second in a series of articles whose ultimate goal is the evaluation of the matrix elements (MEs) of the U(2n) generators in a multishell spin-orbit basis. This extends the existing unitary group approach to spin-dependent configuration interaction (CI) and many-body perturbation theory calculations on molecules to systems where there is a natural partitioning of the electronic orbital space. As a necessary preliminary to obtaining the U(2n) generator MEs in a multishell spin-orbit basis, we must obtain a complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. The zero-shift coefficients were obtained in the first article of the series. in this article, we evaluate the nonzero shift adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis. We then demonstrate that the one-shell versions of these coefficients may be obtained by taking the Gelfand-Tsetlin limit of the two-shell formulas. These coefficients,together with the zero-shift types, then enable us to write down formulas for the U(2n) generator matrix elements in a two-shell spin-orbit basis. Ultimately, the results of the series may be used to determine the many-electron density matrices for a partitioned system. (C) 1998 John Wiley & Sons, Inc.

Matrix elements of U(2n) generators in a multishell spin-orbit basis. III. General formulas

This is the third and final article in a series directed toward the evaluation of the U(2n) generator matrix elements (MEs) in a multishell spin/orbit basis. Such a basis is required for many-electron systems possessing a partitioned orbital space and where spin-dependence is important. The approach taken is based on the transformation properties of the U(2n) generators as an adjoint tensor operator of U(n) x U(2) and application of the Wigner-Eckart theorem. A complete set of adjoint coupling coefficients for the two-shell composite Gelfand-Paldus basis (which is appropriate to the many-electron problem) were obtained in the first and second articles of this series. Ln the first article we defined zero-shift coupling coefficients. These are proportional to the corresponding two-shell del-operator matrix elements. See P. J. Burton and and M. D. Gould, J. Chem. Phys., 104, 5112 (1996), for a discussion of the del-operator and its properties. Ln the second article of the series, the nonzero shift coupling coefficients were derived. Having obtained all the necessary coefficients, we now apply the formalism developed above to obtain the U(2n) generator MEs in a multishell spin-orbit basis. The methods used are based on the work of Gould et al. (see the above reference). (C) 1998 John Wiley & Sons, Inc.

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 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.

Transfer matrix eigenvalues of the anisotropic multiparametric U model

A multiparametric extension of the anisotropic U model is discussed which maintains integrability. The R-matrix solving the Yang-Baxter equation is obtained through a twisting construction applied to the underlying U-q(sl (2/1)) superalgebraic structure which introduces the additional free parameters that arise in the model. Three forms of Bethe ansatz solution for the transfer matrix eigenvalues are given which we show to be equivalent.

The interaction of metal ions and Marimastat with matrix metalloproteinase 9

The effect of a range of metal ions on the ability of Marimastat to inhibit matrix metalloproteinase 9 (MMP-9) was examined in a fluorescence based proteolytic assay. Whilst none of the metals examined significantly affected the inhibitory ability of Marimastat, several metal ions did have a significant effect on MMP-9 activity itself. In the absence of Marimastat, Zn(II) and Fe(II) significantly inhibited MMP-9 activity at metal ion concentrations of 10 and 100 muM, respectively. In both the absence and presence of Marimastat, Cd(II) significantly inhibited MMP-9 at 100 muM. In contrast, 1 mM Co(II) significantly upregulated MMP-9 proteolytic activity. (C) 2003 Elsevier Science Inc. All rights reserved.

An optimized model for rat liver perfusion studies

Conditions which influence the viability, integrity, and extraction efficiency of the isolated perfused rat liver were examined to establish optimal conditions for subsequent work in reperfusion injury studies including the choice of buffer, use of oncotic agents, hematocrit, perfusion flow rate, and pressure. Rat livers were perfused with MOPS-buffered Ringer solution with or without erythrocytes. Perfusates were collected and analyzed for blood gases, electrolytes, enzymes, radioactivity in MID studies, and lignocaine in extraction studies. Liver tissue was sampled for histological examinations, and wet:dry weight of the liver was also determined. MOPS-buffered Ringer solution was found to be superior to Krebs bicarbonate buffer, in terms of pH control and buffering capacity, especially during any prolonged period of liver perfusion. A pH of 7.2 is chosen for perfusion since this is the physiological pH of the portal blood. The presence of albumin was important as an oncotic agent, particularly when erythrocytes were used in the perfusate. Perfusion pressure, resistance, and vascular volume are how-dependent and the inclusion of erythrocytes in the perfusate substantially altered the flow characteristics for perfusion pressure and resistance but not vascular volume. Lignocaine extraction was relatively flow-independent. Perfusion injury as defined by enzyme release and tissue fine structure was closely related to the supply of O-2. The optimal conditions for liver perfusion depend upon an adequate supply of oxygen. This can be achieved by using either erythrocyte-free perfusate at a how rate greater than 6 ml/min/g liver or a 20% erythrocyte-containing perfusate at 2 ml/min/g. (C) 1996 Academic Press, Inc.

Finite element analysis of mechanical properties in discontinuously reinforced metal matrix composites with ultrafine microstructure

This investigation focused on the finite element analyses of elastic and plastic properties of aluminium/alumina composite materials with ultrafine microstructure. The commonly used unit cell model was used to predict the elastic properties. By combining the unit cell model with an indentation model, coupled with experimental indentation measurements, the plastic properties of the composites and the associated strengthening mechanism within the metal matrix material were investigated. The grain size of the matrix material was found to be an important factor influencing the mechanical properties of the composites studied. (C) 1997 Elsevier Science S.A.

Integer matrix diagonalization

We consider algorithms for computing the Smith normal form of integer matrices. A variety of different strategies have been proposed, primarily aimed at avoiding the major obstacle that occurs in such computations-explosive growth in size of intermediate entries. We present a new algorithm with excellent performance. We investigate the complexity of such computations, indicating relationships with NP-complete problems. We also describe new heuristics which perform well in practice. Wie present experimental evidence which shows our algorithm outperforming previous methods. (C) 1997 Academic Press Limited.