Coulomb-hole summations and energies for GW calculations with limited number of empty orbitals: a modified static remainder approach

Ab initio GW calculations are a standard method for computing the spectroscopic properties of many materials. The most computationally expensive part in conventional implementations of the method is the generation and summation over the large number of empty orbitals required to converge the electron self-energy. We propose a scheme to reduce the summation over empty states by the use of a modified static remainder approximation, which is simple to implement and yields accurate self-energies for both bulk and molecular systems requiring a small fraction of the typical number of empty orbitals.

On the Cubicity of Interval Graphs

A k-cube (or ``a unit cube in k dimensions'') is defined as the Cartesian product R-1 x . . . x R-k where R-i (for 1 <= i <= k) is an interval of the form [a(i), a(i) + 1] on the real line. The k-cube representation of a graph G is a mapping of the vertices of G to k-cubes such that the k-cubes corresponding to two vertices in G have a non-empty intersection if and only if the vertices are adjacent. The cubicity of a graph G, denoted as cub(G), is defined as the minimum dimension k such that G has a k-cube representation. An interval graph is a graph that can be represented as the intersection of intervals on the real line - i. e., the vertices of an interval graph can be mapped to intervals on the real line such that two vertices are adjacent if and only if their corresponding intervals overlap. We show that for any interval graph G with maximum degree Delta, cub(G) <= inverted right perpendicular log(2) Delta inverted left perpendicular + 4. This upper bound is shown to be tight up to an additive constant of 4 by demonstrating interval graphs for which cubicity is equal to inverted right perpendicular log(2) Delta inverted left perpendicular.

Effect of human telomerase reverse transcriptase transfection on differentiation in BeWo choriocarcinoma cells

Arrest of proliferation is one of the prerequisites for differentiation of cytotrophoblasts into syncytiotrophoblasts, and thus during differentiation telomerase activity, as well as human telomerase reverse transcriptase (hTERT) expression, is down-regulated. Considering this, it is of interest to investigate whether syncytium formation can be delayed by prolonging the expression of telomerase in cytotrophoblasts. BeWo cells were transfected with pLPC-hTERT retroviral vector and the reverse transcription-polymerase chain reaction analysis for hTERT mRNA concentrations in the transfected cells revealed a several-fold increase in hTERT mRNA compared with the cells transfected with empty vector, and this confirmed that the transfection was successful. An increase in the proliferation, as assessed by bromodeoxyuridine incorporation assay, as well as an increase in mRNA and protein concentration of various cyclins and proliferating cell nuclear antigen, was noticed. The effect of hTERT transfection was also assessed after the addition of forskolin to induce differentiation and it was observed that cell–cell fusion was delayed and differentiation did not occur in hTERT-transfected cells. However, the effects seen were only transient as stable transfection was not possible and the cells were undergoing apoptosis after 72 h, which suggested that apart from hTERT other factors might be important for immortalization of BeWo cells.

State Dependent Attempt Rate Modeling of Single Cell IEEE 802.11 WLANs with Homogeneous Nodes and Poisson Arrivals

Analytical models of IEEE 802.11-based WLANs are invariably based on approximations, such as the well-known mean-field approximations proposed by Bianchi for saturated nodes. In this paper, we provide a new approach for modeling the situation when the nodes are not saturated. We study a State Dependent Attempt Rate (SDAR) approximation to model M queues (one queue per node) served by the CSMA/CA protocol as standardized in the IEEE 802.11 DCF. The approximation is that, when n of the M queues are non-empty, the attempt probability of the n non-empty nodes is given by the long-term attempt probability of n saturated nodes as provided by Bianchi's model. This yields a coupled queue system. When packets arrive to the M queues according to independent Poisson processes, we provide an exact model for the coupled queue system with SDAR service. The main contribution of this paper is to provide an analysis of the coupled queue process by studying a lower dimensional process and by introducing a certain conditional independence approximation. We show that the numerical results obtained from our finite buffer analysis are in excellent agreement with the corresponding results obtained from ns-2 simulations. We replace the CSMA/CA protocol as implemented in the ns-2 simulator with the SDAR service model to show that the SDAR approximation provides an accurate model for the CSMA/CA protocol. We also report the simulation speed-ups thus obtained by our model-based simulation.

