Mechanisms of molecular doping of graphene: A first-principles study

Doping graphene with electron donating or accepting molecules is an interesting approach to introduce carriers into it, analogous to electrochemical doping accomplished in graphene when used in a field-effect transistor. Here, we use first-principles density-functional theory to determine changes in the electronic-structure and vibrational properties of graphene that arise from the adsorption of aromatic molecules such as aniline and nitrobenzene. Identifying the roles of various mechanisms of chemical interaction between graphene and a molecule, we bring out the contrast between electrochemical and molecular doping of graphene. Our estimates of various contributions to shifts in the Raman-active modes of graphene with molecular doping are fundamental to the possible use of Raman spectroscopy in (a) characterization of the nature and concentration of carriers in graphene with molecular doping, and (b) graphene-based chemical sensors.

First-principles electronic structure of the delafossites ABO(2) (A=Cu, Ag, Au; B=Al, Ga, Sc, In, Y): evolution of d(10)-d(10) interactions

While bonding between d(10) atoms and ions in molecular systems has been well studied, less attention has been paid to interactions between such seemingly closed shell species in extended inorganic solids. In this contribution, we present visualizations of the electronic structures of the delafossites ABO(2) (A = Cu, Ag, Au) with particular emphasis on the nature of d(10)-d(10) interactions in the close packed plane of the coinage metal ion. We find that on going from Cu to Ag to Au, the extent of bonding between A and A increases. However, the structures (in terms of distances) of these compounds are largely determined by the strongly ionic 13,11 0 interaction and for the larger B ions Sc, In and Y, the A atoms are sufficiently well-separated that A-A bonding is almost negligible. We also analyze some interesting differences between Ag and Au, including the larger A-O covalency of the Au. The trends in electronic structure suggest that the Ag and Au compounds are not good candidate transparent conducting oxides. (C) 2002 Editions scientifiques et medicales Elsevier SAS. All rights reserved.

Elastic Properties Of Sulphur And Selenium Doped Ternary PbTe Alloys By First Principles

Lead telluride (PbTe) is an established thermoelectric material which can be alloyed with sulphur and selenium to further enhance the thermoelectric properties. Here, a first principles study of ternary alloys PbSxTe(1-x) and PbSexTe(1-x) (0 <= x <= 1) based on the Virtual Crystal Approximation (VCA) is presented for different ratios of the isoelectronic atoms in each series. Equilibrium lattice parameters and elastic constants have been calculated and compared with the reported data. Anisotropy parameter calculated from the stiffness constants showed a slight improvement in anisotropy of elastic properties of the alloys over undoped PbTe. Furthermore, the alloys satisfied the predicted stability criteria from the elastic constants, showing stable structures, which agreed with the previously reported experimental results.

First principles study of metal contacts to monolayer black phosphorous

Atomically thin layered black phosphorous (BP) has recently appeared as an alternative to the transitional metal dichalcogenides for future channel material in a metal-oxide-semiconductor transistor due to its lower carrier effective mass. Investigation of the electronic property of source/drain contact involving metal and two-dimensional material is essential as it impacts the transistor performance. In this paper, we perform a systematic and rigorous study to evaluate the Ohmic nature of the side-contact formed by the monolayer BP (mBP) and metals (gold, titanium, and palladium), which are commonly used in experiments. Employing the Density Functional Theory, we analyse the potential barrier, charge transfer and atomic orbital overlap at the metal-mBP interface in an optimized structure to understand how efficiently carriers could be injected from metal contact to the mBP channel. Our analysis shows that gold forms a Schottky contact with a higher tunnel barrier at the interface in comparison to the titanium and palladium. mBP contact with palladium is found to be purely Ohmic, where as titanium contact demonstrates an intermediate behaviour. (C) 2014 AIP Publishing LLC.

