Linear-scaling density-functional simulations of charged point defects in Al2O3 using hierarchical sparse matrix algebra

We present calculations of formation energies of defects in an ionic solid (Al(2)O(3)) extrapolated to the dilute limit, corresponding to a simulation cell of infinite size. The large-scale calculations required for this extrapolation are enabled by developments in the approach to parallel sparse matrix algebra operations, which are central to linear-scaling density-functional theory calculations. The computational cost of manipulating sparse matrices, whose sizes are determined by the large number of basis functions present, is greatly improved with this new approach. We present details of the sparse algebra scheme implemented in the ONETEP code using hierarchical sparsity patterns, and demonstrate its use in calculations on a wide range of systems, involving thousands of atoms on hundreds to thousands of parallel processes

Electrically induced Dirac fermions in graphene nanoribbons

Graphene nanoribbons are widely regarded as promising building blocks for next-generation carbon-based devices. A critical issue to their prospective applications is whether their electronic structure can be externally controlled. Here, we combine simple model Hamiltonians with extensive first-principles calculations to investigate the response of armchair graphene nanoribbons to transverse electric fields. Such fields can be achieved either upon laterally gating the nanoribbon or incorporating ambipolar chemical codopants along the edges. We reveal that the field induces a semiconductor-to-semimetal transition with the semimetallic phase featuring zero-energy Dirac fermions that propagate along the armchair edges. The transition occurs at critical fields that scale inversely with the width of the nanoribbons. These findings are universal to group-IV honeycomb lattices, including silicene and germanene nanoribbons, irrespective of the type of edge termination. Overall, our results create new opportunities to electrically engineer Dirac semimetallic phases in otherwise semiconducting graphene-like nanoribbons

Prediction of high zT in thermoelectric silicon nanowires with axial germanium heterostructures

