Biophysical and network mechanisms of high frequency extracellular potentials in the rat hippocampus

A fundamental question in neuroscience is how distributed networks of neurons communicate and coordinate dynamically and specifically. Several models propose that oscillating local networks can transiently couple to each other through phase-locked firing. Coherent local field potentials (LFP) between synaptically connected regions is often presented as evidence for such coupling. The physiological correlates of LFP signals depend on many anatomical and physiological factors, however, and how the underlying neural processes collectively generate features of different spatiotemporal scales is poorly understood. High frequency oscillations in the hippocampus, including gamma rhythms (30-100 Hz) that are organized by the theta oscillations (5-10 Hz) during active exploration and REM sleep, as well as sharp wave-ripples (SWRs, 140-200 Hz) during immobility or slow wave sleep, have each been associated with various aspects of learning and memory. Deciphering their physiology and functional consequences is crucial to understanding the operation of the hippocampal network.

We investigated the origins and coordination of high frequency LFPs in the hippocampo-entorhinal network using both biophysical models and analyses of large-scale recordings in behaving and sleeping rats. We found that the synchronization of pyramidal cell spikes substantially shapes, or even dominates, the electrical signature of SWRs in area CA1 of the hippocampus. The precise mechanisms coordinating this synchrony are still unresolved, but they appear to also affect CA1 activity during theta oscillations. The input to CA1, which often arrives in the form of gamma-frequency waves of activity from area CA3 and layer 3 of entorhinal cortex (EC3), did not strongly influence the timing of CA1 pyramidal cells. Rather, our data are more consistent with local network interactions governing pyramidal cells' spike timing during the integration of their inputs. Furthermore, the relative timing of input from EC3 and CA3 during the theta cycle matched that found in previous work to engage mechanisms for synapse modification and active dendritic processes. Our work demonstrates how local networks interact with upstream inputs to generate a coordinated hippocampal output during behavior and sleep, in the form of theta-gamma coupling and SWRs.

Electrocatalytic activity of transition metal substituted heteropolytungstates

The electrochemical and electrocatalytic behavior of a series of heteropolytungstate anions in which a tungsten atom in the well known Keggin structure has been replaced by an iron atom is described. All of the iron substituted ions exhibit a one electron reversible couple associated with the Fe³⁺ center and a pair of two electron waves attributed to electron addition and removal from the tungsten oxo framework. The pH and ionic strength effects upon the various electrochemical processes are discussed and interpreted in terms of a competition between protonation and ion pairing of the highly negatively charged ions.

The anions are efficient catalysts for the electroreduction of hydrogen peroxide. A catalytic mechanism involving a formally Fe(IV) intermediate is proposed. Pulse radiolysis experiments were employed to detect the intermediate and evaluate the rate constants for the reactions in which it is formed and decomposed. A chain mechanism for the catalytic decomposition of hydrogen peroxide in which the Fe center shuttles between the +2, +3, and +4 oxidation states is proposed to explain the non-integral stoichiometry observed for the iron substituted polytungstate catalyzed electroreduction of hydrogen peroxide.

The anions are also efficient electrocatalyst for the electrochemical conversion of nitric oxide to ammonia. The catalyzed reduction does not produce hydroxylamine as an intermediate and appears to depend upon the ability of the multiply reduced heteropolytungstates to deliver several electrons to the bound NO group in a concerted step. A valuable feature of the heteropolytungstates is the ease at which the formal potentials of the several redox couples they exhibit may be shifted by changing the identity of the central heteroatom. Exploitation of this feature provided diagnostic information that was decisive in establishing the mechanism of electrocatalytic reduction.

The iron substituted heteropolytungstates are not degraded by repeated cycling between their oxidized and reduced states. They also show superior activity compared to their unsubstituted analogues, indicating that the Fe center acts as a binding site that facilitates inner-sphere electron transfer processes. The basic electrochemistry of several other transition metal substituted Keggin ions is also described.

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.

DNA-mediated charge transport signaling within the cell

DNA possesses the curious ability to conduct charge longitudinally through the π-stacked base pairs that reside within the interior of the double helix. The rate of charge transport (CT) through DNA has a shallow distance dependence. DNA CT can occur over at least 34 nm, a very long molecular distance. Lastly, DNA CT is exquisitely sensitive to disruptions, such as DNA damage, that affect the dynamics of base-pair stacking. Many DNA repair and DNA-processing enzymes are being found to contain 4Fe-4S clusters. These co-factors have been found in glycosylases, helicases, helicase-nucleases, and even enzymes such as DNA polymerase, RNA polymerase, and primase across the phylogeny. The role of these clusters in these enzymes has remained elusive. Generally, iron-sulfur clusters serve redox roles in nature since, formally, the cluster can exist in multiple oxidation states that can be accessed within a biological context. Taken together, these facts were used as a foundation for the hypothesis that DNA-binding proteins with 4Fe-4S clusters utilize DNA-mediated CT as a means to signal one another to scan the genome as a first step in locating the subtle damage that occurs within a sea of undamaged bases within cells.

