Conditional independence in quantum many-body systems

In this thesis, I will discuss how information-theoretic arguments can be used to produce sharp bounds in the studies of quantum many-body systems. The main advantage of this approach, as opposed to the conventional field-theoretic argument, is that it depends very little on the precise form of the Hamiltonian. The main idea behind this thesis lies on a number of results concerning the structure of quantum states that are conditionally independent. Depending on the application, some of these statements are generalized to quantum states that are approximately conditionally independent. These structures can be readily used in the studies of gapped quantum many-body systems, especially for the ones in two spatial dimensions. A number of rigorous results are derived, including (i) a universal upper bound for a maximal number of topologically protected states that is expressed in terms of the topological entanglement entropy, (ii) a first-order perturbation bound for the topological entanglement entropy that decays superpolynomially with the size of the subsystem, and (iii) a correlation bound between an arbitrary local operator and a topological operator constructed from a set of local reduced density matrices. I also introduce exactly solvable models supported on a three-dimensional lattice that can be used as a reliable quantum memory.

Relating thermodynamics to information theory : the equality of free energy and mutual information

In this thesis we uncover a new relation which links thermodynamics and information theory. We consider time as a channel and the detailed state of a physical system as a message. As the system evolves with time, ever present noise insures that the "message" is corrupted. Thermodynamic free energy measures the approach of the system toward equilibrium. Information theoretical mutual information measures the loss of memory of initial state. We regard the free energy and the mutual information as operators which map probability distributions over state space to real numbers. In the limit of long times, we show how the free energy operator and the mutual information operator asymptotically attain a very simple relationship to one another. This relationship is founded on the common appearance of entropy in the two operators and on an identity between internal energy and conditional entropy. The use of conditional entropy is what distinguishes our approach from previous efforts to relate thermodynamics and information theory.

Engineering nucleic acid mechanisms for regulation and readout of gene expression : conditional Dicer substrate formation and sensitive multiplexed northern blots

Nucleic acids are most commonly associated with the genetic code, transcription and gene expression. Recently, interest has grown in engineering nucleic acids for biological applications such as controlling or detecting gene expression. The natural presence and functionality of nucleic acids within living organisms coupled with their thermodynamic properties of base-pairing make them ideal for interfacing (and possibly altering) biological systems. We use engineered small conditional RNA or DNA (scRNA, scDNA, respectively) molecules to control and detect gene expression. Three novel systems are presented: two for conditional down-regulation of gene expression via RNA interference (RNAi) and a third system for simultaneous sensitive detection of multiple RNAs using labeled scRNAs.

RNAi is a powerful tool to study genetic circuits by knocking down a gene of interest. RNAi executes the logic: If gene Y is detected, silence gene Y. The fact that detection and silencing are restricted to the same gene means that RNAi is constitutively on. This poses a significant limitation when spatiotemporal control is needed. In this work, we engineered small nucleic acid molecules that execute the logic: If mRNA X is detected, form a Dicer substrate that targets independent mRNA Y for silencing. This is a step towards implementing the logic of conditional RNAi: If gene X is detected, silence gene Y. We use scRNAs and scDNAs to engineer signal transduction cascades that produce an RNAi effector molecule in response to hybridization to a nucleic acid target X. The first mechanism is solely based on hybridization cascades and uses scRNAs to produce a double-stranded RNA (dsRNA) Dicer substrate against target gene Y. The second mechanism is based on hybridization of scDNAs to detect a nucleic acid target and produce a template for transcription of a short hairpin RNA (shRNA) Dicer substrate against target gene Y. Test-tube studies for both mechanisms demonstrate that the output Dicer substrate is produced predominantly in the presence of a correct input target and is cleaved by Dicer to produce a small interfering RNA (siRNA). Both output products can lead to gene knockdown in tissue culture. To date, signal transduction is not observed in cells; possible reasons are explored.

Signal transduction cascades are composed of multiple scRNAs (or scDNAs). The need to study multiple molecules simultaneously has motivated the development of a highly sensitive method for multiplexed northern blots. The core technology of our system is the utilization of a hybridization chain reaction (HCR) of scRNAs as the detection signal for a northern blot. To achieve multiplexing (simultaneous detection of multiple genes), we use fluorescently tagged scRNAs. Moreover, by using radioactive labeling of scRNAs, the system exhibits a five-fold increase, compared to the literature, in detection sensitivity. Sensitive multiplexed northern blot detection provides an avenue for exploring the fate of scRNAs and scDNAs in tissue culture.

