Emergence of half-metallicity in suspended NiO chains: Ab initio electronic structure and quantum transport calculations

Contrary to the antiferromagnetic and insulating character of bulk NiO, one-dimensional chains of this material can become half metallic due to the lower coordination of their atoms. Here we present ab initio electronic structure and quantum transport calculations of ideal infinitely long NiO chains and of more realistic short ones suspended between Ni electrodes. While infinite chains are insulating, short suspended chains are half-metallic minority-spin conductors that displays very large magnetoresistance and a spin-valve behavior controlled by a single atom.

Lifshitz transition and modulation of electronic and transport properties of bilayer graphene by sliding and applied normal compressive strain

Using density functional theory (DFT) we investigate the changes in electronic and transport properties of graphene bilayer caused by sliding one of the layers. Change in stacking pattern breaks the lattice symmetry, which results in Lifshitz transition together with the modulation of the electronic structure. Going from AA to AB stacking by sliding along armchair direction leads to a drastic transition in electronic structure from linear to parabolic dispersion. Our transport calculations show a significant change in the overall transmission value for large sliding distances along zigzag direction. The increase in interlayer coupling with normal compressive strain increases the overlapping of conduction and valence band, which leads to further shift in the Dirac points and an enhancement in the Lifshitz transition. The ability to tune the topology of band structure by sliding and/or applying normal compressive strain will open doors for controlled tuning of many physical phenomenon such as Landau levels and quantum Hall effect in graphene. (C) 2015 Elsevier Ltd. All rights reserved.

Ab Initio Calculations of Structural Evolution and Conductance of Benzene-1,4-dithiol on Gold Leads

By performing at) initio density functional theory (DFT) calculations and electronic transport simulations based on the OFT nonequilibrium Green`s functions method we investigate how the conformational changes of a benzene-1,4-dithiol molecule bonded to gold affect the molecular transport as the electrodes are separated from each other. In particular we consider the full evolution of the stretching process until the Junction breaking point and compare results obtained with a standard semilocal exchange and correlation functional to those computed with a self-interaction corrected method. We conclude that the inclusion of self-interaction corrections is fundamental for describing both the molecule conductance and its stability against conformational fluctuations.

Realistic calculations of carbon-based disordered systems

Carbon nanotubes rank amongst potential candidates for a new family of nanoscopic devices, in particular for sensing applications. At the same time that defects in carbon nanotubes act as binding sites for foreign species, our current level of control over the fabrication process does not allow one to specifically choose where these binding sites will actually be positioned. In this work we present a theoretical framework for accurately calculating the electronic and transport properties of long disordered carbon nanotubes containing a large number of binding sites randomly distributed along a sample. This method combines the accuracy and functionality of ab initio density functional theory to determine the electronic structure with a recursive Green`s functions method. We apply this methodology on the problem of nitrogen-rich carbon nanotubes, first considering different types of defects and then demonstrating how our simulations can help in the field of sensor design by allowing one to compute the transport properties of realistic nanotube devices containing a large number of randomly distributed binding sites.

Study of the Influence of Localized Vibrational Modes in Charge Transport Properties at Nanoscale Systems.

In molecular and atomic devices the interaction between electrons and ionic vibrations has an important role in electronic transport. The electron-phonon coupling can cause the loss of the electron's phase coherence, the opening of new conductance channels and the suppression of purely elastic ones. From the technological viewpoint phonons might restrict the efficiency of electronic devices by energy dissipation, causing heating, power loss and instability. The state of the art in electron transport calculations consists in combining ab initio calculations via Density Functional Theory (DFT) with Non-Equilibrium Green's Function formalism (NEGF). In order to include electron-phonon interactions, one needs in principle to include a self-energy scattering term in the open system Hamiltonian which takes into account the effect of the phonons over the electrons and vice versa. Nevertheless this term could be obtained approximately by perturbative methods. In the First Born Approximation one considers only the first order terms of the electronic Green's function expansion. In the Self-Consistent Born Approximation, the interaction self-energy is calculated with the perturbed electronic Green's function in a self-consistent way. In this work we describe how to incorporate the electron-phonon interaction to the SMEAGOL program (Spin and Molecular Electronics in Atomically Generated Orbital Landscapes), an ab initio code for electronic transport based on the combination of DFT + NEGF. This provides a tool for calculating the transport properties of materials' specific system, particularly in molecular electronics. Preliminary results will be presented, showing the effects produced by considering the electron-phonon interaction in nanoscale devices.

Application of reactive fluid transport modeling to hydrothermal systems

