Midlatitude daytime D region ionosphere variations measured from radio atmospherics

We measured the midlatitude daytime ionospheric D region electron density profile height variations in July and August 2005 near Duke University by using radio atmospherics (or sferics for short), which are the high-power, broadband very low frequency (VLF) signals launched by lightning discharges. As expected, the measured daytime D region electron density profile heights showed temporal variations quantitatively correlated with solar zenith angle changes. In the midlatitude geographical regions near Duke University, the observed quiet time heights decreased from ∼80 km near sunrise to ∼71 km near noon when the solar zenith angle was minimum. The measured height quantitative dependence on the solar zenith angle was slightly different from the low-latitude measurement given in a previous work. We also observed unexpected spatial variations not linked to the solar zenith angle on some days, with 15% of days exhibiting regional differences larger than 0.5 km. In these 2 months, 14 days had sudden height drops caused by solar flare X-rays, with a minimum height of 63.4 km observed. The induced height change during a solar flare event was approximately proportional to the logarithm of the X-ray flux. In the long waveband (wavelength, 1-8 Å), an increase in flux by a factor of 10 resulted in 6.3 km decrease of the height at the flux peak time, nearly a perfect agreement with the previous measurement. During the rising and decaying phases of the solar flare, the height changes correlated more consistently with the short, rather than the long, wavelength X-ray flux changes. © 2010 by the American Geophysical Union.

Structural manifestation of the delocalization error of density functional approximations: C(4N+2) rings and C(20) bowl, cage, and ring isomers.

The ground state structure of C(4N+2) rings is believed to exhibit a geometric transition from angle alternation (N < or = 2) to bond alternation (N > 2). All previous density functional theory (DFT) studies on these molecules have failed to reproduce this behavior by predicting either that the transition occurs at too large a ring size, or that the transition leads to a higher symmetry cumulene. Employing the recently proposed perspective of delocalization error within DFT we rationalize this failure of common density functional approximations (DFAs) and present calculations with the rCAM-B3LYP exchange-correlation functional that show an angle-to-bond-alternation transition between C(10) and C(14). The behavior exemplified here manifests itself more generally as the well known tendency of DFAs to bias toward delocalized electron distributions as favored by Huckel aromaticity, of which the C(4N+2) rings provide a quintessential example. Additional examples are the relative energies of the C(20) bowl, cage, and ring isomers; we show that the results from functionals with minimal delocalization error are in good agreement with CCSD(T) results, in contrast to other commonly used DFAs. An unbiased DFT treatment of electron delocalization is a key for reliable prediction of relative stability and hence the structures of complex molecules where many structure stabilization mechanisms exist.

Electron tomography of cryofixed, isometrically contracting insect flight muscle reveals novel actin-myosin interactions.

BACKGROUND: Isometric muscle contraction, where force is generated without muscle shortening, is a molecular traffic jam in which the number of actin-attached motors is maximized and all states of motor action are trapped with consequently high heterogeneity. This heterogeneity is a major limitation to deciphering myosin conformational changes in situ. METHODOLOGY: We used multivariate data analysis to group repeat segments in electron tomograms of isometrically contracting insect flight muscle, mechanically monitored, rapidly frozen, freeze substituted, and thin sectioned. Improved resolution reveals the helical arrangement of F-actin subunits in the thin filament enabling an atomic model to be built into the thin filament density independent of the myosin. Actin-myosin attachments can now be assigned as weak or strong by their motor domain orientation relative to actin. Myosin attachments were quantified everywhere along the thin filament including troponin. Strong binding myosin attachments are found on only four F-actin subunits, the "target zone", situated exactly midway between successive troponin complexes. They show an axial lever arm range of 77°/12.9 nm. The lever arm azimuthal range of strong binding attachments has a highly skewed, 127° range compared with X-ray crystallographic structures. Two types of weak actin attachments are described. One type, found exclusively in the target zone, appears to represent pre-working-stroke intermediates. The other, which contacts tropomyosin rather than actin, is positioned M-ward of the target zone, i.e. the position toward which thin filaments slide during shortening. CONCLUSION: We present a model for the weak to strong transition in the myosin ATPase cycle that incorporates azimuthal movements of the motor domain on actin. Stress/strain in the S2 domain may explain azimuthal lever arm changes in the strong binding attachments. The results support previous conclusions that the weak attachments preceding force generation are very different from strong binding attachments.

Accurate localized resolution of identity approach for linear-scaling hybrid density functionals and for many-body perturbation theory

© 2015 IOP Publishing Ltd and Deutsche Physikalische Gesellschaft.A key component in calculations of exchange and correlation energies is the Coulomb operator, which requires the evaluation of two-electron integrals. For localized basis sets, these four-center integrals are most efficiently evaluated with the resolution of identity (RI) technique, which expands basis-function products in an auxiliary basis. In this work we show the practical applicability of a localized RI-variant ('RI-LVL'), which expands products of basis functions only in the subset of those auxiliary basis functions which are located at the same atoms as the basis functions. We demonstrate the accuracy of RI-LVL for Hartree-Fock calculations, for the PBE0 hybrid density functional, as well as for RPA and MP2 perturbation theory. Molecular test sets used include the S22 set of weakly interacting molecules, the G3 test set, as well as the G2-1 and BH76 test sets, and heavy elements including titanium dioxide, copper and gold clusters. Our RI-LVL implementation paves the way for linear-scaling RI-based hybrid functional calculations for large systems and for all-electron many-body perturbation theory with significantly reduced computational and memory cost.