Efficiency of ab-initio total energy calculations for metals and semiconductors using a plane-wave basis set

Three parallel algorithms for classical molecular dynamics are presented. The first assigns each processor a fixed subset of atoms; the second assigns each a fixed subset of inter-atomic forces to compute; the third assigns each a fixed spatial region. The algorithms are suitable for molecular dynamics models which can be difficult to parallelize efficiently—those with short-range forces where the neighbors of each atom change rapidly. They can be implemented on any distributed-memory parallel machine which allows for message-passing of data between independently executing processors. The algorithms are tested on a standard Lennard-Jones benchmark problem for system sizes ranging from 500 to 100,000,000 atoms on several parallel supercomputers--the nCUBE 2, Intel iPSC/860 and Paragon, and Cray T3D. Comparing the results to the fastest reported vectorized Cray Y-MP and C90 algorithm shows that the current generation of parallel machines is competitive with conventional vector supercomputers even for small problems. For large problems, the spatial algorithm achieves parallel efficiencies of 90% and a 1840-node Intel Paragon performs up to 165 faster than a single Cray C9O processor. Trade-offs between the three algorithms and guidelines for adapting them to more complex molecular dynamics simulations are also discussed.

Fast parallel algorithms for short-range molecular dynamics

QUANTUM ESPRESSO is an integrated suite of computer codes for electronic-structure calculations and materials modeling, based on density-functional theory, plane waves, and pseudopotentials (norm-conserving, ultrasoft, and projector-augmented wave). The acronym ESPRESSO stands for opEn Source Package for Research in Electronic Structure, Simulation, and Optimization. It is freely available to researchers around the world under the terms of the GNU General Public License. QUANTUM ESPRESSO builds upon newly-restructured electronic-structure codes that have been developed and tested by some of the original authors of novel electronic-structure algorithms and applied in the last twenty years by some of the leading materials modeling groups worldwide. Innovation and efficiency are still its main focus, with special attention paid to massively parallel architectures, and a great effort being devoted to user friendliness. QUANTUM ESPRESSO is evolving towards a distribution of independent and interoperable codes in the spirit of an open-source project, where researchers active in the field of electronic-structure calculations are encouraged to participate in the project by contributing their own codes or by implementing their own ideas into existing codes.

/pdf/quantum-espresso-a-modular-and-open-source-software-project-4t6vlupbeh.pdf

QUANTUM ESPRESSO: a modular and open-source software project for quantum simulations of materials

Deposits of clastic carbonate-dominated (calciclastic) sedimentary slope systems in the rock record have been identified mostly as linearly-consistent carbonate apron deposits, even though most ancient clastic carbonate slope deposits fit the submarine fan systems better. Calciclastic submarine fans are consequently rarely described and are poorly understood. Subsequently, very little is known especially in mud-dominated calciclastic submarine fan systems. Presented in this study are a sedimentological core and petrographic characterisation of samples from eleven boreholes from the Lower Carboniferous of Bowland Basin (Northwest England) that reveals a >250 m thick calciturbidite complex deposited in a calciclastic submarine fan setting. Seven facies are recognised from core and thin section characterisation and are grouped into three carbonate turbidite sequences. They include: 1) Calciturbidites, comprising mostly of highto low-density, wavy-laminated bioclast-rich facies; 2) low-density densite mudstones which are characterised by planar laminated and unlaminated muddominated facies; and 3) Calcidebrites which are muddy or hyper-concentrated debrisflow deposits occurring as poorly-sorted, chaotic, mud-supported floatstones. These

Phd by thesis

The surface science of titanium dioxide

We present a sampling method for Brillouin-zone integration in metals which converges exponentially with the number of sampling points, without the loss of precision of normal broadening techniques. The scheme is based on smooth approximants to the \ensuremath{\delta} and step functions which are constructed to give the exact result when integrating polynomials of a prescribed degree. In applications to the simple-cubic tight-binding band as well as to band structures of simple and transition metals, we demonstrate significant improvement over existing methods. The method promises general applicability in the fields of total-energy calculations and many-body physics.

High-precision sampling for Brillouin-zone integration in metals

A stabilization mechanism of zirconia based on oxygen vacancies only

Embrittlement by the segregation of impurity elements to grain boundaries is one of a small number of phenomena that can lead to metallurgical failure by fast fracture. Here we settle a question that has been debated for over a hundred years: how can minute traces of bismuth in copper cause this ductile metal to fail in a brittle manner? Three hypotheses for Bi embrittlement of Cu exist: two assign an electronic effect to either a strengthening or weakening of bonds, the third postulates a simple atomic size effect. Here we report first principles quantum mechanical calculations that allow us to reject the electronic hypotheses, while supporting a size effect. We show that upon segregation to the grain boundary, the large Bi atoms weaken the interatomic bonding by pushing apart the Cu atoms at the interface. The resolution of the mechanism underlying grain boundary weakening should be relevant for all cases of embrittlement by oversize impurities.

Bismuth embrittlement of copper is an atomic size effect

s131d and s231d reconstructions of the (001) SrTiO3 surface were studied using the first-principles fullpotential linear muffin-tin orbital method. Surface energies were calculated as a function of TiO2 chemical potential, oxygen partial pressure and temperature. The s131d unreconstructed surfaces were found to be energetically stable for many of the conditions considered. Under conditions of very low oxygen partial pressure the s231d Ti2O3 reconstruction [Martin R. Castell, Surf. Sci. 505, 1 (2002)] is stable. The question as to why STM images of the s131d surfaces have not been obtained was addressed by calculating charge densities for each surface. These suggest that the s231d reconstructions would be easier to image than the s131d surfaces. The possibility that the presence of oxygen vacancies would destabilise the s131d surfaces was also investigated. If the s131d surfaces are unstable then there exists the further possibility that the s2 31d DL-TiO2 reconstruction [Natasha Erdman et al. Nature (London) 419, 55 (2002)] is stable in a TiO2-rich environment and for pO2 .10−18 atm.

/pdf/srtio3-001-2x1-reconstructions-first-principles-calculations-4l32ek2ec3.pdf

SrTiO3 (001) (2x1) reconstructions: first-principles calculations of surface energy and atomic structure compared with scanning tunnelling microscopy images

We use the local spin-density approximation in density-functional theory to study the electronic structure of ${\mathrm{TiO}}_{2}$ and its (110) stoichiometric surface, as well as ${\mathrm{Ti}}_{2}{\mathrm{O}}_{3}$ and the reduced (110) surface, which are examples of reduced ${\mathrm{TiO}}_{2},$ in order to address the question of how donated electrons are distributed among the previously unoccupied Ti $3d$ states. We conclude that, whereas a split-off band in ${\mathrm{Ti}}_{2}{\mathrm{O}}_{3}$ probably leads to Mott-Hubbard insulating behavior, a picture of itinerant electrons in approximately one-electron states with Stoner exchange splitting is an appropriate description in the case of the reduced surface. Our calculations will be useful in interpreting photoemission experiments as well as measurements of core-level shifts at surfaces.

Anthony Paxton

Papers

High-precision sampling for Brillouin-zone integration in metals

A stabilization mechanism of zirconia based on oxygen vacancies only

Bismuth embrittlement of copper is an atomic size effect

SrTiO3 (001) (2x1) reconstructions: first-principles calculations of surface energy and atomic structure compared with scanning tunnelling microscopy images

Electronic structure of reduced titanium dioxide