On erasure coding for distributed storage and streaming communications

The work presented in this thesis revolves around erasure correction coding, as applied to distributed data storage and real-time streaming communications.

First, we examine the problem of allocating a given storage budget over a set of nodes for maximum reliability. The objective is to find an allocation of the budget that maximizes the probability of successful recovery by a data collector accessing a random subset of the nodes. This optimization problem is challenging in general because of its combinatorial nature, despite its simple formulation. We study several variations of the problem, assuming different allocation models and access models, and determine the optimal allocation and the optimal symmetric allocation (in which all nonempty nodes store the same amount of data) for a variety of cases. Although the optimal allocation can have nonintuitive structure and can be difficult to find in general, our results suggest that, as a simple heuristic, reliable storage can be achieved by spreading the budget maximally over all nodes when the budget is large, and spreading it minimally over a few nodes when it is small. Coding would therefore be beneficial in the former case, while uncoded replication would suffice in the latter case.

Second, we study how distributed storage allocations affect the recovery delay in a mobile setting. Specifically, two recovery delay optimization problems are considered for a network of mobile storage nodes: the maximization of the probability of successful recovery by a given deadline, and the minimization of the expected recovery delay. We show that the first problem is closely related to the earlier allocation problem, and solve the second problem completely for the case of symmetric allocations. It turns out that the optimal allocations for the two problems can be quite different. In a simulation study, we evaluated the performance of a simple data dissemination and storage protocol for mobile delay-tolerant networks, and observed that the choice of allocation can have a significant impact on the recovery delay under a variety of scenarios.

Third, we consider a real-time streaming system where messages created at regular time intervals at a source are encoded for transmission to a receiver over a packet erasure link; the receiver must subsequently decode each message within a given delay from its creation time. For erasure models containing a limited number of erasures per coding window, per sliding window, and containing erasure bursts whose maximum length is sufficiently short or long, we show that a time-invariant intrasession code asymptotically achieves the maximum message size among all codes that allow decoding under all admissible erasure patterns. For the bursty erasure model, we also show that diagonally interleaved codes derived from specific systematic block codes are asymptotically optimal over all codes in certain cases. We also study an i.i.d. erasure model in which each transmitted packet is erased independently with the same probability; the objective is to maximize the decoding probability for a given message size. We derive an upper bound on the decoding probability for any time-invariant code, and show that the gap between this bound and the performance of a family of time-invariant intrasession codes is small when the message size and packet erasure probability are small. In a simulation study, these codes performed well against a family of random time-invariant convolutional codes under a number of scenarios.

Finally, we consider the joint problems of routing and caching for named data networking. We propose a backpressure-based policy that employs virtual interest packets to make routing and caching decisions. In a packet-level simulation, the proposed policy outperformed a basic protocol that combines shortest-path routing with least-recently-used (LRU) cache replacement.

Three-dimensional self-assembly of brush block copolymers to photonic crystals

The development of Ring Opening Metathesis Polymerization has allowed the world of block copolymers to expand into brush block copolymers. Brush block copolymers consist of a polymer backbone with polymeric side chains, forcing the backbone to hold a stretched conformation and giving it a worm-like shape. These brush block copolymers have a number of advantages over tradition block copolymers, including faster self-assembly behavior, larger domain sizes, and much less entanglement. This makes them an ideal candidate in the development of a bottom-up approach to forming photonic crystals. Photonic crystals are periodic nanostructures that transmit and reflect only certain wavelengths of light, forming a band gap. These are used in a number of coatings and other optical uses. One and two dimensional photonic crystals are commercially available, though are often expensive and difficult to manufacture. Previous work has focused on the creation of one dimensional photonic crystals from brush block copolymers. In this thesis, I will focus on the synthesis and characterization of asymmetric brush block copolymers for self-assembly into two and three dimensional photonic crystals. Three series of brush block copolymers were made and characterized by Gel Permeation Chromatography and Nuclear Magnetic Resonance spectroscopy. They were then made into films through compressive thermal annealing and characterized by UV-Vis Spectroscopy and Scanning Electron Microscopy. Evidence of non-lamellar structures were seen, indicating the first reported creation of two or three dimensional photonic crystals from brush block copolymers.

