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.

Matrix Converter-Based Unified Power-Flow Controllers: Advanced Direct Power Control Method

This paper presents a direct power control (DPC) for three-phase matrix converters operating as unified power flow controllers (UPFCs). Matrix converters (MCs) allow the direct ac/ac power conversion without dc energy storage links; therefore, the MC-based UPFC (MC-UPFC) has reduced volume and cost, reduced capacitor power losses, together with higher reliability. Theoretical principles of direct power control (DPC) based on sliding mode control techniques are established for an MC-UPFC dynamic model including the input filter. As a result, line active and reactive power, together with ac supply reactive power, can be directly controlled by selecting an appropriate matrix converter switching state guaranteeing good steady-state and dynamic responses. Experimental results of DPC controllers for MC-UPFC show decoupled active and reactive power control, zero steady-state tracking error, and fast response times. Compared to an MC-UPFC using active and reactive power linear controllers based on a modified Venturini high-frequency PWM modulator, the experimental results of the advanced DPC-MC guarantee faster responses without overshoot and no steady-state error, presenting no cross-coupling in dynamic and steady-state responses.

On the generalized eigenvalue method for energies and matrix elements in lattice field theory

We discuss the generalized eigenvalue problem for computing energies and matrix elements in lattice gauge theory, including effective theories such as HQET. It is analyzed how the extracted effective energies and matrix elements converge when the time separations are made large. This suggests a particularly efficient application of the method for which we can prove that corrections vanish asymptotically as exp(-(E(N+1) - E(n))t). The gap E(N+1) - E(n) can be made large by increasing the number N of interpolating fields in the correlation matrix. We also show how excited state matrix elements can be extracted such that contaminations from all other states disappear exponentially in time. As a demonstration we present numerical results for the extraction of ground state and excited B-meson masses and decay constants in static approximation and to order 1/m(b) in HQET.

Thermal and spectrophotometric studies of new crosslinking method for collagen matrix with glutaraldehyde acetals

Despite the many existing crosslinking procedures, glutaraldehyde (GA) is still the method of choice used in the manufacture of bioprosthesis. The major problems with GA are: (a) uncontrolled reactivity due to the chemical complexity or GA solutions; (b) toxicity due to the release of GA from polymeric crosslinks; and (c) tissue impermeabilization due to polymeric and heterogeneous crosslinks formation, partially responsible for the undesirable calcification of the bioprosthesis. A new method of crosslinking glutaraldehyde acetals has been developed with GA in acid ethanolic solution, and after the distribution inside de matrix, GA is released to crosslinking. Concentrations of hydrochloride acid in ethanolic solutions between 0.1 and 0.001 mol/L with GA concentration between 0.1 and 1.0% were measured in an ultraviolet spectrophotometer to verify the presence of free aldehyde groups and polymeric compounds of GA. After these measurements, the solutions were used to crosslink bovine pericardium. The spectrophotometric results showed that GA was better protected in acetal forms for acid ethanolic solution with HCl at 0.003 mol/L and GA 1.0%(v/v). The shrinkage temperature results of bovine pericardium crosslinked with acetal solutions showed values near 85 C after the exposure to triethylamine vapors.

Pole placement method using the system matrix transfer function and sparsity

This paper presents a pole placement method using both the augmented Jacobian and the corresponding system transfer function matrices. From the manipulation of these matrices a straightforward approach results to get the coefficients of a non-linear system, whose solution gives the parameters of the stabilizers that can provide a pre-specified minimum damping to the system. (C) 2001 Elsevier B.V. Ltd. All rights reserved.

Characterization of silica-carbon mesoporous matrix with embedded nickel nanoparticles synthesized by the polymeric precursor method

Nickel nanoparticles into silica-carbon matrix composites were prepared by using the polymeric precursor method. The effects of the polyester type and the time of pyrolysis on the mesoporosity and nickel particle dispersion into non-aqueous amorphous silica-carbon matrix were investigated by thermogravimetric analysis, adsorption/desorption isotherms and TEM. A well-dispersed metallic phase could be only obtained by using ethylene glycol. Weightier polyesters affected the pyrolysis process due to a combination of more amounts of carbonaceous residues and delaying of pyrolysis process. The post-pyrolyzed composites were successfully cleaned at 200 degrees C for I h in oxygen atmosphere leading to an increase in the surface area and without the occurrence of carbon combustion or nickel nanoparticles oxidation. The matrix composites presented predominantly mesoporous with pore size well defined in 38 angstrom, mainly when tetraethylene glycol was used as polymerizing agent. (C) 2007 Elsevier B.V. All rights reserved.

Measurement of the top quark mass in the lepton plus jets final state with the matrix element method