Geometry of the de-Sitter universe

By making use of the fact that the de-Sitter metric corresponds to a hyperquadric in a five-dimensional flat space, it is shown that the three Robertson-Walker metrics for empty spacetime and positive cosmological constant, corresponding to 3-space of positive, negative and zero curvative, are geometrically equivalent. The 3-spaces correspond to intersections of the hyperquadric by hyperplanes, and the time-like geodesics perpendicular to them correspond to intersections by planes, in all three cases.

Successful recovery of preimplantation embryos by nonsurgical uterine flushing in the bonnet monkey

A major limitation to progress in primate embryology is the lack of an adequate supply of preimplantation embryos. We describe a method for recovering preimplantation-embryos in bonnet monkeys (Macaca radiata ) using a nonsurgical uterine flushing technique similar to the one previously employed in rhesus monkeys. Forty cyclic females were screened for cervical cannulation, and 10% of these had an impassable cervix. Eleven females suitable for cannulation were selected, and 27 menstrual cycles were monitored over a 5-mo period. Seventy-one percent of the cycles showed estrogen peaks, which were observed between Days 9 and 14 of the cycle. Following natural mating, uterine flushings were performed on Days 5 to 8 of pregnancy (Day 0 = the day following the estrogen peak). Of the 27 recovery attempts, 9 (33.3%) resulted in the recovery of ovulation products, including those of an unfertilized oocyte and empty zona (2 cases), retarded cleavage-stage (4 to 8-cell) embryos (4 cases), morula (1 case) and blastocysts (2 cases). These results show, for the first time, that the nonsurgical uterine flushing technique can be successfully performed to recover uterine-stage preimplantation embryos from bonnet monkeys.

Assembly of physalis mottle virus capsid protein in Escherichia coli and the role of amino and carboxy termini in the formation of the icosahedral particles

The coat protein gene of physalis mottle tymovirus (PhMV) was over expressed in Escherichia coli using pET-3d vector. The recombinant protein was found to self assemble into capsids in vivo. The purified recombinant capsids had an apparent s value of 56.5 S and a diameter of 29(±2) nm. In order to establish the role of amino and carboxy-terminal regions in capsid assembly, two amino-terminal deletions clones lacking the first 11 and 26 amino acid residues and two carboxy-terminal deletions lacking the last five and ten amino acid residues were constructed and overexpressed. The proteins lacking N-terminal 11 (PhCPN1) and 26 (PhCPN2) amino acid residues self assembled into T = 3 capsids in vivo, as evident from electron microscopy, ultracentrifugation and agarose gel electrophoresis. The recombinant, PhCPN1 and PhCPN2 capsids were as stable as the empty capsids formed in vivo and encapsidated a small amount of mRNA. The monoclonal antibody PA3B2, which recognizes the epitope within region 22 to 36, failed to react with PhCPN2 capsids while it recognized the recombinant and PhCPN1 capsids. Disassembly of the capsids upon treatment with urea showed that PhCPN2 capsids were most stable. These results demonstrate that the N-terminal 26 amino acid residues are not essential for T = 3 capsid assembly in PhMV. In contrast, both the proteins lacking the C-terminal five and ten amino acid residues were present only in the insoluble fraction and could not assemble into capsids, suggesting that these residues are crucial for folding and assembly of the particles.

Popular Mixed Matchings

