Coupled plasmonic systems and devices : applications in visible metamaterials, nanophotonic circuits, and CMOS imaging

With the size of transistors approaching the sub-nanometer scale and Si-based photonics pinned at the micrometer scale due to the diffraction limit of light, we are unable to easily integrate the high transfer speeds of this comparably bulky technology with the increasingly smaller architecture of state-of-the-art processors. However, we find that we can bridge the gap between these two technologies by directly coupling electrons to photons through the use of dispersive metals in optics. Doing so allows us to access the surface electromagnetic wave excitations that arise at a metal/dielectric interface, a feature which both confines and enhances light in subwavelength dimensions - two promising characteristics for the development of integrated chip technology. This platform is known as plasmonics, and it allows us to design a broad range of complex metal/dielectric systems, all having different nanophotonic responses, but all originating from our ability to engineer the system surface plasmon resonances and interactions. In this thesis, we demonstrate how plasmonics can be used to develop coupled metal-dielectric systems to function as tunable plasmonic hole array color filters for CMOS image sensing, visible metamaterials composed of coupled negative-index plasmonic coaxial waveguides, and programmable plasmonic waveguide network systems to serve as color routers and logic devices at telecommunication wavelengths.

Coupled effects of mechanics, geometry, and chemistry on bio-membrane behavior

Lipid bilayer membranes are models for cell membranes--the structure that helps regulate cell function. Cell membranes are heterogeneous, and the coupling between composition and shape gives rise to complex behaviors that are important to regulation. This thesis seeks to systematically build and analyze complete models to understand the behavior of multi-component membranes.

We propose a model and use it to derive the equilibrium and stability conditions for a general class of closed multi-component biological membranes. Our analysis shows that the critical modes of these membranes have high frequencies, unlike single-component vesicles, and their stability depends on system size, unlike in systems undergoing spinodal decomposition in flat space. An important implication is that small perturbations may nucleate localized but very large deformations. We compare these results with experimental observations.

We also study open membranes to gain insight into long tubular membranes that arise for example in nerve cells. We derive a complete system of equations for open membranes by using the principle of virtual work. Our linear stability analysis predicts that the tubular membranes tend to have coiling shapes if the tension is small, cylindrical shapes if the tension is moderate, and beading shapes if the tension is large. This is consistent with experimental observations reported in the literature in nerve fibers. Further, we provide numerical solutions to the fully nonlinear equilibrium equations in some problems, and show that the observed mode shapes are consistent with those suggested by linear stability. Our work also proves that beadings of nerve fibers can appear purely as a mechanical response of the membrane.

Multiplexed DNA-mediated electrochemistry

The aromatic core of double helical DNA possesses the unique and remarkable ability to form a conduit for electrons to travel over exceptionally long molecular distances. This core of π-stacked nucleobases creates an efficient pathway for charge transfer to proceed that is exquisitely sensitive to even subtle perturbations. Ground state electrochemistry of DNA-modified electrodes has been one of the major techniques used both to investigate and to harness the property of DNA-mediated charge transfer. DNA-modified electrodes have been an essential tool for both gaining insights into the fundamental properties of DNA and, due to the exquisite specificity of DNA-mediated charge transfer for the integrity of the π-stack, for use in next generation diagnostic sensing. Here, multiplexed DNA-modified electrodes are used to (i) gain new insights on the electrochemical coupling of metalloproteins to the DNA π-stack with relevance to the fundaments of in vivo DNA-mediated charge transfer and (ii) enhance the overall sensitivity of DNA-mediated reduction for use in the detection of low abundance diagnostic targets.

First, Methylene Blue (MB′) was covalently attached to DNA through a flexible C12 alkyl linker to yield a new redox reporter for DNA electrochemistry measurements with enhanced sensitivity. Tethered, intercalated MB′ was reduced through DNA-mediated charge transport. The redox signal intensity for MB′-dT-C12-DNA was found to be at least 3 fold larger than that of previously used Nile Blue (NB)-dT-DNA, which is coupled to the base stack via direct conjugation. The signal attenuation, due to an intervening mismatch, and therefore the degree of DNA-mediated reduction, does, however, depend on the DNA film morphology and the backfilling agent used to passivate the surface. These results highlight two possible mechanisms for the reduction of MB′ on the DNA-modified electrode that are distinguishable by their kinetics: reduction mediated by the DNA base pair stack and direct surface reduction of MB′ at the electrode. The extent of direct reduction at the surface can be minimized by overall DNA assembly conditions.