We present a measurement of the top quark mass with the matrix element method in the lepton+jets final state. As the energy scale for calorimeter jets represents the dominant source of systematic uncertainty, the matrix element likelihood is extended by an additional parameter, which is defined as a global multiplicative factor applied to the standard energy scale. The top quark mass is obtained from a fit that yields the combined statistical and systematic jet energy scale uncertainty. Using a data set of 0.4 fb(-1) taken with the D0 experiment at Run II of the Fermilab Tevatron Collider, the mass of the top quark is measured using topological information to be: m(top)(center dot+jets)(topo)=169.2(-7.4)(+5.0)(stat+JES)(-1.4)(+1.5)(syst) GeV, and when information about identified b jets is included: m(top)(center dot+jets)(b-tag)=170.3(-4.5)(+4.1)(stat+ JES)(-1.8)(+1.2)(syst) GeV. The measurements yield a jet energy scale consistent with the reference scale.

Improvements in the score matrix calculation method using parallel score estimating algorithm

The increasing amount of sequences stored in genomic databases has become unfeasible to the sequential analysis. Then, the parallel computing brought its power to the Bioinformatics through parallel algorithms to align and analyze the sequences, providing improvements mainly in the running time of these algorithms. In many situations, the parallel strategy contributes to reducing the computational complexity of the big problems. This work shows some results obtained by an implementation of a parallel score estimating technique for the score matrix calculation stage, which is the first stage of a progressive multiple sequence alignment. The performance and quality of the parallel score estimating are compared with the results of a dynamic programming approach also implemented in parallel. This comparison shows a significant reduction of running time. Moreover, the quality of the final alignment, using the new strategy, is analyzed and compared with the quality of the approach with dynamic programming.

Search for a standard model Higgs boson produced in association with a top-quark pair and decaying to bottom quarks using a matrix element method

Conselho Nacional de Desenvolvimento Científico e Tecnológico (CNPq)

The use of Chebyshev matrix polynomials in the iterative solution of linear equations compared to the method of successive overrelaxtion /

"(This is being submitted in partial fulfillment of the requirements for the degree of Doctor of Philosophy in Mathematics, June 1959.)"

Full S matrix calculation via a single real-symmetric Lanczos recursion: The Lanczos artificial boundary inhomogeneity method

We present an efficient and robust method for the calculation of all S matrix elements (elastic, inelastic, and reactive) over an arbitrary energy range from a single real-symmetric Lanczos recursion. Our new method transforms the fundamental equations associated with Light's artificial boundary inhomogeneity approach [J. Chem. Phys. 102, 3262 (1995)] from the primary representation (original grid or basis representation of the Hamiltonian or its function) into a single tridiagonal Lanczos representation, thereby affording an iterative version of the original algorithm with greatly superior scaling properties. The method has important advantages over existing iterative quantum dynamical scattering methods: (a) the numerically intensive matrix propagation proceeds with real symmetric algebra, which is inherently more stable than its complex symmetric counterpart; (b) no complex absorbing potential or real damping operator is required, saving much of the exterior grid space which is commonly needed to support these operators and also removing the associated parameter dependence. Test calculations are presented for the collinear H+H-2 reaction, revealing excellent performance characteristics. (C) 2004 American Institute of Physics.

Method for the generation and cultivation of functional three-dimensional mammary constructs without exogenous extracellular matrix

During puberty, pregnancy, lactation and postlactation, breast tissue undergoes extensive remodelling and the disruption of these events can lead to cancer. In vitro studies of mammary tissue and its malignant transformation regularly employ mammary epithelial cells cultivated on matrigel or floating collagen rafts. In these cultures, mammary epithelial cells assemble into three-dimensional structures resembling in vivo acini. We present a novel technique for generating functional mammary constructs without the use of matrix substitutes.

A solvent-free matrix application method for matrix-assisted laser desorption/ionization imaging of small molecules

Matrix application continues to be a critical step in sample preparation for matrix-assisted laser desorption/ionization (MALDI) mass spectrometry imaging (MSI). Imaging of small molecules such as drugs and metabolites is particularly problematic because the commonly used washing steps to remove salts are usually omitted as they may also remove the analyte, and analyte spreading is more likely with conventional wet matrix application methods. We have developed a method which uses the application of matrix as a dry, finely divided powder, here referred to as dry matrix application, for the imaging of drug compounds. This appears to offer a complementary method to wet matrix application for the MALDI-MSI of small molecules, with the alternative matrix application techniques producing different ion profiles, and allows the visualization of compounds not observed using wet matrix application methods. We demonstrate its value in imaging clozapine from rat kidney and 4-bromophenyl-1,4-diazabicyclo(3.2.2)nonane-4-carboxylic acid from rat brain. In addition, exposure of the dry matrix coated sample to a saturated moist atmosphere appears to enhance the visualization of a different set of molecules.

An Iterative Method for the Matrix Principal n-th Root

In this paper we give an iterative method to compute the principal n-th root and the principal inverse n-th root of a given matrix. As we shall show this method is locally convergent. This method is analyzed and its numerical stability is investigated.