Topics in randomized numerical linear algebra

This thesis studies three classes of randomized numerical linear algebra algorithms, namely: (i) randomized matrix sparsification algorithms, (ii) low-rank approximation algorithms that use randomized unitary transformations, and (iii) low-rank approximation algorithms for positive-semidefinite (PSD) matrices.

Randomized matrix sparsification algorithms set randomly chosen entries of the input matrix to zero. When the approximant is substituted for the original matrix in computations, its sparsity allows one to employ faster sparsity-exploiting algorithms. This thesis contributes bounds on the approximation error of nonuniform randomized sparsification schemes, measured in the spectral norm and two NP-hard norms that are of interest in computational graph theory and subset selection applications.

Low-rank approximations based on randomized unitary transformations have several desirable properties: they have low communication costs, are amenable to parallel implementation, and exploit the existence of fast transform algorithms. This thesis investigates the tradeoff between the accuracy and cost of generating such approximations. State-of-the-art spectral and Frobenius-norm error bounds are provided.

The last class of algorithms considered are SPSD "sketching" algorithms. Such sketches can be computed faster than approximations based on projecting onto mixtures of the columns of the matrix. The performance of several such sketching schemes is empirically evaluated using a suite of canonical matrices drawn from machine learning and data analysis applications, and a framework is developed for establishing theoretical error bounds.

In addition to studying these algorithms, this thesis extends the Matrix Laplace Transform framework to derive Chernoff and Bernstein inequalities that apply to all the eigenvalues of certain classes of random matrices. These inequalities are used to investigate the behavior of the singular values of a matrix under random sampling, and to derive convergence rates for each individual eigenvalue of a sample covariance matrix.

Regulation of mitochondrial division by the Drp1 receptors

Mitochondria can remodel their membranes by fusing or dividing. These processes are required for the proper development and viability of multicellular organisms. At the cellular level, fusion is important for mitochondrial Ca2+ homeostasis, mitochondrial DNA maintenance, mitochondrial membrane potential, and respiration. Mitochondrial division, which is better known as fission, is important for apoptosis, mitophagy, and for the proper allocation of mitochondria to daughter cells during cellular division.

The functions of proteins involved in fission have been best characterized in the yeast model organism Sarccharomyces cerevisiae. Mitochondrial fission in mammals has some similarities. In both systems, a cytosolic dynamin-like protein, called Dnm1 in yeast and Drp1 in mammals, must be recruited to the mitochondrial surface and polymerized to promote membrane division. Recruitment of yeast Dnm1 requires only one mitochondrial outer membrane protein, named Fis1. Fis1 is conserved in mammals, but its importance for Drp1 recruitment is minor. In mammals, three other receptor proteins—Mff, MiD49, and MiD51—play a major role in recruiting Drp1 to mitochondria. Why mammals require three additional receptors, and whether they function together or separately, are fundamental questions for understanding the mechanism of mitochondrial fission in mammals.

We have determined that Mff, MiD49, or MiD51 can function independently of one another to recruit Drp1 to mitochondria. Fis1 plays a minor role in Drp1 recruitment, suggesting that the emergence of these additional receptors has replaced the system used by yeast. Additionally, we found that Fis1/Mff and the MiDs regulate Drp1 activity differentially. Fis1 and Mff promote constitutive mitochondrial fission, whereas the MiDs activate recruited Drp1 only during loss of respiration.

To better understand the function of the MiDs, we have determined the atomic structure of the cytoplasmic domain of MiD51, and performed a structure-function analysis of MiD49 based on its homology to MiD51. MiD51 adopts a nucleotidyl transferase fold, and binds ADP as a co-factor that is essential for its function. Both MiDs contain a loop segment that is not present in other nucleotidyl transferase proteins, and this loop is used to interact with Drp1 and to recruit it to mitochondria.

Part I. Synchronization of the cell division cycle of HeLa cells in suspension cultures. Part II. Studies on chromosomal proteins of HeLa cells during the cell division cycle

Part I

These studies investigate the potential of single and double treatments with either 5-fluorodeoxyuridine of excess thymidine to induce cell division synchrony in suspension cultures of HeLa cells. The patterns of nucleic acid synthesis and cell proliferation have been analyzed in cultures thus synchronized. Several changes in cell population during long incubation with 5-fluorodeoxyuridine or excess thymidine are also described. These results are subjected to detailed evaluation in terms of the degree and quality of synchrony finally achieved.

Part II

Histones and non-histone proteins associated with interphase and metaphase chromosomes of HeLa cells have been qualitatively and quantitatively analyzed. Histones were fractionated by chromatography on Amberlite CG-50 and further characterized by analytical disc electrophoresis and amino acid analysis of each chromatographic fraction. It is concluded that histones of HeLa cells are comprised of only a small number of major components and that these components are homologous to those of other higher organisms. Of all the histones, arginine-rich histone III alone contains cysteine and can polymerize through formation of intermolecular disulfide bridges between histone III monomers.

A detailed comparison by chromatography and disc electrophoresis established that interphase and metaphase histones are made up of similar components. However, certain quantitative differences in proportions of different histones of interphase and metaphase cells are reported. Indirect evidence indicates that a certain proportion of metaphase histone III is polymerized through intermolecular disulfide links, whereas interphase histone III occurs mainly in the monomeric form.

Metaphase chromosomes are associated with an additional acid-soluble protein fraction which is absent from interphase chromosomes. All of these additional acid-soluble proteins of metaphase chromosomes are shown to be non-histones and it is concluded that the histone/DNA ratio is identical in interphase and metaphase chromosomes. The bulk of acid-soluble non-histone proteins of metaphase chromosomes were found to be polymerized through disulfide bridges; corresponding interphase non-histone proteins displayed no evidence of similar polymerization.