Next, a series of intercalation-based DNA-mediated electrochemical reporters were developed, using a flexible alkane linkage to validate and explore their DNA-mediated reduction. The general mechanism for the reduction of distally bound redox active species, covalently tethered to DNA through flexible alkyl linkages, was established to be an intraduplex DNA-mediated pathway. MB, NB, and anthraquinone were covalently tethered to DNA with three different covalent linkages. The extent of electronic coupling of the reporter was shown to correlate with the DNA binding affinity of the redox active species, supporting an intercalative mechanism. These electrochemical signals were shown to be exceptionally sensitive to a single intervening π-stack perturbation, an AC mismatch, in a densely packed DNA monolayer, which further supports that the reduction is DNA-mediated. Finally, this DNA-mediated reduction of MB occurs primarily via intra- rather than inter duplex intercalation, as probed through varying the proximity and integrity of the neighboring duplex DNA. Further gains to electrochemical sensitivity of our DNA-modified devices were then achieved through the application of electrocatalytic signal amplification using these solvent accessible intercalative reporters, MB-dT-C8, and hemoglobin as a novel electron sink. Electrocatalysis offers an excellent means of electrochemical signal amplification, yet in DNA based sensors, its application has been limited due to strict assembly conditions. We describe the use of hemoglobin as a robust and effective electron sink for electrocatalysis in DNA sensing on low density DNA films. Protein shielding of the heme redox center minimizes direct reduction at the electrode surface and permits assays on low density DNA films. Electrocatalysis of MB that is covalently tethered to the DNA by a flexible alkyl linkage allows for efficient interactions with both the base stack and hemoglobin. Consistent suppression of the redox signal upon incorporation of single CA mismatch in the DNA oligomer demonstrates that both the unamplified and the electrocatalytically amplified redox signals are generated through DNA-mediated charge transport. Electrocatalysis with hemoglobin is robust: it is stable to pH and temperature variations. The utility and applicability of electrocatalysis with hemoglobin is demonstrated through restriction enzyme detection, and an enhancement in sensitivity permits femtomole DNA sampling.

Finally, we expanded the application of our multiplexed DNA-modified electrodes to the electrochemical characterization of DNA-bound proteins containing [4Fe-4S] clusters. DNA-modified electrodes have become an essential tool for the characterization of the redox chemistry of DNA repair proteins that contain redox cofactors. Multiplexed analysis of EndonucleaseIII (EndoIII), a DNA repair protein containing a [4Fe-4S] cluster known to be accessible via DNA-mediated charge transport, elucidated subtle differences in the electrochemical behavior as a function of DNA morphology. DNA-bound EndoIII is seen to have two different electron transfer pathways for reduction, either through the DNA base stack or through direct surface reduction. Closely packed DNA films, where the protein has limited surface accessibility, produce electrochemical signals reflecting electron transfer that is DNA-mediated. The electrochemical comparison of EndoIII mutants, including a new family of mutations altering the electrostatics surrounding the [4Fe-4S] cluster, was able to be quantitatively performed. While little change in the midpoint potential was found for this family of mutants, significant variations in the efficiency of DNA-mediated electron transfer were apparent. Based on the stability of these proteins, examined by circular dichroism, we propose that the electron transfer pathway can be perturbed not only by the removal of aromatic residues, but also through changes in solvation near the cluster.

Chemical-scale studies of G protein-coupled receptors and ligand-gated ion channels

This dissertation describes studies of G protein-coupled receptors (GPCRs) and ligand-gated ion channels (LGICs) using unnatural amino acid mutagenesis to gain high precision insights into the function of these important membrane proteins.

Chapter 2 considers the functional role of highly conserved proline residues within the transmembrane helices of the D2 dopamine GPCR. Through mutagenesis employing unnatural α-hydroxy acids, proline analogs, and N-methyl amino acids, we find that lack of backbone hydrogen bond donor ability is important to proline function. At one proline site we additionally find that a substituent on the proline backbone N is important to receptor function.