Herein we describe a role for 4Fe-4S clusters in DNA-mediated charge transport signaling among EndoIII, MutY, and DinG, which are from distinct repair pathways in E. coli. The DinG helicase is an ATP-dependent helicase that contains a 4Fe-4S cluster. To study the DNA-bound redox properties of DinG, DNA-modified electrochemistry was used to show that the 4Fe-4S cluster of DNA-bound DinG is redox-active at cellular potentials, and shares the 80 mV vs. NHE redox potential of EndoIII and MutY. ATP hydrolysis by DinG increases the DNA-mediated redox signal observed electrochemically, likely reflecting better coupling of the 4Fe-4S cluster to DNA while DinG unwinds DNA, which could have interesting biological implications. Atomic force microscopy experiments demonstrate that DinG and EndoIII cooperate at long range using DNA charge transport to redistribute to regions of DNA damage. Genetics experiments, moreover, reveal that this DNA-mediated signaling among proteins also occurs within the cell and, remarkably, is required for cellular viability under conditions of stress. Knocking out DinG in CC104 cells leads to a decrease in MutY activity that is rescued by EndoIII D138A, but not EndoIII Y82A. DinG, thus, appears to help MutY find its substrate using DNA-mediated CT, but do MutY or EndoIII aid DinG in a similar way? The InvA strain of bacteria was used to observe DinG activity, since DinG activity is required within InvA to maintain normal growth. Silencing the gene encoding EndoIII in InvA results in a significant growth defect that is rescued by the overexpression of RNAseH, a protein that dismantles the substrate of DinG, R-loops. This establishes signaling between DinG and EndoIII. Furthermore, rescue of this growth defect by the expression of EndoIII D138A, the catalytically inactive but CT-proficient mutant of EndoIII, is also observed, but expression of EndoIII Y82A, which is CT-deficient but enzymatically active, does not rescue growth. These results provide strong evidence that DinG and EndoIII utilize DNA-mediated signaling to process DNA damage. This work thus expands the scope of DNA-mediated signaling within the cell, as it indicates that DNA-mediated signaling facilitates the activities of DNA repair enzymes across the genome, even for proteins from distinct repair pathways.

In separate work presented here, it is shown that the UvrC protein from E. coli contains a hitherto undiscovered 4Fe-4S cluster. A broad shoulder at 410 nm, characteristic of 4Fe-4S clusters, is observed in the UV-visible absorbance spectrum of UvrC. Electron paramagnetic resonance spectroscopy of UvrC incubated with sodium dithionite, reveals a spectrum with the signature features of a reduced, [4Fe-4S]+1, cluster. DNA-modified electrodes were used to show that UvrC has the same DNA-bound redox potential, of ~80 mV vs. NHE, as EndoIII, DinG, and MutY. Again, this means that these proteins are capable of performing inter-protein electron transfer reactions. Does UvrC use DNA-mediated signaling to facilitate the repair of its substrates?

UvrC is part of the nucleotide excision repair (NER) pathway in E. coli and is the protein within the pathway that performs the chemistry required to repair bulky DNA lesions, such as cyclopyrimidine dimers, that form as a product of UV irradiation. We tested if UvrC utilizes DNA-mediated signaling to facilitate the efficient repair of UV-induced DNA damage products by helping UvrC locate DNA damage. The UV sensitivity of E. coli cells lacking DinG, a putative signaling partner of UvrC, was examined. Knocking out DinG in E. coli leads to a sensitivity of the cells to UV irradiation. A 5-10 fold reduction in the amount of cells that survive after irradiation with 90 J/m2 of UV light is observed. This is consistent with the hypothesis that UvrC and DinG are signaling partners, but is this signaling due to DNA-mediated CT? Complementing the knockout cells with EndoIII D138A, which can also serve as a DNA CT signaling partner, rescues cells lacking DinG from UV irradiation, while complementing the cells with EndoIII Y82A shows no rescue of viability. These results indicate that there is cross-talk between the NER pathway and DinG via DNA-mediated signaling. Perhaps more importantly, this work also establishes that DinG, EndoIII, MutY, and UvrC comprise a signaling network that seems to be unified by the ability of these proteins to perform long range DNA-mediated CT signaling via their 4Fe-4S clusters.

Electronic structure in five-coordinate complexes of nickel(II) with heavy-donor ligands

I.

Various studies designed to elucidate the electronic structure of the arsenic donor ligand, o-phenylenebisdimethylarsine (diarsine), have been carried out. The electronic spectrum of diarsine has been measured at 300 and 77˚K. Electronic spectra of the molecular complexes of various substituted organoarsines and phosphines with tetracyanoethylene have been measured and used to estimate the relative ionization potentials of these molecules.