Relative mirror symmetry and ramifications of a formula for Gromov-Witten invariants

For a toric Del Pezzo surface S, a new instance of mirror symmetry, said relative, is introduced and developed. On the A-model, this relative mirror symmetry conjecture concerns genus 0 relative Gromov-Witten of maximal tangency of S. These correspond, on the B-model, to relative periods of the mirror to S. Furthermore, for S not necessarily toric, two conjectures for BPS state counts are related. It is proven that the integrality of BPS state counts of the total space of the canonical bundle on S implies the integrality for the relative BPS state counts of S. Finally, a prediction of homological mirror symmetry for the open complement is explored. The B-model prediction is calculated in all cases and matches the known A-model computation for the projective plane.

Reactivity of titanocene methylidene with metal halides, alkene sulfides, and alkene oxides

Titanocene metallacyclobutanes show a wide variety of reactivites with organic and inorganic reagents. Their reactions include methylene transfer to organic carbonyls, formation of enolates, electron transfer from activated alkyl chlorides, olefin metathesis, ring opening polymerization. Recently, preparations of heterobinuclear µ-methylene complexes were reported. In this thesis, mechanistic, synthetic, and structural studies of the heterobinuclear µ-methylene complexes will be described. Also, the reaction of titanocene methylidene trimethylphosphine complex with alkene sulfide and styrene sulfide will be presented.

Heterobinuclear µ-methylene-µ-methyl complexes C_(p2)Ti(µ-CH_2)( µ-CH_3)M(1,5-COD) have been prepared (M = Rh, Ir). X-ray crystallography showed that the methyl group of the complex was bonded to the rhodium and bridges to the titanium through an agostic bond. The ^(1)H,^(13)CNMR, IR spectra along with partial deuteration studies supported the structure in both solution and solid state. Activation of the agostic bond is demonstrated by the equilibration of the µ-CH_3 and µ-CH_2 groups. A nonlinear Arrhenius plot, an unusually large kinetic isotope effect (24(5)), and a large negative activation entropy (-64(3)eu) can be explained by the quantum-mechanical tunneling. Calculated rate constants with Bell-type barrier fitted well with the observed one. This equilibration was best explained by a 4e-4c mechanism (or σ bond metathesis) with the character of quantum-mechanical tunneling.

Heterobinuclear µ-methylene-µ-phenyl complexes were synthesized. Structural study of C_(p2)Ti(µ-CH_(2))(µ-p-Me_(2)NC_(6)H_(4))Rh(l,5-COD) showed that the two metal atoms are bridged by the methylene carbon and the ipso carbon of the p-N,N-dimethylarninophenyl group. The analogous structure of C_(p2))Ti(µ-CH_(2))(µ-o-MeOC_(6)H_(4))Rh(1,5-COD) has been verified by the differential NOE. The aromaticity of the phenyl group observed by ^(1)H NMR, was confirmed by the comparison of the C-C bond lengths in the crystallographic structure. The unusual downfield shifts of the ipso carbon in the ^(13)C NMR are assumed to be an indication of the interaction between the ipso carbon and electron-deficient titanium.

Titanium-platinum heterobinuclear µ-methylene complexes C_(p2)Ti(µ-CH_(2))(µ -X)Pt(Me)(PM_(2)Ph) have been prepared (X= Cl, Me). Structural studies indicate the following:(1) the Ti-CH2 bond possesses residual double bond character, (2) there is a dative Pt→Ti interaction which may be regarded as a π back donation from the platinum atom to the 'Ti=CH_(2)'' group, and (3) the µ-CH_3 group is bound to the titanium atom through a three-center, two-electron agostic bond.

Titanocene (η^(2)-thioformaldehyde)•PMe_3 was prepared from C_(p2)Ti=CH_(2)•PMe_3 and sulfur-containing organic compounds (e.g. alkene sulfide, triphenylphosphine sulfide) including elemental sulfur. Mechanistic studies utilizing trans-styrene sulfide-d_1 suggested the stepwise reaction to explain equimolar mixture of trans- and cis-styrene-d_1 as by-products. The product reacted with methyl iodide to produce cationic titanocene (η_(2)-thiomethoxymethyl) complex. Complexes having less coordinating anion like BF_4 or BPh_4 could be obtained through metathesis. Together with structural analyses, the further reactivities of the complexes have been explored.