In Chapter 3, side chain conformation is probed by mutagenesis of GPCRs and the muscle-type nAChR. Specific side chain rearrangements of highly conserved residues have been proposed to accompany activation of these receptors. These rearrangements were probed using conformationally-biased β-substituted analogs of Trp and Phe and unnatural stereoisomers of Thr and Ile. We also modeled the conformational bias of the unnatural Trp and Phe analogs employed.

Chapters 4 and 5 examine details of ligand binding to nAChRs. Chapter 4 describes a study investigating the importance of hydrogen bonds between ligands and the complementary face of muscle-type and α4β4 nAChRs. A hydrogen bond involving the agonist appears to be important for ligand binding in the muscle-type receptor but not the α4β4 receptor.

Chapter 5 describes a study characterizing the binding of varenicline, an actively prescribed smoking cessation therapeutic, to the α7 nAChR. Additionally, binding interactions to the complementary face of the α7 binding site were examined for a small panel of agonists. We identified side chains important for binding large agonists such as varenicline, but dispensable for binding the small agonist ACh.

Chapter 6 describes efforts to image nAChRs site-specifically modified with a fluorophore by unnatural amino acid mutagenesis. While progress was hampered by high levels of fluorescent background, improvements to sample preparation and alternative strategies for fluorophore incorporation are described.

Chapter 7 describes efforts toward a fluorescence assay for G protein association with a GPCR, with the ultimate goal of probing key protein-protein interactions along the G protein/receptor interface. A wide range of fluorescent protein fusions were generated, expressed in Xenopus oocytes, and evaluated for their ability to associate with each other.

Charge symmetry in ^(13)N and ^(13)C a coupled-channel model

A set of coupled-channel differential equations based on a rotationally distorted optical potential is used to calculate the wave functions required to evaluate the gamma ray transition rate from the first excited state to the ground state in ^(13)C and ^(13)N. The bremsstrahlung differential cross section of low energy protons is also calculated and compared with existing data. The marked similarity between the potentials determined at each resonance level in both nuclei supports the hypothesis of the charge symmetry of nuclear forces by explaining the deviation of the ratios of the experimental E1 transition strengths from unity.

Succinate:ubiquinone oxidoreductase from Paracoccus denitrificans and particulate methane monooxygenase from Methylococcus capsulatus (Bath): experimental and theoretical EPR studies of the metal cofactors

This thesis summarizes the application of conventional and modern electron paramagnetic resonance (EPR) techniques to establish proximity relationships between paramagnetic metal centers in metalloproteins and between metal centers and magnetic ligand nuclei in two important and timely membrane proteins: succinate:ubiquinone oxidoreductase (SQR) from Paracoccus denitrificans and particulate methane monooxygenase (pMMO) from Methylococcus capsulatus. Such proximity relationships are thought to be critical to the biological function and the associated biochemistry mediated by the metal centers in these proteins. A mechanistic understanding of biological function relies heavily on structure-function relationships and the knowledge of how molecular structure and electronic properties of the metal centers influence the reactivity in metalloenzymes. EPR spectroscopy has proven to be one of the most powerful techniques towards obtaining information about interactions between metal centers as well as defining ligand structures. SQR is an electron transport enzyme wherein the substrates, organic and metallic cofactors are held relatively far apart. Here, the proximity relationships of the metallic cofactors were studied through their weak spin-spin interactions by means of EPR power saturation and electron spin-lattice (T_1) measurements, when the enzyme was poised at designated reduction levels. Analysis of the electron T_1 measurements for the S-3 center when the b-heme is paramagnetic led to a detailed analysis of the dipolar interactions and distance determination between two interacting metal centers. Studies of ligand environment of the metal centers by electron spin echo envelope modulation (ESEEM) spectroscopy resulted in the identication of peptide nitrogens as coupled nuclei in the environment of the S-1 and S-3 centers.

Finally, an EPR model was developed to describe the ferromagnetically coupled trinuclear copper clusters in pMMO when the enzyme is oxidized. The Cu(II) ions in these clusters appear to be strongly exchange coupled, and the EPR is consistent with equilateral triangular arrangements of type 2 copper ions. These results offer the first glimpse of the magneto-structural correlations for a trinuclear copper cluster of this type, which, until the work on pMMO, has had no precedent in the metalloprotein literature. Such trinuclear copper clusters are even rare in synthetic models.