Uv photolysis of arsines in frozen solution (96˚K) has yielded thermally labile, paramagnetic products. These include the molecular cations of the photolyzed compounds. The species (diars)⁺ exhibits hyper-fine splitting due to two equivalent ⁷⁵As(I=3/2) nuclei. Resonances due to secondary products are reported and assignments discussed.

Evidence is presented for the involvement of d-orbitals in the bonding of arsines. In (diars)⁺ there is mixing of arsenic “lone-pair” orbitals with benzene ring π-orbitals.

II.

Detailed electronic spectral measurements at 300 and 77˚K have been carried out on five-coordinate complexes of low-spin nickel(II), including complexes of both trigonal bipyramidal (TBP) and square pyramidal (SPY) geometry. TBP complexes are of the form NiLX⁺ (X=halide or cyanide,

L = Qƭ(CH₂)₃As(CH₃)₂]₃ or

P [hexagon - Q'CH₃] , Q = P, As,

Q’=S, Se).

The electronic spectra of these compounds exhibit a novel feature at low temperature. The first ligand field band, which is asymmetric in the room temperature solution spectrum, is considerably more symmetrical at 77˚K. This effect is interpreted in terms of changes in the structure of the complex.

The SPY complexes are of the form Ni(diars)₂X^z (X=CL, Br, CNS, CN, thiourea, NO₂, As). On the basis of the spectral results, the d-level ordering is concluded to be xy ˂ xz, yz ˂ z² ˂˂ x² - y². Central to this interpretation is identification of the symmetry-allowed ¹A₁ → ¹E (xz, yz → x² - y²) transition. This assignment was facilitated by the low temperature measurements.

An assignment of the charge-transfer spectra of the five-coordinate complexes is reported, and electronic spectral criteria for distinguishing the two limiting geometries are discussed.

I. Numerical solution of exact pair equations. II. Numerical solution of first-order pair equations

Part I

Numerical solutions to the S-limit equations for the helium ground state and excited triplet state and the hydride ion ground state are obtained with the second and fourth difference approximations. The results for the ground states are superior to previously reported values. The coupled equations resulting from the partial wave expansion of the exact helium atom wavefunction were solved giving accurate S-, P-, D-, F-, and G-limits. The G-limit is -2.90351 a.u. compared to the exact value of the energy of -2.90372 a.u.

Part II

The pair functions which determine the exact first-order wavefunction for the ground state of the three-electron atom are found with the matrix finite difference method. The second- and third-order energies for the (1s1s)¹S, (1s2s)³S, and (1s2s)¹S states of the two-electron atom are presented along with contour and perspective plots of the pair functions. The total energy for the three-electron atom with a nuclear charge Z is found to be E(Z) = -1.125•Z² +1.022805•Z-0.408138-0.025515•(1/Z)+O(1/Z²)a.u.

Supernovae explosions induced by pair production instability

Stars with a core mass greater than about 30 M_⊙ become dynamically unstable due to electron-positron pair production when their central temperature reaches 1.5-2.0 x 10^{9 0}K. The collapse and subsequent explosion of stars with core masses of 45, 52, and 60 M_⊙ is calculated. The range of the final velocity of expansion (3,400 – 8,500 km/sec) and of the mass ejected (1 – 40 M_⊙) is comparable to that observed for type II supernovae.

An implicit scheme of hydrodynamic difference equations (stable for large time steps) used for the calculation of the evolution is described.

For fast evolution the turbulence caused by convective instability does not produce the zero entropy gradient and perfect mixing found for slower evolution. A dynamical model of the convection is derived from the equations of motion and then incorporated into the difference equations.

Experimental study study of antiproton-proton annihilation into a pair of charged Pi-Mesons or K-Mesons for incident antiproton mementa in the range from 0.72 GeV/c to 2.62 GeV/c

The cross sections for the two antiproton-proton annihilation-in-flight modes,

ˉp + p → π⁺ + π^-

ˉp + p → k⁺ + k^-

were measured for fifteen laboratory antiproton beam momenta ranging from 0.72 to 2.62 GeV/c. No magnets were used to determine the charges in the final state. As a result, the angular distributions were obtained in the form [dσ/dΩ (Θ_C.M.) + dσ/dΩ (π – Θ_C.M.)] for 45 ≲ Θ_C.M. ≲ 135°.

A hodoscope-counter system was used to discriminate against events with final states having more than two particles and antiproton-proton elastic scattering events. One spark chamber was used to record the track of each of the two charged final particles. A total of about 40,000 pictures were taken. The events were analyzed by measuring the laboratory angle of the track in each chamber. The value of the square of the mass of the final particles was calculated for each event assuming the reaction

ˉp + p → a pair of particles with equal masses.

