Selected configuration interaction with truncation energy error and application to the Ne atom

Selected configuration interaction (SCI) for atomic and molecular electronic structure calculations is reformulated in a general framework encompassing all CI methods. The linked cluster expansion is used as an intermediate device to approximate CI coefficients BK of disconnected configurations (those that can be expressed as products of combinations of singly and doubly excited ones) in terms of CI coefficients of lower-excited configurations where each K is a linear combination of configuration-state-functions (CSFs) over all degenerate elements of K. Disconnected configurations up to sextuply excited ones are selected by Brown's energy formula, ΔEK=(E-HKK)BK2/(1-BK2), with BK determined from coefficients of singly and doubly excited configurations. The truncation energy error from disconnected configurations, Δdis, is approximated by the sum of ΔEKS of all discarded Ks. The remaining (connected) configurations are selected by thresholds based on natural orbital concepts. Given a model CI space M, a usual upper bound ES is computed by CI in a selected space S, and EM=E S+ΔEdis+δE, where δE is a residual error which can be calculated by well-defined sensitivity analyses. An SCI calculation on Ne ground state featuring 1077 orbitals is presented. Convergence to within near spectroscopic accuracy (0.5 cm-1) is achieved in a model space M of 1.4× 109 CSFs (1.1 × 1012 determinants) containing up to quadruply excited CSFs. Accurate energy contributions of quintuples and sextuples in a model space of 6.5 × 1012 CSFs are obtained. The impact of SCI on various orbital methods is discussed. Since ΔEdis can readily be calculated for very large basis sets without the need of a CI calculation, it can be used to estimate the orbital basis incompleteness error. A method for precise and efficient evaluation of ES is taken up in a companion paper

Select-divide-and-conquer method for large-scale configuration interaction

A select-divide-and-conquer variational method to approximate configuration interaction (CI) is presented. Given an orthonormal set made up of occupied orbitals (Hartree-Fock or similar) and suitable correlation orbitals (natural or localized orbitals), a large N-electron target space S is split into subspaces S0,S1,S2,...,SR. S0, of dimension d0, contains all configurations K with attributes (energy contributions, etc.) above thresholds T0={T0egy, T0etc.}; the CI coefficients in S0 remain always free to vary. S1 accommodates KS with attributes above T1≤T0. An eigenproblem of dimension d0+d1 for S0+S 1 is solved first, after which the last d1 rows and columns are contracted into a single row and column, thus freezing the last d1 CI coefficients hereinafter. The process is repeated with successive Sj(j≥2) chosen so that corresponding CI matrices fit random access memory (RAM). Davidson's eigensolver is used R times. The final energy eigenvalue (lowest or excited one) is always above the corresponding exact eigenvalue in S. Threshold values {Tj;j=0, 1, 2,...,R} regulate accuracy; for large-dimensional S, high accuracy requires S 0+S1 to be solved outside RAM. From there on, however, usually a few Davidson iterations in RAM are needed for each step, so that Hamiltonian matrix-element evaluation becomes rate determining. One μhartree accuracy is achieved for an eigenproblem of order 24 × 106, involving 1.2 × 1012 nonzero matrix elements, and 8.4×109 Slater determinants

An assessment of gene prediction accuracy in large DNA sequences

One of the first useful products from the human genome will be a set of predicted genes. Besides its intrinsic scientific interest, the accuracy and completeness of this data set is of considerable importance for human health and medicine. Though progress has been made on computational gene identification in terms of both methods and accuracy evaluation measures, most of the sequence sets in which the programs are tested are short genomic sequences, and there is concern that these accuracy measures may not extrapolate well to larger, more challenging data sets. Given the absence of experimentally verified large genomic data sets, we constructed a semiartificial test set comprising a number of short single-gene genomic sequences with randomly generated intergenic regions. This test set, which should still present an easier problem than real human genomic sequence, mimics the approximately 200kb long BACs being sequenced. In our experiments with these longer genomic sequences, the accuracy of GENSCAN, one of the most accurate ab initio gene prediction programs, dropped significantly, although its sensitivity remained high. Conversely, the accuracy of similarity-based programs, such as GENEWISE, PROCRUSTES, and BLASTX was not affected significantly by the presence of random intergenic sequence, but depended on the strength of the similarity to the protein homolog. As expected, the accuracy dropped if the models were built using more distant homologs, and we were able to quantitatively estimate this decline. However, the specificities of these techniques are still rather good even when the similarity is weak, which is a desirable characteristic for driving expensive follow-up experiments. Our experiments suggest that though gene prediction will improve with every new protein that is discovered and through improvements in the current set of tools, we still have a long way to go before we can decipher the precise exonic structure of every gene in the human genome using purely computational methodology.

Impaired fast-spiking, suppressed cortical inhibition and increased susceptibility to seizures in mice lacking Kv3.2 K+ channel proteins