Iridium and rhodium analogues of the Shilov cycle catalyst; and the investigation and applications of the Reduction-Coupled Oxo Activation (ROA) mechanistic motif towards alkane upgrading

This dissertation will cover several disparate topics, with the overarching theme centering on the investigation of organometallic C-H activation and hydrocarbon transformation and upgrading. Chapters 2 and 3 discuss iridium and rhodium analogues of the Shilov cycle catalyst for methane to methanol oxidation, and Chapter 4 on the recently discovered ROA mechanistic motif in catalysts for various alkane partial oxidation reactions. In addition, Chapter 5 discusses the mechanism of nickel pyridine bisoxazoline Negishi catalysts for asymmetric and stereoconvergent C-C coupling, and the appendices discuss smaller projects on rhodium H/D exchange catalysts and DFT method benchmarking.

Speculation-Aware Resource Allocation for Cluster Schedulers

Real-time demand response is essential for handling the uncertainties of renewable generation. Traditionally, demand response has been focused on large industrial and commercial loads, however it is expected that a large number of small residential loads such as air conditioners, dish washers, and electric vehicles will also participate in the coming years. The electricity consumption of these smaller loads, which we call deferrable loads, can be shifted over time, and thus be used (in aggregate) to compensate for the random fluctuations in renewable generation.

In this thesis, we propose a real-time distributed deferrable load control algorithm to reduce the variance of aggregate load (load minus renewable generation) by shifting the power consumption of deferrable loads to periods with high renewable generation. The algorithm is model predictive in nature, i.e., at every time step, the algorithm minimizes the expected variance to go with updated predictions. We prove that suboptimality of this model predictive algorithm vanishes as time horizon expands in the average case analysis. Further, we prove strong concentration results on the distribution of the load variance obtained by model predictive deferrable load control. These concentration results highlight that the typical performance of model predictive deferrable load control is tightly concentrated around the average-case performance. Finally, we evaluate the algorithm via trace-based simulations.

Direct simulation of proton-coupled electron transfer reaction dynamics and mechanisms

Proton-coupled electron transfer (PCET) reactions are ubiquitous throughout chemistry and biology. However, challenges arise in both the the experimental and theoretical investigation of PCET reactions; the rare-event nature of the reactions and the coupling between quantum mechanical electron- and proton-transfer with the slower classical dynamics of the surrounding environment necessitates the development of robust simulation methodology. In the following dissertation, novel path-integral based methods are developed and employed for the direct simulation of the reaction dynamics and mechanisms of condensed-phase PCET.

Structure Prediction of G-Protein Coupled Receptors

G-protein coupled receptors (GPCRs) form a large family of proteins and are very important drug targets. They are membrane proteins, which makes computational prediction of their structure challenging. Homology modeling is further complicated by low sequence similarly of the GPCR superfamily.

In this dissertation, we analyze the conserved inter-helical contacts of recently solved crystal structures, and we develop a unified sequence-structural alignment of the GPCR superfamily. We use this method to align 817 human GPCRs, 399 of which are nonolfactory. This alignment can be used to generate high quality homology models for the 817 GPCRs.

To refine the provided GPCR homology models we developed the Trihelix sampling method. We use a multi-scale approach to simplify the problem by treating the transmembrane helices as rigid bodies. In contrast to Monte Carlo structure prediction methods, the Trihelix method does a complete local sampling using discretized coordinates for the transmembrane helices. We validate the method on existing structures and apply it to predict the structure of the lactate receptor, HCAR1. For this receptor, we also build extracellular loops by taking into account constraints from three disulfide bonds. Docking of lactate and 3,5-dihydroxybenzoic acid shows likely involvement of three Arg residues on different transmembrane helices in binding a single ligand molecule.

Protein structure prediction relies on accurate force fields. We next present an effort to improve the quality of charge assignment for large atomic models. In particular, we introduce the formalism of the polarizable charge equilibration scheme (PQEQ) and we describe its implementation in the molecular simulation package Lammps. PQEQ allows fast on the fly charge assignment even for reactive force fields.