We calculate the thermoelectric figure of merit, zT=S^2GT/(\kappa_l+\kappa_e), for p-type Si nanowires with axial Ge heterostructures using a combination of first-principles density-functional theory, interatomic potentials, and Landauer-Buttiker transport theory. We consider nanowires with up to 8400 atoms and twelve Ge axial heterostructures along their length. We find that introducing heterostructures always reduces S^2G, and that our calculated increases in zT are predominantly driven by associated decreases in \kappa_l. Of the systems considered, nanowires with a regular distribution of Ge heterostructures have the highest figure-of-merit: zT=3, an order of magnitude larger than the equivalent pristine nanowire. Even in the presence of realistic structural disorder, in the form of small variations in length of the heterostructures, zT remains several times larger than that of the pristine case, suggesting that axial heterostructuring is a promising route to high-zT thermoelectric nanowires.Comment: In print with Europhysics Letters (2011

Maximally localized Wannier functions in LaMnO3 within PBE+U, hybrid functionals, and partially self-consistent GW: an efficient route to construct ab-initio tight-binding parameters for e_g perovskites

Using the newly developed VASP2WANNIER90 interface we have constructed maximally localized Wannier functions (MLWFs) for the e_g states of the prototypical Jahn-Teller magnetic perovskite LaMnO3 at different levels of approximation for the exchange-correlation kernel. These include conventional density functional theory (DFT) with and without additional on-site Hubbard U term, hybrid-DFT, and partially self-consistent GW. By suitably mapping the MLWFs onto an effective e_g tight-binding (TB) Hamiltonian we have computed a complete set of TB parameters which should serve as guidance for more elaborate treatments of correlation effects in effective Hamiltonian-based approaches. The method-dependent changes of the calculated TB parameters and their interplay with the electron-electron (el-el) interaction term are discussed and interpreted. We discuss two alternative model parameterizations: one in which the effects of the el-el interaction are implicitly incorporated in the otherwise "noninteracting" TB parameters, and a second where we include an explicit mean-field el-el interaction term in the TB Hamiltonian. Both models yield a set of tabulated TB parameters which provide the band dispersion in excellent agreement with the underlying ab initio and MLWF bands.Comment: 30 pages, 7 figure

Robotic Wireless Sensor Networks

In this chapter, we present a literature survey of an emerging, cutting-edge, and multi-disciplinary field of research at the intersection of Robotics and Wireless Sensor Networks (WSN) which we refer to as Robotic Wireless Sensor Networks (RWSN). We define a RWSN as an autonomous networked multi-robot system that aims to achieve certain sensing goals while meeting and maintaining certain communication performance requirements, through cooperative control, learning and adaptation. While both of the component areas, i.e., Robotics and WSN, are very well-known and well-explored, there exist a whole set of new opportunities and research directions at the intersection of these two fields which are relatively or even completely unexplored. One such example would be the use of a set of robotic routers to set up a temporary communication path between a sender and a receiver that uses the controlled mobility to the advantage of packet routing. We find that there exist only a limited number of articles to be directly categorized as RWSN related works whereas there exist a range of articles in the robotics and the WSN literature that are also relevant to this new field of research. To connect the dots, we first identify the core problems and research trends related to RWSN such as connectivity, localization, routing, and robust flow of information. Next, we classify the existing research on RWSN as well as the relevant state-of-the-arts from robotics and WSN community according to the problems and trends identified in the first step. Lastly, we analyze what is missing in the existing literature, and identify topics that require more research attention in the future

A converse approach to the calculation of NMR shielding tensors

We introduce an alternative approach to the first-principles calculation of NMR shielding tensors. These are obtained from the derivative of the orbital magnetization with respect to the application of a microscopic, localized magnetic dipole. The approach is simple, general, and can be applied to either isolated or periodic systems. Calculated results for simple hydrocarbons, crystalline diamond, and liquid water show very good agreement with established methods and experimental results.Comment: submitted to Phys. Rev. Lett; second half completely rewritten to address referee's comment

wannier90: A tool for obtaining maximally-localised Wannier functions

We present Wannier90, a program for calculating maximally-localised Wannier functions (MLWF) from a set of Bloch energy bands that may or may not be attached to or mixed with other bands. The formalism works by minimising the total spread of the MLWF in real space. This done in the space of unitary matrices that describe rotations of the Bloch bands at each k-point. As a result, Wannier90 is independent of the basis set used in the underlying calculation to obtain the Bloch states. Therefore, it may be interfaced straightforwardly to any electronic structure code. The locality of MLWF can be exploited to compute band-structure, density of states and Fermi surfaces at modest computational cost. Furthermore, Wannier90 is able to output MLWF for visualisation and other post-processing purposes. Wannier functions are already used in a wide variety of applications. These include analysis of chemical bonding in real space; calculation of dielectric properties via the modern theory of polarisation; and as an accurate and minimal basis set in the construction of model Hamiltonians for large-scale systems, in linear-scaling quantum Monte Carlo calculations, and for efficient computation of material properties, such as the anomalous Hall coefficient. Wannier90 is freely available under the GNU General Public License from http://www.wannier.org/Comment: Submitted to Computer Physics Communication

Effect of teeth micro-geometrical form modification on contact kinematics and efficiency of high performance transmissions

Light weight, compactness and efficiency are key objectives in high performance vehicular transmission systems, which are subject to large variations in torque and power. Pitch line velocities of up to 52 m/s and teeth pair contact pressures of up to 3 GPa are routinely encountered under race conditions. Contact patch asymmetry due to angular misalignments between input and output shafts leads to the generation of high edge stress discontinuities on gear flanks, inducing fatigue spalling which affects system durability. Crowning is widely used as a palliative measure to mitigate these undesired effects. These problems can be further exacerbated by contact footprint truncation. The paper presents a new approach to modelling the kinematics and contact micro-geometry of meshing conjunctions of involute spur gears with profile and lead modifications. A time-efficient analytical method is presented to accurately determine the contact footprint and kinematics, leading to the solution of highly loaded non-Newtonian mixed thermo-elastohydrodynamic contact under the extreme prevalent conditions of high performance vehicular transmissions. The effect of tooth form modification on contact footprint truncation, contact kinematics and generated frictional power loss is investigated. This approach has not hitherto been reported in literature

Metastatic primary adenocarcinoma of the bladder in a twenty-five years old woman

Primary adenocarcinoma of the bladder is a rare tumor. The classification between primary vesical and urachal is debated. We present the case of a young female who presented clinicopathological features of a metastatic urachal adenocarcinoma, but the histological result revealed primary adenocarcinoma of the bladder contrary to expectancy. To the best of our knowledge this is the first reported case of a metastatic adenocarcinoma of the bladder in a 25 years old female. This case emphasizes the challenge for urologists to recognize and manage this aggressive tumor in the setting described

Low-Energy Effective Hamiltonian and the Surface States of Ca_3PbO

The band structure of Ca_3PbO, which possesses a three-dimensional massive Dirac electron at the Fermi energy, is investigated in detail. Analysis of the orbital weight distributions on the bands obtained in the first-principles calculation reveals that the bands crossing the Fermi energy originate from the three Pb-p orbitals and three Ca-dx2y2 orbitals. Taking these Pb-p and Ca-dx2y2 orbitals as basis wave functions, a tight-binding model is constructed. With the appropriate choice of the hopping integrals and the strength of the spin-orbit coupling, the constructed model sucessfully captures important features of the band structure around the Fermi energy obtained in the first-principles calculation. By applying the suitable basis transformation and expanding the matrix elements in the series of the momentum measured from a Dirac point, the low-energy effective Hamiltonian of this model is explicitely derived and proved to be a Dirac Hamiltonain. The origin of the mass term is also discussed. It is shown that the spin-orbit coupling and the orbitals other than Pb-p and Ca-dx2y2 orbitals play important roles in making the mass term finite. Finally, the surface band structures of Ca_3PbO for several types of surfaces are investigated using the constructed tight-binding model. We find that there appear nontrivial surface states that cannot be explained as the bulk bands projected on the surface Brillouin zone. The relation to the topological insulator is also discussed.Comment: 11 page