The factors responsible for the condensed configuration and metabolic inactivity of metaphase chromosomes are discussed in light of these findings.

The relationship between histone and DNA synthesis in nondividing differentiated chicken erythrocyte cells and in rapidly dividing undifferentiated HeLa cells is also investigated. Of all the histones, only arginine-rich histones are synthesized in mature erythrocytes. Histone synthesis in HeLa cells was studied in both unsynchronized and synchronized cultures. In HeLa cells, only part of the synthesis of all histone fractions is dependent on concurrent DNA synthesis, whereas all histones are synthesized in varying degrees even in the absence of DNA synthesis.

Structure theorems for noncommutative complete local rings

If R is a ring with identity, let N(R) denote the Jacobson radical of R. R is local if R/N(R) is an artinian simple ring and ∩N(R)ⁱ = 0. It is known that if R is complete in the N(R)-adic topology then R is equal to (B)_n, the full n by n matrix ring over B where E/N(E) is a division ring. The main results of the thesis deal with the structure of such rings B. In fact we have the following.

If B is a complete local algebra over F where B/N(B) is a finite dimensional normal extension of F and N(B) is finitely generated as a left ideal by k elements, then there exist automorphisms g_i,...,g_k of B/N(B) over F such that B is a homomorphic image of B/N[[x₁,…,x_k;g₁,…,g_k]] the power series ring over B/N(B) in noncommuting indeterminates x_i, where x_ib = g_i(b)x_i for all b ϵ B/N.

Another theorem generalizes this result to complete local rings which have suitable commutative subrings. As a corollary of this we have the following. Let B be a complete local ring with B/N(B) a finite field. If N(B) is finitely generated as a left ideal by k elements then there exist automorphisms g₁,…,g_k of a v-ring V such that B is a homomorphic image of V [[x₁,…,x_k;g₁,…,g_k]].

In both these results it is essential to know the structure of N(B) as a two sided module over a suitable subring of B.

The Riesz space structure of an Abelian W*-algebra

Let M be an Abelian W*-algebra of operators on a Hilbert space H. Let M₀ be the set of all linear, closed, densely defined transformations in H which commute with every unitary operator in the commutant M’ of M. A well known result of R. Pallu de Barriere states that if ɸ is a normal positive linear functional on M, then ɸ is of the form T → (Tx, x) for some x in H, where T is in M. An elementary proof of this result is given, using only those properties which are consequences of the fact that ReM is a Dedekind complete Riesz space with plenty of normal integrals. The techniques used lead to a natural construction of the class M₀, and an elementary proof is given of the fact that a positive self-adjoint transformation in M₀ has a unique positive square root in M₀. It is then shown that when the algebraic operations are suitably defined, then M₀ becomes a commutative algebra. If ReM₀ denotes the set of all self-adjoint elements of M₀, then it is proved that ReM₀ is Dedekind complete, universally complete Riesz spaces which contains ReM as an order dense ideal. A generalization of the result of R. Pallu de la Barriere is obtained for the Riesz space ReM₀ which characterizes the normal integrals on the order dense ideals of ReM₀. It is then shown that ReM₀ may be identified with the extended order dual of ReM, and that ReM₀ is perfect in the extended sense.

Some secondary questions related to the Riesz space ReM are also studied. In particular it is shown that ReM is a perfect Riesz space, and that every integral is normal under the assumption that every decomposition of the identity operator has non-measurable cardinal. The presence of atoms in ReM is examined briefly, and it is shown that ReM is finite dimensional if and only if every order bounded linear functional on ReM is a normal integral.

Relativistic quark models and the current algebra

In this thesis we are concerned with finding representations of the algebra of SU(3) vector and axial-vector charge densities at infinite momentum (the "current algebra") to describe the mesons, idealizing the real continua of multiparticle states as a series of discrete resonances of zero width. Such representations would describe the masses and quantum numbers of the mesons, the shapes of their Regge trajectories, their electromagnetic and weak form factors, and (approximately, through the PCAC hypothesis) pion emission or absorption amplitudes.

We assume that the mesons have internal degrees of freedom equivalent to being made of two quarks (one an antiquark) and look for models in which the mass is SU(3)-independent and the current is a sum of contributions from the individual quarks. Requiring that the current algebra, as well as conditions of relativistic invariance, be satisfied turns out to be very restrictive, and, in fact, no model has been found which satisfies all requirements and gives a reasonable mass spectrum. We show that using more general mass and current operators but keeping the same internal degrees of freedom will not make the problem any more solvable. In particular, in order for any two-quark solution to exist it must be possible to solve the "factorized SU(2) problem," in which the currents are isospin currents and are carried by only one of the component quarks (as in the K meson and its excited states).

In the free-quark model the currents at infinite momentum are found using a manifestly covariant formalism and are shown to satisfy the current algebra, but the mass spectrum is unrealistic. We then consider a pair of quarks bound by a potential, finding the current as a power series in 1/m where m is the quark mass. Here it is found impossible to satisfy the algebra and relativistic invariance with the type of potential tried, because the current contributions from the two quarks do not commute with each other to order 1/m³. However, it may be possible to solve the factorized SU(2) problem with this model.

The factorized problem can be solved exactly in the case where all mesons have the same mass, using a covariant formulation in terms of an internal Lorentz group. For a more realistic, nondegenerate mass there is difficulty in covariantly solving even the factorized problem; one model is described which almost works but appears to require particles of spacelike 4-momentum, which seem unphysical.

Although the search for a completely satisfactory model has been unsuccessful, the techniques used here might eventually reveal a working model. There is also a possibility of satisfying a weaker form of the current algebra with existing models.