A one-dimensional multi-component reactive fluid transport algorithm, 1DREACT (Steefel, 1993) was used to investigate different fluid-rock interaction systems. A major short coming of mass transport calculations which include mineral reactions is that solid solutions occurring in many minerals are not treated adequately. Since many thermodynamic models of solid solutions are highly non-linear, this can seriously impact on the stability and efficiency of the solution algorithms used. Phase petrology community saw itself faced with a similar predicament 10 years ago. To improve performance and reliability, phase equilibrium calculations have been using pseudo compounds. The same approach is used here in the first, using the complex plagioclase solid solution as an example. Thermodynamic properties of a varying number of intermediate plagioclase phases were calculated using ideal molecular, Al-avoidance, and non-ideal mixing models. These different mixing models can easily be incorporated into the simulations without modification of the transport code. Simulation results show that as few as nine intermediate compositions are sufficient to characterize the diffusional profile between albite and anorthite. Hence this approach is very efficient, and can be used with little effort. A subsequent chapter reports the results of reactive fluid transport modeling designed to constrain the hydrothermal alteration of Paleoproterozoic sediments of the Southern Lake Superior region. Field observations reveal that quartz-pyrophyllite (or kaolinite) bearing assemblages have been transformed into muscovite-pyrophyllite-diaspore bearing assemblages due to action of fluids migrating along permeable flow channels. Fluid-rock interaction modeling with an initial qtz-prl assemblage and a K-rich fluid simulates the formation of observed mineralogical transformation. The bulk composition of the system evolves from an SiO2-rich one to an Al2O3+K2O-rich one. Simulations show that the fluid flow was up-temperature (e.g. recharge) and that fluid was K-rich. Pseudo compound approach to include solid solutions in reactive transport models was tested in modeling hydrothermal alteration of Icelandic basalts. Solid solutions of chlorites, amphiboles and plagioclase were included as the secondary mineral phases. Saline and fresh water compositions of geothermal fluids were used to investigate the effect of salinity on alteration. Fluid-rock interaction simulations produce the observed mineral transformations. They show that roughly the same alteration minerals are formed due to reactions with both types of fluid which is in agreement with the field observations. A final application is directed towards the remediation of nitrate rich groundwaters. Removal of excess nitrate from groundwater by pyrite oxidation was modeled using the reactive fluid transport algorithm. Model results show that, when a pyrite-bearing, permeable zone is placed in the flow path, nitrate concentration in infiltrating water can be significantly lowered, in agreement with proposals from the literature. This is due to nitrogen reduction. Several simulations investigate the efficiency of systems with different mineral reactive surface areas, reactive barrier zone widths, and flow rates to identify the optimum setup.

Optimization of multidimensional cross-section tables for few-group core calculations

Multigroup diffusion codes for three dimensional LWR core analysis use as input data pre-generated homogenized few group cross sections and discontinuity factors for certain combinations of state variables, such as temperatures or densities. The simplest way of compiling those data are tabulated libraries, where a grid covering the domain of state variables is defined and the homogenized cross sections are computed at the grid points. Then, during the core calculation, an interpolation algorithm is used to compute the cross sections from the table values. Since interpolation errors depend on the distance between the grid points, a determined refinement of the mesh is required to reach a target accuracy, which could lead to large data storage volume and a large number of lattice transport calculations. In this paper, a simple and effective procedure to optimize the distribution of grid points for tabulated libraries is presented. Optimality is considered in the sense of building a non-uniform point distribution with the minimum number of grid points for each state variable satisfying a given target accuracy in k-effective. The procedure consists of determining the sensitivity coefficients of k-effective to cross sections using perturbation theory; and estimating the interpolation errors committed with different mesh steps for each state variable. These results allow evaluating the influence of interpolation errors of each cross section on k-effective for any combination of state variables, and estimating the optimal distance between grid points.

Electronic structure and transport properties of atomic NiO spinvalves

Ab initio quantum transport calculations show that short NiO chains suspended in Ni nanocontacts present a very strong spin-polarization of the conductance.The generalized gradient approximation we use here predicts a similar polarization of the conductance as the one previously computed with non-local exchange, confirming the robustness of the result. Their use as nanoscopic spinvalves is proposed.

Single crystalline ultrathin gold nanowires: Promising nanoscale interconnects

Using first principles based density functional calculation we study the mechanical, electronic and transport properties of single crystalline gold nanowires. While nanowires with the diameter less than 2 nm retain hexagonal cross-section, the larger diameter wires show a structural smoothening leading to circular cross-section. These structural changes significantly affect the mechanical properties of the wires, however, strength remains comparable to the bulk. The transport calculations reveal that the conductivity of these wires are in good agreement with experiments. The combination of good mechanical, electronic and transport properties make these wires promising as interconnects for nano devices. Copyright 2013 Author(s). This article is distributed under a Creative Commons Attribution 3.0 Unported License. http://dx.doi.org/10.1063/1.4796188]

Thermoelectric properties of beta-FeSi2

We investigate the thermoelectric properties of beta-FeSi2 using first principles electronic structure and Boltzmann transport calculations. We report a high thermopower for both p- and n-type beta-FeSi2 over a wide range of carrier concentration and in addition find the performance for n-type to be higher than for the p-type. Our results indicate that, depending upon temperature, a doping level of 3 x 10(20) to 2 x 10(21) cm(-3) may optimize the thermoelectric performance. (C) 2013 AIP Publishing LLC.

