Lanczos subspace filter diagonalization: Homogeneous recursive filtering and a low-storage method for the calculation of matrix elements

We develop a new iterative filter diagonalization (FD) scheme based on Lanczos subspaces and demonstrate its application to the calculation of bound-state and resonance eigenvalues. The new scheme combines the Lanczos three-term vector recursion for the generation of a tridiagonal representation of the Hamiltonian with a three-term scalar recursion to generate filtered states within the Lanczos representation. Eigenstates in the energy windows of interest can then be obtained by solving a small generalized eigenvalue problem in the subspace spanned by the filtered states. The scalar filtering recursion is based on the homogeneous eigenvalue equation of the tridiagonal representation of the Hamiltonian, and is simpler and more efficient than our previous quasi-minimum-residual filter diagonalization (QMRFD) scheme (H. G. Yu and S. C. Smith, Chem. Phys. Lett., 1998, 283, 69), which was based on solving for the action of the Green operator via an inhomogeneous equation. A low-storage method for the construction of Hamiltonian and overlap matrix elements in the filtered-basis representation is devised, in which contributions to the matrix elements are computed simultaneously as the recursion proceeds, allowing coefficients of the filtered states to be discarded once their contribution has been evaluated. Application to the HO2 system shows that the new scheme is highly efficient and can generate eigenvalues with the same numerical accuracy as the basic Lanczos algorithm.

A review of the use of the basic cation saturation ratio and the "ideal" soil

The use of 'balanced' Ca, Mg, and K ratios, as prescribed by the basic cation saturation ratio (BCSR) concept, is still used by some private soil-testing laboratories for the interpretation of soil analytical data. This review aims to examine the suitability of the BCSR concept as a method for the interpretation of soil analytical data. According to the BCSR concept, maximum plant growth will be achieved only when the soil’s exchangeable Ca, Mg, and K concentrations are approximately 65 % Ca, 10 % Mg, and 5 % K (termed the ‘ideal soil’). This ‘ideal soil’ was originally proposed by Firman Bear and co-workers in New Jersey (USA) during the 1940s as a method of reducing luxury K uptake by alfalfa (Medicago sativa L.). At about the same time, William Albrecht, working in Missouri (USA), concluded through his own investigations that plants require a soil with a high Ca saturation for optimal growth. Whilst it now appears that several of Albrecht’s experiments were fundamentally flawed, the BCSR (‘balanced soil’) concept has been widely promoted, suggesting that the prescribed cationic ratios provide optimum chemical, physical, and biological soil properties. Our examination of data from numerous studies (particularly those of Albrecht and Bear, themselves) would suggest that, within the ranges commonly found in soils, the chemical, physical, and biological fertility of a soil is generally not influenced by the ratios of Ca, Mg, and K. The data do not support the claims of the BCSR, and continued promotion of the BCSR will result in the inefficient use of resources in agriculture and horticulture.

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.

hnRNP A2 selectively binds the cytoplasmic transport sequence of myelin basic protein mRNA

Segregation of mRNAs in the cytoplasm of polar cells has been demonstrated for proteins involved in Xenopus and Drosophila oogenesis, and for some proteins in somatic cells. It is assumed that vectorial transport of the messages is generally responsible for this localization. The mRNA encoding the basic protein of central nervous system myelin is selectively transported to the distal ends of the processes of oligodendrocytes, where it is anchored to the myelin membrane and translated. This transport is dependent on a 21-nucleotide cis-acting segment of the 3'-untranslated region (RTS). Proteins that bind to this cis-acting segment have now been isolated from extracts of rat brain. A group of six 35-42-kDa proteins bind to a 35-base oligoribonucleotide incorporating the RTS, but not to several oligoribonucleotides with the same composition but randomized sequences, thus establishing specificity for the base sequence in the RTS. The most abundant of these proteins has been identified, by Edman sequencing of tryptic peptides and mass spectroscopy, as heterogeneous nuclear ribonucleoprotein (hnRNP) A2, a 36-kDa member of a family of proteins that are primarily, but not solely, intranuclear. This protein was most abundant in samples from rat brain and testis, with lower amounts in other tissues. It was separated from the other polypeptides by using reverse-phase HPLC and shown to retain preferential association with the RTS. In cultured oligodendrocytes, hnRNP A2 was demonstrated by confocal microscopy to be distributed throughout the nucleus, cell soma, and processes.

Analysis of pore-fluid pressure gradient and effective vertical-stress gradient distribution in layered hydrodynamic systems