We study the problem of matching applicants to jobs under one-sided preferences: that is, each applicant ranks a non-empty subset of jobs under an order of preference, possibly involving ties. A matching M is said to be rnore popular than T if the applicants that prefer M to T outnumber those that prefer T to M. A matching is said to be popular if there is no matching more popular than it. Equivalently, a matching M is popular if phi(M,T) >= phi(T, M) for all matchings T, where phi(X, Y) is the number of applicants that prefer X to Y. Previously studied solution concepts based oil the popularity criterion are either not guaranteed to exist for every instance (e.g., popular matchings) or are NP-hard to compute (e.g., least unpopular matchings). This paper addresses this issue by considering mixed matchings. A mixed matching is simply a probability distributions over matchings in the input graph. The function phi that compares two matchings generalizes in a natural manner to mixed matchings by taking expectation. A mixed matching P is popular if phi(P,Q) >= phi(Q,P) for all mixed matchings Q. We show that popular mixed matchings always exist. and we design polynomial time algorithms for finding them. Then we study their efficiency and give tight bounds on the price of anarchy and price of stability of the popular matching problem.

Electronic Structure and Bonding in Neutral and Dianionic Boradiphospholes: R ' BC2P2R2 (R=H, tBu, R '=H, Ph)

Classical and non-classical isomers of both neutral and dianionic BC2P2H3 species, which are isolobal to Cp+ and Cp-, are studied at both B3LYP/6-311++G(d,p) and G3B3 levels of theory. The global minimum structure given by B3LYP/6-311+ + G(d,p) for BC2P2H3 is based on a vinylcyclopropenyl-type structure, whereas BC2P2H32- has a planar aromatic cyclopentadienyl-ion-like structure. However, at the G3B3 level, there are three low-energy isomers for BC2P2H3: 1)tricyclopentane, 2) nido and 3) vinylcyclopropenyl-type structures, all within 1.7 kcal mol(-1) of each other. On the contrary, for the dianionic species the cyclic planar structure is still the minimum. In comparison to the isolobal Cp+ and HnCnP5-n+ isomers, BC2P2H3 shows a competition between pi-delocalised vinylcyclopropenyl- and cluster-type structures (nido and tricyclopentane). Substitution of H on C by tBu, and H on B by Ph, in BC2P2H3 increases the energy difference between the low-lying isomers, giving the lowest energy structure as a tricyclopentane type. Similar substitution in BC2P2H32- merely favours different positional isomers of the cyclic planar geometry, as observed in 1) isoelectronic neutral heterodiphospholes EtBu2C2P2 (E=S, Se, Te), 2) monoanionic heterophospholyl rings EtBu2C2P2 (E=P-, As-, Sb-) and 3) polyphospholyl rings anions tBu(5-n)C(n)P(5-n) (n=0-5). The principal factors that affect the stability of three-, four-, and five-membered ring and acyclic geometrical and positional isomers of neutral and dianionic BC2P2H3 isomers appear to be: 1) relative bond strengths, 2) availability of electrons for the empty 2p boron orbital and 3) steric effects of the tBu groups in the HBC(2)P(2)tBu(2) systems.

On the Cubicity of AT-Free Graphs and Circular-Arc Graphs

A unit cube in k dimensions (k-cube) is defined as the Cartesian product R-1 x R-2 x ... x R-k where R-i (for 1 <= i <= k) is a closed interval of the form [a(i), a(i) + 1] on the real line. A graph G on n nodes is said to be representable as the intersection of k-cubes (cube representation in k dimensions) if each vertex of C can be mapped to a k-cube such that two vertices are adjacent in G if and only if their corresponding k-cubes have a non-empty intersection. The cubicity of G denoted as cub(G) is the minimum k for which G can be represented as the intersection of k-cubes. An interesting aspect about cubicity is that many problems known to be NP-complete for general graphs have polynomial time deterministic algorithms or have good approximation ratios in graphs of low cubicity. In most of these algorithms, computing a low dimensional cube representation of the given graph is usually the first step. We give an O(bw . n) algorithm to compute the cube representation of a general graph G in bw + 1 dimensions given a bandwidth ordering of the vertices of G, where bw is the bandwidth of G. As a consequence, we get O(Delta) upper bounds on the cubicity of many well-known graph classes such as AT-free graphs, circular-arc graphs and cocomparability graphs which have O(Delta) bandwidth. Thus we have: 1. cub(G) <= 3 Delta - 1, if G is an AT-free graph. 2. cub(G) <= 2 Delta + 1, if G is a circular-arc graph. 3. cub(G) <= 2 Delta, if G is a cocomparability graph. Also for these graph classes, there axe constant factor approximation algorithms for bandwidth computation that generate orderings of vertices with O(Delta) width. We can thus generate the cube representation of such graphs in O(Delta) dimensions in polynomial time.