The complex C_(p2)TiOCH_(2)CH(Ph)CH_2 was prepared from the compound C_(p2)Ti=CH_(2)-PMe_3 and styrene oxide. The product was characterized with ^(1)H-^(1)H correlated 2-dimensional NMR, selective decoupling of ^(1)H NMR, and differential NOE. Stereospecificity of deuterium in the product was lost when trans-styrene oxide-d_1 was allowed to react. Relative rates of the reaction were measured with varying substituents on the phenyl ring. Better linearity (r = -0.98, p^(+) = -0.79) was observed with σ_(p)^(+)than σ(r = -0.87, p = -1.26). The small magnitude of p^+ value and stereospecificity loss during the formation of product were best explained by the generation of biradicals, but partial generation of charge cannot be excluded. Carbonylation of the product followed by exposure to iodine yields the corresponding β-phenyl γ-lactone.

Towards conditional RNAi : shape and sequence transduction with small conditional RNAs

RNA interference (RNAi) is a powerful biological pathway allowing for sequence-specific knockdown of any gene of interest. While RNAi is a proven tool for probing gene function in biological circuits, it is limited by being constitutively ON and executes the logical operation: silence gene Y. To provide greater control over post-transcriptional gene silencing, we propose engineering a biological logic gate to implement “conditional RNAi.” Such a logic gate would silence gene Y only upon the expression of gene X, a completely unrelated gene, executing the logic: if gene X is transcribed, silence independent gene Y. Silencing of gene Y could be confined to a specific time and/or tissue by appropriately selecting gene X.

To implement the logic of conditional RNAi, we present the design and experimental validation of three nucleic acid self-assembly mechanisms which detect a sub-sequence of mRNA X and produce a Dicer substrate specific to gene Y. We introduce small conditional RNAs (scRNAs) to execute the signal transduction under isothermal conditions. scRNAs are small RNAs which change conformation, leading to both shape and sequence signal transduction, in response to hybridization to an input nucleic acid target. While all three conditional RNAi mechanisms execute the same logical operation, they explore various design alternatives for nucleic acid self-assembly pathways, including the use of duplex and monomer scRNAs, stable versus metastable reactants, multiple methods of nucleation, and 3-way and 4-way branch migration.

We demonstrate the isothermal execution of the conditional RNAi mechanisms in a test tube with recombinant Dicer. These mechanisms execute the logic: if mRNA X is detected, produce a Dicer substrate targeting independent mRNA Y. Only the final Dicer substrate, not the scRNA reactants or intermediates, is efficiently processed by Dicer. Additional work in human whole-cell extracts and a model tissue-culture system delves into both the promise and challenge of implementing conditional RNAi in vivo.

The relative abundances of Sn, Te, Xe, Ba, and Ce in the cosmic radiation

Elements with even atomic number (Z) in the interval 50 ≤ Z ≤ 58 have been resolved in the cosmic radiation using the Heavy Nuclei Experiment on the HEAO-3 satellite. Their relative abundances have been compared with the results expected from pure r-process material, pure s-process material, and solar system material, both with and without a modification due to possible first ionization potential effects. Such effects may be the result of the preferential acceleration, and hence enhancement in the cosmic rays, or those elements having low first ionization potentials. We find that our measurements are inconsistent with pure r-process material at the greater than 98% confidence level whether or not the first ionization potential adjustments are made.

In addition, we have compared our results with mixtures having varying ratios of pure r-process material to pure s-process material. We find that, if no first ionization potential effects are included,

(r/s)CRS/(r/s)SS = 0.20^+0.18_-0.14

where CRS refers to the cosmic ray source and SS refers to the solar system, consistent with having an almost pure s-process source. If the first ionization potential adjustments are applied

(r/s)CRS/(r/s)SS = 1.5^+1.1_-0.7

consistent with a solar system mixture.

Electron transfer in chemically and genetically modified myoglobins