Finite differences and a coupled analytic technique with applications to explosions and earthquakes

An analytic technique is developed that couples to finite difference calculations to extend the results to arbitrary distance. Finite differences and the analytic result, a boundary integral called two-dimensional Kirchhoff, are applied to simple models and three seismological problems dealing with data. The simple models include a thorough investigation of the seismologic effects of a deep continental basin. The first problem is explosions at Yucca Flat, in the Nevada test site. By modeling both near-field strong-motion records and teleseismic P-waves simultaneously, it is shown that scattered surface waves are responsible for teleseismic complexity. The second problem deals with explosions at Amchitka Island, Alaska. The near-field seismograms are investigated using a variety of complex structures and sources. The third problem involves regional seismograms of Imperial Valley, California earthquakes recorded at Pasadena, California. The data are shown to contain evidence of deterministic structure, but lack of more direct measurements of the structure and possible three-dimensional effects make two-dimensional modeling of these data difficult.

Computational Predictions of G Protein-Coupled Receptor Structures and Binding Sites

G protein-coupled receptors (GPCRs) are the largest family of proteins within the human genome. They consist of seven transmembrane (TM) helices, with a N-terminal region of varying length and structure on the extracellular side, and a C-terminus on the intracellular side. GPCRs are involved in transmitting extracellular signals to cells, and as such are crucial drug targets. Designing pharmaceuticals to target GPCRs is greatly aided by full-atom structural information of the proteins. In particular, the TM region of GPCRs is where small molecule ligands (much more bioavailable than peptide ligands) typically bind to the receptors. In recent years nearly thirty distinct GPCR TM regions have been crystallized. However, there are more than 1,000 GPCRs, leaving the vast majority of GPCRs with limited structural information. Additionally, GPCRs are known to exist in a myriad of conformational states in the body, rendering the static x-ray crystal structures an incomplete reflection of GPCR structures. In order to obtain an ensemble of GPCR structures, we have developed the GEnSeMBLE procedure to rapidly sample a large number of variations of GPCR helix rotations and tilts. The lowest energy GEnSeMBLE structures are then docked to small molecule ligands and optimized. The GPCR family consists of five subfamilies with little to no sequence homology between them: class A, B1, B2, C, and Frizzled/Taste2. Almost all of the GPCR crystal structures have been of class A GPCRs, and much is known about their conserved interactions and binding sites. In this work we particularly focus on class B1 GPCRs, and aim to understand that family’s interactions and binding sites both to small molecules and their native peptide ligands. Specifically, we predict the full atom structure and peptide binding site of the glucagon-like peptide receptor and the TM region and small molecule binding sites for eight other class B1 GPCRs: CALRL, CRFR1, GIPR, GLR, PACR, PTH1R, VIPR1, and VIPR2. Our class B1 work reveals multiple conserved interactions across the B1 subfamily as well as a consistent small molecule binding site centrally located in the TM bundle. Both the interactions and the binding sites are distinct from those seen in the more well-characterized class A GPCRs, and as such our work provides a strong starting point for drug design targeting class B1 proteins. We also predict the full structure of CXCR4 bound to a small molecule, a class A GPCR that was not closely related to any of the class A GPCRs at the time of the work.

Design of antenna-coupled lumped-element titanium nitride KIDs for long-wavelength multi-band continuum imaging

Many applications in cosmology and astrophysics at millimeter wavelengths including CMB polarization, studies of galaxy clusters using the Sunyaev-Zeldovich effect (SZE), and studies of star formation at high redshift and in our local universe and our galaxy, require large-format arrays of millimeter-wave detectors. Feedhorn and phased-array antenna architectures for receiving mm-wave light present numerous advantages for control of systematics, for simultaneous coverage of both polarizations and/or multiple spectral bands, and for preserving the coherent nature of the incoming light. This enables the application of many traditional "RF" structures such as hybrids, switches, and lumped-element or microstrip band-defining filters.