Semiconductor-like Sensitivity in Metallic Ultrathin Gold Nanowire-Based Sensors

Due to the ease of modification of electronic structure upon analyte adsorption, semiconductors have been the preferred materials as chemical sensors. At reduced dimension, however, the sensitivity of semiconductor-based sensors deteriorates significantly due to passivation, and often by increased band gap caused by quantum confinement. Using first-principles density functional theory combined with Boltzmann transport calculations, we demonstrate semiconductor-like sensitivity toward chemical species in ultrathin gold nanowires (AuNWs). The sensing mechanism is governed by the modification of the electronic structure of the AuNW as well as scattering of the charge carriers by analyte adsorption. Most importantly, the sensitivity exhibits a linear relationship with the electron affinities of the respective analytes. Based on this relationship, we propose an empirical parameter, which can predict an analyte-specific sensitivity of a AuNW, rendering them as effective sensors for a wide range of chemical an alytes.

Strain-induced electronic phase transition and strong enhancement of thermopower of TiS2

Using first-principles density functional theory calculations, we show a semimetal to semiconducting electronic phase transition for bulk TiS2 by applying uniform biaxial tensile strain. This electronic phase transition is triggered by charge transfer from Ti to S, which eventually reduces the overlap between Ti-(d) and S-(p) orbitals. The electronic transport calculations show a large anisotropy in electrical conductivity and thermopower, which is due to the difference in the effective masses along the in-plane and out-of-plane directions. Strain-induced opening of band gap together with changes in dispersion of bands lead to threefold enhancement in thermopower for both p-and n-type TiS2. We further demonstrate that the uniform tensile strain, which enhances the thermoelectric performance, can be achieved by doping TiS2 with larger iso-electronic elements such as Zr or Hf at Ti sites.

High thermopower and ultra low thermal conductivity in Cd-based Zintl phase compounds

By combining first principles density functional theory and electronic as well as lattice Boltzmann transport calculations, we unravel the excellent thermoelectric properties of Zintl phase compounds ACd(2)Sb(2) (where, A = Ca, Ba, Sr). The calculated electronic structures of these compounds show charge carrier pockets and heavy light bands near the band edge, which lead to a large power factor. Furthermore, we report large Gruneisen parameters and low phonon group velocity indicating essential strong anharmonicity in these compounds, which resulted in low lattice thermal conductivity. The combination of low thermal conductivity and the excellent transport properties give a high ZT value of similar to 1.4-1.9 in CaCd2Sb2 and BaCd2Sb2 at moderate p and n-type doping. Our results indicate that well optimized Cd-based Zintl phase compounds have the potential to match the performance of conventional thermoelectric materials.

Métodos espectronodais para cálculos de transporte de partículas neutras com fonte fixa na formulação de ordenadas discretas e multigrupo de energia

Um método espectronodal é desenvolvido para problemas de transporte de partículas neutras de fonte fixa, multigrupo de energia em geometria cartesiana na formulação de ordenadas discretas (SN). Para geometria unidimensional o método espectronodal multigrupo denomina-se método spectral Greens function (SGF) com o esquema de inversão nodal (NBI) que converge solução numérica para problemas SN multigrupo em geometria unidimensional, que são completamente livre de erros de truncamento espacial para ordem L de anisotropia de espalhamento desde que L < N. Para geometria X; Y o método espectronodal multigrupo baseia-se em integrações transversais das equações SN no interior dos nodos de discretização espacial, separadamente nas direções coordenadas x e y. Já que os termos de fuga transversal são aproximados por constantes, o método nodal resultante denomina-se SGF-constant nodal (SGF-CN), que é aplicado a problemas SN multigrupo de fonte fixa em geometria X; Y com espalhamento isotrópico. Resultados numéricos são apresentados para ilustrar a eficiência dos códigos SGF e SGF-CN e a precisão das soluções numéricas convergidas em cálculos de malha grossa.

Método numérico de Matriz Resposta acoplado a um esquema de reconstrução espacial analítica para cálculos unidimensionais de transporte de nêutrons na formulação de ordenadas discretas multigrupo de energia com fonte fixa

Um método de Matriz Resposta (MR) é descrito para gerar soluções numéricas livres de erros de truncamento espacial para problemas multigrupo de transporte de nêutrons com fonte fixa e em geometria unidimensional na formulação de ordenadas discretas (SN). Portanto, o método multigrupo MR com esquema iterativo de inversão nodal parcial (NBI) converge valores numéricos para os fluxos angulares nas fronteiras das regiões que coincidem com os valores da solução analítica das equações multigrupo SN, afora os erros de arredondamento da aritmética finita computacional. É também desenvolvido um esquema numérico de reconstrução espacial, que fornece a saída para os fluxos escalares de nêutrons em cada grupo de energia em um intervalo qualquer do domínio definido pelo usuário, com um passo de avanço também escolhido pelo usuário. Resultados numéricos são apresentados para ilustrar a precisão do presente método em cálculos de malha grossa.