The temperature dependences of the reduction potentials (E^o') of wildtype human myoglobin (Mb) and three site-directed mutants have been measured by using thin-layer spectroelectrochemistry. Residue Val68, which is in van der Waals contact with the heme in Mb, has been replaced by Glu, Asp, and Asn. At pH 7.0, reduction of the heme iron (III) in the former two proteins is accompanied by uptake of a proton by the protein. The changes in E^o', and the standard entropy (ΔS^o') and enthalpy (ΔH^o') of reduction in the mutant proteins were determined relative to values for wild-type; the change in E^o' at 25°C was about -200 millivolts for the Glu and Asp mutants, and about -80 millivolts for the Asn mutant. Reduction of Fe(III) to Fe(II) in the Glu and Asp mutants is accompanied by uptake of a proton. These studies demonstrate that Mb can tolerate substitution of a buried hydrophobic group by potentially charged and polar residues, and that such amino acid replacements can lead to substantial changes in the redox thermodynamics of the protein.

Through analysis of the temperature dependence and shapes of NMR dispersion signals, it is determined that a water molecule is bound to the sixth coordination site of the ferric heme in the Val68Asp and in the Val68Asn recombinant proteins while the carboxyl group of the sidechain of Glu68 occupies this position in Val68Glu. The relative rhombic distortions in the ESR spectra of these mutant proteins combined with H₂¹⁷O and spin interconversion experiments performed on them confirm the conclusions of the NMRD study.

The rates of intramolecular electron transfer (ET) of (NH₃)₅Ru-His48 (Val68Asp, His81GIn, Cys110AIa)Mb and (NH₃)₅Ru-His48 (Val68GIu,His81GIn,Cys110Ala)Mb were measured to be .85(3)s^-1 and .30(2)s^-1, respectively. This data supports the hypothesis that entropy of 111 reduction and reorganization energy of ET are inversely related. The rates of forward and reverse ET for (NH₃)₅ Ru-His48 (Val68GIu, His81 GIn, Cys110AIa)ZnMb -7.2(5)•10⁴s^-1and 1.4(2)•10⁵s^-1, respectively- demonstrate that the placement of a highly polar residue nearby does not significantly change the reorganization energy of the photoactive Zn porphyrin.

The distal histidine imidazoles of (NH₃)₄isnRu-His48 SWMb and (NH₃)₅Ru-His48 SWMb were cyanated with BrCN. The intramolecular ET rates of these BrCN-modified Mb derivatives are 5.5(6)s^-1 and 3.2(5)s^-1, respectively. These respective rates are 20 and 10 times faster than those of their noncyanated counterparts after the differences in ET rate from driving force are scaled according to the Marcus equation. This increase in ET rate of the cyanated Mb derivatives is attributed to lower reorganization energy since the cyanated Mb heme is pentacoordinate in both oxidation states; whereas, the native Mb heme loses a water molecule upon reduction so that it changes from six to five coordinate. The reorganization energy from Fe-OH₂ dissociation is estimated to be .2eV. This conclusion is used to reconcile data from previous experiments in our lab. ET in photoactive porphyrin-substituted myoglobins proceed faster than predicted by Marcus Theory when it is assumed that the only difference in ET parameters between photoactive porphyrins and native heme systems is driving force. However, the data can be consistently fit to Marcus Theory if one corrects for the smaller reorganization in the photoactive porphyrin systems since they do not undergo a coordination change upon ET.

Finally, the intramolecular ET rate of (NH₃)₄isnRu-His48 SWMb was measured to be 3.0(4)s^-1. This rate is within experimental error of that for (NH₃)₄pyrRu-His48 SWMb even though the former has 80mV more driving force. One likely possibility for this observation is that the tetraamminepyridineruthenium group undergoes less reorganization upon ET than the tetraammineisonicotinamideruthenium group. Moreover, analysis of the (NH₃)₄isnRu-His48 SWMb experimental system gives a likely explanation of why ET was not observed previously in (NH₃)₄isnRu-Cytochrome C.

Quantifying Earthquake Collapse Risk of Tall Steel Moment Frame Buildings Using Rupture-to-Rafters Simulations

There is a sparse number of credible source models available from large-magnitude past earthquakes. A stochastic source model generation algorithm thus becomes necessary for robust risk quantification using scenario earthquakes. We present an algorithm that combines the physics of fault ruptures as imaged in laboratory earthquakes with stress estimates on the fault constrained by field observations to generate stochastic source models for large-magnitude (Mw 6.0-8.0) strike-slip earthquakes. The algorithm is validated through a statistical comparison of synthetic ground motion histories from a stochastically generated source model for a magnitude 7.90 earthquake and a kinematic finite-source inversion of an equivalent magnitude past earthquake on a geometrically similar fault. The synthetic dataset comprises of three-component ground motion waveforms, computed at 636 sites in southern California, for ten hypothetical rupture scenarios (five hypocenters, each with two rupture directions) on the southern San Andreas fault. A similar validation exercise is conducted for a magnitude 6.0 earthquake, the lower magnitude limit for the algorithm. Additionally, ground motions from the Mw7.9 earthquake simulations are compared against predictions by the Campbell-Bozorgnia NGA relation as well as the ShakeOut scenario earthquake. The algorithm is then applied to generate fifty source models for a hypothetical magnitude 7.9 earthquake originating at Parkfield, with rupture propagating from north to south (towards Wrightwood), similar to the 1857 Fort Tejon earthquake. Using the spectral element method, three-component ground motion waveforms are computed in the Los Angeles basin for each scenario earthquake and the sensitivity of ground shaking intensity to seismic source parameters (such as the percentage of asperity area relative to the fault area, rupture speed, and risetime) is studied.

Under plausible San Andreas fault earthquakes in the next 30 years, modeled using the stochastic source algorithm, the performance of two 18-story steel moment frame buildings (UBC 1982 and 1997 designs) in southern California is quantified. The approach integrates rupture-to-rafters simulations into the PEER performance based earthquake engineering (PBEE) framework. Using stochastic sources and computational seismic wave propagation, three-component ground motion histories at 636 sites in southern California are generated for sixty scenario earthquakes on the San Andreas fault. The ruptures, with moment magnitudes in the range of 6.0-8.0, are assumed to occur at five locations on the southern section of the fault. Two unilateral rupture propagation directions are considered. The 30-year probabilities of all plausible ruptures in this magnitude range and in that section of the fault, as forecast by the United States Geological Survey, are distributed among these 60 earthquakes based on proximity and moment release. The response of the two 18-story buildings hypothetically located at each of the 636 sites under 3-component shaking from all 60 events is computed using 3-D nonlinear time-history analysis. Using these results, the probability of the structural response exceeding Immediate Occupancy (IO), Life-Safety (LS), and Collapse Prevention (CP) performance levels under San Andreas fault earthquakes over the next thirty years is evaluated.

Furthermore, the conditional and marginal probability distributions of peak ground velocity (PGV) and displacement (PGD) in Los Angeles and surrounding basins due to earthquakes occurring primarily on the mid-section of southern San Andreas fault are determined using Bayesian model class identification. Simulated ground motions at sites within 55-75km from the source from a suite of 60 earthquakes (Mw 6.0 − 8.0) primarily rupturing mid-section of San Andreas fault are considered for PGV and PGD data.

Algebraic Techniques in Coding Theory: Entropy Vectors, Frames, and Constrained Coding

The study of codes, classically motivated by the need to communicate information reliably in the presence of error, has found new life in fields as diverse as network communication, distributed storage of data, and even has connections to the design of linear measurements used in compressive sensing. But in all contexts, a code typically involves exploiting the algebraic or geometric structure underlying an application. In this thesis, we examine several problems in coding theory, and try to gain some insight into the algebraic structure behind them.

The first is the study of the entropy region - the space of all possible vectors of joint entropies which can arise from a set of discrete random variables. Understanding this region is essentially the key to optimizing network codes for a given network. To this end, we employ a group-theoretic method of constructing random variables producing so-called "group-characterizable" entropy vectors, which are capable of approximating any point in the entropy region. We show how small groups can be used to produce entropy vectors which violate the Ingleton inequality, a fundamental bound on entropy vectors arising from the random variables involved in linear network codes. We discuss the suitability of these groups to design codes for networks which could potentially outperform linear coding.

The second topic we discuss is the design of frames with low coherence, closely related to finding spherical codes in which the codewords are unit vectors spaced out around the unit sphere so as to minimize the magnitudes of their mutual inner products. We show how to build frames by selecting a cleverly chosen set of representations of a finite group to produce a "group code" as described by Slepian decades ago. We go on to reinterpret our method as selecting a subset of rows of a group Fourier matrix, allowing us to study and bound our frames' coherences using character theory. We discuss the usefulness of our frames in sparse signal recovery using linear measurements.

The final problem we investigate is that of coding with constraints, most recently motivated by the demand for ways to encode large amounts of data using error-correcting codes so that any small loss can be recovered from a small set of surviving data. Most often, this involves using a systematic linear error-correcting code in which each parity symbol is constrained to be a function of some subset of the message symbols. We derive bounds on the minimum distance of such a code based on its constraints, and characterize when these bounds can be achieved using subcodes of Reed-Solomon codes.