Magnetic behavior of basic iron(III) compounds. I. Magnetic behavior of the basic iron sulphates and basic iron chromate. II. Magnetic behavior of delta ferric oxide hydroxide

I. The thermomagnetic behavior and infrared spectroscopic features of KFe₃(SO₄)₂(OH)₆ (jarosite), (H₃O)Fe₃(SO₄)₂ (OH)₆ (hydronium jarosite), KFe₃(CrO₄)₂ (OH)₆, Fe(OH)SO₄ (basic iron sulfate), and Fe(OH)CrO₄ (basic iron chromate) are reported. Fe(OH)CrO₄ and KFe₃(CrO₄)₂ (OH)₆ are shown to be weak ferro magnets with Curie temperatures of 73 and 71 °K, respectively. This unusual magnetic behavior is rationalized in terms of the ionic spin configurations of the phases. Exchange coupling through chromate bridging groups is shown to be weak.

II. The magnetic behavior and the influence of preparative history on the magnetic behavior of δFeO(OH) is reported. δFeO(OH) is shown to be a fine-particulate, uniaxial, magnetic species. Magnetization data for this species are shown to be consistent with the existence of magnetically inactive layers surrounding magnetic particles.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

The interactions of basic proteins and DNA

I. ELECTROPHORESIS OF THE NUCLEIC ACIDS

A zone electrophoresis apparatus using ultraviolet optics has been constructed to study nucleic acids at concentrations less than 0.004%. Native DNA has a mobility about 15% higher than denatured DNA over a range of conditions. Otherwise, the electrophoretic mobility is independent of molecular weight, base composition or source. DNA mobilities change in the expected way with pH but the fractional change in mobility is less than the calculated change in charge. A small decrease in mobility accompanies an increase in ionic strength. RNA’s from various sources have mobilities slightly lower than denatured DNA except for s-RNA which travels slightly faster. The important considerations governing the mobility of nucleic acids appear to be the nature of the hydrodynamic segment, and the binding of counterions. The differences between electrophoresis and sedimentation stem from the fact that all random coil polyelectrolytes are fundamentally free draining in electrophoresis.

II. THE CYTOCHROME C/DNA COMPLEX

The basic protein, cytochrome c, has been complexed to DNA. Up to a cytochrome:DNA mass ratio of 2, a single type of complex is formed. Dissociation of this complex occurs between 0.05F and 0.1F NaCl. The complexing of cytochrome to DNA causes a slight increase in the melting temperature of the DNA, and a reduction of the electrophoretic mobility proportional to the decrease in net charge. Above a cytochrome:DNA mass ratio of 2.5, a different type of complex is formed. The results suggest that complexes such as are formed in the Kleinschmidt technique of electron microscopy would not exist in bulk solution and are exclusively film phenomena.

III. STUDIES OF THE ELECTROPHORESIS AND MELTING BEHAVIOUR OF NUCLEOHISTONES

Electrophoresis studies on reconstituted nucleohistones indicate that the electrophoretic mobility for these complexes is a function of the net charge of the complex. The mobility is therefore dependent on the charge density of the histone complexing the DNA, as well as on the histone/DNA ratio. It is found that the different histones affect the transition from native to denatured DNA in different ways. It appears that histone I is exchanging quite rapidly between DNA molecules in 0.01 F salt, while histone II is irreversibly bound. Histone III-IV enhances the capacity of non-strand separated denatured DNA to reanneal. Studies on native nucleoproteins indicate that there are no gene-sized uncomplexed DNA regions in any preparations studied.

IV. THE DISSOCIATION OF HISTONE FROM CALF THYMUS CROMATIN

Calf thymus nucleoprotein was treated with varying concentrations of NaCl. The identity of the histones associated and dissociated from the DNA at each salt concentration was determined by gel electrophoresis. It was found that there is no appreciable histone dissociation below 0.4 F NaCl. The lysine rich histones dissociate between 0.4 and 0.5 F NaCl. Their dissociation is accompanies by a marked increase in the solubility of the chromatin. The moderately lysine rich histones dissociate mainly between 0.8 and 1.1 F NaCl. There are two arginine rich histone components: the first dissociates between 0.8 F and 1.1 F NaCl, but the second class is the very last to be dissociated from the DNA (dissociation beginning at 1.0 F NaCl). By 2.0 F NaCl, essentially all the histones are dissociated.