Hydrophobic Peptide Channels and Encapsulated Water Wires

Peptide nanotubes with filled and empty pores and close-packed structures are formed in closely related pentapeptides. Enantiomorphic sequences, Boc-(D)Pro-Aib-Xxx-Aib-Val-OMe (Xxx = Leu, 1; Val, 2; Ala, 3; Phe, 4) and Boc-Pro-Aib-(D)Xxx-Aib-(D)Val-OMe ((XXX)-X-D = (D)Leu, 5; (D)Val, 6; (D)Ala, 7; (D)Phe, 8), yield molecular structures with a very similar backbone conformation but varied packing patterns in crystals. Peptides 1, 2, 5, and 6 show tubular structures with the molecules self-assembling along the crystallographic six-fold axis (c-axis) and revealing a honeycomb arrangement laterally (ab plane). Two forms of entrapped water wires have been characterized in 2: 2a with d(O center dot center dot center dot O) = 2.6 angstrom and 2b with d(O center dot center dot center dot O) = 3.5 angstrom. The latter is observed in 6 (6a) also. A polymorphic form of 6 (6b), grown from a solution of methanol-water, was observed to crystallize in a monoclinic system as a close-packed structure. Single-file water wire arrangements encapsulated inside hydrophobic channels formed by peptide nanotubes could be established by modeling the published structures in the cases of a cyclic peptide and a dipeptide. In all the entrapped water wires, each water molecule is involved in a hydrogen bond with a previous and succeeding water molecule. The O-H group of the water not involved in any hydrogen bond does not seem to be involved in an energetically significant interaction with the nanotube interior, a general feature of the one-dimensional water wires encapsulated in hydrophobic environements. Water wires in hydrophobic channels are contrasted with the single-file arrangements in amphipathic channels formed by aquaporins.

The core and shapley value analysis for cooperative formation of procurement networks

Formation of high value procurement networks involves a bottom-up assembly of complex production, assembly, and exchange relationships through supplier selection and contracting decisions, where suppliers are intelligent and rational agents who act strategically. In this paper we address the problem of forming procurement networks for items with value adding stages that are linearly arranged We model the problem of Procurement Network Formation (PNF) for multiple units of a single item as a cooperative game where agents cooperate to form a surplus maximizing procurement network and then share the surplus in a stable and fair manner We first investigate the stability of such networks by examining the conditions under which the core of the game is non-empty. We then present a protocol, based on the extensive form game realization of the core, for forming such networks so that the resulting network is stable. We also mention a key result when the Shapley value is applied as a solution concept.

Structure and assembly of Sesbania mosaic virus