Voltage-gated K+ channels of the Kv3 subfamily have unusual electrophysiological properties, including activation at very depolarized voltages (positive to −10 mV) and very fast deactivation rates, suggesting special roles in neuronal excitability. In the brain, Kv3 channels are prominently expressed in select neuronal populations, which include fast-spiking (FS) GABAergic interneurons of the neocortex, hippocampus, and caudate, as well as other high-frequency firing neurons. Although evidence points to a key role in high-frequency firing, a definitive understanding of the function of these channels has been hampered by a lack of selective pharmacological tools. We therefore generated mouse lines in which one of the Kv3 genes, Kv3.2, was disrupted by gene-targeting methods. Whole-cell electrophysiological recording showed that the ability to fire spikes at high frequencies was impaired in immunocytochemically identified FS interneurons of deep cortical layers (5-6) in which Kv3.2 proteins are normally prominent. No such impairment was found for FS neurons of superficial layers (2-4) in which Kv3.2 proteins are normally only weakly expressed. These data directly support the hypothesis that Kv3 channels are necessary for high-frequency firing. Moreover, we found that Kv3.2 −/− mice showed specific alterations in their cortical EEG patterns and an increased susceptibility to epileptic seizures consistent with an impairment of cortical inhibitory mechanisms. This implies that, rather than producing hyperexcitability of the inhibitory interneurons, Kv3.2 channel elimination suppresses their activity. These data suggest that normal cortical operations depend on the ability of inhibitory interneurons to generate high-frequency firing.

K(+) channel expression distinguishes subpopulations of parvalbumin- and somatostatin-containing neocortical interneurons

Kv3.1 and Kv3.2 K+ channel proteins form similar voltage-gated K+ channels with unusual properties, including fast activation at voltages positive to −10 mV and very fast deactivation rates. These properties are thought to facilitate sustained high-frequency firing. Kv3.1 subunits are specifically found in fast-spiking, parvalbumin (PV)-containing cortical interneurons, and recent studies have provided support for a crucial role in the generation of the fast-spiking phenotype. Kv3.2 mRNAs are also found in a small subset of neocortical neurons, although the distribution of these neurons is different. We raised antibodies directed against Kv3.2 proteins and used dual-labeling methods to identify the neocortical neurons expressing Kv3.2 proteins and to determine their subcellular localization. Kv3.2 proteins are prominently expressed in patches in somatic and proximal dendritic membrane as well as in axons and presynaptic terminals of GABAergic interneurons. Kv3.2 subunits are found in all PV-containing neurons in deep cortical layers where they probably form heteromultimeric channels with Kv3.1 subunits. In contrast, in superficial layer PV-positive neurons Kv3.2 immunoreactivity is low, but Kv3.1 is still prominently expressed. Because Kv3.1 and Kv3.2 channels are differentially modulated by protein kinases, these results raise the possibility that the fast-spiking properties of superficial- and deep-layer PV neurons are differentially regulated by neuromodulators. Interestingly, Kv3.2 but not Kv3.1 proteins are also prominent in a subset of seemingly non-fast-spiking, somatostatin- and calbindin-containing interneurons, suggesting that the Kv3.1–Kv3.2 current type can have functions other than facilitating high-frequency firing.

Motivational effects of cannabinoids are mediated by μ-opioid and k-opioid receptor

Repeated THC administration produces motivational and somaticadaptive changes leading to dependence in rodents. Toinvestigate the molecular basis for cannabinoid dependenceand its possible relationship with the endogenous opioid system,we explored Δ9-tetrahydrocannabinol (THC) activity in mice lacking μ-, δ- or κ-opioid receptor genes. Acute THCinduced hypothermia, antinociception, and ypolocomotion remained unaffected in these mice, whereas THC tolerance and withdrawal were minimally modified in mutant animals. In contrast, profound phenotypic changes are observed in several place conditioning protocols that reveal both THC rewarding and aversive properties. Absence of μ receptors abolishes THC place preference. Deletion of κ receptors ablates THC place aversion and furthermore unmasks THC place preference. Thus, an opposing activity of μ- and κ-opioid receptors in modulating reward pathways forms the basis for the dual euphoric–dysphoric activity of THC.

K. Iran i la proliferació nuclear

Conté: ¿Irán, potencia nuclear?; El reto nuclear de Irán; Dos respuestas al reto nuclear; Con la vista puesta en Irán; La provocación nuclear norcoreana; Irán y la oferta de Obama

Geographical effects on the accuracy of textile trade data:an international approach for 1913

Foreign trade statistics are the main data source to the study of international trade.However its accuracy has been under suspicion since Morgernstern published hisfamous work in 1963. Federico and Tena (1991) have resumed the question arguing thatthey can be useful in an adequate level of aggregation. But the geographical assignmentproblem remains unsolved. This article focuses on the spatial variable through theanalysis of the reliability of textile international data for 1913. A geographical biasarises between export and import series, but because of its quantitative importance it canbe negligible in an international scale.

