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.

Enamel matrix derivative induces matrix synthesis by cultured human periodontal fibroblast cells

Background: Periodontal wound healing and regeneration require that new matrix be synthesized, creating an environment into which cells can migrate. One agent which has been described as promoting periodontal regeneration is an enamel matrix protein derivative (EMD). Since no specific growth factors have been identified in EMD preparations, it is postulated that EMD acts as a matrix enhancement factor. This study was designed to investigate the effect of EMD in vitro on matrix synthesis by cultured periodontal fibroblasts. Methods: The matrix response of the cells was evaluated by determination of the total proteoglycan synthesis, glycosaminoglycan profile, and hyaluronan synthesis by the uptake of radiolabeled precursors. The response of the individual proteoglycans, versican, decorin, and biglycan were examined at the mRNA level by Northern blot analysis. Hyaluronan synthesis was probed by identifying the isotypes of hyaluronan synthase (HAS) expressed in periodontal fibroblasts as HAS-2 and HAS-3 and the effect of EMD on the levels of mRNA for each enzyme was monitored by reverse transcription polymerase chain reaction (RTPCR). Comparisons were made between gingival fibroblast (GF) cells and periodontal ligament (PDLF) cells. Results: EMD was found to significantly affect the synthesis of the mRNAs for the matrix proteoglycans versican, biglycan, and decorin, producing a response similar to, but potentially greater than, mitogenic cytokines. EMD also stimulated hyaluronan synthesis in both GF and PDLF cells. Although mRNA for HAS-2 was elevated in GF after exposure to EMD, the PDLF did not show a similar response. Therefore, the point at which the stimulation of hyaluronan becomes effective may not be at the level of stimulation of the mRNA for hyaluronan synthase, but, rather, at a later point in the pathway of regulation of hyaluronan synthesis. In all cases, GF cells appeared to be more responsive to EMD than PDLF cells in vitro. Conclusions: EMD has the potential to significantly modulate matrix synthesis in a manner consistent with early regenerative events.

Matrix metalloproteinases and their inhibitors in oral lichen planus