Sesbania mosaic virus (SeMV) is a ss-RNA (4149 nt) plant sobemovirus isolated from farmer's field around Tirupathi, Andhra Pradesh. The viral capsid (30 nm diameter) consists of 180 copies of protein subunits (MW 29 kDa) organized with icosahedral symmetry. In order to understand the mechanism of assembly of SeMV, a large number of deletion and substitution mutants of the coat protein (CP) were constructed. Recombinant SeMV CP (rCP) as well as the N-terminal rCP deletion mutant Delta N22 were found to assemble in E. coli into virus-like particles (VLPs). Delta N36 and Delta N65 mostly formed smaller particles consisting of 60 protein subunits. Although particlem assembly was not affected due to the substitution of aspartates (D14 and D149) that coordinate calcium ions by asparagines, the stability of the resulting capsids was drastically reduced. Deletion of residues forming a characteristic beta-annulus at the icosahedral 3-folds did not affect the assembly of VLPs. Mutation of a single tryptophan, which occurs near the icosahedral fivefold axis to glutamate or lysine, resulted in the disruption of the capsid leading to soluble dimers that resembled the quasi-dimer structure of the native virus. Replacement of positively charged residues in the amino terminal segment of CP resulted in the formation of empty shells. Based on these observations, a plausible mechanism of assembly is proposed.

Terahertz wave characteristics of a single-walled carbon nanotube containing a fluid flow using the nonlocal Timoshenko beam model

This paper presents the effect of nonlocal scaling parameter on the terahertz wave propagation in fluid filled single walled carbon nanotubes (SWCNTs). The SWCNT is modeled as a Timoshenko beam,including rotary inertia and transverse shear deformation by considering the nonlocal scale effects. A uniform fluid velocity of 1000 m/s is assumed. The analysis shows that, for a fluid filled SWCNT, the wavenumbers of flexural and shear waves will increase and the corresponding wave speeds will decrease as compared to an empty SWCNT. The nonlocal scale parameter introduces certain band gap region in both flexural and shear wave mode where no wave propagation occurs. This is manifested in the wavenumber plots as the region where the wavenumber tends to infinite (or wave speed tends to zero). The frequency at which this phenomenon occurs is called the ``escape frequency''. The effect of fluid density on the terahertz wave propagation in SWCNT is also studied and the analysis shows that as the fluid becomes denser, the wave speeds will decrease. The escape frequency decreases with increase in nonlocal scaling parameter, for both wave modes. We also show that the effect of fluid density and velocity are negligible on the escape frequencies of flexural and shear wave modes. (C) 2010 Elsevier B.V. All rights reserved.

Modular design of synthetic protein mimics. Crystal structures, assembly, and hydration of two 15- and 16-residue apolar, leucyl-rich helical peptides

Two IS- and 16-residue peptides containing a-aminoisobutyric acid (Aib) have been synthesized, as part of a strategy to construct stereochemically rigid peptide helices, in a modular approach to design of protein mimics. The peptides Boc-(Val-Ala-Leu-Aib),-OMe ( I ) and Boc-Val-Ala-Leu-Aib-Val-Ala-Leu-(Val-Ala-Leu-Aib()11z)- OhaMvee been crystallized.Both crystals are stable only in the presence of mother liquor or water. The crystal data are as follows. I: C78H140N16019~2H20,P2,, a = 16.391 (3) A, b = 16.860 (3) A, c = 18.428 (3) A, p = 103.02 (I)O, Z = 2, R = 9.6% for 3445 data with lFol >30(F), resolution 0.93 A. 11: C7,Hl,,N,S018.7.5H,0, C2221, a = 18.348 ( 5 ) A, b = 47.382 (1 1) A, c = 24.157 ( 5 ) A, Z =8, R = l0,6%, for 3147 data with lFol > 3a(F), resolution 1.00 A. The 15-residue peptide (11) is entirely a helical, while the 16-residue peptide ( I ) has a short segment of 310 helix at the N terminus. The packing of the helices in the crystals is rather incfficicnt with no particular attractions between Leu-Leu side chains, or any other pair. Both crystals have fairly large voids, which are filled with water molecules in a disordered fashion. Water molecule sites near the polar head-to-tail regions are well detcrmined, those closer to the hydrophobic side chains less so and a number of possible water sites in the remaining "empty" space are not determined. No interdigitation of Leu side chains is observed in the crystal as is hypothesized in the "leucine zipper" class of DNA binding proteins.