On the accuracy of Latin American trade statistics: A nonparametric test for 1925

This paper proposes a nonparametric test in order to establish the level of accuracy of theforeign trade statistics of 17 Latin American countries when contrasted with the trade statistics of the main partners in 1925. The Wilcoxon Matched-Pairs Ranks test is used to determine whether the differences between the data registered by exporters and importers are meaningful, and if so, whether the differences are systematic in any direction. The paper tests for the reliability of the data registered for two homogeneous products, petroleum and coal, both in volume and value. The conclusion of the several exercises performed is that we cannot accept the existence of statistically significant differences between the data provided by the exporters and the registered by the importing countries in most cases. The qualitative historiography of Latin American describes its foreign trade statistics as mostly unusable. Our quantitative results contest this view.

Prediction of subcutaneous fatty acid composition from ham of CLA fed boars and gilts using a NIRS fiber optic probe

The aim of this work was the use of NIR technology by direct application of a fiber optic probe on back fat to analyze the fatty acid composition of CLA fed boars and gilts. 265 animals were fed 3 different diets and the fatty acid profile of back fat from Gluteus medius was analyzed using gas chromatography and FT-NIR. Spectra were acquired using a Bruker Optics Matrix-F duplex spectrometer equipped with a fiber optic probe (IN-268-2). Oleic and stearic fatty acids were predicted accurately; myristic, vaccenic and linoleic fatty acids were predicted with lower accuracy, while palmitic and α-linolenic fatty acids were poorly predicted. The relative percentage of fatty acids and NIR spectra showed differences in fatty acid composition of back fat from pigs fed CLA which increased the relative percentage of SFA and PUFA while MUFA decreased. Results suggest that a NIR fiber optic probe can be used to predict total saturated and unsaturated fatty acid composition, as well as the percentage of stearic and oleic. NIR showed potential as a rapid and easily implemented method to discriminate carcasses from animals fed different diets.

Comparing Random-based and k-Anonymity-Based Algorithms for Graph Anonymization

Recently, several anonymization algorithms have appeared for privacy preservation on graphs. Some of them are based on random-ization techniques and on k-anonymity concepts. We can use both of them to obtain an anonymized graph with a given k-anonymity value. In this paper we compare algorithms based on both techniques in orderto obtain an anonymized graph with a desired k-anonymity value. We want to analyze the complexity of these methods to generate anonymized graphs and the quality of the resulting graphs.

Planar homography: accuracy analysis and applications

Projective homography sits at the heart of many problems in image registration. In addition to many methods for estimating the homography parameters (R.I. Hartley and A. Zisserman, 2000), analytical expressions to assess the accuracy of the transformation parameters have been proposed (A. Criminisi et al., 1999). We show that these expressions provide less accurate bounds than those based on the earlier results of Weng et al. (1989). The discrepancy becomes more critical in applications involving the integration of frame-to-frame homographies and their uncertainties, as in the reconstruction of terrain mosaics and the camera trajectory from flyover imagery. We demonstrate these issues through selected examples

Multivalent display of the antimicrobial peptides BP100 and BP143

Carbohydrates are considered as promising templates for the display of multiple copies of antimicrobial peptides. Herein, wedescribe the design and synthesis of chimeric structures containing two or four copies of the antimicrobial peptidesKKLFKKILKYL-NH2 (BP100) and KKLfKKILKYL-NH2 (BP143) attached to the carbohydrate template cyclodithioerythritol(cDTE) or α-D-galactopyranoside (Galp). The synthesis involved the preparation of the corresponding peptide aldehyde followedby coupling to an aminooxy-functionalized carbohydrate template. After purification, the multivalent display systems were obtainedin high purities (90–98%) and in good yields (42–64%). These compounds were tested against plant and human pathogenic bacteriaand screened for their cytotoxicity on eukaryotic cells. They showed lower MIC values than the parent peptides against the bacteriaanalyzed. In particular, the carbopeptides derived from cDTE and Galp, which contained two or four copies of BP100, respectively,were 2- to 8-fold more active than the monomeric peptide against the phytopathogenic bacteria. These results suggest thatpreassembling antimicrobial peptides to multimeric structures is not always associated with a significant improvement of theactivity. In contrast, the carbopeptides synthesized were active against human red blood cells pointing out that peptide preassemblyis critical for the hemolytic activity. Notably, peptide preassembly resulted in an enhanced bactericidal effect

UBVRI photometry of G, K, M Hipparcos stars

UBVRI data are presented for a set of 229 late-type stars, most of them being high proper motion stars. All these data are part of the Input Catalog planned observations for the Hipparcos mission.