The properties of the extracted nucleoprotein were studied. The electrophoretic mobility increases and the melting temperature decreases as more histones are dissociated from the DNA. A comparison with the dissociation of histones from DNA in NaClO₄ shows that to dissociate the same class of histones, the concentration of NaCl required is twice that of NaClO₄.

λ-designs and related combinatorial configurations

This thesis deals with two problems. The first is the determination of λ-designs, combinatorial configurations which are essentially symmetric block designs with the condition that each subset be of the same cardinality negated. We construct an infinite family of such designs from symmetric block designs and obtain some basic results about their structure. These results enable us to solve the problem for λ = 3 and λ = 4. The second problem deals with configurations related to both λ -designs and (ѵ, k, λ)-configurations. We have (n-1) k-subsets of {1, 2, ..., n}, S₁, ..., S_n-1 such that S_i ∩ S_j is a λ-set for i ≠ j. We obtain specifically the replication numbers of such a design in terms of n, k, and λ with one exceptional class which we determine explicitly. In certain special cases we settle the problem entirely.

On the Cesari fixed point method in a Banach space

In a paper published in 1961, L. Cesari [1] introduces a method which extends certain earlier existence theorems of Cesari and Hale ([2] to [6]) for perturbation problems to strictly nonlinear problems. Various authors ([1], [7] to [15]) have now applied this method to nonlinear ordinary and partial differential equations. The basic idea of the method is to use the contraction principle to reduce an infinite-dimensional fixed point problem to a finite-dimensional problem which may be attacked using the methods of fixed point indexes.

The following is my formulation of the Cesari fixed point method:

Let B be a Banach space and let S be a finite-dimensional linear subspace of B. Let P be a projection of B onto S and suppose Г≤B such that pГ is compact and such that for every x in PГ, P^-1x∩Г is closed. Let W be a continuous mapping from Г into B. The Cesari method gives sufficient conditions for the existence of a fixed point of W in Г.

Let I denote the identity mapping in B. Clearly y = Wy for some y in Г if and only if both of the following conditions hold:

(i) Py = PWy.

(ii) y = (P + (I - P)W)y.

Definition. The Cesari fixed paint method applies to (Г, W, P) if and only if the following three conditions are satisfied:

(1) For each x in PГ, P + (I - P)W is a contraction from P^-1x∩Г into itself. Let y(x) be that element (uniqueness follows from the contraction principle) of P^-1x∩Г which satisfies the equation y(x) = Py(x) + (I-P)Wy(x).

(2) The function y just defined is continuous from PГ into B.

(3) There are no fixed points of PWy on the boundary of PГ, so that the (finite- dimensional) fixed point index i(PWy, int PГ) is defined.

Definition. If the Cesari fixed point method applies to (Г, W, P) then define i(Г, W, P) to be the index i(PWy, int PГ).

The three theorems of this thesis can now be easily stated.

Theorem 1 (Cesari). If i(Г, W, P) is defined and i(Г, W, P) ≠0, then there is a fixed point of W in Г.

Theorem 2. Let the Cesari fixed point method apply to both (Г, W, P₁) and (Г, W, P₂). Assume that P₂P₁=P₁P₂=P₁ and assume that either of the following two conditions holds:

(1) For every b in B and every z in the range of P₂, we have that ‖b=P₂b‖ ≤ ‖b-z‖

(2)P₂Г is convex.

Then i(Г, W, P₁) = i(Г, W, P₂).

Theorem 3. If Ω is a bounded open set and W is a compact operator defined on Ω so that the (infinite-dimensional) Leray-Schauder index i_LS(W, Ω) is defined, and if the Cesari fixed point method applies to (Ω, W, P), then i(Ω, W, P) = i_LS(W, Ω).

Theorems 2 and 3 are proved using mainly a homotopy theorem and a reduction theorem for the finite-dimensional and the Leray-Schauder indexes. These and other properties of indexes will be listed before the theorem in which they are used.