Elastic effects on relaxation volume tensor calculations

Relaxation volume tensors quantify the effect of stress on diffusion of crystal defects. Continuum linear elasticity predicts that calculations of these parameters using periodic boundary conditions do not suffer from systematic deviations due to elastic image effects and should be independent of supercell size or symmetry. In practice, however, calculations of formation volume tensors of the interstitial in Stillinger-Weber silicon demonstrate that changes in bonding at the defect affect the elastic moduli and result in system-size dependent relaxation volumes. These vary with the inverse of the system size. Knowing the rate of convergence permits accurate estimates of these quantities from modestly sized calculations. Furthermore, within the continuum linear elasticity assumptions the average stress can be used to estimate the relaxation volume tensor from constant volume calculations.Comment: 31 pages, 6 figures, submitted to Phys. Rev.

Non-Periodic Finite-Element Formulation of Orbital-Free Density Functional Theory

We propose an approach to perform orbital-free density functional theory calculations in a non-periodic setting using the finite-element method. We consider this a step towards constructing a seamless multi-scale approach for studying defects like vacancies, dislocations and cracks that require quantum mechanical resolution at the core and are sensitive to long range continuum stresses. In this paper, we describe a local real space variational formulation for orbital-free density functional theory, including the electrostatic terms and prove existence results. We prove the convergence of the finite-element approximation including numerical quadratures for our variational formulation. Finally, we demonstrate our method using examples

Two groups of phenylalanine biosynthetic operon leader peptides genes: a high level of apparently incidental frameshifting in decoding Escherichia coli pheL

The bacterial pheL gene encodes the leader peptide for the phenylalanine biosynthetic operon. Translation of pheL mRNA controls transcription attenuation and, consequently, expression of the downstream pheA gene. Fifty-three unique pheL genes have been identified in sequenced genomes of the gamma subdivision. There are two groups of pheL genes, both of which are short and contain a run(s) of phenylalanine codons at an internal position. One group is somewhat diverse and features different termination and 5′-flanking codons. The other group, mostly restricted to Enterobacteria and including Escherichia coli pheL, has a conserved nucleotide sequence that ends with UUC_CCC_UGA. When these three codons in E. coli pheL mRNA are in the ribosomal E-, P- and A-sites, there is an unusually high level, 15%, of +1 ribosomal frameshifting due to features of the nascent peptide sequence that include the penultimate phenylalanine. This level increases to 60% with a natural, heterologous, nascent peptide stimulator. Nevertheless, studies with different tRNAPro mutants in Salmonella enterica suggest that frameshifting at the end of pheL does not influence expression of the downstream pheA. This finding of incidental, rather than utilized, frameshifting is cautionary for other studies of programmed frameshifting

Roadmap on Electronic Structure Codes in the Exascale Era

Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing

Roadmap on Electronic Structure Codes in the Exascale Era

Electronic structure calculations have been instrumental in providing many important insights into a range of physical and chemical properties of various molecular and solid-state systems. Their importance to various fields, including materials science, chemical sciences, computational chemistry and device physics, is underscored by the large fraction of available public supercomputing resources devoted to these calculations. As we enter the exascale era, exciting new opportunities to increase simulation numbers, sizes, and accuracies present themselves. In order to realize these promises, the community of electronic structure software developers will however first have to tackle a number of challenges pertaining to the efficient use of new architectures that will rely heavily on massive parallelism and hardware accelerators. This roadmap provides a broad overview of the state-of-the-art in electronic structure calculations and of the various new directions being pursued by the community. It covers 14 electronic structure codes, presenting their current status, their development priorities over the next five years, and their plans towards tackling the challenges and leveraging the opportunities presented by the advent of exascale computing

4D Imaging and Diffraction Dynamics of Single-Particle Phase Transition in Heterogeneous Ensembles

In this Letter, we introduce conical-scanning dark-field imaging in four-dimensional (4D) ultrafast electron microscopy to visualize single-particle dynamics of a polycrystalline ensemble undergoing phase transitions. Specifically, the ultrafast metal–insulator phase transition of vanadium dioxide is induced using laser excitation and followed by taking electron-pulsed, time-resolved images and diffraction patterns. The single-particle selectivity is achieved by identifying the origin of all constituent Bragg spots on Debye–Scherrer rings from the ensemble. Orientation mapping and dynamic scattering simulation of the electron diffraction patterns in the monoclinic and tetragonal phase during the transition confirm the observed behavior of Bragg spots change with time. We found that the threshold temperature for phase recovery increases with increasing particle sizes and we quantified the observation through a theoretical model developed for single-particle phase transitions. The reported methodology of conical scanning, orientation mapping in 4D imaging promises to be powerful for heterogeneous ensemble, as it enables imaging and diffraction at a given time with a full archive of structural information for each particle, for example, size, morphology, and orientation while minimizing radiation damage to the specimen

Classification of Plant Associated Bacteria Using RIF, a Computationally Derived DNA Marker

A DNA marker that distinguishes plant associated bacteria at the species level and below was derived by comparing six sequenced genomes of Xanthomonas, a genus that contains many important phytopathogens. This DNA marker comprises a portion of the dnaA replication initiation factor (RIF). Unlike the rRNA genes, dnaA is a single copy gene in the vast majority of sequenced bacterial genomes, and amplification of RIF requires genus-specific primers. In silico analysis revealed that RIF has equal or greater ability to differentiate closely related species of Xanthomonas than the widely used ribosomal intergenic spacer region (ITS). Furthermore, in a set of 263 Xanthomonas, Ralstonia and Clavibacter strains, the RIF marker was directly sequenced in both directions with a success rate approximately 16% higher than that for ITS. RIF frameworks for Xanthomonas, Ralstonia and Clavibacter were constructed using 682 reference strains representing different species, subspecies, pathovars, races, hosts and geographic regions, and contain a total of 109 different RIF sequences. RIF sequences showed subspecific groupings but did not place strains of X. campestris or X. axonopodis into currently named pathovars nor R. solanacearum strains into their respective races, confirming previous conclusions that pathovar and race designations do not necessarily reflect genetic relationships. The RIF marker also was sequenced for 24 reference strains from three genera in the Enterobacteriaceae: Pectobacterium, Pantoea and Dickeya. RIF sequences of 70 previously uncharacterized strains of Ralstonia, Clavibacter, Pectobacterium and Dickeya matched, or were similar to, those of known reference strains, illustrating the utility of the frameworks to classify bacteria below the species level and rapidly match unknown isolates to reference strains. The RIF sequence frameworks are available at the online RIF database, RIFdb, and can be queried for diagnostic purposes with RIF sequences obtained from unknown strains in both chromatogram and FASTA format