A theoretical analysis is carried out to investigate the pore-fluid pressure gradient and effective vertical-stress gradient distribution in fluid saturated porous rock masses in layered hydrodynamic systems. Three important concepts, namely the critical porosity of a porous medium, the intrinsic Fore-fluid pressure and the intrinsic effective vertical stress of the solid matrix, are presented and discussed. Using some basic scientific principles, we derive analytical solutions and explore the conditions under which either the intrinsic pore-fluid pressure gradient or the intrinsic effective vertical-stress gradient can be maintained at the value of the lithostatic pressure gradient. Even though the intrinsic pore-fluid pressure gradient can be maintained at the value of the lithostatic pressure gradient in a single layer, it is impossible to maintain it at this value in all layers in a layered hydrodynamic system, unless all layers have the same permeability and porosity simultaneously. However, the intrinsic effective vertical-stress gradient of the solid matrix can be maintained at a value close to the lithostatic pressure gradient in all layers in any layered hydrodynamic system within the scope of this study.

Type 1 and Type 2 cytokines: from basic science to fungal infections

At the present time, it is clear that Th1 responses afford protection against the fungi; however, the development, maintenance and function of the protective immune responses are complex mechanisms and are influenced by multiple factors. The route of infection has been shown to affect initial cytokine production and, consequently, the induction of protective Th1 responses. The ability of different isolates of the same fungal agent to induce and sustain a protective response has also been emphasized. Protective immune responses have been shown to vary in genetically different mouse strains after infection. In addition, these protective responses, such as cellular influx and cytokine production, also vary within the same animal depending on the tissue infected. The functional dominance of certain cytokines over others in influencing development and maintenance of protective responses has been discussed. Certain cytokines may act differently in hosts lacking important components of their innate or immune repertoire. It is evident from these presentations that a more comprehensive understanding of the protective mechanisms against different fungal agents is emerging. However, there is still much to learn before cytokine modulatory therapy can be used effectively without risk in the human host.

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.

TFEC is a macrophage-restricted member of the microphthalmia-TFE subfamily of basic helix-loop-helix leucine zipper transcription factors

The murine homologue of the TFEC was cloned as part of an analysis of the expression of the microphthalmia-TFE (MiT) subfamily of transcription factors in macrophages. TFEC, which most likely acts as a transcriptional repressor in heterodimers with other MiT family members, was identified in cells of the mononuclear phagocyte lineage, coexpressed,vith all other known MiT subfamily members (Mitf, TFE3, TFEB), Northern blot analysis of several different cell lineages indicated that the expression of murine TFEC (mTFEC) was restricted to macrophages. A 600-bp fragment of the TATA-less putative proximal promoter of TFEC shares features with many known macrophage-specific promoters and preferentially directs luciferase expression in the RAW264.7 macrophage cell line in transient transfection assays. Five of six putative Ets motifs identified in the TFEC promoter bind the macrophage-restricted transcription factor PU,I under in vitro conditions and in transfected 3T3 fibroblasts; the minimal luciferase activity of the TFEC promoter could be induced by coexpression of PU.1 or the related transcription factor Ets-2. The functional importance of the tissue-restricted expression of TFEC and a possible role in macrophage-specific gene regulation require further investigation, but are likely to be linked to the role of the other MiT family members in this lineage.

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.

Identification of the glycosaminoglycan keratan sulfate in the prostatic secretory cell

BACKGROUND. Prostate secretory granules (PSG) represent the basic secretory unit of the prostate gland, containing many of its exocrine proteases. Recent analysis of intraluminal corpora amylacea, a proposed by-product of PSG secretion, detected sulfated glycosaminoglycans (GAG) possibly keratan sulfate (KS),indicating a secretory mechanism for GAG in the human prostate surface epithelial cell. METHODS. Immunostains using anti-KS and anti-prostate-specific antigen (PSA) were evaluated on 10 sequential radical prostatectomy specimens, three of which had received neoadjuvant antiandrogen therapy. Extracts of normal secretory tissue as well as a sample composed almost entirely of prostatic stroma were subjected to Western blot analysis, using the same antibody panel. RESULTS. Keratan sulfate secretion from the normal prostate epithelial cell has been confirmed and correlates, as does PSA, with the presence of cytoplasmic PSG. No such correlation exists in most adenocarcinomas or in benign epithelium after androgen ablation. Western blot analyses confirmed tissue immunostains and demonstrated a secretory proteoglycan of 70-95 kDa. CONCLUSIONS. Recognition of PSG heralds a novel secretory mechanism within the human prostate gland that is linked to the secretion of KS. The role of KS in normal prostate secretion remains unknown, although it appears downregulated in neoplastic and androgen-ablated cells. (C) 2000 Wiley-Liss, Inc.