The effect of gamma-gamma ` interface on the tensile and shear strengths of nickel-based superalloys: A first-principles study

First-principles density functional theory has been used to evaluate the shear and cleavage strength in terms of Griffith work and generalized stacking fault energy (GSF) of (001) plane for gamma, gamma' and gamma-gamma' system as a function of distance from the gamma/gamma' interface. Calculation of Griffith work suggests higher cleavage energy for bulk gamma as compared to gamma' while the GSF calculation suggests higher shear strength for bulk gamma' as compared to gamma. It has been found that the shear strength of the cubic plane of the gamma/gamma' interface is marginally lower than those of bulk gamma and gamma' phases. (C) 2014 Elsevier B.V. All rights reserved.

First-Principles Study on Structural Stability of Belite

A first-principles study was carried out to investigate the stability of the crystal structure of beta-form belite (beta-C2S) substituted by Sr atoms as trace impurities for Ca atoms in CaOx polyhedra. The effect of the connection types of CaOx polyhedral, in the form of common-edge bond and common-face bond, upon the crystal stability is described. The Ca-Ca interatomic distance closely relates to the hydraulic activity of beta-C2S. The beta-C2S substituted by an Sr atom for Ca(1) atoms having seven Ca-O bonds is energetically more stable than that substituted by an Sr atom for Ca(2) atoms having eight Ca-O bonds. The Sr-doped beta-C2S having a common face bond with SrOx polyhedra is energetically more favorable and results in structural stability compared with that having a common edge bond with SrOx polyhedra.

Design of symmetric TIM barrel proteins from first principles

Background: Computational protein design is a rapidly maturing field within structural biology, with the goal of designing proteins with custom structures and functions. Such proteins could find widespread medical and industrial applications. Here, we have adapted algorithms from the Rosetta software suite to design much larger proteins, based on ideal geometric and topological criteria. Furthermore, we have developed techniques to incorporate symmetry into designed structures. For our first design attempt, we targeted the (alpha/beta)(8) TIM barrel scaffold. We gained novel insights into TIM barrel folding mechanisms from studying natural TIM barrel structures, and from analyzing previous TIM barrel design attempts. Methods: Computational protein design and analysis was performed using the Rosetta software suite and custom scripts. Genes encoding all designed proteins were synthesized and cloned on the pET20-b vector. Standard circular dichroism and gel chromatographic experiments were performed to determine protein biophysical characteristics. 1D NMR and 2D HSQC experiments were performed to determine protein structural characteristics. Results: Extensive protein design simulations coupled with ab initio modeling yielded several all-atom models of ideal, 4-fold symmetric TIM barrels. Four such models were experimentally characterized. The best designed structure (Symmetrin-1) contained a polar, histidine-rich pore, forming an extensive hydrogen bonding network. Symmetrin-1 was easily expressed and readily soluble. It showed circular dichroism spectra characteristic of well-folded alpha/beta proteins. Temperature melting experiments revealed cooperative and reversible unfolding, with a T-m of 44 degrees C and a Gibbs free energy of unfolding (Delta G degrees) of 8.0 kJ/mol. Urea denaturing experiments confirmed these observations, revealing a C-m of 1.6 M and a Delta G degrees of 8.3 kJ/mol. Symmetrin-1 adopted a monomeric conformation, with an apparent molecular weight of 32.12 kDa, and displayed well resolved 1D-NMR spectra. However, the HSQC spectrum revealed somewhat molten characteristics. Conclusions: Despite the detection of molten characteristics, the creation of a soluble, cooperatively folding protein represents an advancement over previous attempts at TIM barrel design. Strategies to further improve Symmetrin-1 are elaborated. Our techniques may be used to create other large, internally symmetric proteins.

First principles DFT investigation Of Yttrium-decorated Boron-Nitride Nanotube: Electronic Structure and Hydrogen Storage

The electronic structure and hydrogen storage capability of Yttrium-doped BNNTs has been theoretically investigated using first principles density functional theory (DFT). Yttrium atom prefers the hollow site in the center of the hexagonal ring with a binding energy of 0.8048eV. Decorating by Y makes the system half-metallic and magnetic with a magnetic moment of 1.0 mu(B). Y decorated Boron-Nitride (8,0) nanotube can adsorb up to five hydrogen molecules whose average binding energy is computed as 0.5044eV. All the hydrogen molecules are adsorbed with an average desorption temperature of 644.708 K. Taking that the Y atoms can be placed only in alternate hexagons, the implied wt% comes out to be 5.31%, a relatively acceptable value for hydrogen storage materials. Thus, this system can serve as potential hydrogen storage medium.

First-principles study of the electronic structure and dynamical electronic excitations of strong spin-orbit metal Pb: bulk, surface and nanosized films

En la presente tesis se ha realizado el estudio de primeros principios (esto es, sinhacer uso de parámetros ajustables) de la estructura electrónica y la dinámica deexcitaciones electrónicas en plomo, tanto en volumen como en superficie y en formade películas de espesor nanométrico. Al presentar el plomo un número atómico alto(82), deben tenerse en cuenta los efectos relativistas. Con este fin, el doctorando haimplementado el acoplo espín-órbita en los códigos computacionales que hanrepresentado la principal herramienta de trabajo.En volumen, se han encontrado fuertes efectos relativistas asi como de lalocalización de los electrones, tanto en la respuesta dieléctrica (excitacioneselectrónicas colectivas) como en el tiempo de vida de electrones excitados. Lacomparación de nuestros resultados con medidas experimentales ha ayudado aprofundizar en dichos efectos.En el estudio de las películas a escala nanométrica se han hallado fuertes efectoscuánticos debido al confinamiento de los estados electrónicos. Dichos efectos semanifiestan tanto en el estado fundamental (en acuerdo con estudiosexperimentales), como en la respuesta dieléctrica a través de la aparición y dinámicade plasmones de diversas características. Los efectos relativistas, a pesar de no serimportantes en la estructura electrónica de las películas, son los responsables de ladesaparación del plasmón de baja energía en nuestros resultados.

First principles based multiparadigm modeling of electronic structures and dynamics

Electronic structures and dynamics are the key to linking the material composition and structure to functionality and performance.

An essential issue in developing semiconductor devices for photovoltaics is to design materials with optimal band gaps and relative positioning of band levels. Approximate DFT methods have been justified to predict band gaps from KS/GKS eigenvalues, but the accuracy is decisively dependent on the choice of XC functionals. We show here for CuInSe₂ and CuGaSe₂, the parent compounds of the promising CIGS solar cells, conventional LDA and GGA obtain gaps of 0.0-0.01 and 0.02-0.24 eV (versus experimental values of 1.04 and 1.67 eV), while the historically first global hybrid functional, B3PW91, is surprisingly the best, with band gaps of 1.07 and 1.58 eV. Furthermore, we show that for 27 related binary and ternary semiconductors, B3PW91 predicts gaps with a MAD of only 0.09 eV, which is substantially better than all modern hybrid functionals, including B3LYP (MAD of 0.19 eV) and screened hybrid functional HSE06 (MAD of 0.18 eV).

The laboratory performance of CIGS solar cells (> 20% efficiency) makes them promising candidate photovoltaic devices. However, there remains little understanding of how defects at the CIGS/CdS interface affect the band offsets and interfacial energies, and hence the performance of manufactured devices. To determine these relationships, we use the B3PW91 hybrid functional of DFT with the AEP method that we validate to provide very accurate descriptions of both band gaps and band offsets. This confirms the weak dependence of band offsets on surface orientation observed experimentally. We predict that the CBO of perfect CuInSe₂/CdS interface is large, 0.79 eV, which would dramatically degrade performance. Moreover we show that band gap widening induced by Ga adjusts only the VBO, and we find that Cd impurities do not significantly affect the CBO. Thus we show that Cu vacancies at the interface play the key role in enabling the tunability of CBO. We predict that Na further improves the CBO through electrostatically elevating the valence levels to decrease the CBO, explaining the observed essential role of Na for high performance. Moreover we find that K leads to a dramatic decrease in the CBO to 0.05 eV, much better than Na. We suggest that the efficiency of CIGS devices might be improved substantially by tuning the ratio of Na to K, with the improved phase stability of Na balancing phase instability from K. All these defects reduce interfacial stability slightly, but not significantly.

A number of exotic structures have been formed through high pressure chemistry, but applications have been hindered by difficulties in recovering the high pressure phase to ambient conditions (i.e., one atmosphere and room temperature). Here we use dispersion-corrected DFT (PBE-ulg flavor) to predict that above 60 GPa the most stable form of N₂O (the laughing gas in its molecular form) is a 1D polymer with an all-nitrogen backbone analogous to cis-polyacetylene in which alternate N are bonded (ionic covalent) to O. The analogous trans-polymer is only 0.03-0.10 eV/molecular unit less stable. Upon relaxation to ambient conditions both polymers relax below 14 GPa to the same stable non-planar trans-polymer, accompanied by possible electronic structure transitions. The predicted phonon spectrum and dissociation kinetics validate the stability of this trans-poly-NNO at ambient conditions, which has potential applications as a new type of conducting polymer with all-nitrogen chains and as a high-energy oxidizer for rocket propulsion. This work illustrates in silico materials discovery particularly in the realm of extreme conditions.

Modeling non-adiabatic electron dynamics has been a long-standing challenge for computational chemistry and materials science, and the eFF method presents a cost-efficient alternative. However, due to the deficiency of FSG representation, eFF is limited to low-Z elements with electrons of predominant s-character. To overcome this, we introduce a formal set of ECP extensions that enable accurate description of p-block elements. The extensions consist of a model representing the core electrons with the nucleus as a single pseudo particle represented by FSG, interacting with valence electrons through ECPs. We demonstrate and validate the ECP extensions for complex bonding structures, geometries, and energetics of systems with p-block character (C, O, Al, Si) and apply them to study materials under extreme mechanical loading conditions.

Despite its success, the eFF framework has some limitations, originated from both the design of Pauli potentials and the FSG representation. To overcome these, we develop a new framework of two-level hierarchy that is a more rigorous and accurate successor to the eFF method. The fundamental level, GHA-QM, is based on a new set of Pauli potentials that renders exact QM level of accuracy for any FSG represented electron systems. To achieve this, we start with using exactly derived energy expressions for the same spin electron pair, and fitting a simple functional form, inspired by DFT, against open singlet electron pair curves (H₂ systems). Symmetric and asymmetric scaling factors are then introduced at this level to recover the QM total energies of multiple electron pair systems from the sum of local interactions. To complement the imperfect FSG representation, the AMPERE extension is implemented, and aims at embedding the interactions associated with both the cusp condition and explicit nodal structures. The whole GHA-QM+AMPERE framework is tested on H element, and the preliminary results are promising.

First-principles study of the electronic structures and absorption spectra for the PbMoO4 crystal with lead vacancy

The electronic structures and absorption spectra for the perfect PbMoO4 crystal and the crystal containing lead vacancy V-Pb(2-) with lattice structure optimized are calculated using density functional theory code CASTEP. The calculated absorption spectra of the PbMoO4 crystal containing V-Pb(2-) exhibit three absorption bands peaking at 2.0 eV (620 nm), 3.0 eV (413 run) and 3.3 eV (375 nm), which are in good agreement with experimental values. The theory predicts that the 390 nm, 430 nm and 580 run absorption bands are related to the existence of V-Pb(2-) in the PbMoO4 crystal.

Gaussian Approximation Potential: an interatomic potential derived from first principles Quantum Mechanics

Simulation of materials at the atomistic level is an important tool in studying microscopic structure and processes. The atomic interactions necessary for the simulation are correctly described by Quantum Mechanics. However, the computational resources required to solve the quantum mechanical equations limits the use of Quantum Mechanics at most to a few hundreds of atoms and only to a small fraction of the available configurational space. This thesis presents the results of my research on the development of a new interatomic potential generation scheme, which we refer to as Gaussian Approximation Potentials. In our framework, the quantum mechanical potential energy surface is interpolated between a set of predetermined values at different points in atomic configurational space by a non-linear, non-parametric regression method, the Gaussian Process. To perform the fitting, we represent the atomic environments by the bispectrum, which is invariant to permutations of the atoms in the neighbourhood and to global rotations. The result is a general scheme, that allows one to generate interatomic potentials based on arbitrary quantum mechanical data. We built a series of Gaussian Approximation Potentials using data obtained from Density Functional Theory and tested the capabilities of the method. We showed that our models reproduce the quantum mechanical potential energy surface remarkably well for the group IV semiconductors, iron and gallium nitride. Our potentials, while maintaining quantum mechanical accuracy, are several orders of magnitude faster than Quantum Mechanical methods.

First-principles energetics of water: a many-body analysis

Standard forms of density-functional theory (DFT) have good predictive power for many materials, but are not yet fully satisfactory for solid, liquid and cluster forms of water. We use a many-body separation of the total energy into its 1-body, 2-body (2B) and beyond-2-body (B2B) components to analyze the deficiencies of two popular DFT approximations. We show how machine-learning methods make this analysis possible for ice structures as well as for water clusters. We find that the crucial energy balance between compact and extended geometries can be distorted by 2B and B2B errors, and that both types of first-principles error are important.