About 20,000 events were found to be either annihilation into π ^±-pair or k ^±-pair events. The two different charged meson pair modes were also distinctly separated.

The average differential cross section of ˉp + p → π⁺ + π^- varied from ~ 25 µb/sr at antiproton beam momentum 0.72 GeV/c (total energy in center-of-mass system, √s = 2.0 GeV) to ~ 2 µb/sr at beam momentum 2.62 GeV/c (√s = 2.64 GeV). The most striking feature in the angular distribution was a peak at Θ_C.M. = 90° (cos Θ_C.M. = 0) which increased with √s and reached a maximum at √s ~ 2.1 GeV (beam momentum ~ 1.1 GeV/c). Then it diminished and seemed to disappear completely at √s ~ 2.5 GeV (beam momentum ~ 2.13 GeV/c). A valley in the angular distribution occurred at cos Θ_C.M. ≈ 0.4. The differential cross section then increased as cos Θ_C.M. approached 1.

The average differential cross section for ˉp + p → k⁺ + k^- was about one third of that of the π^±-pair mode throughout the energy range of this experiment. At the lower energies, the angular distribution, unlike that of the π^±-pair mode, was quite isotropic. However, a peak at Θ_C.M. = 90° seemed to develop at √s ~ 2.37 GeV (antiproton beam momentum ~ 1.82 GeV/c). No observable change was seen at that energy in the π^±-pair cross section.

The possible connection of these features with the observed meson resonances at 2.2 GeV and 2.38 GeV, and its implications, were discussed.

I. Site effects and exciton structure in molecular crystals - Benzene. II. The first and second triplet states of solid benzene

The effect of intermolecular coupling in molecular energy levels (electronic and vibrational) has been investigated in neat and isotopic mixed crystals of benzene. In the isotopic mixed crystals of C₆H₆, C₆H₅D, m-C₆H₄D₂, p-C₆H₄D₂, sym-C₆H₃D₃, C₆D₅H, and C₆D₆ in either a C₆H₆ or C₆D₆ host, the following phenomena have been observed and interpreted in terms of a refined Frenkel exciton theory: a) Site shifts; b) site group splittings of the degenerate ground state vibrations of C₆H₆, C₆D₆, and sym-C₆H₃D₃; c) the orientational effect for the isotopes without a trigonal axis in both the ¹B_2u electronic state and the ground state vibrations; d) intrasite Fermi resonance between molecular fundamentals due to the reduced symmetry of the crystal site; and e) intermolecular or intersite Fermi resonance between nearly degenerate states of the host and guest molecules. In the neat crystal experiments on the ground state vibrations it was possible to observe many of these phenomena in conjunction with and in addition to the exciton structure.

To theoretically interpret these diverse experimental data, the concepts of interchange symmetry, the ideal mixed crystal, and site wave functions have been developed and are presented in detail. In the interpretation of the exciton data the relative signs of the intermolecular coupling constants have been emphasized, and in the limit of the ideal mixed crystal a technique is discussed for locating the exciton band center or unobserved exciton components. A differentiation between static and dynamic interactions is made in the Frenkel limit which enables the concepts of site effects and exciton coupling to be sharpened. It is thus possible to treat the crystal induced effects in such a fashion as to make their similarities and differences quite apparent.

A calculation of the ground state vibrational phenomena (site shifts and splittings, orientational effects, and exciton structure) and of the crystal lattice modes has been carried out for these systems. This calculation serves as a test of the approximations of first order Frenkel theory and the atom-atom, pair wise interaction model for the intermolecular potentials. The general form of the potential employed was V(r) = Be^-Cr - A/r⁶ ; the force constants were obtained from the potential by assuming the atoms were undergoing simple harmonic motion.

In part II the location and identification of the benzene first and second triplet states (³B_1u and ³E_1u) is given.

Rydberg states of diatomic and polyatomic molecules using model potentials

A simple model potential is used to calculate Rydberg series for the molecules: nitrogen, oxygen, nitric oxide, carbon monoxide, carbon dioxide, nitrogen dioxide, nitrous oxide, acetylene, formaldehyde, formic acid, diazomethane, ketene, ethylene, allene, acetaldehyde, propyne, acrolein, dimethyl ether, 1, 3-butadiene, 2-butene, and benzene. The model potential for a molecule is taken as the sum of atomic potentials, which are calibrated to atomic data and contain no further parameters. Our results agree with experimentally measured values to within 5-10% in all cases. The results of these calculations are applied to many unresolved problems connected with the above molecules. Some of the more notable of these problems are the reassignment of states in carbon monoxide, the first ionization potential of nitrogen dioxide, the interpretation of the V state in ethylene, and the mystery bands in substituted ethylenes, the identification of the R and R’ series in benzene and the determination of the orbital scheme in benzene from electron impact data.