Background: Oral lichen planus (OLP) is characterized by a subepithelial lymphocytic infiltrate, basement membrane (BM) disruption, intra-epithelial T-cell migration and apoptosis of basal keratinocytes. BM damage and T-cell migration in OLP may be mediated by matrix metalloproteinases (MMPs). Methods: We examined the distribution, activation and cellular sources of MMPs and their inhibitors (TIMPs) in OLP using immunohistochemistry, ELISA, RT-PCR and zymography. Results: MMP-2 and -3 were present in the epithelium while MMP-9 was associated with the inflammatory infiltrate. MMP-9 and TIMP-1 secretion by OLP lesional T cells was greater than OLP patient (p

Matrix-product codes over Fq

Codes C-1,...,C-M of length it over F-q and an M x N matrix A over F-q define a matrix-product code C = [C-1 (...) C-M] (.) A consisting of all matrix products [c(1) (...) c(M)] (.) A. This generalizes the (u/u + v)-, (u + v + w/2u + v/u)-, (a + x/b + x/a + b + x)-, (u + v/u - v)- etc. constructions. We study matrix-product codes using Linear Algebra. This provides a basis for a unified analysis of /C/, d(C), the minimum Hamming distance of C, and C-perpendicular to. It also reveals an interesting connection with MDS codes. We determine /C/ when A is non-singular. To underbound d(C), we need A to be 'non-singular by columns (NSC)'. We investigate NSC matrices. We show that Generalized Reed-Muller codes are iterative NSC matrix-product codes, generalizing the construction of Reed-Muller codes, as are the ternary 'Main Sequence codes'. We obtain a simpler proof of the minimum Hamming distance of such families of codes. If A is square and NSC, C-perpendicular to can be described using C-1(perpendicular to),...,C-M(perpendicular to) and a transformation of A. This yields d(C-perpendicular to). Finally we show that an NSC matrix-product code is a generalized concatenated code.

Thermal characterization of copolymers obtained by radiation grafting of styrene onto poly(tetrafluoroethyleneperfluoropropylvinylether) substrates: thermal decomposition, melting behavior and low-temperature transitions

Differential scanning calorimetric (DSC) and thermogravimetric analysis (TGA) have been used to study the thermal decomposition, the melting behavior and low-temperature transitions of copolymers obtained by radiation-induced grafting of styrene onto poly (tetrafluoroethylene- perfluoropropylvinylether) (PFA) substrates. PFA with different contents of perfluoropropylvinylether (PPVE) as a comonomer have been investigated. A two step degradation pattern was observed from TGA thermograms of all the grafted copolymers, which was attributed to degradation of PSTY followed by the degradation of the PFA backbone at higher temperature. One broad melting peak can be identified for all copolymers, which has two components in the samples with higher PPVE content. The melting peak, crystal-crystal transition and the degree of crystallinity of the grafted copolymers increases with radiation grafting up to 50 kGy, followed by a decrease at higher doses. No such decrease was observed in the ungrafted PFA samples after irradiation. This indicated that the changes in the heats of transitions and crystallinity at low doses are due to the radiation effects on the microstructure of PFA (chain scission), whereas at higher doses the grafted PSTY is the driving force behind these changes. (C) 2001 Elsevier Science Ltd. All rights reserved.

Development and implementation of the population Fisher information matrix for the evaluation of population pharmacokinetic designs.

In population pharmacokinetic studies, the precision of parameter estimates is dependent on the population design. Methods based on the Fisher information matrix have been developed and extended to population studies to evaluate and optimize designs. In this paper we propose simple programming tools to evaluate population pharmacokinetic designs. This involved the development of an expression for the Fisher information matrix for nonlinear mixed-effects models, including estimation of the variance of the residual error. We implemented this expression as a generic function for two software applications: S-PLUS and MATLAB. The evaluation of population designs based on two pharmacokinetic examples from the literature is shown to illustrate the efficiency and the simplicity of this theoretic approach. Although no optimization method of the design is provided, these functions can be used to select and compare population designs among a large set of possible designs, avoiding a lot of simulations.

Relationship of glycosaminoglycan and matrix changes to vascular smooth cell phenotype modulation in rabbit arteries after acute injury

Purpose: The phenotype of vascular smooth muscle cells (SMCs) is altered in several arterial pathologies, including the neointima formed after acute arterial injury. This study examined the time course of this phenotypic change in relation to changes in the amount and distribution of matrix glycosaminoglycans. Methods: The immunochemical staining of heparan sulphates (HS) and chondroitin sulphates (CS) in the extracellular matrix of the arterial wall was examined at early points after balloon catheter injury of the rabbit carotid artery. SMC phenotype was assessed by means of ultrastructural morphometry of the cytoplasmic volume fraction of myofilaments. The proportions of cell and matrix components in the media were analyzed with similar morphometric techniques. Results: HS and CS were shown in close association with SMCs of the uninjured arterial media as well as being more widespread within the matrix. Within 6 hours after arterial injury, there was loss of the regular pericellular distribution of both HS and CS, which was associated with a significant expansion in the extracellular space. This preceded the change in ultrastructural phenotype of the SMCs. The glycosaminoglycan loss was most exaggerated at 4 days, after which time the HS and CS reappeared around the medial SMCs. SMCs of the recovering media were able to rapidly replace their glycosaminoglycans, whereas SMCs of the developing neointima failed to produce HS as readily as they produced CS. Conclusions: These studies indicate that changes in glycosaminoglycans of the extracellular matrix precede changes in SMC phenotype after acute arterial injury. In the recovering arterial media, SMCs replace their matrix glycosaminoglycans rapidly, whereas the newly established neointima fails to produce similar amounts of heparan sulphates.

Time-dependent master equation simulation of complex elementary reactions in combustion: Application to the reaction of (CH2)-C-1 with C2H2 from 300-2000 K

Computational simulations of the title reaction are presented, covering a temperature range from 300 to 2000 K. At lower temperatures we find that initial formation of the cyclopropene complex by addition of methylene to acetylene is irreversible, as is the stabilisation process via collisional energy transfer. Product branching between propargyl and the stable isomers is predicted at 300 K as a function of pressure for the first time. At intermediate temperatures (1200 K), complex temporal evolution involving multiple steady states begins to emerge. At high temperatures (2000 K) the timescale for subsequent unimolecular decay of thermalized intermediates begins to impinge on the timescale for reaction of methylene, such that the rate of formation of propargyl product does not admit a simple analysis in terms of a single time-independent rate constant until the methylene supply becomes depleted. Likewise, at the elevated temperatures the thermalized intermediates cannot be regarded as irreversible product channels. Our solution algorithm involves spectral propagation of a symmetrised version of the discretized master equation matrix, and is implemented in a high precision environment which makes hitherto unachievable low-temperature modelling a reality.

Time evolution in the unimolecular master equation at low temperatures: full spectral solution with scalable iterative methods and high precision

A new method is presented to determine an accurate eigendecomposition of difficult low temperature unimolecular master equation problems. Based on a generalisation of the Nesbet method, the new method is capable of achieving complete spectral resolution of the master equation matrix with relative accuracy in the eigenvectors. The method is applied to a test case of the decomposition of ethane at 300 K from a microcanonical initial population with energy transfer modelled by both Ergodic Collision Theory and the exponential-down model. The fact that quadruple precision (16-byte) arithmetic is required irrespective of the eigensolution method used is demonstrated. (C) 2001 Elsevier Science B.V. All rights reserved.

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.

Freeform fabrication of aluminum metal-matrix composites

A series of metal-matrix composites were formed by extrusion freeform, fabrication of a sinterable aluminum alloy in combination with silicon carbide particles and whiskers, carbon fibers, alumina particles, and hollow flyash cenospheres. Silicon carbide particles were most successful in that the composites retained high density with up to 20 vol% of reinforcement and the strength approximately doubles over the strength of the metal matrix alone. Comparison with simple models suggests that this unexpectedly high degree of reinforcement can be attributed to the concentration of small silicon carbide particles around the larger metal powder. This fabrication method also allows composites to be formed with hollow spheres that cannot be formed by other powder or melt methods.

Effects of moist chilling and solid matrix priming on germination of loblolly pine (Pinus taeda L.) seeds

Loblolly pine ( Pinus taeda L.) seeds from sources with a mild climate under maritime influence (North Carolina) required shorter moist chilling to achieve maximum germination vigor than seeds from sources with a harsher continental climate (Oklahoma). Solid matrix priming (SMP) for 6 d achieved as much as 60 d of moist chilling to improve rapidity, synchrony and completeness of germination for three of the four families studied. SMP after moist chilling increased the rapidity, synchrony and completeness of germination. The benefit of SMP was greatest for non-stratified seeds and the benefit decreased with length of moist chilling. In general, delaying planting for one week after SMP had minor effects on germination when seeds were kept in the SMP matrix at 4 degreesC. Delayed planting after SMP can increase germination rapidity and synchrony of seeds that have received long moist chilling and reduce the benefit of SMP in non-moist-chilled seeds.