Conjugate-gradient method for computing the Moore-Penrose inverse and rank of a matrix

A conjugate-gradient method is developed for computing the Moore-Penrose generalized inverseA of a matrix and the associated projectors, by using the least-square characteristics of both the method and the inverseA . Two dual algorithms are introduced for computing the least-square and the minimum-norm generalized inverses, as well asA . It is shown that (i) these algorithms converge for any starting approximation; (ii) if they are started from the zero matrix, they converge toA ; and (iii) the trace of a sequence of approximations multiplied byA is a monotone increasing function converging to the rank ofA. A practical way of compensating the self-correcting feature in the computation ofA is devised by using the duality of the algorithms. Comparison with Ben-Israel's method is made through numerical examples. The conjugate-gradient method has an advantage over Ben-Israel's method.After having completed the present paper, the author received from Professor M. R. Hestenes his paper entitledPseudo Inverses and Conjugate Gradients. This paper treated the same subject and appeared in Communications of the ACM, Vol. 18, pp. 40–43, 1975.     

Some model theory for monotone quantifiers

A generalization of the existential and universal quantifier, the monotone quantifiers, are studied. It is shown that the model theory for monotone quantifiers behaves very much like classical model theory. Completeness theorems, definability theorems and preservation theorems are given. Ultraproducts, reduced products and Back and Forth arguments are studied.This work is part of the research program of the G.N.S.A.G. group of the Italian C.N.R.     

Crown compounds for anions. The nature of chemical bonds in the complexes of halide anions with cyclic trimeric perfluoro-o-phenylenemercury and some of its analogs

Geometries and electronic structures of the complexes of halide anions with cyclic trimerico-phenylenemercury, (o-C₆H₄Hg)₃, perfluoro-o-phenylenemercury, (o-C₆F₄Hg)₃, vinylenemercury, (C₂H₂Hg)₃, and perfluorovinylenemercury, (C₂F₂Hg)₃, were modelled by the MNDO method. Calculations were performed for [L-X]^– semisandwich complexes, [X-L-X]^2– bipyramidal complexes, and [L-X-L]^– sandwich complexes (where X=Hal,L is a mercury-containing macrocycle). Based on the results of calculations, we concluded that it was advantageous to describe the chemical bonding between halide anions and mercury-containing macrocycles in terms of generalized chemical bonds, which were successfully used for -complexes of transition metals. In the [L-X]^– semisandwich complexes, the halide anion and the metallacycle are involved in the formation of three generalized chemical bonds: one headlight-shaped -bond and two two-lobe -bonds. In the [X-L-X]^2– bipyramidal complexes, each halide anion forms three generalized chemical bonds with the macrocycle. It is possible because the macrocycleL has unoccupied molecular orbitals suitable for the formation of such bonds; these MOs consist mainly of the orbitals of mercury atoms directed toward both the upper and lower halogen atoms. In the [L-X-L]^– sandwich complexes, the halide anion cannot be bonded to each ringvia three bonds, and, hence, an unsymmetrical structure is formed, in which the rings are located at different distances from the central atom: the [L-X]^– semisandwich complex solvated by macrocycleL.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 6, pp. 1035–1042, June, 1995.The authors are grateful to V. I. Faustov for valuable remarks.This work was supported by the Russian Foundation for Basic Research (Project No. 93-03-18342).     

Net information measures for modified Yukawa and Hulthén potentials

The dimensional analyses of the position and momentum variances‐based quantum mechanical Heisenberg uncertainty measure, as well as the entropic information measures given by the Shannon information entropy sum and the product of Fisher information measures are carried out for two widely used nonrelativistic isotropic exponential‐cosine screened Coulomb potentials generated by multiplying the superpositions of (i) Yukawa‐like, ?Z(e^?μr/r), and (ii) Hulthén‐like, ?Zμ(1/(e^μr ? 1)), potentials by cos(bμr) followed by addition of the term a/r², where a and b ≥ 0, μ are the screening parameters and Z, in case of atoms, denotes the nuclear charge. Under the spherical symmetry, all the information measures considered are shown to be independent of the scaling of the set [μ, Z] at a fixed value of μ/Z, a, and b and the other parameters defining the superpositions of the potentials. Numerical results are presented, which support the validity of the scaling properties. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

An open optimization problem in statistical quality control

We ask the experts in global optimization if there is an efficient solution to an optimization problem in acceptance sampling: Here, one often has incomplete prior information about the quality of incoming lots. Given a cost model, a decision rule for the inspection of a lot may then be designed that minimizes the maximum loss compatible with the available information. The resulting minimax problem is sometimes hard to solve, as the loss functions may have several local maxima which vary in an unpredictable way with the parameters of the decision rule.     

Fractional noise

Fractional noiseN(t),t 0, is a stochastic process for every , and is defined as the fractional derivative or fractional integral of white noise. For = 1 we recover Brownian motion and for = 1/2 we findf ^–1-noise. For 1/2 1, a superposition of fractional noise is related to the fractional diffusion equation.     

Wiener distributions and white noise analysis

The paper describes the structure of a new space of generalized Wiener functionals, , called the Wiener algebra, or space of Wiener distributions, and demonstrates its use in the white noise analysis. The concepts of derivatives and integrals for multi-time parameter generalized stochastic process:^N are introduced, and a derivative version of Itô's lemma is proved. The algebraic structure of and its lattice of subspaces is elaborated, and within this framework a generalized version of the Malliavin calculus is presented.     

Assembly maps,K-theory, and hyperbolic groups

Following Connes and Moscovici, we show that the Baum-Connes assembly map forK _*(C*v) is rationally injective when is word-hyperbolic, implying the Equivariant Novikov conjecture for such groups. Using this result in topologicalK-theory and Borel-Karoubi regulators, we also show that the corresponding generalized assembly map in algebraicK-theory is rationally injective.     

Application of torsion angle molecular dynamics for efficient sampling of protein conformations

We investigate the application of torsion angle molecular dynamics (TAMD) to augment conformational sampling of peptides and proteins. Interesting conformational changes in proteins mainly involve torsional degrees of freedom. Carrying out molecular dynamics in torsion space does not only explicitly sample the most relevant degrees of freedom, but also allows larger integration time steps with elimination of the bond and angle degrees of freedom. However, the covalent geometry needs to be fixed during internal coordinate dynamics, which can introduce severe distortions to the underlying potential surface in the extensively parameterized modern Cartesian-based protein force fields. A "projection" approach (Katritch et al. J Comput Chem 2003, 24, 254-265) is extended to construct an accurate internal coordinate force field (ICFF) from a source Cartesian force field. Torsion crossterm corrections constructed from local molecular fragments, together with softened van der Waals and electrostatic interactions, are used to recover the potential surface and incorporate implicit bond and angle flexibility. MD simulations of dipeptide models demonstrate that full flexibility in both the backbone phi/psi and side chain chi1 angles are virtually restored. The efficacy of TAMD in enhancing conformational sampling is then further examined by folding simulations of small peptides and refinement experiments of protein NMR structures. The results show that an increase of several fold in conformational sampling efficiency can be reliably achieved. The current study also reveals some complicated intrinsic properties of internal coordinate dynamics, beyond energy conservation, that can limit the maximum size of the integration time step and thus the achievable gain in sampling efficiency.     

Classification of extensions of torus algebra II

We use extension theory and algebraic methods to give a complete characterization of extensions of torus algebra by stable Cuntz algebras,and prove certain classification theorems of these extension algebras.