Simultaneously, kinetic inductance detectors (KIDs) using high-resistivity materials like titanium nitride are an attractive sensor option for large-format arrays because they are highly multiplexable and because they can have sensitivities reaching the condition of background-limited detection. A KID is a LC resonator. Its inductance includes the geometric inductance and kinetic inductance of the inductor in the superconducting phase. A photon absorbed by the superconductor breaks a Cooper pair into normal-state electrons and perturbs its kinetic inductance, rendering it a detector of light. The responsivity of KID is given by the fractional frequency shift of the LC resonator per unit optical power.

However, coupling these types of optical reception elements to KIDs is a challenge because of the impedance mismatch between the microstrip transmission line exiting these architectures and the high resistivity of titanium nitride. Mitigating direct absorption of light through free space coupling to the inductor of KID is another challenge. We present a detailed titanium nitride KID design that addresses these challenges. The KID inductor is capacitively coupled to the microstrip in such a way as to form a lossy termination without creating an impedance mismatch. A parallel plate capacitor design mitigates direct absorption, uses hydrogenated amorphous silicon, and yields acceptable noise. We show that the optimized design can yield expected sensitivities very close to the fundamental limit for a long wavelength imager (LWCam) that covers six spectral bands from 90 to 400 GHz for SZE studies.

Excess phase (frequency) noise has been observed in KID and is very likely caused by two-level systems (TLS) in dielectric materials. The TLS hypothesis is supported by the measured dependence of the noise on resonator internal power and temperature. However, there is still a lack of a unified microscopic theory which can quantitatively model the properties of the TLS noise. In this thesis we derive the noise power spectral density due to the coupling of TLS with phonon bath based on an existing model and compare the theoretical predictions about power and temperature dependences with experimental data. We discuss the limitation of such a model and propose the direction for future study.

Proton-Coupled Reduction of N2 Facilitated by Molecular Fe Complexes

The activation of Fe-coordinated N₂ via the formal addition of hydrogen atom equivalents is explored in this thesis. These reactions may occur in nitrogenase enzymes during the biological conversion of N₂ to NH₃. To understand these reactions, the N₂ reactivity of a series of molecular Fe(N₂) platforms is investigated. A trigonal pyramidal, carbon-ligated Fe^I complex was prepared that displays a similar geometry to that of the resting state 'belt' Fe atoms of nitrogenase. Upon reduction, this species was shown to coordinate N₂, concomitant with significant weakening of the C-Fe interaction. This hemilability of the axial ligand may play a critical role in mediating the interconversion of Fe(N_xH_y) species during N₂ conversion to NH₃. In fact, a trigonal pyramidal borane-ligated Fe complex was shown to catalyze this transformation, generating up to 8.49 equivalents of NH₃. To shed light on the mechanistic details of this reaction, protonation of a borane-ligated Fe(N₂) complex was investigated and found to give rise to a mixture of species that contains an iron hydrazido(2-) [Fe(NNH₂)] complex. The identification of this species is suggestive of an early N-N bond cleavage event en route to NH₃ production, but the highly-reactive nature of this complex frustrated direct attempts to probe this possibility. A structurally-analogous silyl-ligated Fe(N₂) complex was found to react productively with hydrogen atom equivalents, giving rise to an isolable Fe(NNH₂) species. Spectroscopic and crystallographic studies benefited from the enhanced stability of this complex relative to the borane analogue. One-electron reduction of this species initiates a spontaneous disproportionation reaction with an iron hydrazine [Fe(NH₂NH₂)] complex as the predominant reaction product. This transformation provides support for an Fe-mediated N₂ activation mechanism that proceeds via a late N-N bond cleavage. In hopes of gaining more fundamental insight into these reactions, a series of Fe(CN) complexes were prepared and reacted with hydrogen-atom equivalents. Significant quantities of CH₄ and NH₃ are generated in these reactions as a result of complete C-N bond activation. A series of Fe(CNH_x) were found to be exceptionally stable and may be intermediates in these reactions. The stability of these compounds permitted collection of thermodynamic parameters pertinent to the unique N-H bonds. This data is comparatively discussed with the theoretically-predicted data of the N₂-derived Fe(NNH_x) species. Exceptionally-weak N-H bond enthalpies are found for many of these compounds, and sheds light on their short-lived nature and tendency to evolve H₂. As a whole, these works both establish and provide a means to understand Fe-mediated N₂ activation via the addition of hydrogen atom equivalents.