Block Tridiagonal Matrices in Electronic Structure Calculations
DEFF Research Database (Denmark)
Petersen, Dan Erik
This thesis focuses on some of the numerical aspects of the treatment of the electronic structure problem, in particular that of determining the ground state electronic density for the non–equilibrium Green’s function formulation of two–probe systems and the calculation of transmission in the Lan...
Cobalamins uncovered by modern electronic structure calculations
DEFF Research Database (Denmark)
Kepp, Kasper Planeta; Ryde, Ulf
electronic-structure calculations, in particular density functional methods, the understanding of the molecular mechanism of cobalamins has changed dramatically, going from a dominating view of trans-steric strain effects to a much more complex view involving an arsenal of catalytic strategies. Among these...
Multilevel domain decomposition for electronic structure calculations
Barrault, M; Hager, W W; Le Bris, C
2005-01-01
We introduce a new multilevel domain decomposition method (MDD) for electronic structure calculations within semi-empirical and Density Functional Theory (DFT) frameworks. This method iterates between local fine solvers and global coarse solvers, in the spirit of domain decomposition methods.
Multigrid Methods in Electronic Structure Calculations
Briggs, E L; Bernholc, J
1996-01-01
We describe a set of techniques for performing large scale ab initio calculations using multigrid accelerations and a real-space grid as a basis. The multigrid methods provide effective convergence acceleration and preconditioning on all length scales, thereby permitting efficient calculations for ill-conditioned systems with long length scales or high energy cut-offs. We discuss specific implementations of multigrid and real-space algorithms for electronic structure calculations, including an efficient multigrid-accelerated solver for Kohn-Sham equations, compact yet accurate discretization schemes for the Kohn-Sham and Poisson equations, optimized pseudo\\-potentials for real-space calculations, efficacious computation of ionic forces, and a complex-wavefunction implementation for arbitrary sampling of the Brillioun zone. A particular strength of a real-space multigrid approach is its ready adaptability to massively parallel computer architectures, and we present an implementation for the Cray-T3D with essen...
Isogeometric analysis in electronic structure calculations
Cimrman, Robert; Kolman, Radek; Tůma, Miroslav; Vackář, Jiří
2016-01-01
In electronic structure calculations, various material properties can be obtained by means of computing the total energy of a system as well as derivatives of the total energy w.r.t. atomic positions. The derivatives, also known as Hellman-Feynman forces, require, because of practical computational reasons, the discretized charge density and wave functions having continuous second derivatives in the whole solution domain. We describe an application of isogeometric analysis (IGA), a spline modification of finite element method (FEM), to achieve the required continuity. The novelty of our approach is in employing the technique of B\\'ezier extraction to add the IGA capabilities to our FEM based code for ab-initio calculations of electronic states of non-periodic systems within the density-functional framework, built upon the open source finite element package SfePy. We compare FEM and IGA in benchmark problems and several numerical results are presented.
Isogeometric analysis in electronic structure calculations
Czech Academy of Sciences Publication Activity Database
Cimrman, R.; Novák, M.; Kolman, Radek; Tůma, Miroslav; Vackář, Jiří
Ostrava: Ústav geoniky AV ČR, 2014 - (Blaheta, R.; Starý, J.; Sysalová, D.). s. 49-49 ISBN 978-80-86407-47-0. [Modelling 2014. 02.06.2014-06.06.2014, Rožnov pod Radhoštěm] R&D Projects: GA ČR(CZ) GAP101/12/2315; GA ČR(CZ) GAP108/11/0853 Institutional support: RVO:61388998 ; RVO:68378271 ; RVO:67985807 Keywords : isogeometric analysis * electronic structure calculations * density functional theory Subject RIV: JJ - Other Materials
Multilevel domain decomposition for electronic structure calculations
International Nuclear Information System (INIS)
We introduce a new multilevel domain decomposition method (MDD) for electronic structure calculations within semi-empirical and density functional theory (DFT) frameworks. This method iterates between local fine solvers and global coarse solvers, in the spirit of domain decomposition methods. Using this approach, calculations have been successfully performed on several linear polymer chains containing up to 40,000 atoms and 200,000 atomic orbitals. Both the computational cost and the memory requirement scale linearly with the number of atoms. Additional speed-up can easily be obtained by parallelization. We show that this domain decomposition method outperforms the density matrix minimization (DMM) method for poor initial guesses. Our method provides an efficient preconditioner for DMM and other linear scaling methods, variational in nature, such as the orbital minimization (OM) procedure
METHODS OF CALCULATING THE ELECTRONIC AND ATOMIC STRUCTURES OF INTERFACES
Sutton, A
1985-01-01
Methods of calculating the electronic and atomic structures of interfaces are described. An introduction to pseudopotentials and LCAO methods is given. Methods of calculating the electronic structure of an interface with a given atomic structure are considered. The feasibility of total energy calculations, in which the atomic and electronic structures are calculated simultaneously, is discussed.
Real-time feedback from iterative electronic structure calculations
Vaucher, Alain C; Reiher, Markus
2015-01-01
Real-time feedback from iterative electronic structure calculations requires to mediate between the inherently unpredictable execution times of the iterative algorithm employed and the necessity to provide data in fixed and short time intervals for real-time rendering. We introduce the concept of a mediator as a component able to deal with infrequent and unpredictable reference data to generate reliable feedback. In the context of real-time quantum chemistry, the mediator takes the form of a surrogate potential that has the same local shape as the first-principles potential and can be evaluated efficiently to deliver atomic forces as real-time feedback. The surrogate potential is updated continuously by electronic structure calculations and guarantees to provide a reliable response to the operator for any molecular structure. To demonstrate the application of iterative electronic structure methods in real-time reactivity exploration, we implement self-consistent semi-empirical methods as the data source and a...
Elongation method for electronic structure calculations of random DNA sequences.
Orimoto, Yuuichi; Liu, Kai; Aoki, Yuriko
2015-10-30
We applied ab initio order-N elongation (ELG) method to calculate electronic structures of various deoxyribonucleic acid (DNA) models. We aim to test potential application of the method for building a database of DNA electronic structures. The ELG method mimics polymerization reactions on a computer and meets the requirements for linear scaling computational efficiency and high accuracy, even for huge systems. As a benchmark test, we applied the method for calculations of various types of random sequenced A- and B-type DNA models with and without counterions. In each case, the ELG method maintained high accuracy with small errors in energy on the order of 10(-8) hartree/atom compared with conventional calculations. We demonstrate that the ELG method can provide valuable information such as stabilization energies and local densities of states for each DNA sequence. In addition, we discuss the "restarting" feature of the ELG method for constructing a database that exhaustively covers DNA species. PMID:26337429
Electronic structure of crystalline uranium nitride: LCAO DFT calculations
International Nuclear Information System (INIS)
The results of electronic structure calculations performed for the first time for crystalline uranium nitride and using a LCAO basis are discussed. For calculations we used the density functional method with the PW91 exchange correlation potential and a variety of relativistic core potentials for the uranium atom. The calculated atomization energy of the crystal agrees well with the experimental data and with the results of calculations with the plane wave basis. It is shown that a chemical bond in crystalline uranium nitride is a metal covalent bond. The metal component of the bond is due to the 5f electrons localized on the uranium atom and having energies near the Fermi level and the bottom of the conduction band. The covalent component of the chemical bond results from an overlap between the uranium 6d and 7s valence orbitals and the nitrogen 2p atomic orbitals. Inclusion of the 5f electrons in the core of the uranium atom introduces relatively minor changes in the calculated binding energy and electron density distribution
Electronic structure of crystalline uranium nitride: LCAO DFT calculations
International Nuclear Information System (INIS)
The results of the first LCAO DFT calculations of cohesive energy, band structure and charge distribution in uranium nitride (UN) crystal are presented and discussed. The calculations are made with the uranium atom relativistic effective core potentials, including 60, 78 and 81 electrons in the core. It is demonstrated that the chemical bonding in UN crystal has a metallic-covalent nature. Three 5f-electrons are localized on the U atom and occupy the states near the Fermi level. The metallic nature of the crystal is due to the f-character of both the valence-band top and the conduction-band bottom. The covalent bonds are formed by the interaction of 7s- and 6d-states of the uranium atom with the 2p-states of the nitrogen atom. It is shown that the inclusion of 5f-electrons in the atomic core introduces small changes in the calculated cohesive energy of UN crystal and electron charge distribution. However, the inclusion of 5s-, 5p-, 5d-electrons in the valence shell allows the better agreement with the calculated and experimental cohesive-energy value. (copyright 2008 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim) (orig.)
Electronic Structure and Molecular Dynamics Calculations for KBH4
Papaconstantopoulos, Dimitrios; Shabaev, Andrew; Hoang, Khang; Mehl, Michael; Kioussis, Nicholas
2012-02-01
In the search for hydrogen storage materials, alkali borohydrides MBH4 (M=Li, Na, K) are especially interesting because of their light weight and the high number of hydrogen atoms per metal atom. Electronic structure calculations can give insights into the properties of these complex hydrides and provide understanding of the structural properties and of the bonding of hydrogen. We have performed first-principles density-functional theory (DFT) and tight-binding (TB) calculations for KBH4 in both the high temperature (HT) and low temperature (LT) phases to understand its electronic and structural properties. Our DFT calculations were carried out using the VASP code. The results were then used as a database to develop a tight-binding Hamiltonian using the NRL-TB method. This approach allowed for computationally efficient calculations of phonon frequencies and elastic constants using the static module of the NRL-TB, and also using the molecular dynamics module to calculate mean-square displacements and formation energies of hydrogen vacancies.
Thick-Restart Lanczos Method for Electronic Structure Calculations
International Nuclear Information System (INIS)
This paper describes two recent innovations related to the classic Lanczos method for eigenvalue problems, namely the thick-restart technique and dynamic restarting schemes. Combining these two new techniques we are able to implement an efficient eigenvalue problem solver. This paper will demonstrate its effectiveness on one particular class of problems for which this method is well suited: linear eigenvalue problems generated from non-self-consistent electronic structure calculations
Statistical learning for alloy design from electronic structure calculations
Broderick, Scott R.
The objective of this thesis is to explore how statistical learning methods can contribute to the interpretation and efficacy of electronic structure calculations. This study develops new applications of statistical learning and data mining methods to both semi-empirical and density functional theory (DFT) calculations. Each of these classes of electronic structure calculations serves as templates for different data driven discovery strategies for materials science applications. In our study of semi-empirical methods, we take advantage of the ability of data mining methods to quantitatively assess high dimensional parameterization schemes. The impact of this work includes the development of accelerated computational schemes for developing reduced order models. Another application is the use of these informatics based techniques to serve as a means for estimating parameters when data for such calculations are not available. Using density of states (DOS) spectra derived from DFT calculations we have demonstrated the classification power of singular value decomposition methods to accurately develop structural and stoichiometric classifications of compounds. Building on this work we have extended this analytical strategy to apply the predictive capacity of informatics methods to develop a new and far more robust modeling approach for DOS spectra, addressing an issue that has gone relatively unchallenged over two decades. By exploring a diverse array of materials systems (metals, ceramics, different crystal structures) this work has laid the foundations for expanding the linkages between statistical learning and statistical thermodynamics. The results of this work provide exciting new opportunities in computational based design of materials that have not been explored before.
Electronic Structure of Silicon Nanowires Matrix from Ab Initio Calculations.
Monastyrskii, Liubomyr S; Boyko, Yaroslav V; Sokolovskii, Bogdan S; Potashnyk, Vasylyna Ya
2016-12-01
An investigation of the model of porous silicon in the form of periodic set of silicon nanowires has been carried out. The electronic energy structure was studied using a first-principle band method-the method of pseudopotentials (ultrasoft potentials in the basis of plane waves) and linearized mode of the method of combined pseudopotentials. Due to the use of hybrid exchange-correlation potentials (B3LYP), the quantitative agreement of the calculated value of band gap in the bulk material with experimental data is achieved. The obtained results show that passivation of dangling bonds with hydrogen atoms leads to substantial transformation of electronic energy structure. At complete passivation of the dangling silicon bonds by hydrogen atoms, the band gap value takes the magnitude which substantially exceeds that for bulk silicon. The incomplete passivation gives rise to opposite effect when the band gap value decreases down the semimetallic range. PMID:26768147
Multi-million atom electronic structure calculations for quantum dots
Usman, Muhammad
Quantum dots grown by self-assembly process are typically constructed by 50,000 to 5,000,000 structural atoms which confine a small, countable number of extra electrons or holes in a space that is comparable in size to the electron wavelength. Under such conditions quantum dots can be interpreted as artificial atoms with the potential to be custom tailored to new functionality. In the past decade or so, these nanostructures have attracted significant experimental and theoretical attention in the field of nanoscience. The new and tunable optical and electrical properties of these artificial atoms have been proposed in a variety of different fields, for example in communication and computing systems, medical and quantum computing applications. Predictive and quantitative modeling and simulation of these structures can help to narrow down the vast design space to a range that is experimentally affordable and move this part of nanoscience to nano-Technology. Modeling of such quantum dots pose a formidable challenge to theoretical physicists because: (1) Strain originating from the lattice mismatch of the materials penetrates deep inside the buffer surrounding the quantum dots and require large scale (multi-million atom) simulations to correctly capture its effect on the electronic structure, (2) The interface roughness, the alloy randomness, and the atomistic granularity require the calculation of electronic structure at the atomistic scale. Most of the current or past theoretical calculations are based on continuum approach such as effective mass approximation or k.p modeling capturing either no or one of the above mentioned effects, thus missing some of the essential physics. The Objectives of this thesis are: (1) to model and simulate the experimental quantum dot topologies at the atomistic scale; (2) to theoretically explore the essential physics i.e. long range strain, linear and quadratic piezoelectricity, interband optical transition strengths, quantum confined
An Extensive Database of Electronic Structure Calculations between Transition Metals
Sayed, Shereef; Papaconstantopoulos, Dimitrios
Density Functional Theory and its derived application methods, such as the Augmented Plane Wave (APW) method, have shown great success in predicting the fundamental properties of materials. In this work, we apply the APW method to explore the properties of diatomic pairs of transition metals in the CsCl structure, for all possible combinations. A total of 435 compounds have been studied. The predicted Density of States, and Band Structures are presented, along with predicted electron-phonon coupling and Stoner Criterion, in order to identify potential new superconducting or ferromagnetic materials. This work is performed to demonstrate the concept of ``high-throughput'' calculations at the crossing-point of ``Big Data'' and materials science. Us Dept of Energy.
Electronic structure calculations toward new potentially AChE inhibitors
de Paula, A. A. N.; Martins, J. B. L.; Gargano, R.; dos Santos, M. L.; Romeiro, L. A. S.
2007-10-01
The main purpose of this study was the use of natural non-isoprenoid phenolic lipid of cashew nut shell liquid from Anacardium occidentale as lead material for generating new potentially candidates of acetylcholinesterase inhibitors. Therefore, we studied the electronic structure of 15 molecules derivatives from the cardanol using the following groups: methyl, acetyl, N, N-dimethylcarbamoyl, N, N-dimethylamine, N, N-diethylamine, piperidine, pyrrolidine, and N-benzylamine. The calculations were performed at RHF level using 6-31G, 6-31G(d), 6-31+G(d) and 6-311G(d,p) basis functions. Among the proposed compounds we found that the structures with substitution by acetyl, N, N-dimethylcarbamoyl, N, N-dimethylamine, and pyrrolidine groups were better correlated to rivastigmine indicating possible activity.
Gradient type optimization methods for electronic structure calculations
Zhang, Xin; Wen, Zaiwen; Zhou, Aihui
2013-01-01
The density functional theory (DFT) in electronic structure calculations can be formulated as either a nonlinear eigenvalue or direct minimization problem. The most widely used approach for solving the former is the so-called self-consistent field (SCF) iteration. A common observation is that the convergence of SCF is not clear theoretically while approaches with convergence guarantee for solving the latter are often not competitive to SCF numerically. In this paper, we study gradient type methods for solving the direct minimization problem by constructing new iterations along the gradient on the Stiefel manifold. Global convergence (i.e., convergence to a stationary point from any initial solution) as well as local convergence rate follows from the standard theory for optimization on manifold directly. A major computational advantage is that the computation of linear eigenvalue problems is no longer needed. The main costs of our approaches arise from the assembling of the total energy functional and its grad...
Adaptations in Electronic Structure Calculations in Heterogeneous Environments
Energy Technology Data Exchange (ETDEWEB)
Talamudupula, Sai [Iowa State Univ., Ames, IA (United States)
2011-01-01
Modern quantum chemistry deals with electronic structure calculations of unprecedented complexity and accuracy. They demand full power of high-performance computing and must be in tune with the given architecture for superior e ciency. To make such applications resourceaware, it is desirable to enable their static and dynamic adaptations using some external software (middleware), which may monitor both system availability and application needs, rather than mix science with system-related calls inside the application. The present work investigates scienti c application interlinking with middleware based on the example of the computational chemistry package GAMESS and middleware NICAN. The existing synchronous model is limited by the possible delays due to the middleware processing time under the sustainable runtime system conditions. Proposed asynchronous and hybrid models aim at overcoming this limitation. When linked with NICAN, the fragment molecular orbital (FMO) method is capable of adapting statically and dynamically its fragment scheduling policy based on the computing platform conditions. Signi cant execution time and throughput gains have been obtained due to such static adaptations when the compute nodes have very di erent core counts. Dynamic adaptations are based on the main memory availability at run time. NICAN prompts FMO to postpone scheduling certain fragments, if there is not enough memory for their immediate execution. Hence, FMO may be able to complete the calculations whereas without such adaptations it aborts.
Large Scale Electronic Structure Calculations using Quantum Chemistry Methods
Scuseria, Gustavo E.
1998-03-01
This talk will address our recent efforts in developing fast, linear scaling electronic structure methods for large scale applications. Of special importance is our fast multipole method( M. C. Strain, G. E. Scuseria, and M. J. Frisch, Science 271), 51 (1996). (FMM) for achieving linear scaling for the quantum Coulomb problem (GvFMM), the traditional bottleneck in quantum chemistry calculations based on Gaussian orbitals. Fast quadratures(R. E. Stratmann, G. E. Scuseria, and M. J. Frisch, Chem. Phys. Lett. 257), 213 (1996). combined with methods that avoid the Hamiltonian diagonalization( J. M. Millam and G. E. Scuseria, J. Chem. Phys. 106), 5569 (1997) have resulted in density functional theory (DFT) programs that can be applied to systems containing many hundreds of atoms and ---depending on computational resources or level of theory-- to many thousands of atoms.( A. D. Daniels, J. M. Millam and G. E. Scuseria, J. Chem. Phys. 107), 425 (1997). Three solutions for the diagonalization bottleneck will be analyzed and compared: a conjugate gradient density matrix search (CGDMS), a Hamiltonian polynomial expansion of the density matrix, and a pseudo-diagonalization method. Besides DFT, our near-field exchange method( J. C. Burant, G. E. Scuseria, and M. J. Frisch, J. Chem. Phys. 105), 8969 (1996). for linear scaling Hartree-Fock calculations will be discussed. Based on these improved capabilities, we have also developed programs to obtain vibrational frequencies (via analytic energy second derivatives) and excitation energies (through time-dependent DFT) of large molecules like porphyn or C_70. Our GvFMM has been extended to periodic systems( K. N. Kudin and G. E. Scuseria, Chem. Phys. Lett., in press.) and progress towards a Gaussian-based DFT and HF program for polymers and solids will be reported. Last, we will discuss our progress on a Laplace-transformed \\cal O(N^2) second-order pertubation theory (MP2) method.
Quasiparticle GW calculations within the GPAW electronic structure code
DEFF Research Database (Denmark)
Hüser, Falco
The GPAW electronic structure code, developed at the physics department at the Technical University of Denmark, is used today by researchers all over the world to model the structural, electronic, optical and chemical properties of materials. They address fundamental questions in material science...... properties are to a large extent governed by the physics on the atomic scale, that means pure quantum mechanics. For many decades, Density Functional Theory has been the computational method of choice, since it provides a fairly easy and yet accurate way of determining electronic structures and related...... respect to the system one wants to investigate by choosing a certain functional or by tuning parameters. A succesful alternative is the so-called GW approximation. It is mathematically precise and gives a physically well-founded description of the complicated electron interactions in terms of screening...
Electronic structure calculations on defects and impurities in semiconductors
International Nuclear Information System (INIS)
Self-consistent tight-binding methods are developed and used to investigate a number of defects in silicon and in four different lll-V compound semiconductors, GaP, GaAs, GaSb, and InP. The wave functions of defect states are calculated with the use of the 'largest weight method'. The hyperfine interaction parameters are in turn derived from the calculated wave functions of gap states. Substitutional phosphorus-vacancy (V-P), interstitial hydrogen-vacancy (V-H), and substitutional phosphorus-interstitial hydrogen-vacancy (V-H-P) complexes in silicon are studied in detail. We demonstrate that in the V-H and V-H-P complexes the electrical activity of the silicon dangling bonds are well passivated by phosphorus atoms at substitutional position through Coulomb attractions and by hydrogen atoms at bonding positions through strong orbital interactions, and all the remaining electrical activity in these complexes can be accounted for by those silicon dangling bonds which have neither been attached by hydrogen atoms nor replaced by phosphorus dangling bonds. In the V-P complexes, the phosphorus-dangling-bond states are found to interact with the silicon-dangling-bond states and, therefore, significant contributions from phosphorus dangling bonds to the electrically active gap states of these complexes are found. Detailed investigations on neutral and charged vacancies and divacancies in GaP, GaAs, GaSb, and InP are carried out. We find that with the Fermi level at an energy around the midpoint of the fundamental band gap, an isolated cation and an isolated anion vacancy in each of the four compound semiconductors have opposite charge states and may thus attract each other. We demonstrate that a divacancy in these compounds can have many charge states, and can introduce many energy levels into the fundamental band gap, revealing the complex nature of the defect. A simple one-electron molecular-orbital model accounting for the basic feature of the electronic structure of a
2007-01-01
Density functional calculations of electronic structure, total energy, structural distortions, and magnetism for hydrogenated single-layer, bilayer, and multi-layer graphene are performed. It is found that hydrogen-induced magnetism can survives only at very low concentrations of hydrogen (single-atom regime) whereas hydrogen pairs with optimized structure are usually nonmagnetic. Chemisorption energy as a function of hydrogen concentration is calculated, as well as energy barriers for hydrog...
Boukhvalov, D W; Katsnelson, M. I.; Lichtenstein, A. I.
2008-01-01
Density functional calculations of electronic structure, total energy, structural distortions, and magnetism for hydrogenated single-layer, bilayer, and multi-layer graphene are performed. It is found that hydrogen-induced magnetism can survives only at very low concentrations of hydrogen (single-atom regime) whereas hydrogen pairs with optimized structure are usually nonmagnetic. Chemisorption energy as a function of hydrogen concentration is calculated, as well as energy barriers for hydrog...
Unfolding method for first-principles LCAO electronic structure calculations
Lee, Chi-Cheng; Yamada-Takamura, Yukiko; Ozaki, Taisuke
2013-08-01
Unfolding the band structure of a supercell to a normal cell enables us to investigate how symmetry breakers such as surfaces and impurities perturb the band structure of the normal cell. We generalize the unfolding method, originally developed based on Wannier functions, to the linear combination of atomic orbitals (LCAO) method, and present a general formula to calculate the unfolded spectral weight. The LCAO basis set is ideal for the unfolding method because the basis functions allocated to each atomic species are invariant regardless of the existence of surface and impurity. The unfolded spectral weight is well defined by the property of the LCAO basis functions. In exchange for the property, the non-orthogonality of the LCAO basis functions has to be taken into account. We show how the non-orthogonality can be properly incorporated in the general formula. As an illustration of the method, we calculate the dispersive quantized spectral weight of a ZrB2 slab and show strong spectral broadening in the out-of-plane direction, demonstrating the usefulness of the unfolding method.
Unfolding method for first-principles LCAO electronic structure calculations
International Nuclear Information System (INIS)
Unfolding the band structure of a supercell to a normal cell enables us to investigate how symmetry breakers such as surfaces and impurities perturb the band structure of the normal cell. We generalize the unfolding method, originally developed based on Wannier functions, to the linear combination of atomic orbitals (LCAO) method, and present a general formula to calculate the unfolded spectral weight. The LCAO basis set is ideal for the unfolding method because the basis functions allocated to each atomic species are invariant regardless of the existence of surface and impurity. The unfolded spectral weight is well defined by the property of the LCAO basis functions. In exchange for the property, the non-orthogonality of the LCAO basis functions has to be taken into account. We show how the non-orthogonality can be properly incorporated in the general formula. As an illustration of the method, we calculate the dispersive quantized spectral weight of a ZrB2 slab and show strong spectral broadening in the out-of-plane direction, demonstrating the usefulness of the unfolding method. (paper)
Accelerating VASP electronic structure calculations using graphic processing units
Hacene, Mohamed
2012-08-20
We present a way to improve the performance of the electronic structure Vienna Ab initio Simulation Package (VASP) program. We show that high-performance computers equipped with graphics processing units (GPUs) as accelerators may reduce drastically the computation time when offloading these sections to the graphic chips. The procedure consists of (i) profiling the performance of the code to isolate the time-consuming parts, (ii) rewriting these so that the algorithms become better-suited for the chosen graphic accelerator, and (iii) optimizing memory traffic between the host computer and the GPU accelerator. We chose to accelerate VASP with NVIDIA GPU using CUDA. We compare the GPU and original versions of VASP by evaluating the Davidson and RMM-DIIS algorithms on chemical systems of up to 1100 atoms. In these tests, the total time is reduced by a factor between 3 and 8 when running on n (CPU core + GPU) compared to n CPU cores only, without any accuracy loss. © 2012 Wiley Periodicals, Inc.
Electronic structure calculations on lithium battery electrolyte salts.
Johansson, Patrik
2007-03-28
New lithium salts for non-aqueous liquid, gel and polymeric electrolytes are crucial due to the limiting role of the electrolyte in modern lithium batteries. The solvation of any lithium salt to form an electrolyte solution ultimately depends on the strength of the cation-solvent vs. the cation-anion interaction. Here, the latter is probed via HF, B3LYP and G3 theory gas-phase calculations for the dissociation reaction: LiX Li(+) + X(-). Furthermore, a continuum solvation method (C-PCM) has been applied to mimic solvent effects. Anion volumes were also calculated to facilitate a discussion on ion conductivities and cation transport numbers. Judging from the present results, synthesis efforts should target heterocyclic anions with a size of ca. 150 A(3) molecule(-1) to render new highly dissociative lithium salts that result in electrolytes with high cation transport numbers. PMID:17356757
Ab Initio factorized LCAO calculation of the electronic structure of α-SiO2
International Nuclear Information System (INIS)
The authors report on the results of calculations of the electronic structure of α-quartz that were made using first principles, factorized linear combination of atomic orbitals method. Results were obtained for the primitive 9-atom, and orthorhombic 18- and 72-atom unit cells. Application of this method to the calculation of the electronic structure of the neutral oxygen vacancy in α-quartz is discussed and results obtained using a 72-atom unit cell are given
High Resolution Measurements and Electronic Structure Calculations of a Diazanaphthalene
Gruet, Sébastien; Goubet, Manuel; Pirali, Olivier
2014-06-01
Polycyclic Aromatic Hydrocarbons (PAHs) have long been suspected to be the carriers of so called Unidentified Infrared Bands (UIBs). Most of the results published in the literature report rotationally unresolved spectra of pure carbon as well as heteroatom-containing PAHs species. To date for this class of molecules, the principal source of rotational informations is ruled by microwave (MW) spectroscopy while high resolution measurements reporting rotational structure of the infrared (IR) vibrational bands are very scarce. Recently, some high resolution techniques provided interesting new results to rotationally resolve the IR and far-IR bands of these large carbonated molecules of astrophysical interest. One of them is to use the bright synchrotron radiation as IR continuum source of a high resolution Fourier transform (FTIR) spectrometer. We report the very complementary analysis of the [1,6] naphthyridine (a N-bearing PAH) for which we recorded the microwave spectrum at the PhLAM laboratory (Lille) and the high resolution far-infrared spectrum on the AILES beamline at synchrotron facility SOLEIL. MW spectroscopy provided highly accurate rotational constants in the ground state to perform Ground State Combinations Differences (GSCD) allowing the analysis of the two most intense FT-FIR bands in the 50-900 wn range. Moreover, during this presentation the negative value of the inertial defect in the GS of the molecule will be discussed. A. Leger, J. L. Puget, Astron. Astrophys. 137, L5-L8 (1984) L. J. Allamandola et al. Astrophys. J. 290, L25-L28 (1985). Z. Kisiel et al. J. Mol. Spectrosc. 217, 115 (2003) S. Thorwirth et al. Astrophys. J. 662, 1309 (2007) D. McNaughton et al. J. Chem. Phys. 124, 154305 (2011). S. Albert et al. Faraday Discuss. 150, 71-99 (2011) B. E. Brumfield et al. Phys. Chem. Lett. 3, 1985-1988 (2012) O. Pirali et al. Phys. Chem. Chem. Phys. 15, 10141 (2013).
Energy Technology Data Exchange (ETDEWEB)
Larsen, Ross E.
2016-05-12
We introduce two simple tight-binding models, which we call fragment frontier orbital extrapolations (FFOE), to extrapolate important electronic properties to the polymer limit using electronic structure calculations on only a few small oligomers. In particular, we demonstrate by comparison to explicit density functional theory calculations that for long oligomers the energies of the highest occupied molecular orbital (HOMO), the lowest unoccupied molecular orbital (LUMO), and of the first electronic excited state are accurately described as a function of number of repeat units by a simple effective Hamiltonian parameterized from electronic structure calculations on monomers, dimers and, optionally, tetramers. For the alternating copolymer materials that currently comprise some of the most efficient polymer organic photovoltaic devices one can use these simple but rigorous models to extrapolate computed properties to the polymer limit based on calculations on a small number of low-molecular-weight oligomers.
Calculation of the valence electron structures of alloying cementite and its biphase interface
Institute of Scientific and Technical Information of China (English)
无
2001-01-01
The valence electron structures of alloying cementite θ-(Fe, M)3C and ε-(Fe, M)3C andthose of the biphase interfaces between them and α-Fe are calculated with Yu's empirical electrontheory of solid and molecules. The calculation results accord with the actual behavior of alloys.
Electronic structure of the heavy fermion superconductor Ce2PdIn8: Experiment and calculations
International Nuclear Information System (INIS)
The electronic structure of a heavy-fermion superconductor Ce2PdIn8 was investigated by means of X-ray photoelectron spectroscopy (XPS) and ab initio density functional band structure calculations. The Ce 3d core-level XPS spectra point to stable trivalent configuration of Ce atoms that is also reproduced in the band structure calculations within the generalized gradient approximation GGA+U approach. Analysis of the 3d9f2 weight in the 3d XPS spectra within the Gunnarsson-Schönhammer model suggests that the onsite hybridization energy between Ce 4f and the conduction band states, Δfs, is ∼120 meV, which is about 30 meV larger than Δfs in isostructural Ce2TIn8 compounds with T = Co, Rh, and Ir. Taking into account a Coulomb repulsion U on both the Ce 4f and Pd 4d states in electronic band structure calculations, a satisfactory agreement was found between the calculated density of states (DOS) and the measured valence band XPS spectra. - Highlights: • XPS data validated strong electronic correlations in superconducting Ce2PdIn8. • DFT calculations reproduced XPS spectra measured for Ce2PdIn8. • Crucial role of Pd d electrons in the HF behavior of Ce2PdIn8 was established
Ab-initio calculations of electronic structure and optical properties of TiAl alloy
Energy Technology Data Exchange (ETDEWEB)
Hussain, Altaf [Department of Physics, The Islamia University of Bahawalpur, Bahawalpur 63120 (Pakistan); Sikandar Hayat, Sardar, E-mail: sikandariub@yahoo.co [Department of Physics, Hazara University, Mansehra 21300 (Pakistan); Choudhry, M.A. [Department of Physics, The Islamia University of Bahawalpur, Bahawalpur 63120 (Pakistan)
2011-05-01
The electronic structures and optical properties of TiAl intermetallic alloy system are studied by the first-principle orthogonalized linear combination of atomic orbitals method. Results on the band structure, total and partial density of states, localization index, effective atomic charges, and optical conductivity are presented and discussed in detail. Total density of states spectra reveal that (near the Fermi level) the majority of the contribution is from Ti-3d states. The effective charge calculations show an average charge transfer of 0.52 electrons from Ti to Al in primitive cell calculations of TiAl alloy. On the other hand, calculations using supercell approach reveal an average charge transfer of 0.48 electrons from Ti to Al. The localization index calculations, of primitive cell as well as of supercell, show the presence of relatively localized states even above the Fermi level for this alloy. The calculated optical conductivity spectra of TiAl alloy are rich in structures, showing the highest peak at 5.73 eV for supercell calculations. Calculations of the imaginary part of the linear dielectric function show a prominent peak at 5.71 eV and a plateau in the range 1.1-3.5 eV.
Atomic and Electronic Structures of C_60+BN Nanopeapods from ab initio Pseudopotential Calculations
Trave, Andrea; Ribeiro, Filipe; Louie, Steven G.; Cohen, Marvin L.
2004-03-01
Nanopeapods are structures of nanometric size consisting of an external carbon nanotube encapsulating a chain or complex array of fullerenes. Recent calculations and experiments have proven that nanopeapods can be obtained assembling fullerenes within boron nitride nanotubes, creating novel materials of possible interest for electronic transport applications. To improve the understanding of the properties of these composite systems, as compared to empty nanotubes and carbon nanopeapods, ab-initio total energy calculations have been performed within the pseudopotential Density Functional Theory in local density approximation. Results of these calculations on the energetics and geometrical deformations involved in the encapsulation will be presented, followed by a discussion of the consequences on the electronic structures of these systems, with particular focus on aspects relevant to electronic transport phenomena. This work is supported by NFS (Grant DMR00-87088) and DOE (Contract DE-AC03-76SF00098), using computational resources at NERSC and NPACI.
Self-consistent cluster-embedding calculation method and the calculated electronic structure of NiO
International Nuclear Information System (INIS)
The self-consistent cluster-embedding method is discussed theoretically. A definition of the total energy for an embedded cluster has been introduced. The method has two advantages. (i) It can describe both localized and band properties, including their excitations. (ii) It can give a good description of the magnetic properties for both spin-ordered and spin-disordered states. The electronic structure of NiO is studied using a high-quality basis set to calculate the electronic structure of a small embedded cluster and an antiferromagnetic insulating ground state is obtained. The picture has both localized and band properties. A small energy gap separates the unoccupied and occupied nickel 3d orbitals which are well localized. Each 3d orbital is attached to a particular nickel ion. Below the 3d levels are two diffuse oxygen 2p bands, and above the 3d levels are oxygen 3s, nickel 4s, and oxygen 3p bands. Experimental data concerning photoemission and optical absorption can be interpreted naturally. The spin magnetic moment of the nickel ion is calculated correctly. The simulation of the spin-disordered state shows that NiO remains as an insulator in the paramagnetic state. The Neel temperature of NiO is calculated directly to give a reasonable result. The Hubbard U parameter for nickel 3d electrons is estimated. The calculation shows that the excited nickel 3d electrons are also well localized and the overlaps are less than 4.5%. We propose the following: The overlap of the excited 3d electrons is too small to form a metallic band, but the overlap is sufficient for the ''hole'' to migrate through the crystal. In this sense, NiO is a charge-transfer insulator with a gap of about 4 eV (mostly from oxygen to nickel)
Ab initio calculations of electronic structure of anatase TiO2
Institute of Scientific and Technical Information of China (English)
Chen Qiang; Cao Hong-Hong
2004-01-01
This paper presents the results of the self-consistent calculations on the electronic structure of anatase phase of TiO2. The calculations were performed using the full potential-linearized augmented plane wave method (FP-LAPW)in the framework of the density functional theory (DFT) with the generalized gradient approximation (GGA). The fully optimized structure, obtained by minimizing the total energy and atomic forces, is in good agreement with experiment.We also calculated the band structure and the density of states. In particular, the calculated band structure prefers an indirect transition between wlence and conduction bands of anatase TiO2, which may be helpful for clarifying the ambiguity in other theoretical works.
International Nuclear Information System (INIS)
The recoil proton polarization for the quasielastic electron-proton scattering is represented as a contraction of the electron structure and the hard part of the polarization dependent contribution into cross-section. The calculation of the hard part with first order radiative correction is performed. The obtained representation includes the leading radiative corrections in all orders of perturbation theory and the main part of the second order next-to-leading ones
International Nuclear Information System (INIS)
The structural, electronic and magnetic properties of free standing Au-Pd bimetallic atomic chain is studied using ab-initio method. It is found that electronic and magnetic properties of chains depend on position of atoms and number of atoms. Spin polarization factor for different atomic configuration of atomic chain is calculated predicting a half metallic behavior. It suggests a total spin polarised transport in these chains
GPAW - massively parallel electronic structure calculations with Python-based software
DEFF Research Database (Denmark)
Enkovaara, Jussi; Romero, Nichols A.; Shende, Sameer;
2011-01-01
Electronic structure calculations are a widely used tool in materials science and large consumer of supercomputing resources. Traditionally, the software packages for these kind of simulations have been implemented in compiled languages, where Fortran in its different versions has been the most...
Electronic structure of Co-phthalocyanine calculated by GGA+U and hybrid functional methods
International Nuclear Information System (INIS)
Graphical abstract: Electronic structure of Co-phthalocyanine molecule has been calculated using GGA+U and B3LYP methods. The results are in good agreement with experimental observations. Abstract: Electronic structure calculations have been performed for the Co-phthalocyanine molecule using density functional theory (DFT) within the framework of Generalized Gradient Approximation (GGA). The electronic correlation in Co 3d orbitals is treated in terms of the GGA+U method in the framework of the Hubbard model. We find that for U = 6 eV, the calculated structural parameters as well as the spectral features are in good agreement with the experimental findings. From our calculation both the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) are dominated by the pyrrole carbon, with a HOMO-LUMO gap of about 1.4 eV. The GGA+U results obtained with U = 6 eV compare reasonably well with the calculations performed using Gaussian basis set and hybrid functionals in terms of ground state geometry, spin state and spectral features. The calculated valence band photoemission spectrum is in quite good agreement with the recently published experimental results.
Full potential calculation of structural, electronic and optical properties of KMgF3
International Nuclear Information System (INIS)
A theoretical study of the structural, electronic and optical properties of KMgF3 is presented using the full-potential linearized augmented plane wave method (FP-LAPW). In this approach, the local density approximation was used for the exchange-correlation potentials. First, we present the main features of the structural and electronic properties of this compound, where the electronic band structure shows that the fundamental energy gap is indirect. The contribution of the different bands was analysed from the total and partial density of states curves. The different interband transitions have been determined from the imaginary part of the dielectric function. The results are compared with previous calculations and with experimental measurements. The present work also deals with the behaviour of electronic properties, namely, the energy band gaps, and the valence bandwidth of KMgF3 subject of hydrostatic pressures up to 30 GPa
Large scale electronic structure calculations in the study of the condensed phase
Dam, H.J.J. van; Guest, M.F.; Sherwood, P.; Thomas, J.M.H.; van Lenthe, J.H.; van Lingen, J.N.J.; Bailey, C. L.; Bush, I. J.
2006-01-01
We consider the role that large-scale electronic structure computations can now play in the modelling of the condensed phase. To structure our analysis, we consider four distict ways in which today's scientific targets can be re-scoped to take advantage of advances in computing resources: 1. time to solution-performing the same calculation, with delivery of the simulation in shorter elapsed time; 2. Size-applying today's methods to a more extensive problem; 3. Accuracy-replacing current physi...
Flocke, N; Lotrich, V
2008-12-01
For the new parallel implementation of electronic structure methods in ACES III (Lotrich et al., in preparation) the present state-of-the-art algorithms for the evaluation of electronic integrals and their generalized derivatives were implemented in new object oriented codes with attention paid to efficient execution on modern processors with a deep hierarchy of data storage including multiple caches and memory banks. Particular attention has been paid to define proper integral blocks as basic building objects. These objects are stand-alone units and are no longer tied to any specific software. They can hence be used by any quantum chemistry code without modification. The integral blocks can be called at any time and in any sequence during the execution of an electronic structure program. Evaluation efficiency of these integral objects has been carefully tested and it compares well with other fast integral programs in the community. Correctness of the objects has been demonstrated by several application runs on real systems using the ACES III program. PMID:18496792
Electronic Structure of KFe2Se2 from First-Principles Calculations
International Nuclear Information System (INIS)
Electronic structures and magnetic properties for iron-selenide KFe2Se2 are studied by first-principles calculations. The ground state is collinear antiferromagnetic with calculated 2.26μB magnetic moment on Fe atoms; and the J1 and J2 coupling strengths are calculated to be 0.038eV and 0.029eV. The states around EF are dominated by the Fe 3d orbitals which hybridize noticeably to the Se 4p orbitals. While the band structure of KFe2Se2 is similar to a heavily electron-doped BaFe2AS2 or FeSe system, the Fermi surface of KFe2Se2 is much closer to the FeSe system since the electron sheets around M are symmetric with respect to x—y exchange. These features, as well as the absence of Fermi surface nesting, suggest that the parent KFe2Se2 could be regarded as an electron doped FeSe system with possible local moment magnetism. (condensed matter: electronic structure, electrical, magnetic, and optical properties)
Electronic structure calculations for Zn S xSe1-x
International Nuclear Information System (INIS)
Energy band gaps and electron effective mass as well as their composition dependence are the most critical parameters for band structure calculations of semiconductor alloys. Therefore, an accurate knowledge of these parameters is very important. Unfortunately, there is a limited experimental and theoretical information in the literature regarding the electronic band parameters for zinc blende Zn S xSe1-x. This has incited US to carry out such calculations. For this purpose, we have used the empirical pseudo potential method within the Virtual Crystal Approximation and the effect of compositional disorder is treated as an effective potential. The band gap variation versus sulfur concentration x shows two different behaviors: clear diminution of gap energy for low concentrations, and quasi-linear behavior with a small bowing for large values of x. Furthermore, the calculated effective mass shows that the disorder is not only compositional but also structural
Landau, Arie; Kaprálová-Žďánská, Petra Ruth; Moiseyev, Nimrod
2015-01-01
Complex eigenvalues, resonances, play an important role in large variety of fields in physics and chemistry. For example, in cold molecular collision experiments and electron scattering experiments, autoionizing and pre-dissociative metastable resonances are generated. However, the computation of complex resonance eigenvalues is difficult, since it requires severe modifications of standard electronic structure codes and methods. Here we show how resonance eigenvalues, positions and widths, can be calculated using the standard, widely used, electronic-structure packages. Our method enables the calculations of the complex resonance eigenvalues by using analytical continuation procedures (such as Pad\\'{e}). The key point in our approach is the existence of narrow analytical passages from the real axis to the complex energy plane. In fact, the existence of these analytical passages relies on using finite basis sets. These passages become narrower as the basis set becomes more complete, whereas in the exact limit,...
First-principle calculation of GaAs electronic structure by doping with Mn and P
International Nuclear Information System (INIS)
The geometry structure of Mn, P doped 64-atom supercell of GaAs (Ga1-xMnxAs1-yPy) was optimized and studied by the self-consistent full-potential linearized augmented plane wave method (FPLAPW) based on the density functional theory (DFT). Cell parameters of both doped and undoped were calculated theoretically. Band structure, binding energies, partial density of states, mulliken charges, and electron density different of doped GaAs crystals were calculated and analyzed in detail. The results revealed that the both doped compounds are true half-metallic ferromagnets and the near distance doped one is a stable ground state.
Structural, Elastic, and Electronic Properties of ReB2: A First-Principles Calculation
Directory of Open Access Journals (Sweden)
Run Long
2008-02-01
Full Text Available The structural, elastic, and electronic properties of the hard material ReB2 have been investigated by means of density functional theory. The calculated equilibrium structural parameters of ReB2 are in agreement with the experimental results. Our result of bulk modulus shows that it is a low compressible material. Furthermore, the elastic anisotropy is discussed by investigating the elastic stiffness constants. The charge density and the electronic properties indicate that the covalent bonding of Re-B and B-B plays an important role in formation of a hard material. The good metallicity and hardness of ReB2 might serve as hard conductors.
Xiao, Ling-Ping; Zeng, Zhi; Chen, Xiao-Jia
2016-06-01
The pressure effect on the geometrical and electronic structures of crystalline naphthalene is calculated up to 30 GPa by performing density functional calculations. The lattice parameters a, b, and c, decrease by 1.77 Å (-20.4%), 0.85 Å (-14.1%), and 0.91 Å (-8.2%), respectively, while the monoclinic angle β increases by 3.95° in this pressure region. At the highest pressure of 30 GPa the unit cell volume decreases by 62.7%. The detailed analysis of the molecular arrangement within crystal structure reveals that the molecular motion becomes more and more localized, and hints towards the evolution of intermolecular interaction with pressure. Moreover, the electronic structure of naphthalene under high pressure is also discussed. A pressure induced decrease of the band gap is observed.
Suleiman, Mohammed S. H.; Joubert, Daniel P.
2015-11-01
In the present work, the atomic and the electronic structures of Au3N, AuN and AuN2 are investigated using first-principles density-functional theory (DFT). We studied cohesive energy vs. volume data for a wide range of possible structures of these nitrides. Obtained data were fitted to a Birch-Murnaghan third-order equation of state (EOS) so as to identify the most likely candidates for the true crystal structure in this subset of the infinite parameter space, and to determine their equilibrium structural parameters. The analysis of the electronic properties was achieved by the calculations of the band structure and the total and partial density of states (DOS). Some possible pressure-induced structural phase transitions have been pointed out. Further, we carried out GW0 calculations within the random-phase approximation (RPA) to the dielectric tensor to investigate the optical spectra of the experimentally suggested modification: Au3N(D09). Obtained results are compared with experiment and with some available previous calculations.
Minimal parameter implicit solvent model for ab initio electronic structure calculations
Dziedzic, Jacek; Skylaris, Chris-Kriton; Mostofi, Arash A; Payne, Mike C
2011-01-01
We present an implicit solvent model for ab initio electronic structure calculations which is fully self-consistent and is based on direct solution of the nonhomogeneous Poisson equation. The solute cavity is naturally defined in terms of an isosurface of the electronic density according to the formula of Fattebert and Gygi (J. Comp. Chem. 23, 6 (2002)). While this model depends on only two parameters, we demonstrate that by using appropriate boundary conditions and dispersion-repulsion contributions, solvation energies obtained for an extensive test set including neutral and charged molecules show dramatic improvement compared to existing models. Our approach is implemented in, but not restricted to, a linear-scaling density functional theory (DFT) framework, opening the path for self-consistent implicit solvent DFT calculations on systems of unprecedented size, which we demonstrate with calculations on a 2615-atom protein-ligand complex.
Density functional calculation of the structural and electronic properties of germanium quantum dots
International Nuclear Information System (INIS)
We apply first principles density functional computational methods to study the structures, densities of states (DOS), and higher occupied molecular orbital (HOMO) – lowest unoccupied molecular orbital (LUMO) gaps of selected free-standing Ge semiconductor quantum dots up to 1.8nm. Our calculations are performed using numerical atomic orbital approach where linear combination of atomic orbital was applied. The surfaces of the quantum dots was passivized by hydrogen atoms. We find that surface passivation does affect the electronic properties associated with the changes of surface state, electron localization, and the energy gaps of germanium nanocrystals as well as the confinement of electrons inside the quantum dots (QDs). Our study shows that the energy gaps of germanium quantum dots decreases with the increasing dot diameter. The size-dependent variations of the computed HOMO-LUMO gaps in our quantum dots model were found to be consistent with the effects of quantum confinement reported in others theoretical and experimental calculation
The structural and electronic properties of amorphous HgCdTe from first-principles calculations
International Nuclear Information System (INIS)
Amorphous mercury cadmium telluride (a-MCT) model structures, with x being 0.125 and 0.25, are obtained from first-principles calculations. We generate initial structures by computation alchemy method. It is found that most atoms in the network of amorphous structures tend to be fourfold and form tetrahedral structures, implying that the chemical ordered continuous random network with some coordination defects is the ideal structure for a-MCT. The electronic structure is also concerned. The gap is found to be 0.30 and 0.26 eV for a-Hg0.875Cd0.125Te and a-Hg0.75Cd0.25Te model structures, independent of the composition. By comparing with the properties of crystalline MCT with the same composition, we observe a blue-shift of energy band gap. The localization of tail states and its atomic origin are also discussed. (paper)
Density functional calculation of equilibrium geometry and electronic structure of pyrite
Institute of Scientific and Technical Information of China (English)
邱冠周; 肖奇; 胡岳华; 徐竞
2001-01-01
The equilibrium geometry and electronic structure of pyrite has been studied using self-consistent density-functional theory within the local density approximation (LDA). The optimum bulk geometry is in good agreement with crystallographic data. The calculated band structure and density of states in the region around the Fermi energy show that valence-band maximum (VBM) is at X (100), and the conduction-band minimum (CBM) is at G (000). The indirect and direct band gaps are 0.6eV and 0.74eV, respectively. The calculated contour map of difference of charge density shows excess charge in nonbonding d electron states on the Fe sites. The density increases between sulfur nuclei and between iron and sulfur nuclei qualitatively reveal that S-S bond and Fe-S bond are covalent binding.
An approach to first principles electronic structure calculation by symbolic-numeric computation
Kikuchi, Akihito
2013-01-01
This article is an introduction to a new approach to first principles electronic structure calculation. The starting point is the Hartree-Fock-Roothaan equation, in which molecular integrals are approximated by polynomials by way of Taylor expansion with respect to atomic coordinates and other variables. It leads to a set of polynomial equations whose solutions are eigenstate, which is designated as algebraic molecular orbital equation. Symbolic computation, especially, Gr\\"obner bases theory, enables us to rewrite the polynomial equations into more trimmed and tractable forms with identical roots, from which we can unravel the relationship between physical parameters (wave function, atomic coordinates, and others) and numerically evaluate them one by one in order. Furthermore, this method is a unified way to solve the electronic structure calculation, the optimization of physical parameters, and the inverse problem as a forward problem.
A proposal to first principles electronic structure calculation: Symbolic-Numeric method
Kikuchi, Akihito
2012-01-01
This study proposes an approach toward the first principles electronic structure calculation with the aid of symbolic-numeric solving. The symbolic computation enables us to express the Hartree-Fock-Roothaan equation in an analytic form and approximate it as a set of polynomial equations. By use of the Grobner basis technique, the polynomial equations are transformed into other ones which have identical roots. The converted equations take more convenient forms which will simplify numerical procedures, from which we can derive necessary physical properties in order, in an a la carte way. This method enables us to solve the electronic structure calculation, the optimization of any kind, or the inverse problem as a forward problem in a unified way, in which there is no need for iterative self-consistent procedures with trials and errors.
Electronic structure of KFe2Se2 from first-principles calculations
International Nuclear Information System (INIS)
Electronic structures and magnetic properties for iron-selenide KFe2Se2 are studied by first-principles calculations. The ground state is collinear antiferromagnetic with calculated 2.26μB magnetic moment on Fe atoms; and the J1 and J2 coupling strengths are calculated to be 0.038 eV and 0.029 eV. The states around EF are dominated by the Fe 3d orbitals which hybridize noticeably to the Se 4p orbitals. While the band structure of KFe2Se2 is similar to a heavily electron-doped BaFe2AS2 or FeSe system, the Fermi surface of KFe2Se2 is much closer to the FeSe system since the electron sheets around M are symmetric with respect to x-y exchange. These features, as well as the absence of Fermi surface nesting, suggest that the parent KFe2Se2 could be regarded as an electron doped FeSe system with possible local moment magnetism. (authors)
Multi-Center Electronic Structure Calculations for Plasma Equation of State
Energy Technology Data Exchange (ETDEWEB)
Wilson, B G; Johnson, D D; Alam, A
2010-12-14
We report on an approach for computing electronic structure utilizing solid-state multi-center scattering techniques, but generalized to finite temperatures to model plasmas. This approach has the advantage of handling mixtures at a fundamental level without the imposition of ad hoc continuum lowering models, and incorporates bonding and charge exchange, as well as multi-center effects in the calculation of the continuum density of states.
Hao, Yajiang; Inhester, Ludger; Hanasaki, Kota; Son, Sang-Kil; Santra, Robin
2015-07-01
We present the implementation of an electronic-structure approach dedicated to ionization dynamics of molecules interacting with x-ray free-electron laser (XFEL) pulses. In our scheme, molecular orbitals for molecular core-hole states are represented by linear combination of numerical atomic orbitals that are solutions of corresponding atomic core-hole states. We demonstrate that our scheme efficiently calculates all possible multiple-hole configurations of molecules formed during XFEL pulses. The present method is suitable to investigate x-ray multiphoton multiple ionization dynamics and accompanying nuclear dynamics, providing essential information on the chemical dynamics relevant for high-intensity x-ray imaging. PMID:26798806
Efficient electronic structure calculation for molecular ionization dynamics at high x-ray intensity
Hao, Yajiang; Hanasaki, Kota; Son, Sang-Kil; Santra, Robin
2015-01-01
We present the implementation of an electronic-structure approach dedicated to ionization dynamics of molecules interacting with x-ray free-electron laser (XFEL) pulses. In our scheme, molecular orbitals for molecular core-hole states are represented by linear combination of numerical atomic orbitals that are solutions of corresponding atomic core-hole states. We demonstrate that our scheme efficiently calculates all possible multiple-hole configurations of molecules formed during XFEL pulses. The present method is suitable to investigate x-ray multiphoton multiple ionization dynamics and accompanying nuclear dynamics, providing essential information on the chemical dynamics relevant for high-intensity x-ray imaging.
Efficient electronic structure calculation for molecular ionization dynamics at high x-ray intensity
Directory of Open Access Journals (Sweden)
Yajiang Hao
2015-07-01
Full Text Available We present the implementation of an electronic-structure approach dedicated to ionization dynamics of molecules interacting with x-ray free-electron laser (XFEL pulses. In our scheme, molecular orbitals for molecular core-hole states are represented by linear combination of numerical atomic orbitals that are solutions of corresponding atomic core-hole states. We demonstrate that our scheme efficiently calculates all possible multiple-hole configurations of molecules formed during XFEL pulses. The present method is suitable to investigate x-ray multiphoton multiple ionization dynamics and accompanying nuclear dynamics, providing essential information on the chemical dynamics relevant for high-intensity x-ray imaging.
Energy Technology Data Exchange (ETDEWEB)
Wills, John M [Los Alamos National Laboratory; Mattsson, Ann E [Sandia National Laboratories
2012-06-06
Brooks, Johansson, and Skriver, using the LMTO-ASA method and considerable insight, were able to explain many of the ground state properties of the actinides. In the many years since this work was done, electronic structure calculations of increasing sophistication have been applied to actinide elements and compounds, attempting to quantify the applicability of DFT to actinides and actinide compounds and to try to incorporate other methodologies (i.e. DMFT) into DFT calculations. Through these calculations, the limits of both available density functionals and ad hoc methodologies are starting to become clear. However, it has also become clear that approximations used to incorporate relativity are not adequate to provide rigorous tests of the underlying equations of DFT, not to mention ad hoc additions. In this talk, we describe the result of full-potential LMTO calculations for the elemental actinides, comparing results obtained with a full Dirac basis with those obtained from scalar-relativistic bases, with and without variational spin-orbit. This comparison shows that the scalar relativistic treatment of actinides does not have sufficient accuracy to provide a rigorous test of theory and that variational spin-orbit introduces uncontrolled errors in the results of electronic structure calculations on actinide elements.
A novel Gaussian-Sinc mixed basis set for electronic structure calculations
International Nuclear Information System (INIS)
A Gaussian-Sinc basis set methodology is presented for the calculation of the electronic structure of atoms and molecules at the Hartree–Fock level of theory. This methodology has several advantages over previous methods. The all-electron electronic structure in a Gaussian-Sinc mixed basis spans both the “localized” and “delocalized” regions. A basis set for each region is combined to make a new basis methodology—a lattice of orthonormal sinc functions is used to represent the “delocalized” regions and the atom-centered Gaussian functions are used to represent the “localized” regions to any desired accuracy. For this mixed basis, all the Coulomb integrals are definable and can be computed in a dimensional separated methodology. Additionally, the Sinc basis is translationally invariant, which allows for the Coulomb singularity to be placed anywhere including on lattice sites. Finally, boundary conditions are always satisfied with this basis. To demonstrate the utility of this method, we calculated the ground state Hartree–Fock energies for atoms up to neon, the diatomic systems H2, O2, and N2, and the multi-atom system benzene. Together, it is shown that the Gaussian-Sinc mixed basis set is a flexible and accurate method for solving the electronic structure of atomic and molecular species
International Nuclear Information System (INIS)
We investigate the cohesive energy, heat of formation, elastic constant and electronic band structure of transition metal diborides TMB2 (TM = Hf, Ta, W, Re, Os and Ir, Pt) in the Pmmn space group using the ab initio pseudopotential total energy method. Our calculations indicate that there is a relationship between elastic constant and valence electron concentration (VEC): the bulk modulus and shear modulus achieve their maximum when the VEC is in the range of 6.8-7.2. In addition, trends in the elastic constant are well explained in terms of electronic band structure analysis, e.g., occupation of valence electrons in states near the Fermi level, which determines the cohesive energy and elastic properties. The maximum in bulk modulus and shear modulus is attributed to the nearly complete filling of TM d-B p bonding states without filling the antibonding states. On the basis of the observed relationship, we predict that alloying W and Re in the orthorhombic structure OsB2 might be harder than alloying the Ir element. Indeed, the further calculations confirmed this expectation
International Nuclear Information System (INIS)
Electronic-structure calculations of elemental praseodymium are presented. Several approximations are used to describe the Pr f electrons. It is found that the low-pressure, trivalent phase is well described using either the self-interaction corrected (SIC) local-spin-density (LSD) approximation or the generalized-gradient approximation (GGA) with spin and orbital polarization (OP). In the SIC-LSD approach the Pr f electrons are treated explicitly as localized with a localization energy given by the self-interaction of the f orbital. In the GGA+OP scheme the f-electron localization is described by the onset of spin and orbital polarization, the energetics of which is described by spin-moment formation energy and a term proportional to the total orbital moment, Lz2. The high-pressure phase is well described with the f electrons treated as band electrons, in either the LSD or the GGA approximations, of which the latter describes more accurately the experimental equation of state. The calculated pressure of the transition from localized to delocalized behavior is 280 kbar in the SIC-LSD approximation and 156 kbar in the GGA+OP approach, both comparing favorably with the experimentally observed transition pressure of 210 kbar. copyright 1997 The American Physical Society
International Nuclear Information System (INIS)
The defect states and optical absorption enhancement induced by twin boundaries in silicon are investigated by first-principle calculation. The defect states in the forbidden bands are identified and based on the established electronic structures, the dielectric functions and absorption coefficients are derived. An important result of our calculations is that visible light absorption by the twinning configuration is enhanced significantly, indicating that twinning structures possibly play an important role in silicon-based photovoltaic devices. - Highlights: • Defect states and optical absorption enhancement induced by twin boundaries in silicon are investigated theoretically. • Dielectric functions and absorption coefficients are derived. • Enhanced visible light absorption by the twinning configuration is demonstrated. • Twinning structures play an important role in silicon-based photovoltaic devices
Energy Technology Data Exchange (ETDEWEB)
Liu, X.X.; Liu, L.Z. [Key Laboratory of Modern Acoustics, MOE, Institute of Acoustics and Collaborative Innovation Center of Advanced Microstructures, National Laboratory of Solid State Microstructures, Nanjing University, Nanjing 210093 (China); Wu, X.L., E-mail: hkxlwu@nju.edu.cn [Key Laboratory of Modern Acoustics, MOE, Institute of Acoustics and Collaborative Innovation Center of Advanced Microstructures, National Laboratory of Solid State Microstructures, Nanjing University, Nanjing 210093 (China); Department of Physics, NingBo University, NingBo 315301 (China); Chu, Paul K. [Department of Physics and Materials Science, City University of Hong Kong, Tat Chee Avenue, Kowloon, Hong Kong (China)
2015-07-03
The defect states and optical absorption enhancement induced by twin boundaries in silicon are investigated by first-principle calculation. The defect states in the forbidden bands are identified and based on the established electronic structures, the dielectric functions and absorption coefficients are derived. An important result of our calculations is that visible light absorption by the twinning configuration is enhanced significantly, indicating that twinning structures possibly play an important role in silicon-based photovoltaic devices. - Highlights: • Defect states and optical absorption enhancement induced by twin boundaries in silicon are investigated theoretically. • Dielectric functions and absorption coefficients are derived. • Enhanced visible light absorption by the twinning configuration is demonstrated. • Twinning structures play an important role in silicon-based photovoltaic devices.
Structural, electronic, and optical properties of NiAl3: First-principles calculations
Saniz, R; Ye, Lin-Hui; Shishidou, Tatsuya; Freeman, A. J.
2006-01-01
We report ab initio density-functional calculations of the structural, electronic, and optical properties of NiAl3, using the full-potential linearized augmented plane wave method within the generalized gradient approximation to the exchange-correlation potential. The D011 structure is found to be energetically favorable over both the cubic L12 and A15 phases. The density of states around the Fermi energy, including a pseudogap just above it, is dominated by strongly hybridized Ni d and Al p ...
Landau, Arie; Haritan, Idan; Kaprálová-Žd'ánská, Petra Ruth; Moiseyev, Nimrod
2016-05-19
Complex eigenvalues, resonances, play an important role in a large variety of fields in physics and chemistry. For example, in cold molecular collision experiments and electron scattering experiments, autoionizing and predissociative metastable resonances are generated. However, the computation of complex resonance requires modifications of standard electronic structure codes and methods, which are not always straightforward, in addition, application of complex codes requires more computational efforts. Here we show how resonance eigenvalues, positions and widths, can be calculated using the standard, widely used, electronic-structure packages. Our method enables the calculations of the complex resonance eigenvalues by using analytical continuation procedures (such as Padé). The key point in our approach is the existence of narrow analytical passages from the real axis to the complex energy plane. In fact, the existence of these analytical passages relies on using finite basis sets. These passages become narrower as the basis set becomes more complete, whereas in the exact limit, these passages to the complex plane are closed. As illustrative numerical examples we calculated the autoionization Feshbach resonances of helium, hydrogen anion, and hydrogen molecule. We show that our results are in an excellent agreement with the results obtained by other theoretical methods and with available experimental results. PMID:26677725
Electronics Environmental Benefits Calculator
U.S. Environmental Protection Agency — The Electronics Environmental Benefits Calculator (EEBC) was developed to assist organizations in estimating the environmental benefits of greening their purchase,...
Directory of Open Access Journals (Sweden)
Yu Wang
2002-01-01
Full Text Available Abstract:We investigate a theoretical model of molecular metalwire constructed from linear polynuclear metal complexes. In particular we study the linear Crn metal complex and Cr molecular metalwire. The electron density distributions of the model nanowire and the linear Crn metal complexes, with n = 3, 5, and 7, are calculated by employing CRYSTAL98 package with topological analysis. The preliminary results indicate that the bonding types between any two neighboring Cr are all the same, namely the polarized open-shell interaction. The pattern of electron density distribution in metal complexes resembles that of the model Cr nanowire as the number of metal ions increases. The conductivity of the model Cr nanowire is also tested by performing the band structure calculation.
Ab initio calculations on twisted graphene/hBN: Electronic structure and STM image simulation
Correa, J. D.; Cisternas, E.
2016-09-01
By performing ab initio calculations we obtained theoretical scanning tunneling microscopy (STM) images and studied the electronic properties of graphene on a hexagonal boron-nitrite (hBN) layer. Three different stack configurations and four twisted angles were considered. All calculations were performed using density functional theory, including van der Waals interactions as implemented in the SIESTA ab initio package. Our results show that the electronic structure of graphene is preserved, although some small changes are induced by the interaction with the hBN layer, particularly in the total density of states at 1.5 eV under the Fermi level. When layers present a twisted angle, the density of states shows several van Hove singularities under the Fermi level, which are associated to moiré patterns observed in theoretical STM images.
Electronic structure calculations of europium chalcogenides EuS and EuSe
Energy Technology Data Exchange (ETDEWEB)
Rached, D.; Ameri, M.; Rabah, M.; Benkhettou, N.; Dine el Hannani, M. [Laboratoire des Materiaux Appliques, Centre de Recherche, Route de Mascara, Universite de Sidi-Bel-Abbes, Sidi Bel Abbes 22000 (Algeria); Khenata, R. [Laboratoire des Materiaux Appliques, Centre de Recherche, Route de Mascara, Universite de Sidi-Bel-Abbes, Sidi Bel Abbes 22000 (Algeria); Laboratoire de Physique Quantique et de Modelisation Mathematique de la Matiere (LPQ3M), Centre Universitaire de Mascara, Mascara 29000 (Algeria); Bouhemadou, A. [Departement de Physique, Faculte des Sciences, Universite Ferhat Abbes, 19000 Setif (Algeria)
2007-06-15
We have performed ab-initio self-consistent calculations on the full-potential linear muffin-tin orbital method with the local-density approximation and local spin-density approximation to investigate the structural and electronic properties of EuS and EuSe in its stable (NaCl-B1) and high-pressure phases. The magnetic phase stability was determined from the total energy calculations for both the nonmagnetic (NM) and magnetic (M) phases. These theoretical calculations clearly indicate that both at ambient and high pressures, the magnetic phase is more stable than the nonmagnetic phase. The transition pressure at which these compounds undergo the structural phase transition from NaCl-B1 to CsCl-B2 phase is calculated. The elastic constants at equilibrium in both NaCl-B1 and CsCl-B2 structures are also determined. (copyright 2007 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim) (orig.)
Electronic structure calculations of europium chalcogenides EuS and EuSe
International Nuclear Information System (INIS)
We have performed ab-initio self-consistent calculations on the full-potential linear muffin-tin orbital method with the local-density approximation and local spin-density approximation to investigate the structural and electronic properties of EuS and EuSe in its stable (NaCl-B1) and high-pressure phases. The magnetic phase stability was determined from the total energy calculations for both the nonmagnetic (NM) and magnetic (M) phases. These theoretical calculations clearly indicate that both at ambient and high pressures, the magnetic phase is more stable than the nonmagnetic phase. The transition pressure at which these compounds undergo the structural phase transition from NaCl-B1 to CsCl-B2 phase is calculated. The elastic constants at equilibrium in both NaCl-B1 and CsCl-B2 structures are also determined. (copyright 2007 WILEY-VCH Verlag GmbH and Co. KGaA, Weinheim) (orig.)
Theoretical calculations on the atomic and electronic structure of β-SiC(110) surface
Institute of Scientific and Technical Information of China (English)
无
2002-01-01
We present a theoretical calculation of the atomic and electronic structure of β-SiC and its non-polar (110) surface using the full potential linear augmented plane wave (FPLAPW) approach. The calculated lattice constant and bulk modulus of β-SiC crystal are in excellent agreement with experimental data. The atomic and electronic structure of β-SiC(110) surface has been calculated by employing the slab and supercell model. It is found that the surface is characterized by a top-layer bond-length-contracting rotation relaxation in which the Si-surface atom moves closer towards the substrate while the C-surface atom moves outward. This relaxation is analogous to that of Ⅲ-Ⅴ semi-conductor surface. The driving mechanism for this atomic rearrangement is that the Si atom tends to a planar sp2-like bonding situation with its three N neighbors and the N atom tends to a p3-like bonding with its three Si neighbors. Furthermore, surface relaxation induces the change from metallic to semiconducting characterization.
Genovese, Luigi; Deutsch, Thierry
2015-12-21
Discretizing an analytic function on a uniform real-space grid is often done via a straightforward collocation method. This is ubiquitous in all areas of computational physics and quantum chemistry. An example in density functional theory (DFT) is given by the external potential or the pseudo-potential describing the interaction between ions and electrons. The accuracy of the collocation method used is therefore very important for the reliability of subsequent treatments like self-consistent field solutions of the electronic structure problems. By construction, the collocation method introduces numerical artifacts typical of real-space treatments, like the so-called egg-box error, which may spoil the numerical stability of the description when the real-space grid is too coarse. As the external potential is an input of the problem, even a highly precise computational treatment cannot cope this inconvenience. We present in this paper a new quadrature scheme that is able to exactly preserve the moments of a given analytic function even for large grid spacings, while reconciling with the traditional collocation method when the grid spacing is small enough. In the context of real-space electronic structure calculations, we show that this method improves considerably the stability of the results for large grid spacings, opening up the path towards reliable low-accuracy DFT calculations with a reduced number of degrees of freedom. PMID:26372293
Analytic methods for the calculation of the electronic structure of solids
International Nuclear Information System (INIS)
Andersen's atomic-sphere approximation has been utilized with approximations based upon linear-combination of atomic orbitals (LCAO) theory to obtain approximate energy-band parameters for solids. Simple analytic expressions for the bandwidth and position of the band center have been derived that require only free-atom wave functions evaluated at the Wigner-Seitz atomic-sphere radius. For convenience, the method has been named the atomic surface method (ASM). The following simple analytic expressions for the band parameters have been derived from the ASM: (i) The bandwidth is equal to the product of h2/m, the gradient of the electron density at the atomic-sphere radius, and the surface area of the sphere; (ii) the average band energy is shifted from the atomic-term-value energy by an amount given by the product of the bandwidth, electron density at the atomic-sphere radius, and atomic-sphere volume. The theory has been applied without adjustable parameters to the transition metals and f-shell metals with use of tabulated Hartree-Fock wave functions and is in reasonable agreement with full band-structure calculations. The same analysis is applied to atomic core states under compression and is also in reasonable agreement with complete band-structure calculations. The 2s and 2p states of Na and Al have been calculated to the point where they merge with the conduction band as free-electron states. These bandwidths and shifts are also written in terms of the atomic term values by using the asymptotic form of the radial wave function. Finally, the LCAO energy bands of Ni are calculated with use of the ASM parameters
Real-space ab-initio electronic structure calculations using SfePy
Czech Academy of Sciences Publication Activity Database
Cimrman, R.; Novák, Matyáš; Kolman, Radek; Vackář, Jiří
Plzeň: University of West Bohemia, 2015 - (Adámek, V.). s. 21-22 ISBN 978-80-261-0568-8. [Computational Mechanics 2015 /31./- conference with international participation /31./. 09.11.2015-11.11.2015, Špičák] R&D Projects: GA ČR(CZ) GAP108/11/0853; GA ČR(CZ) GAP101/12/2315 Institutional support: RVO:61388998 ; RVO:68378271 Keywords : real-space ab-initio electronic structure calculations * finite element method * isogeometric analysis Subject RIV: BE - Theoretical Physics
Tatiya Chokbunpiam; Patchanita Thamyongkit; Oraphan Saengsawang; Supot Hannongbua
2010-01-01
This study aimed to design a new series of compounds consisting of a porphyrin macrocycle linked to a perylene unit via a thiophenic bridge. The structural and electronic properties of the molecules, and the effects of mono- and di-substituents R on C3 and R′ on C4 of the thiophene ring were investigated using a quantum calculation approach. The results from the method validation revealed that using the density functional theory approach at B3LYP/6–31G(d) data set was the optimal one, conside...
Cimrman, Robert; Kolman, Radek; Tůma, Miroslav; Vackář, Jiří
2015-01-01
We compare convergence of isogeometric analysis (IGA), a spline modification of finite element method (FEM), with FEM in the context of our real space code for ab-initio electronic structure calculations of non-periodic systems. The convergence is studied on simple sub-problems that appear within the density functional theory approximation to the Schr\\"odinger equation: the Poisson problem and the generalized eigenvalue problem. We also outline the complete iterative algorithm seeking a fixed point of the charge density of a system of atoms or molecules, and study IGA/FEM convergence on a benchmark problem of nitrogen atom.
Calculated electronic and magnetic structure of screw dislocations in alpha iron
Energy Technology Data Exchange (ETDEWEB)
Odbadrakh, K.; Rusanu, A.; Stocks, G. Malcolm; Samolyuk, G. D.; Eisenbach, M.; Wang, Yang; Nicholson, D. M.
2011-01-01
Local atomic magnetic moments in crystalline Fe are perturbed by the presence of dislocations. The effects are most pronounced near the dislocation core and decay slowly as the strain field of the dislocation decreases with distance. We have calculated local moments using the locally self-consistent multiple scattering (LSMS) method for a supercell containing a screw-dislocation quadrupole. Finite size effects are found to be significant indicating that dislocation cores affect the electronic structure and magnetic moments of neighboring dislocations. The influence of neighboring dislocations points to a need to study individual dislocations from first principles just as they appear amid surrounding atoms in large-scale classical force field simulations. An approach for the use of the LSMS to calculate local moments in subvolumes of large atomic configurations generated in the course of classical molecular dynamics simulation of dislocationdynamics is discussed.
Calculated electronic and magnetic structure of screw dislocations in alpha iron
Energy Technology Data Exchange (ETDEWEB)
Odbadrakh, Khorgolkhuu [ORNL; Rusanu, Aurelian [ORNL; Stocks, George Malcolm [ORNL; Samolyuk, German D [ORNL; Eisenbach, Markus [ORNL; Wang, Yang Nmn [ORNL; Nicholson, Don M [ORNL
2011-01-01
Local atomic magnetic moments in crystalline Fe are perturbed by the presence of dislocations. The effects are most pronounced near the dislocation core and decay slowly as the strain field of the dislocation decreases with distance. We have calculated local moments using the locally self-consistent multiple scattering (LSMS) method for a supercell containing a screw-dislocation quadrupole. Finite size effects are found to be significant indicating that dislocation cores affect the electronic structure and magnetic moments of neighboring dislocations. The influence of neighboring dislocations points to a need to study individual dislocations from first principles just as they appear amid surrounding atoms in large-scale classical force field simulations. An approach for the use of the LSMS to calculate local moments in subvolumes of large atomic configurations generated in the course of classical molecular dynamics simulation of dislocation dynamics is discussed. VC2011 American Institute of Physics. [doi:10.1063/1.3562217
Cao, Jun; Xie, Zhi-Zhong; Yu, Xiaodong
2016-08-01
In the present work, the combined electronic structure calculations and surface hopping simulations have been performed to investigate the excited-state decay of the parent oxazole in the gas phase. Our calculations show that the S2 state decay of oxazole is an ultrafast process characterized by the ring-opening and ring-closure of the five-membered oxazole ring, in which the triplet contribution is minor. The ring-opening involves the Osbnd C bond cleavage affording the nitrile ylide and airine intermediates, while the ring-closure gives rise to a bicyclic species through a 2sbnd 5 bond formation. The azirine and bicyclic intermediates in the S0 state are very likely involved in the phototranspositions of oxazoles. This is different from the previous mechanism in which these intermediates in the T1 state have been proposed for these phototranspositions.
Electronic structure of alkali metal hydrides on data of cluster calculations by LCAO MO SCF CNDO
International Nuclear Information System (INIS)
The results of quantum-chemical study in where by M = Li, Na, K, Rb and Cs, are presented. The calculation expresses the expected electron density distributions in hydrides on the hydrogen and metal atoms as well as the energy characteristics: M-H, M-M and compounds binding energies. The latter ones qualitatively correlate with the binding energies of LiH-CsH compounds. The calculated values for the Fermi energy and the width of the forbidden zone at the Fermi level make it possible to suppose that the ideally formed lithium hydride crystal will be characterized by the highest electrical resistance. It is established that quantum-chemical characteristics of the MH hydrides structure change nonmonotonously by transfer from Li to Cs
Polfus, Jonathan M; Bjørheim, Tor S; Norby, Truls; Haugsrud, Reidar
2012-09-01
The nitrogen related defect chemistry and electronic structure of wide band gap oxides are investigated by density functional theory defect calculations of N(O)(q), NH(O)(×), and (NH2)(O)(·) as well as V(O)(··) and OH(O)(·) in MgO, CaO, SrO, Al(2)O(3), In(2)O(3), Sc(2)O(3), Y(2)O(3), La(2)O(3), TiO(2), SnO(2), ZrO(2), BaZrO(3), and SrZrO(3). The N(O)(q) acceptor level is found to be deep and the binding energy of NH(O)(×) with respect to N(O)' and (OH(O)(·) is found to be significantly negative, i.e. binding, in all of the investigated oxides. The defect structure of the oxides was found to be remarkably similar under reducing and nitriding conditions (1 bar N(2), 1 bar H(2) and 1 × 10(-7) bar H(2)O): NH(O)(×) predominates at low temperatures and [N(O)'] = 2[V(O)(··) predominates at higher temperatures (>900 K for most of the oxides). Furthermore, we evaluate how the defect structure is affected by non-equilibrium conditions such as doping and quenching. In terms of electronic structure, N(O)' is found to introduce isolated N-2p states within the band gap, while the N-2p states of NH(O)(×) are shifted towards, or overlap with the VBM. Finally, we assess the effect of nitrogen incorporation on the proton conducting properties of oxides and comment on their corrosion resistance in nitriding atmospheres in light of the calculated defect structures. PMID:22828729
Song, Xiaowei; Fagiani, Matias R.; Gewinner, Sandy; Schöllkopf, Wieland; Asmis, Knut R.; Bischoff, Florian A.; Berger, Fabian; Sauer, Joachim
2016-06-01
We use cryogenic ion trap vibrational spectroscopy in combination with quantum chemical calculations to study the structure of mono- and dialuminum oxide anions. The infrared photodissociation spectra of D2-tagged AlO1-4- and Al2O3-6- are measured in the region from 400 to 1200 cm-1. Structures are assigned based on a comparison to simulated harmonic and anharmonic IR spectra derived from electronic structure calculations. The monoaluminum anions contain an even number of electrons and exhibit an electronic closed-shell ground state. The Al2O3-6- anions are oxygen-centered radicals. As a result of a delicate balance between localization and delocalization of the unpaired electron, only the BHLYP functional is able to qualitatively describe the observed IR spectra of all species with the exception of AlO3-. Terminal Al-O stretching modes are found between 1140 and 960 cm-1. Superoxo and peroxo stretching modes are found at higher (1120-1010 cm-1) and lower energies (850-570 cm-1), respectively. Four modes in-between 910 and 530 cm-1 represent the IR fingerprint of the common structural motif of dialuminum oxide anions, an asymmetric four-member Al-(O)2-Al ring.
Energy Technology Data Exchange (ETDEWEB)
Singh, David J [ORNL; Safa-Sefat, Athena [ORNL; McGuire, Michael A [ORNL; Sales, Brian C [ORNL; Mandrus, David [ORNL; VanBebber, L. H. [University of Tennessee, Knoxville (UTK); Keppens, Veerle [University of Tennessee, Knoxville (UTK)
2009-01-01
We report single crystal synthesis, specific heat and resistivity measurements and electronic structure calculations for BaCr2As2. This material is a metal with itinerant antiferromagnetism, similar to the parent phases of Fe-based high temperature superconductors, but differs in magnetic order. Comparison of bare band structure density of states and the low temperature specific heat implies a mass renormalization of 2. BaCr2As2 shows stronger transition metal - pnictogen covalency than the Fe compounds, and in this respect is more similar to BaMn2As2. This provides an explanation for the observation that Ni and Co doping is effective in the Fe-based superconductors, but Cr or Mn doping is not.
Seiler, Christian
2016-01-01
A formalism for electronic-structure calculations is presented that is based on the functional renormalization group (FRG). The traditional FRG has been formulated for systems that exhibit a translational symmetry with an associated Fermi surface, which can provide the organization principle for the renormalization group (RG) procedure. We here advance an alternative formulation, where the RG-flow is organized in the energy-domain rather than in k-space. This has the advantage that it can also be applied to inhomogeneous matter lacking a band-structure, such as disordered metals or molecules. The energy-domain FRG ({\\epsilon}FRG) presented here accounts for Fermi-liquid corrections to quasi-particle energies and particle-hole excitations. It goes beyond the state of the art GW-BSE, because in {\\epsilon}FRG the Bethe-Salpeter equation (BSE) is solved in a self-consistent manner. An efficient implementation of the approach that has been tested against exact diagonalization calculations and calculations based on...
Time domain numerical calculations of the short electron bunch wakefields in resistive structures
Energy Technology Data Exchange (ETDEWEB)
Tsakanian, Andranik
2010-10-15
The acceleration of electron bunches with very small longitudinal and transverse phase space volume is one of the most actual challenges for the future International Linear Collider and high brightness X-Ray Free Electron Lasers. The exact knowledge on the wake fields generated by the ultra-short electron bunches during its interaction with surrounding structures is a very important issue to prevent the beam quality degradation and to optimize the facility performance. The high accuracy time domain numerical calculations play the decisive role in correct evaluation of the wake fields in advanced accelerators. The thesis is devoted to the development of a new longitudinally dispersion-free 3D hybrid numerical scheme in time domain for wake field calculation of ultra short bunches in structures with walls of finite conductivity. The basic approaches used in the thesis to solve the problem are the following. For materials with high but finite conductivity the model of the plane wave reflection from a conducting half-space is used. It is shown that in the conductive half-space the field components perpendicular to the interface can be neglected. The electric tangential component on the surface contributes to the tangential magnetic field in the lossless area just before the boundary layer. For high conducting media, the task is reduced to 1D electromagnetic problem in metal and the so-called 1D conducting line model can be applied instead of a full 3D space description. Further, a TE/TM (''transverse electric - transverse magnetic'') splitting implicit numerical scheme along with 1D conducting line model is applied to develop a new longitudinally dispersion-free hybrid numerical scheme in the time domain. The stability of the new hybrid numerical scheme in vacuum, conductor and bound cell is studied. The convergence of the new scheme is analyzed by comparison with the well-known analytical solutions. The wakefield calculations for a number of
International Nuclear Information System (INIS)
Uranium dioxide UO2 is the standard nuclear fuel used in pressurized water reactors. During in-reactor operation, the fission of uranium atoms yields a wide variety of fission products (FP) which create numerous point defects while slowing down in the material. Point defects and FP govern in turn the evolution of the fuel physical properties under irradiation. In this study, we use electronic structure calculations in order to better understand the fuel behavior under irradiation. In particular, we investigate point defect behavior, as well as the stability of three volatile FP: iodine, krypton and xenon. In order to take into account the strong correlations of uranium 5f electrons in UO2, we use the DFT+U approximation, based on the density functional theory. This approximation, however, creates numerous metastable states which trap the system and induce discrepancies in the results reported in the literature. To solve this issue and to ensure the ground state is systematically approached as much as possible, we use a method based on electronic occupancy control of the correlated orbitals. We show that the DFT+U approximation, when used with electronic occupancy control, can describe accurately point defect and fission product behavior in UO2 and provide quantitative information regarding point defect transport properties in the oxide fuel. (author)
International Nuclear Information System (INIS)
It was recently shown that the energy resolution of Ce-doped LaBr3 scintillator radiation detectors can be crucially improved by co-doping with Sr, Ca, or Ba. Here, we outline a mechanism for this enhancement on the basis of electronic structure calculations. We show that (i) Br vacancies are the primary electron traps during the initial stage of thermalization of hot carriers, prior to hole capture by Ce dopants; (ii) isolated Br vacancies are associated with deep levels; (iii) Sr doping increases the Br vacancy concentration by several orders of magnitude; (iv) SrLa binds to VBr resulting in a stable neutral complex; and (v) association with Sr causes the deep vacancy level to move toward the conduction band edge. The latter is essential for reducing the effective carrier density available for Auger quenching during thermalization of hot carriers. Subsequent de-trapping of electrons from SrLa–VBr complexes can activate Ce dopants that have previously captured a hole leading to luminescence. This mechanism implies an overall reduction of Auger quenching of free carriers, which is expected to improve the linearity of the photon light yield with respect to the energy of incident electron or photon
International Nuclear Information System (INIS)
An accurate and efficient method is described for the evaluation of electrostatic contributions in LCAO electronic structure calculations. The charge density rho(r) is decomposed into rho/sup(1)(r), a component whose rapid variation near any nucleus reproduces that of rho(r) to a very good approximation, and a remainder density deltarho(r)equivalentrho(r)-rho/sup(1)(r), which is thereby guaranteed to be slowly varying in space. The power of the decomposition resides in the fact that rho/sup(1)(r) can be expressed exactly as a sum of one-center densities, without the use of any fit procedure. Because rho/sup(1)(r) is a sum of one-center multipolar densities, the Hartree potential is a function with a simple one-dimensional integral representation, and its matrix elements can be obtained by performing one-dimensional integrals over it. Since deltarho(r) is spatially slowly varying, the Hartree potential to which it corresponds and the matrix elements of this potential can accurately be evaluated on a relatively coarse coordinate space mesh, using fast Fourier transforms. The method is illustrated via molecular structure calculations for N2 and NH3. The calculations are accurate to a few percent when the required integrals over deltarho(r) and deltaV(r) are performed on a mesh of spacing 0.4 a.u. The N--N bond length and stretch frequency are found to equal 2.10 a.u. and 2.3 x 103 cm-1, respectively. The equilibrium N--H bond length and H--N--H angle are calculated to be 1.93 a.u. and 1050, respectively, while the NH3 inversion barrier turns out to equal 0.25 eV. These results are in good agreement with earlier calculations
Energy Technology Data Exchange (ETDEWEB)
Chauvin, C
2005-11-15
This thesis is devoted to the definition and the implementation of a multi-resolution method to determine the fundamental state of a system composed of nuclei and electrons. In this work, we are interested in the Density Functional Theory (DFT), which allows to express the Hamiltonian operator with the electronic density only, by a Coulomb potential and a non-linear potential. This operator acts on orbitals, which are solutions of the so-called Kohn-Sham equations. Their resolution needs to express orbitals and density on a set of functions owing both physical and numerical properties, as explained in the second chapter. One can hardly satisfy these two properties simultaneously, that is why we are interested in orthogonal and bi-orthogonal wavelets basis, whose properties of interpolation are presented in the third chapter. We present in the fourth chapter three dimensional solvers for the Coulomb's potential, using not only the preconditioning property of wavelets, but also a multigrid algorithm. Determining this potential allows us to solve the self-consistent Kohn-Sham equations, by an algorithm presented in chapter five. The originality of our method consists in the construction of the stiffness matrix, combining a Galerkin formulation and a collocation scheme. We analyse the approximation properties of this method in case of linear Hamiltonian, such as harmonic oscillator and hydrogen, and present convergence results of the DFT for small electrons. Finally we show how orbital compression reduces considerably the number of coefficients to keep, while preserving a good accuracy of the fundamental energy. (author)
International Nuclear Information System (INIS)
The methods that are actively used for electronic structure calculations of low-lying states of heavy- and superheavy-element compounds are briefly described. The advantages and disadvantages of the Dirac-Coulomb-Breit Hamiltonian, Huzinaga-type potential, shape-consistent Relativistic Effective Core Potential (RECP), and Generalized RECP are discussed. The nonvariational technique of the electron-structure restoration in atomic cores after the RECP calculation of a molecule is presented. The features of the approaches accounting for electron correlation, the configuration interaction and coupled cluster methods, are also described. The results of calculations on E113, E114, U, and other heavy-atom systems are presented
International Nuclear Information System (INIS)
Systematic electronic structure calculations have been performed for (CH3CN)n-(n=2-10) anion clusters with the hybrid B3LYP and non-hybrid PW91 density-functional methods in order to understand the stabilization mechanism of an acetonitrile dimer radical anion core by solvent molecules. Since the excess negative charge is mainly localized on N atoms in the dimer anion core, solvent acetonitrile molecules are bound to the N atoms by C-H...Nδ- hydrogen-bond-like attractive interaction with the binding energy per bond being about 10-13kcal/mol. Due to this stabilization mechanism, the anion cluster for n>=4-6 is stable with respect to the electron autodetachment. Geometry optimization was also carried out for the (CH3CN)6- anion cluster where an excess electron was internally trapped. The size dependence of the stabilization energy and vertical detachment energy for the (CH3CN)n- anion clusters is discussed
Electronic Structure Calculations and Adaptation Scheme in Multi-core Computing Environments
Energy Technology Data Exchange (ETDEWEB)
Seshagiri, Lakshminarasimhan; Sosonkina, Masha; Zhang, Zhao
2009-05-20
Multi-core processing environments have become the norm in the generic computing environment and are being considered for adding an extra dimension to the execution of any application. The T2 Niagara processor is a very unique environment where it consists of eight cores having a capability of running eight threads simultaneously in each of the cores. Applications like General Atomic and Molecular Electronic Structure (GAMESS), used for ab-initio molecular quantum chemistry calculations, can be good indicators of the performance of such machines and would be a guideline for both hardware designers and application programmers. In this paper we try to benchmark the GAMESS performance on a T2 Niagara processor for a couple of molecules. We also show the suitability of using a middleware based adaptation algorithm on GAMESS on such a multi-core environment.
Energy Technology Data Exchange (ETDEWEB)
Hinsche, Nicki; Yavorski, Bogdan; Zahn, Peter; Mertig, Ingrid [Martin-Luther-Universitaet, Institut fuer Physik, Halle/S. (Germany)
2010-07-01
Starting from bulk silicon, we studied the valley splitting due to symmetry breaking that occurs in rolled-up Si. Valley splitting in Si was studied recently because of tetragonal distortion and quantum well effects in heterostructures. The new aspect in nowadays experimentally accessible rolled-up Si tubes is that symmetry breaking occurs in all spatial directions. As a result, splitting of the six-fold degenerate conduction-band minimum is expected to be lifted. This has a strong influence on the transport properties as well. In detail, the anisotropy of the effective masses of charge carriers contributing to the conductivity in different directions are studied in dependence on the applied strain. The electronic structure is calculated self consistently within the framework of density functional theory. The transport properties of the promising thermoelectric material are studied in the diffusive limit of transport applying the Boltzmann theory in relaxation time approximation.
Is C50 a superaromat? Evidence from electronic structure and ring current calculations.
Matías, Ana Sanz; Havenith, Remco W A; Alcamí, Manuel; Ceulemans, Arnout
2016-04-28
The fullerene-50 is a 'magic number' cage according to the 2(N + 1)(2) rule. For the three lowest isomers of C50 with trigonal and pentagonal symmetries, we calculate the sphericity index, the spherical parentage of the occupied π-orbitals, and the current density in an applied magnetic field. The minimal energy isomer, with D3 symmetry, comes closest to a spherical aromat or a superaromat. In the D5h bond-stretch isomers the electronic structure shows larger deviations from the ideal spherical shells, with hybridisation or even reversal of spherical parentages. It is shown that relative stabilities of fullerene cages do not correlate well with aromaticity, unlike the magnetic properties which are very sensitive indicators of spherical aromaticity. Superaromatic diamagnetism in the D3 cage is characterized by global diatropic currents, which encircle the whole cage. The breakdown of sphericity in the D5h cages gives rise to local paratropic countercurrents. PMID:26444568
International Nuclear Information System (INIS)
Graphical abstract: Gas-phase UV photoelectron spectrum of the thermal decomposition of 5-aminotetrazole (5ATZ), obtained at 245 oC, and mechanism underlying the thermal dissociation of 2H-5ATZ. Research highlights: → Electronic structure of 5ATZ studied by photoelectron spectroscopy. → Gas-phase 5-ATZ exists mainly as the 2H-tautomer. → Thermal decomposition of 5ATZ gives N2, NH2CN, HN3 and HCN, at 245 oC. → HCN can be originated from a carbene intermediate. - Abstract: The electronic properties and thermal decomposition of 5-aminotetrazole (5ATZ) are investigated using UV photoelectron spectroscopy (UVPES) and theoretical calculations. Simulated spectra of both 1H- and 2H-5ATZ, based on electron propagator methods, are produced in order to study the relative gas-phase tautomer population. The thermal decomposition results are rationalized in terms of intrinsic reaction coordinate (IRC) calculations. 5ATZ yields a HOMO ionization energy of 9.44 ± 0.04 eV and the gas-phase 5ATZ assumes mainly the 2H-form. The thermal decomposition of 5ATZ leads to the formation of N2, HN3 and NH2CN as the primary products, and HCN from the decomposition of a intermediate CH3N3 compound. The reaction barriers for the formation of HN3 and N2 from 2H-5ATZ are predicted to be ∼228 and ∼150 kJ/mol, at the G2(MP2) level, respectively. The formation of HCN and HNNH from the thermal decomposition of a CH3N3 carbene intermediate is also investigated.
Institute of Scientific and Technical Information of China (English)
WU WenXia; XUE ZhiYong; HONG Xing; LI XiuMei; GUO YongQuan
2009-01-01
The valence electronic structures of Fe, Co and Ni have been investigated with Empirical Electron Theory of Solids and Molecules. The magnetic moments, Curie temperature, cohesive energy and melting point have been calculated according to the valence electronic structure. These calculations fit the experimental data very well. Based on the calculations, the magnetic moments are proportional to the number of 3d magnetic electrons. Curie temperatures are related to the magnetic electrons and the bond lengths between magnetic atoms. Cohesive energies increase with the increase of the number of covalent electrons, and the decrease of the number of magnetic and dumb pair electrons. The melting point is mainly related to the number of covalent electron pairs distributed in the strongest bond. The contribution from the lattice electrons is very small, the dumb pair electrons weaken the melting point; however, the contribution to melting point of the magnetic electrons can be neglected. It reveals that the magnetic and thermal properties are closely related to the valence electronic structures, and the changes or transitions between the electrons obviously affect the physical properties.
Institute of Scientific and Technical Information of China (English)
无
2009-01-01
The valence electronic structures of Fe, Co and Ni have been investigated with Empirical Electron Theory of Solids and Molecules. The magnetic moments, Curie temperature, cohesive energy and melting point have been calculated according to the valence electronic structure. These calculations fit the experimental data very well. Based on the calculations, the magnetic moments are proportional to the number of 3d magnetic electrons. Curie temperatures are related to the magnetic electrons and the bond lengths between magnetic atoms. Cohesive energies increase with the increase of the number of covalent electrons, and the decrease of the number of magnetic and dumb pair electrons. The melting point is mainly related to the number of covalent electron pairs distributed in the strongest bond. The contribution from the lattice electrons is very small, the dumb pair electrons weaken the melting point; however, the contribution to melting point of the magnetic electrons can be neglected. It reveals that the magnetic and thermal properties are closely related to the valence electronic structures, and the changes or transitions between the electrons obviously affect the physical properties.
An approach to first principles electronic structure calculation by symbolic-numeric computation
Directory of Open Access Journals (Sweden)
Akihito Kikuchi
2013-04-01
Full Text Available There is a wide variety of electronic structure calculation cooperating with symbolic computation. The main purpose of the latter is to play an auxiliary role (but not without importance to the former. In the field of quantum physics [1-9], researchers sometimes have to handle complicated mathematical expressions, whose derivation seems almost beyond human power. Thus one resorts to the intensive use of computers, namely, symbolic computation [10-16]. Examples of this can be seen in various topics: atomic energy levels, molecular dynamics, molecular energy and spectra, collision and scattering, lattice spin models and so on [16]. How to obtain molecular integrals analytically or how to manipulate complex formulas in many body interactions, is one such problem. In the former, when one uses special atomic basis for a specific purpose, to express the integrals by the combination of already known analytic functions, may sometimes be very difficult. In the latter, one must rearrange a number of creation and annihilation operators in a suitable order and calculate the analytical expectation value. It is usual that a quantitative and massive computation follows a symbolic one; for the convenience of the numerical computation, it is necessary to reduce a complicated analytic expression into a tractable and computable form. This is the main motive for the introduction of the symbolic computation as a forerunner of the numerical one and their collaboration has won considerable successes. The present work should be classified as one such trial. Meanwhile, the use of symbolic computation in the present work is not limited to indirect and auxiliary part to the numerical computation. The present work can be applicable to a direct and quantitative estimation of the electronic structure, skipping conventional computational methods.
International Nuclear Information System (INIS)
Quantum mechanical ab initio calculation constitutes the biggest portion of the computer time in material science and chemical science simulations. As a computer center like NERSC, to better serve these communities, it will be very useful to have a prediction for the future trends of ab initio calculations in these areas. Such prediction can help us to decide what future computer architecture can be most useful for these communities, and what should be emphasized on in future supercomputer procurement. As the size of the computer and the size of the simulated physical systems increase, there is a renewed interest in using the real space grid method in electronic structure calculations. This is fueled by two factors. First, it is generally assumed that the real space grid method is more suitable for parallel computation for its limited communication requirement, compared with spectrum method where a global FFT is required. Second, as the size N of the calculated system increases together with the computer power, O(N) scaling approaches become more favorable than the traditional direct O(N3) scaling methods. These O(N) methods are usually based on localized orbital in real space, which can be described more naturally by the real space basis. In this report, the author compares the real space methods versus the traditional plane wave (PW) spectrum methods, for their technical pros and cons, and the possible of future trends. For the real space method, the author focuses on the regular grid finite different (FD) method and the finite element (FE) method. These are the methods used mostly in material science simulation. As for chemical science, the predominant methods are still Gaussian basis method, and sometime the atomic orbital basis method. These two basis sets are localized in real space, and there is no indication that their roles in quantum chemical simulation will change anytime soon. The author focuses on the density functional theory (DFT), which is the most
Magnetic state and electronic structure of plutonium from "first principles" calculations
Czech Academy of Sciences Publication Activity Database
Anisimov, V.I.; Shorikov, A.O.; Kuneš, Jan
444-445, - (2007), s. 42-49. ISSN 0925-8388 Institutional research plan: CEZ:AV0Z10100521 Keywords : metals * electron-electron interactions * electronic band structure Subject RIV: BM - Solid Matter Physics ; Magnetism Impact factor: 1.455, year: 2007
Hydrogen trapping in δ-Pu: insights from electronic structure calculations
International Nuclear Information System (INIS)
Density functional theory calculations have been performed to provide details of the structural and charge-transfer details related to the solid solution of hydrogen in (δ)-plutonium. We follow the Flanagan model that outlines the process by which hydrogen interacts with a metal to produce hydride phases, via a sequence of surface, interstitial and defect-bound (trapped) states. Due to the complexities of the electronic structure in plutonium solid-state systems, we take the pragmatic approach of adopting the ‘special quasirandom structure’ to disperse the atomic magnetic moments. We find that this approach produces sound structural and thermodynamic properties in agreement with the available experimental data. In δ-Pu, hydrogen has an exothermic binding energy to all of the states relevant in the Flanagan model, and, furthermore, is anionic in all these states. The charge transfer is maximized (i.e. most negative for hydrogen) in the hydride phase. The pathway from surface to hydride is sequentially exothermic, in the order surface < interstitial < grain boundary < vacancy < hydride (hydride being the most exothermic state). Thus, we find that there is no intermediate state that involves an endothermic increase in energy, consistent with the general experimental observations that the hydriding reaction in plutonium metal can proceed with zero apparent activation barrier. (paper)
Calculated Electronic and Magnetic Structure of Screw Dislocations in Alpha Iron
Energy Technology Data Exchange (ETDEWEB)
Odbadrakh, Khorgolkhuu [ORNL; Rusanu, Aurelian [ORNL; Stocks, George Malcolm [ORNL; Samolyuk, German D [ORNL; Eisenbach, Markus [ORNL; Wang, Yang [Pittsburgh Supercomputing Center; Nicholson, Don M [ORNL
2011-01-01
Local atomic magnetic moments in crystalline Fe are perturbed by the presence of dislocations. The effects are most pronounced near the dislocation core and decay slowly as the strain field of the dislocation decreases with distance. We have calculated the local moments using the Locally Self-consistent Multiple Scattering (LSMS) method for an 1848 atom supercell containing a screw- dislocation quadrupole. The atomic positions were determined by relaxation with an embedded atom force field. Finite size effects are found to be significant for this small cell size indicating that dislocation cores affect the electronic structure and magnetic moments of neighboring dislocations. The influence of neighboring dislocations point to a need to study individual dislocations from first principles just as they appear amidst surrounding atoms in large scale classical force field simulations. An approach for the use of the LSMS to calculate local moments in sub-volumes of large atomic configurations generated in the course of classical MD simulation of dislocation dynamics is discussed.
Institute of Scientific and Technical Information of China (English)
LIU Zhilin; LIN Cheng; LIU Yan; GUO Yanchang
2005-01-01
Combined with the phase transformations in rolling, the phase configuration, the tensile strength, and the yield strength with different terminal rolling grain sizes in Q235 strip steel have been theoretically calculated using the covalent electron number (nA) of the strongest bond in phase cells and the interface electron density difference (Ap) in alloys. The calculated results agree well with the results of real production. Therefore, the calculation method of terminal rolling tensile and yield strength in the non-quenched-tempered steel containing pearlite is given by the alloying electron structure parameters.
Lihua Xiao; Yuchang Su; Hongyang Chen; Min Jiang; Sainan Liu; Zexing Hu; Ruifeng Liu; Ping Peng; Yuanlong Mu; Diya Zhu
2011-01-01
The electronic structure and the optical performance of YB6 were investigated by first-principles calculations within the framework of density functional theory. It was found that the calculated results are in agreement with the relevant experimental data. Our theoretical studies showed that YB6 is a promising solar radiation shielding material for windows.
Directory of Open Access Journals (Sweden)
Lihua Xiao
2011-06-01
Full Text Available The electronic structure and the optical performance of YB6 were investigated by first-principles calculations within the framework of density functional theory. It was found that the calculated results are in agreement with the relevant experimental data. Our theoretical studies showed that YB6 is a promising solar radiation shielding material for windows.
International Nuclear Information System (INIS)
Density functional theory (DFT) is the most widely used ab initio method in material simulations. It accounts for 75% of the NERSC allocation time in the material science category. The DFT can be used to calculate the electronic structure, the charge density, the total energy and the atomic forces of a material system. With the advance of the HPC power and new algorithms, DFT can now be used to study thousand atom systems in some limited ways (e.g, a single selfconsistent calculation without atomic relaxation). But there are many problems which either requires much larger systems (e.g, >100,000 atoms), or many total energy calculation steps (e.g. for molecular dynamics or atomic relaxations). Examples include: grain boundary, dislocation energies and atomic structures, impurity transport and clustering in semiconductors, nanostructure growth, electronic structures of nanostructures and their internal electric fields. Due to the O(N3) scaling of the conventional DFT algorithms (as implemented in codes like Qbox, Paratec, Petots), these problems are beyond the reach even for petascale computers. As the proposed petascale computers might have millions of processors, new computational paradigms and algorithms are needed to solve the above large scale problems. In particular, O(N) scaling algorithms with parallelization capability up to millions of processors are needed. For a large material science problem, a natural approach to achieve this goal is by divide-and-conquer method: to spatially divide the system into many small pieces, and solve each piece by a small local group of processors. This solves the O(N) scaling and the parallelization problem at the same time. However, the challenge of this approach is for how to divide the system into small pieces and how to patch them up without the trace of the spatial division. Here, we present a linear scaling 3 dimensional fragment (LS3DF) method which uses a novel division-patching scheme that cancels out the artificial
Energy Technology Data Exchange (ETDEWEB)
Lin, Lin; Yang, Chao; Lu, Jiangfeng; Ying, Lexing; E, Weinan
2009-09-25
We present an efficient parallel algorithm and its implementation for computing the diagonal of $H^-1$ where $H$ is a 2D Kohn-Sham Hamiltonian discretized on a rectangular domain using a standard second order finite difference scheme. This type of calculation can be used to obtain an accurate approximation to the diagonal of a Fermi-Dirac function of $H$ through a recently developed pole-expansion technique \\cite{LinLuYingE2009}. The diagonal elements are needed in electronic structure calculations for quantum mechanical systems \\citeHohenbergKohn1964, KohnSham 1965,DreizlerGross1990. We show how elimination tree is used to organize the parallel computation and how synchronization overhead is reduced by passing data level by level along this tree using the technique of local buffers and relative indices. We analyze the performance of our implementation by examining its load balance and communication overhead. We show that our implementation exhibits an excellent weak scaling on a large-scale high performance distributed parallel machine. When compared with standard approach for evaluating the diagonal a Fermi-Dirac function of a Kohn-Sham Hamiltonian associated a 2D electron quantum dot, the new pole-expansion technique that uses our algorithm to compute the diagonal of $(H-z_i I)^-1$ for a small number of poles $z_i$ is much faster, especially when the quantum dot contains many electrons.
International Nuclear Information System (INIS)
We present an efficient parallel algorithm and its implementation for computing the diagonal of H-1 where H is a 2D Kohn-Sham Hamiltonian discretized on a rectangular domain using a standard second order finite difference scheme. This type of calculation can be used to obtain an accurate approximation to the diagonal of a Fermi-Dirac function of H through a recently developed pole-expansion technique LinLuYingE2009. The diagonal elements are needed in electronic structure calculations for quantum mechanical systems HohenbergKohn1964, KohnSham 1965,DreizlerGross1990. We show how elimination tree is used to organize the parallel computation and how synchronization overhead is reduced by passing data level by level along this tree using the technique of local buffers and relative indices. We analyze the performance of our implementation by examining its load balance and communication overhead. We show that our implementation exhibits an excellent weak scaling on a large-scale high performance distributed parallel machine. When compared with standard approach for evaluating the diagonal a Fermi-Dirac function of a Kohn-Sham Hamiltonian associated a 2D electron quantum dot, the new pole-expansion technique that uses our algorithm to compute the diagonal of (H-zi I)-1 for a small number of poles zi is much faster, especially when the quantum dot contains many electrons.
Electronic structure and optical properties of B/P-doped amorphous Si calculated by first-principles
International Nuclear Information System (INIS)
Highlights: • Short-range order in a-Si lead to the similar electronic structure and optical properties with c-Si. • Long-range disorder of a-Si lead to the different electronic structure and optical properties. • Localized states predominately determine the optical properties in visible-light region of a-Si. • B/P-doping have no obvious effects for the electronic structure and optical properties of a-Si. - Abstract: In order to understand the electronic structures, optical properties, and explain the experimental observations of B/P-doped amorphous Si, the relevant micro-structure and properties have been calculated by simulated annealing and DFT+U methods. Based on the calculated results, the short-range order features of micro-structure in amorphous Si lead to the similar electronic structure and optical properties with crystalline Si, owing to the short-range order reflects the nature of atomic chemical bonding and plays a major role in the decision of fundamental characteristics of amorphous Si. What is important, the long-range disorder features of micro-structure lead to the different electronic structure and optical properties of amorphous Si, in compared with crystalline Si. Especially, the localized states caused by structural defects predominately determined the optical properties in visible-light region. The findings in the present work could well explain the experimental observations in literatures, and are helpful for the development of amorphous Si based functional materials
Calculation of electronic structure of YBa2Cu3O7-δ in LCAO MO approximation
International Nuclear Information System (INIS)
On the basis of calculation by CNDO methods study of copper oxidation degree and valence, oxygen bond character, instability of charge states, their causes is carried out. An attempt of studying electron state density in YBa2Cu3O7-δ compound and also electron behaviour near Fermi surface depending on oxygen and copper atoms state is made using cluster calculations. It is supposed that at 0.5 2Cu3O7-δ copper atoms are in a state close to 3d10 state, therewith oxygen atoms are in a state close to 0-1 one
International Nuclear Information System (INIS)
Highlights: • The band gaps for CaB6, SrB6 and BaB6 depend sensitively on the values of lattice constant a and positional parameter z. • The order in elastic anisotropy is CaB6 > SrB6 > BaB6. • There are LO/TO splitting lines in the range of 5–10 THz at G point. - Abstract: The electronic structures, mechanical and thermodynamic properties of alkaline-earth hexaborides MB6 (M = Ca, Sr or Ba) are calculated from first principles using density functional theory combined with the quasi-harmonic approximation. These three alkaline-earth hexaborides are semiconductors with a slightly increased trend for their band gaps as M orders from Ca to Ba. Their band gaps depend sensitively on the values of lattice constant a and internal parameter z. The polycrystalline values of the elastic constants and bulk, shear and Young’s moduli are consistent with those determined experimentally. All alkaline-earth hexaborides have strongly anisotropic elastic properties in the order of CaB6 > SrB6 > BaB6. By using the phonon calculations, the thermodynamic properties are investigated. The obtained phonon dispersion relations for CaB6, SrB6, and BaB6 show similar features and there are LO/TO splitting lines in the range of 5–10 THz. Finally, the thermal conductivities of CaB6, SrB6 and BaB6 are evaluated via Clarke’s model and Cahill’s model
Corsetti, Fabiano
2014-01-01
The implementation of the orbital minimization method (OMM) for solving the self-consistent Kohn-Sham (KS) problem for electronic structure calculations in a basis of non-orthogonal numerical atomic orbitals of finite-range is reported. We explore the possibilities for using the OMM as an exact cubic-scaling solver for the KS problem, and compare its performance with that of explicit diagonalization in realistic systems. We analyze the efficiency of the method depending on the choice of line search algorithm and on two free parameters, the scale of the kinetic energy preconditioning and the eigenspectrum shift. The results of several timing tests are then discussed, showing that the OMM can achieve a noticeable speedup with respect to diagonalization even for minimal basis sets for which the number of occupied eigenstates represents a significant fraction of the total basis size (>15%). We investigate the hard and soft parallel scaling of the method on multiple cores, finding a performance equal to or better ...
Liu, Ming-Yang; Huang, Yang; Chen, Qing-Yuan; Cao, Chao; He, Yao
2016-01-01
We study the equilibrium geometry and electronic structure of alloyed and doped arsenene sheets based on the density functional theory calculations. AsN, AsP and SbAs alloys possess indirect band gap and BiAs is direct band gap. Although AsP, SbAs and BiAs alloyed arsenene sheets maintain the semiconducting character of pure arsenene, they have indirect-direct and semiconducting-metallic transitions by applying biaxial strain. We find that B- and N-doped arsenene render p-type semiconducting character, while C- and O-doped arsenene are metallic character. Especially, the C-doped arsenene is spin-polarization asymmetric and can be tuned into the bipolar spin-gapless semiconductor by the external electric field. Moreover, the doping concentration can effectively affect the magnetism of the C-doped system. Finally, we briefly study the chemical molecule adsorbed arsenene. Our results may be valuable for alloyed and doped arsenene sheets applications in mechanical sensors and spintronic devices in the future. PMID:27373712
Liu, Ming-Yang; Huang, Yang; Chen, Qing-Yuan; Cao, Chao; He, Yao
2016-07-01
We study the equilibrium geometry and electronic structure of alloyed and doped arsenene sheets based on the density functional theory calculations. AsN, AsP and SbAs alloys possess indirect band gap and BiAs is direct band gap. Although AsP, SbAs and BiAs alloyed arsenene sheets maintain the semiconducting character of pure arsenene, they have indirect-direct and semiconducting-metallic transitions by applying biaxial strain. We find that B- and N-doped arsenene render p-type semiconducting character, while C- and O-doped arsenene are metallic character. Especially, the C-doped arsenene is spin-polarization asymmetric and can be tuned into the bipolar spin-gapless semiconductor by the external electric field. Moreover, the doping concentration can effectively affect the magnetism of the C-doped system. Finally, we briefly study the chemical molecule adsorbed arsenene. Our results may be valuable for alloyed and doped arsenene sheets applications in mechanical sensors and spintronic devices in the future.
Theoretical study of the electronic structure with dipole moment calculations of barium monofluoride
Tohme, Samir N.; Korek, Mahmoud
2015-12-01
The potential energy curves have been investigated for the 41 lowest doublet and quartet electronic states in the 2s+1Λ± representation below 55,000 cm-1 of the molecule BaF via CASSCF and MRCI (single and double excitations with Davidson correction) calculations. Twenty-five electronic states have been studied here theoretically for the first time. The crossing and avoided crossing of 20 doublet electronic states have been studied in the region 30,000-50,000 cm-1. The harmonic frequency ωe, the internuclear distance Re, the rotational constant Be, the electronic energy with respect to the ground state Te, and the permanent and transition dipole moments have been calculated in addition to static dipole polarizability of the ground state. By using the canonical functions approach, the eigenvalue Ev, the rotational constant Bv, and the abscissas of the turning points Rmin and Rmax have been calculated for the electronic states up to the vibrational level v=98. The comparison of these values with the theoretical results available in the literature shows a very good agreement.
International Nuclear Information System (INIS)
One calculated four certain modifications of SrZrO3 crystal of various symmetry: a cubic one, a tetragonal one and two orthorhombic ones, by the density functional technique in the basis of linear combination of atomic orbitals. One carried out comparison analysis of electron properties of the investigated crystals based on the calculated band structures and distribution densities of electron states (the complete ones and atomic state designed ones). The calculation base relative stability of different modifications correlates adequately with the experimental data on phase transitions in SrZrO3 crystal: less symmetric low-temperature modifications are more stable ones
International Nuclear Information System (INIS)
We present first-principles study of the electronic and the optical properties for the intermetallic trialuminides ScAl3 compound using the full-potential linear augmented plane wave method within density-functional theory. We have employed the generalized gradient approximation (GGA), which is based on exchange-correlation energy optimization to calculate the total energy. Also we have used the Engel-Vosko GGA formalism, which optimizes the corresponding potential for calculating the electronic band structure and optical properties. The electronic specific heat coefficient (γ), which is a function of density of states, can be calculated from the density of states at Fermi energy N(EF). The N(EF) of the phase L12 is found to be lower than that of D022 structure which confirms the stability of L12 structure. We found that the dispersion of the band structure of D022 is denser than L12 phase. The linear optical properties were calculated. The evaluations are based on calculations of the energy band structure. - Graphical abstract: Crystal structure of L12 (a) and D022 (b) phases of ScAl3 compound.
Energy Technology Data Exchange (ETDEWEB)
Werwiński, M. [Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego 17, 60-179 Poznań (Poland); Szajek, A. [Institute of Molecular Physics, Polish Academy of Sciences, Smoluchowskiego 17, 60-179 Poznań (Poland); Centre for Advanced Materials and Smart Structures, Polish Academy of Sciences, Okólna 2, 50-950 Wrocław (Poland); Ślebarski, A. [Institute of Physics, University of Silesia, Uniwersytecka 4, 40-007 Katowice (Poland); Centre for Advanced Materials and Smart Structures, Polish Academy of Sciences, Okólna 2, 50-950 Wrocław (Poland); Kaczorowski, D., E-mail: D.Kaczorowski@int.pan.wroc.pl [Institute of Low Temperature and Structure Research, Polish Academy of Sciences, P. O. Box 1410, 50-950 Wrocław (Poland); Centre for Advanced Materials and Smart Structures, Polish Academy of Sciences, Okólna 2, 50-950 Wrocław (Poland)
2015-10-25
The electronic structure of a heavy-fermion superconductor Ce{sub 2}PdIn{sub 8} was investigated by means of X-ray photoelectron spectroscopy (XPS) and ab initio density functional band structure calculations. The Ce 3d core-level XPS spectra point to stable trivalent configuration of Ce atoms that is also reproduced in the band structure calculations within the generalized gradient approximation GGA+U approach. Analysis of the 3d{sup 9}f{sup 2} weight in the 3d XPS spectra within the Gunnarsson-Schönhammer model suggests that the onsite hybridization energy between Ce 4f and the conduction band states, Δ{sub fs}, is ∼120 meV, which is about 30 meV larger than Δ{sub fs} in isostructural Ce{sub 2}TIn{sub 8} compounds with T = Co, Rh, and Ir. Taking into account a Coulomb repulsion U on both the Ce 4f and Pd 4d states in electronic band structure calculations, a satisfactory agreement was found between the calculated density of states (DOS) and the measured valence band XPS spectra. - Highlights: • XPS data validated strong electronic correlations in superconducting Ce{sub 2}PdIn{sub 8}. • DFT calculations reproduced XPS spectra measured for Ce{sub 2}PdIn{sub 8}. • Crucial role of Pd d electrons in the HF behavior of Ce{sub 2}PdIn{sub 8} was established.
The calculation of the electronic structure of deep traps in silicon
International Nuclear Information System (INIS)
In the present thesis the LMTO-ASA Green's function method is applied to a real system (Silicon crystal with trap) for the first time. The electronic structure of the point defects for H and He on the tetrahedral (Tp) sites in Si shows agreement with experimental results. Moreover, the theory can distinguish the two Tp sites for S, Se and Te. Then the diffusion of the chalcogen atoms in Si is investigated and the electronic structure of the 3d-transition element ions in Silicon is examined. (BHO)
Ab-initio calculations of electronic, transport, and structural properties of boron phosphide
Energy Technology Data Exchange (ETDEWEB)
Ejembi, J. I.; Nwigboji, I. H.; Franklin, L.; Malozovsky, Y.; Zhao, G. L.; Bagayoko, D., E-mail: diola-bagayoko@subr.edu [Department of Physics, Southern University and A and M College, Baton Rouge, Louisiana 70813 (United States)
2014-09-14
We present results from ab-initio, self-consistent density functional theory calculations of electronic and related properties of zinc blende boron phosphide (zb-BP). We employed a local density approximation potential and implemented the linear combination of atomic orbitals formalism. This technique follows the Bagayoko, Zhao, and Williams method, as enhanced by the work of Ekuma and Franklin. The results include electronic energy bands, densities of states, and effective masses. The calculated band gap of 2.02 eV, for the room temperature lattice constant of a=4.5383 Å, is in excellent agreement with the experimental value of 2.02±0.05 eV. Our result for the bulk modulus, 155.7 GPa, agrees with experiment (152–155 GPa). Our predictions for the equilibrium lattice constant and the corresponding band gap, for very low temperatures, are 4.5269 Å and 2.01 eV, respectively.
Electronic structure of cubic ScF$_3$ from first-principles calculations
Zhgun, P.; Bocharov, D.; Piskunov, S.; Kuzmin, A; Purans, J.
2012-01-01
The first-principles calculations have been performed to investigate the ground state properties of cubic scandium trifluoride (ScF$_3$) perovskite. Using modified hybrid exchange-correlation functionals within the density functional theory (DFT) we have comprehensively compared the electronic properties of ScF$_3$ obtained by means of the linear combination of atomic orbitals (LCAO) and projector augmented-waves (PAW) methods. Both methods allowed us to reproduce the lattice constant experim...
First principle calculations of structural phase transition and electronic properties in AmTe
Energy Technology Data Exchange (ETDEWEB)
Pataiya, Jagdeesh, E-mail: jagdish-pet@yahoo.co.in; Makode, C. [Sagar Institute of Research & Technology, Bhopal, 462041 (India); Aynyas, Mahendra [Department of Physics, C.S.A. Govt. P.G. College, Sehore, 466001 (India); Singh, A.; Sanyal, S. P. [Department of Physics, Barkatullah University, Bhopal, 462026 (India)
2015-06-24
The tight-binding linear muffin-tin orbital (TB-LMTO) with in the local density approximation is used to calculate total energy, lattice parameters, bulk modulus, density of states and energy band structure of americium telluride at ambient as well as at high pressure. It is found that AmTe is stable in NaCl – type structure under ambient pressure. The phase transition pressure was found to be 15.0 GPa from NaCl-type (B{sub 1}-phase) structure to CsCl-type (B{sub 2}-phase) structure for this compound. From energy band diagram it is observed that AmTe exhibit metallic behaviour. The calculated ground state properties such as lattice parameters and bulk modulus are in general good agreement with the available results.
First principle calculations of structural phase transition and electronic properties in AmTe
International Nuclear Information System (INIS)
The tight-binding linear muffin-tin orbital (TB-LMTO) with in the local density approximation is used to calculate total energy, lattice parameters, bulk modulus, density of states and energy band structure of americium telluride at ambient as well as at high pressure. It is found that AmTe is stable in NaCl – type structure under ambient pressure. The phase transition pressure was found to be 15.0 GPa from NaCl-type (B1-phase) structure to CsCl-type (B2-phase) structure for this compound. From energy band diagram it is observed that AmTe exhibit metallic behaviour. The calculated ground state properties such as lattice parameters and bulk modulus are in general good agreement with the available results
Unfolding method for the first-principles LCAO electronic structure calculations
Lee, Chi-Cheng; Yamada-Takamura, Yukiko; Ozaki, Taisuke
2012-01-01
Unfolding the band structure of a supercell to a normal cell enables us to investigate how symmetry breakers such as surfaces and impurities perturb the band structure of the normal cell. We generalize the unfolding method, originally developed based on Wannier functions, to the linear combination of atomic orbitals (LCAO) method, and present a general formula to calculate the unfolded spectral weight. The LCAO basis set is ideal for the unfolding method because of the invariance that basis f...
Nonrelativistic structure calculations of two-electron ions in a strongly coupled plasma environment
Energy Technology Data Exchange (ETDEWEB)
Bhattacharyya, S.; Saha, J. K.; Mukherjee, T. K.
2015-04-01
In this work, the controversy between the interpretations of recent measurements on dense aluminum plasma created with the Linac coherent light source (LCLS) x-ray free electron laser (FEL) and the Orion laser has been addressed. In both kinds of experiments, heliumlike and hydrogenlike spectral lines are used for plasma diagnostics. However, there exist no precise theoretical calculations for He-like ions within a dense plasma environment. The strong need for an accurate theoretical estimate for spectral properties of He-like ions in a strongly coupled plasma environment leads us to perform ab initio calculations in the framework of the Rayleigh-Ritz variation principle in Hylleraas coordinates where an ion-sphere potential is used. An approach to resolve the long-drawn problem of numerical instability for evaluating two-electron integrals with an extended basis inside a finite domain is presented here. The present values of electron densities corresponding to the disappearance of different spectral lines obtained within the framework of an ion-sphere potential show excellent agreement with Orion laser experiments in Al plasma and with recent theories. Moreover, this method is extended to predict the critical plasma densities at which the spectral lines of H-like and He-like carbon and argon ions disappear. Incidental degeneracy and level-crossing phenomena are being reported for two-electron ions embedded in strongly coupled plasma. Thermodynamic pressure experienced by the ions in their respective ground states inside the ion spheres is also reported.
Graphene allotropes: stability, structural and electronic properties from DF-TB calculations
ENYASHIN A.N.; Ivanovskii, A. L.
2010-01-01
Using the density-functional-based tight-binding method we performed a systematic comparative study of stability, structural and electronic properties for 12 various types of graphene allotropes, which are likely candidates for engineering of novel graphene-like materials.
Electronic fine structure calculation of [Gd(DOTA)(H2O)]- using LF-DFT: The zero field splitting
Senn, Florian; Helm, Lothar; Borel, Alain; Daul, Claude A.
2012-01-01
Zerfo field splitting plays an important role in determining the electron spin relaxation of Gd(III) in solution. We understand the ZFS as an effect depending on the f electron structure and treat it in the framework of ligand field-density functional theory (LF-DFT). We apply this theory to calculate the ZFS of [Gd(DOTA)(H2O)]- from first principles, having an insight concerning the contributions determining the ZFS.
Electronic fine structure calculation of [Gd(DOTA)(H₂O)]⁻ using LF-DFT: The zero field splitting
Senn, Florian; Helm, Lothar; Borel, Alain; Daul, Claude A.
2012-01-01
Zerfo field splitting plays an important role in determining the electron spin relaxation of Gd(III) in solution. We understand the ZFS as an effect depending on the f electron structure and treat it in the framework of ligand field-density functional theory (LF-DFT). We apply this theory to calculate the ZFS of [Gd(DOTA)(H₂O)]⁻ from first principles, having an insight concerning the contributions determining the ZFS.
All-electron Bethe-Salpeter calculations for shallow-core x-ray absorption near-edge structures
Olovsson, W.; Tanaka, I.; Mizoguchi, T.; Puschnig, P.; Ambrosch-Draxl, C.
2009-01-01
X-ray absorption near-edge structure spectra are calculated by fully solving the electron/core-hole Bethe-Salpeter equation (BSE) in an all-electron framework. We study transitions from shallow core states, including the Mg L2,3 edge in MgO, the Li K edge in the Li halides LiF, LiCl, LiBr, and LiI, as well as Li2O. We illustrate the advantage of the many-body approach over a core-hole supercell calculation. Both schemes lead to strongly bound excitons, but the nonlocal treatment of the electr...
Energy Technology Data Exchange (ETDEWEB)
Hegde, Ganesh, E-mail: ganesh.h@ssi.samsung.com; Bowen, R. Chris [Advanced Logic Lab, Samsung Semiconductor Inc., Austin, TX 78754 (United States)
2015-10-15
The accuracy of a single s-orbital representation of Cu towards enabling multi-thousand atom ab initio calculations of electronic structure is evaluated in this work. If an electrostatic compensation charge of 0.3 electron per atom is used in this basis representation, the electronic transmission in bulk and nanocrystalline Cu can be made to compare accurately to that obtained with a Double Zeta Polarized basis set. The use of this representation is analogous to the use of single band effective mass representation for semiconductor electronic structure. With a basis of just one s-orbital per Cu atom, the representation is extremely computationally efficient and can be used to provide much needed ab initio insight into electronic transport in nanocrystalline Cu interconnects at realistic dimensions of several thousand atoms.
International Nuclear Information System (INIS)
The accuracy of a single s-orbital representation of Cu towards enabling multi-thousand atom ab initio calculations of electronic structure is evaluated in this work. If an electrostatic compensation charge of 0.3 electron per atom is used in this basis representation, the electronic transmission in bulk and nanocrystalline Cu can be made to compare accurately to that obtained with a Double Zeta Polarized basis set. The use of this representation is analogous to the use of single band effective mass representation for semiconductor electronic structure. With a basis of just one s-orbital per Cu atom, the representation is extremely computationally efficient and can be used to provide much needed ab initio insight into electronic transport in nanocrystalline Cu interconnects at realistic dimensions of several thousand atoms
Photophysics of Auramine-O: electronic structure calculations and nonadiabatic dynamics simulations.
Xie, Bin-Bin; Xia, Shu-Hua; Chang, Xue-Ping; Cui, Ganglong
2016-01-01
Diphenylmethane dyes are very useful photoinduced molecular rotors; however, their photophysical mechanisms are still elusive until now. In this work, we adopted combined static electronic structure calculations (MS-CASPT2//CASSCF) and trajectory-based surface-hopping dynamics simulations (OM2/MRCI) to study the S1 excited-state relaxation mechanism of a representative diphenylmethane dye Auramine-O. On the basis of the optimized S1 minima and the computed emission bands, we have for the first time assigned experimentally proposed three transient states (i.e. S1-LE, S1-I1 or S1-I2, and S1-II). Mechanistically, upon irradiation to the S1 state, the system first relaxes to the locally excited S1 minimum (S1-LE). Starting from this point, there exist two kinds of relaxation paths to S1-II. In the sequential path, the system first evolves into S1-I1 or S1-I2 and then runs into S1-II; in the concerted one, the system, bypassing S1-I1 and S1-I2, directly runs into S1-II. In addition, the system can decay to the S0 state in the vicinity of three S1/S0 conical intersections i.e. S1S0-I1, S1S0-I2, and S1S0-II. In the S1 dynamic simulations, 54% trajectories decay to the S0 state via S1S0-II; the remaining trajectories are de-excited to the S0 state via S1S0-I1 (11%) and S1S0-I2 (35%). Our present theoretical investigation does not support the experimentally proposed S1 excited-state hypothesis that the intramolecular rotation of the two dimethyl groups around the C-N bond is responsible for the rapid decay of the emission band at about 500 nm; instead, it should be heavily interrelated with the rotation of the two dimethylanilino groups. Finally, this work provides important mechanistic insights into similar diphenylmethane dyes. PMID:26615798
Structural, electronic and optical properties of CdxZn1−xS alloys from first-principles calculations
International Nuclear Information System (INIS)
Structural, electronic and optical properties as well as structural phase transitions of ternary alloy CdxZn1−xS have been investigated using the first-principles calculations based on the density functional theory. We found that the crystal structure of CdxZn1−xS alloys transforms from wurtzite to zinc blende as Cd content of x=0.83. Effect of Cd content on electronic structures of CdxZn1−xS alloys has been studied. The bandgaps of CdxZn1−xS alloys with wurtzite and zinc blende structures decrease with the increase of Cd content. Furthermore, dielectric constant and absorption coefficient also have been discussed in detail. - Highlights: • The systematic calculation on structural, electronic and optical properties of ternary CdxZn1−xS alloy are performed. • We build two structures of ternary CdxZn1−xS: wurtzite and zinc blende. • The fundamental bandgap decreases, accompanying this structural phase transition. • We report the structural phase transforms of CdxZn1−xS alloys from wurtzite to zinc blende
Gidofalvi, Gergely
2014-01-01
Molecule-optimized basis sets, based on approximate natural orbitals, are developed for accelerating the convergence of quantum calculations with strongly correlated (multi-referenced) electrons. We use a low-cost approximate solution of the anti-Hermitian contracted Schr{\\"o}dinger equation (ACSE) for the one- and two-electron reduced density matrices (RDMs) to generate an approximate set of natural orbitals for strongly correlated quantum systems. The natural-orbital basis set is truncated to generate a molecule-optimized basis set whose rank matches that of a standard correlation-consistent basis set optimized for the atoms. We show that basis-set truncation by approximate natural orbitals can be viewed as a one-electron unitary transformation of the Hamiltonian operator and suggest an extension of approximate natural-orbital truncations through two-electron unitary transformations of the Hamiltonian operator, such as those employed in the solution of the ACSE. The molecule-optimized basis set from the ACS...
Calculation of the electronic and magnetic structures of 3d impurities in the Hcp Fe matrix
International Nuclear Information System (INIS)
In this work we investigate the local magnetic properties and the electronic structure of HCP Fe, as well introducing transition metals atoms 3d (Cs, Ti, Cr, Mn, Co, Ni, Cu, Zn) in HCP iron matrix. We employed the discrete variational method (DVM), which is an orbital molecular method which incorporate the Hartree-Fock-Slater theory and the linear combination of atomic orbitals (LCAO), in the self-consistent charge approximation and the local density approximation of Von Barth and Hedin to the exchange-correlation potential. We used the embedded cluster model to investigate the electronic structure and the local magnetic properties for the central atom of a cluster of 27 atoms immersed in the microcrystal representing the HCP Fe. (author)
Electronic Structure and Elastic Properties of Ti3AlC from First-Principles Calculations
Institute of Scientific and Technical Information of China (English)
DU Yu-Lei
2009-01-01
We perform a first-principles study on the electronic structure and elastic properties of Ti3AlC with an antiper-ovskite structure. The absence of band gap at the Fermi level and the finite value of the density of states at the Fermi energy reveal the metallic behavior of this compound. The elastic constants of Ti_3AlC are derived yielding c_(11)=356 GPa, c_(12)= 55 GPa, c_(44)=157 GPa. The bulk modulus B, shear modulus G and Young's modulus E are determined to be 156, 151 and 342 GPa, respectively. These properties are compared with those of Ti_3AlC_2 and Ti_2AlC with a layered structure in the Ti-Al-C system and Fe_3AlC with the same antiperovskite structure.
Ab-initio calculation of electronic structure of partially inverted manganese ferrite
Czech Academy of Sciences Publication Activity Database
Chlan, V.; Novák, Pavel
2010-01-01
Roč. 322, 9-12 (2010), s. 1056-1058. ISSN 0304-8853 R&D Projects: GA ČR GA202/08/0541; GA ČR GA202/06/0051 Institutional research plan: CEZ:AV0Z10100521 Keywords : manganese ferrite * electronic structure Subject RIV: BM - Solid Matter Physics ; Magnetism Impact factor: 1.689, year: 2010
DEFF Research Database (Denmark)
Romero, N. A.; Glinsvad, Christian; Larsen, Ask Hjorth; Enkovaara, J.; Shende, S.; Morozov, V. A.; Mortensen, Jens Jørgen
2013-01-01
Density function theory (DFT) is the most widely employed electronic structure method because of its favorable scaling with system size and accuracy for a broad range of molecular and condensed-phase systems. The advent of massively parallel supercomputers has enhanced the scientific community...... GPAW code was ported an optimized for the Blue Gene/P architecture. We present our algorithmic parallelization strategy and interpret the results for a number of benchmark test cases....
Banerjee, Amartya S; Hu, Wei; Yang, Chao; Pask, John E
2016-01-01
The Discontinuous Galerkin (DG) electronic structure method employs an adaptive local basis set to solve the equations of density functional theory in a discontinuous Galerkin framework. The methodology is implemented in the Discontinuous Galerkin Density Functional Theory (DGDFT) code for large-scale parallel electronic structure calculations. In DGDFT, the basis is generated on-the-fly to capture the local material physics, and can systematically attain chemical accuracy with only a few tens of degrees of freedom per atom. Hence, DGDFT combines the key advantage of planewave basis sets in terms of systematic improvability with that of localized basis sets in reducing basis size. A central issue for large-scale calculations, however, is the computation of the electron density from the discretized Hamiltonian in an efficient and scalable manner. We show in this work how Chebyshev polynomial filtered subspace iteration (CheFSI) can be used to address this issue and push the envelope in large-scale materials si...
International Nuclear Information System (INIS)
A scheme of calculation of the electronic structure of a solid state surface and chemisorbed molecules is discussed. The method of the Green's function and MO LCAO approximation are used which permits to perform calculations, taking into account the whole crystal but not its fragment only, with the accuracy adopted by quantum chemistry. Results of model calculations are presented: chemisorption of hydrogen-like atom on the (100) face of the one-band crystal model and dispersion curves for the density of states of nickel (100) face. (Auth.)
Egawa, Toru; Kameyama, Akiyo; Takeuchi, Hiroshi
2006-08-01
The molecular structures of vanillin (4-hydroxy-3-methoxybenzaldehyde), isovanillin (3-hydroxy-4-methoxybenzaldehyde) and ethylvanillin (3-ethoxy-4-hydroxybenzaldehyde) were determined by means of gas electron diffraction. Among them, vanillin and ethylvanillin have a vanilla odor but isovanillin smells differently. The nozzle temperatures were 125, 173 and 146 °C, for vanillin, isovanillin and ethylvanillin, respectively. The results of MP2 and B3LYP calculations with the 6-31G** basis set were used as supporting information. The MP2 calculations predicted that vanillin and isovanillin have two stable conformers and ethylvanillin has four stable conformers. The electron diffraction data were found to be consistent with these conformational compositions. The determined structural parameters ( rg and ∠ α) of vanillin are as follows: =1.397(4) Å; r(C 1-C aldehyde)=1.471(←) Å; r(C 3-O Me)=1.374(9) Å; r(C 4-O H)=1.361(←) Å; r(O-C Me)=1.428(←) Å; r(C dbnd6 O)=1.214(8) Å; =1.110(11) Å; r(O-H)=0.991(←) Å; ∠C 6-C 1-C 2=120.6(2)°; ∠C 1-C 2-C 3=118.8(←)°; ∠C 1-C 6-C 5=120.1(←)°; ∠C 2-C 1-C aldehyde=122.7(18)°; ∠C 1-C dbnd6 O=119.4(16)°; ∠C 4-C 3-O Me=112.2(12)°; ∠C 3-C 4-O H=119.1(←)°; ∠C 3-O-C=121.7(29)°. Those of isovanillin are as follows: =1.402(4) Å; r(C 1-C aldehyde)=1.479(←) Å; r(C 4-O Me)=1.369(9) Å; r(C 3-O H)=1.357(←) Å; r(O-C Me)=1.422(←) Å; r(C dbnd6 O)=1.221(9) Å; =1.114(14) Å; r(O-H)=0.995(←) Å; ∠C 6-C 1-C 2=120.2(3)°; ∠C 1-C 2-C 3=119.0(←)°; ∠C 1-C 6-C 5=119.9(←)°; ∠C 2-C 1-C aldehyde=124.6(25)°; ∠C 1-C dbnd6 O=121.3(24)°; ∠C 3-C 4-O Me=114.4(12)°; ∠C 4-C 3-O H=121.2(←)°; ∠C 4-O-C=123.8(26)°. Those of ethylvanillin are as follows: =1.397(6) Å; r(C 1-C aldehyde)=1.471(←) Å; r(C 3-O Et)=1.365(13) Å; r(C 4-O H)=1.352(←) Å; r(O-C Et)=1.427(←) Å; r(C-C Et)=1.494(21) Å; r(C dbnd6 O)=1.206(9) Å; =1.109(10) Å; r(O-H)=0.990(←) Å; ∠C 6-C 1-C 2=120.2(3)°;
International Nuclear Information System (INIS)
We present an approach to solid-state electronic-structure calculations based on the finite-element method. In this method, the basis functions are strictly local, piecewise polynomials. Because the basis is composed of polynomials, the method is completely general and its convergence can be controlled systematically. Because the basis functions are strictly local in real space, the method allows for variable resolution in real space; produces sparse, structured matrices, enabling the effective use of iterative solution methods; and is well suited to parallel implementation. The method thus combines the significant advantages of both real-space-grid and basis-oriented approaches and so promises to be particularly well suited for large, accurate ab initio calculations. We develop the theory of our approach in detail, discuss advantages and disadvantages, and report initial results, including electronic band structures and details of the convergence of the method. copyright 1999 The American Physical Society
First-principle calculations on the structural and electronic properties of hard C{sub 11}N{sub 4}
Energy Technology Data Exchange (ETDEWEB)
Li, Dongxu, E-mail: lidongxu@hqu.edu.cn [College of Materials Science and Engineering, Huaqiao University, Xiamen 361021 (China); Shi, Jiancheng; Lai, Mengling; Li, Rongkai [College of Materials Science and Engineering, Huaqiao University, Xiamen 361021 (China); Yu, Dongli [State Key Laboratory of Metastable Materials Science and Technology, Yanshan University, Qinhuangdao 066004 (China)
2014-09-15
A graphite-like C{sub 11}N{sub 4} model was built by stacking graphene and a C{sub 3}N{sub 4} triazine layer and simulated by first principle calculations, which transfers to a diamond-like structure under high pressure. The structural, mechanical, and electronic properties of both materials were calculated. The elastic constants of both materials satisfy the Born-criterion. Furthermore, no imaginary frequencies were observed in phonon calculations. The diamond-like C{sub 11}N{sub 4} is semiconducting and consists of polyhedral and hollow C–N cages. The Vickers hardness of diamond-like C{sub 11}N{sub 4} was calculated to be 58 GPa. The phase transformation from graphite-like to diamond-like C{sub 11}N{sub 4} is proposed to occur at approximately 27.2 GPa based on the pressure-dependent enthalpy.
Term structure of 4d-electron configurations and calculated spectrum in Sn-isonuclear sequence
International Nuclear Information System (INIS)
Theoretical calculations of term structure are carried out for the ground configurations 4dw, of atomic ions in the Sn isonuclear sequence. Atomic computations are performed to give a detailed account of the transitions in Sn+6 to Sn+13 ions. The spectrum is calculated for the most important excited configurations 4p5 4dn+1, 4dn-1 4f1, and 4dn-1 5p1 with respect to the ground configuration 4dn, with n=8-1, respectively. The importance of 4p-4d, 4d-4f, and 4d-5p transitions is stressed, as well as the need for the configuration-interaction CI treatment of the Δn=0 transitions. In the region of importance for extreme ultraviolet (EUV) lithography around 13.4nm, the strongest lines were expected to be 4dn-4p5 4dn+1 and 4dn-4dn-1 4f1
Jónsson, Elvar Ö; Puska, Martti; Jónsson, Hannes
2016-01-01
An implementation of the generalized Pipek-Mezey method [Lehtola, S.; J\\'onsson, H. J. Chem. Theory Comput. 2014, 10, 642] for generating localized orbitals in periodic systems, i.e. Wannier functions, is described. The projector augmented wave (PAW) formalism for the representation of atomic core electrons is included in the implementation, which has been developed within the atomic simulation environment (ASE) software library. The implementation supports several different kinds of representations for the wave function, including real-space grids, plane waves or a linear combination of atomic orbitals. The implementation is tailored to the GPAW program but can easily be adapted to use output from various other electronic structure software packages such as ABINIT, NWChem, or VASP through interfaces in ASE. Generalized Pipek-Mezey Wannier functions (PMWF) are presented for both isolated molecules, as well as systems with periodicity in one, two and three dimensions. The method gives a set of highly localized...
Ogitsu, Tadashi; Gygi, Francois; Reed, John; Schwegler, Eric; Galli, Giulia
2007-03-01
Boron exhibits the most complex structure of all elemental solids, with more than 300 atoms per unit cell arranged in interconnecting icosahedra, and some crystallographic positions occupied with a probability of less than one. The precise determination of the ground state geometry of boron---the so-called β-boron structure--has been elusive and its electronic and bonding properties have been difficult to rationalize. Using lattice model Monte Carlo optimization techniques and ab-initio simulations, we have shown that a defective, quasi-ordered β solid is the most stable structure at zero as well as finite T. In the absence of partially occupied sites (POS), the perfect β-boron crystal is unstable; the presence of POS lower its internal energy below that of an ordered α-phase, not mere an entropic effect. We present a picture of the intricate and unique bonding in boron based on maximally localized Wannier (MLWF) functions, which indicates that the presence of POS provides a subtle, yet essential spatial balance between electron deficient and fully saturated bonds. This work was performed under the auspices of the U.S. Dept. of Energy at the University of California/ LLNL under contract no. W-7405-Eng-48.
Sun, Feng; Wang, Li; Stoumpos, Constantinos C.
2016-08-01
The synthesis, structure, and characterization of a new centrosymmetric borate Pb2O[BO2(OH)] based on anion-centered OPb4 tetrahedra are reported. Pb2O[BO2(OH)] crystallizes in monoclinic space group C2/m with a=12.725(7) Å, b=5.698(3) Å, c=7.344(4) Å, β=116.277(6)°. The electronic band structure and density of states of Pb2O[BO2(OH)] have been calculated via the density functional theory (DFT). Electron density difference calculation indicates that lone-pair electrons of Pb2+ cation should be stereoactive.
Many-body electronic structure calculations of Eu-doped ZnO
Lorke, M.; Frauenheim, T.; da Rosa, A. L.
2016-03-01
The formation energies and electronic structure of europium-doped zinc oxide has been determined using DFT and many-body G W methods. In the absence of intrisic defects, we find that the europium-f states are located in the ZnO band gap with europium possessing a formal charge of 2+. On the other hand, the presence of intrinsic defects in ZnO allows intraband f -f transitions otherwise forbidden in atomic europium. This result corroborates with recently observed photoluminescence in the visible red region S. Geburt et al. [Nano Lett. 14, 4523 (2014), 10.1021/nl5015553].
International Nuclear Information System (INIS)
First-principles calculations of structural, electronic, optical, elastic, mechanical properties, and Born effective charges of monoclinic HfO2 are performed with the plane-wave pseudopotential technique based on the density-functional theory. The calculated structural properties are consistent with the previous theoretical and experimental results. The electronic structure reveals that monoclinic HfO2 has an indirect band gap. The analyses of density of states and Mulliken charges show mainly covalent nature in Hf-O bonds. Optical properties, including the dielectric function, refractive index, extinction coefficient, reflectivity, absorption coefficient, loss function, and optical conductivity each as a function of photon energy are calculated and show an optical anisotropy. Moreover, the independent elastic constants, bulk modulus, shear modulus, Young's modulus, Poisson's ratio, compressibility, Lamé constant, sound velocity, Debye temperature, and Born effective charges of monoclinic HfO2 are obtained, which may help to understand monoclinic HfO2 for future work. (condensed matter: electronic structure, electrical, magnetic, and optical properties)
Sarkar, Kanchan; Sharma, Rahul; Bhattacharyya, S P
2010-03-01
A density matrix based soft-computing solution to the quantum mechanical problem of computing the molecular electronic structure of fairly long polythiophene (PT) chains is proposed. The soft-computing solution is based on a "random mutation hill climbing" scheme which is modified by blending it with a deterministic method based on a trial single-particle density matrix [P((0))(R)] for the guessed structural parameters (R), which is allowed to evolve under a unitary transformation generated by the Hamiltonian H(R). The Hamiltonian itself changes as the geometrical parameters (R) defining the polythiophene chain undergo mutation. The scale (λ) of the transformation is optimized by making the energy [E(λ)] stationary with respect to λ. The robustness and the performance levels of variants of the algorithm are analyzed and compared with those of other derivative free methods. The method is further tested successfully with optimization of the geometry of bipolaron-doped long PT chains. PMID:26613302
A Initio Lcao Electronic Structure Calculations of Layered Transition Metal Compounds.
Dawson, William G.
1987-09-01
Available from UMI in association with The British Library. In this work the electronic structure of three systems of layered transition metal compounds are examined using an ab initio tight binding (LCAO) method using the Xalpha exchange/correlation approximation: group VI ditellurides, group IV trichalcogenides and quaternary copper oxide defect-perovskites. A chemical pseudopotential argument is presented in order to justify the use of a small basis set of atomic orbitals. The group VI transition metal compounds MoTe_2 and WTe _2 show strong metal-metal interactions and MoTe_2 undergoes an unusual phase transition with the lattice parameter perpendicular to the layers decreasing with increasing temperature. The group IV transition metal trichalcogenides provide a useful series for study due to their quasi-1-dimensional character and the occurrence of two closely related structural variants. The atypical compound ZrTe_3 is given special attention because of its apparent semimetallic nature. The final group of compounds studied are the high Tc superconducting ceramics Ba-La-Cu-O and Ba-Y-Cu-O. The technological importance of compounds with zero resistance and showing the Meissner effect (expelling magnetic fields) above liquid nitrogen temperatures and the, as yet, undefined nature of the mechanism of superconductivity stresses the need to carefully examine the electronic structure of these materials. The role of oxygen vacancies, the charge state of the copper ions and the possibility of structural phase transitions are some of the topics considered here. The use of an atomic-orbital basis allows a comparatively straightforward description of the chemical bonding in a crystal--especially useful when the unit cell contains a large number of atoms.
Institute of Scientific and Technical Information of China (English)
无
2009-01-01
Calculations have been made for single-walled zigzag(n,0) carbon nanotubes containing substitutional boron impurity atoms using ab initio density functional theory.It is found that the formation energies of these nanotubes depend on the tube diameter,as do the electronic properties,and show periodic fea-ture that results from their different π bonding structures compared to those of perfect zigzag carbon nanotubes.When more boron atoms are incorporated into a single-walled zigzag carbon nanotube,the substitutional boron atoms tend to come together to form structure of BC3 nanodomains,and B-doped tubes have striking acceptor states above the top of the valence bands.For the structure of BC3,there are two kinds of configurations with different electronic structures.
Czech Academy of Sciences Publication Activity Database
Kuneš, Jan; Novák, Pavel; Schmid, R.; Blaha, P.; Schwarz, K.
2001-01-01
Roč. 64, - (2001), s. 153102-1-153102-3. ISSN 0163-1829 Grant ostatní: CZ-AT(XX) project No.KONTAKT 1999/21 Institutional research plan: CEZ:A02/98:Z1-010-914 Keywords : actinide compounds * ab initio electronic structure calculations Subject RIV: BM - Solid Matter Physics ; Magnetism Impact factor: 3.070, year: 2001
Mn-doped Ga(As,P) and (Al,Ga)As ferromagnetic semiconductors: electronic structure calculations
Czech Academy of Sciences Publication Activity Database
Mašek, Jan; Kudrnovský, Josef; Máca, František; Sinova, J.; MacDonald, A. H.; Champion, R.P.; Gallagher, B. L.; Jungwirth, Tomáš
2007-01-01
Roč. 75, č. 4 (2007), 045202/1-045202/6. ISSN 1098-0121 R&D Projects: GA ČR GA202/05/0575; GA ČR GA202/04/0583 Institutional research plan: CEZ:AV0Z10100521; CEZ:AV0Z10100520 Keywords : ferromagnetic semiconductors * electronic structure calculations Subject RIV: BM - Solid Matter Physics ; Magnetism Impact factor: 3.172, year: 2007
Energy Technology Data Exchange (ETDEWEB)
Tohme, Samir N.; Korek, Mahmoud, E-mail: mahmoud.korek@bau.edu.lb, E-mail: fkorek@yahoo.com; Awad, Ramadan [Faculty of Science, Beirut Arab University, P.O. Box 11-5020 Riad El Solh, Beirut 1107 2809 (Lebanon)
2015-03-21
Ab initio techniques have been applied to investigate the electronic structure of the LiYb molecule. The potential energy curves have been computed in the Born–Oppenheimer approximation for the ground and 29 low-lying doublet and quartet excited electronic states. Complete active space self-consistent field, multi-reference configuration interaction, and Rayleigh Schrödinger perturbation theory to second order calculations have been utilized to investigate these states. The spectroscopic constants, ω{sub e}, R{sub e}, B{sub e}, …, and the static dipole moment, μ, have been investigated by using the two different techniques of calculation with five different types of basis. The eigenvalues, E{sub v}, the rotational constant, B{sub v}, the centrifugal distortion constant, D{sub v}, and the abscissas of the turning points, R{sub min} and R{sub max}, have been calculated by using the canonical functions approach. The comparison between the values of the present work, calculated by different techniques, and those available in the literature for several electronic states shows a very good agreement. Twenty-one new electronic states have been studied here for the first time.
Tohme, Samir N.; Korek, Mahmoud; Awad, Ramadan
2015-03-01
Ab initio techniques have been applied to investigate the electronic structure of the LiYb molecule. The potential energy curves have been computed in the Born-Oppenheimer approximation for the ground and 29 low-lying doublet and quartet excited electronic states. Complete active space self-consistent field, multi-reference configuration interaction, and Rayleigh Schrödinger perturbation theory to second order calculations have been utilized to investigate these states. The spectroscopic constants, ωe, Re, Be, …, and the static dipole moment, μ, have been investigated by using the two different techniques of calculation with five different types of basis. The eigenvalues, Ev, the rotational constant, Bv, the centrifugal distortion constant, Dv, and the abscissas of the turning points, Rmin and Rmax, have been calculated by using the canonical functions approach. The comparison between the values of the present work, calculated by different techniques, and those available in the literature for several electronic states shows a very good agreement. Twenty-one new electronic states have been studied here for the first time.
Pan, Yong; Guan, Weiming
2016-09-01
MoS3 has attracted considerable attention as potential hydrogen storage material due to the interaction between the hydrogen and unsaturated sulfur atoms. However, its structure and physical properties are unknown. By means of first-principles approach and Inorganic crystal structure Database (ISCD), we systematically investigated the structure, relevant physical and thermodynamic properties of MoS3. Phonon dispersion, electronic structure, band structure and heat capacity are calculated in detail. We predicted the orthorhombic B2ab (SrS3-type) and tetragonal P-421m (BaS3-type) structures of MoS3, which prefers to form the SrS3-type (Space group: B2ab, No.41) structure at the ground state. High pressure results in structural transition from SrS3-type structure to BaS3-type structure. This sulfide exhibits a degree of metallic behavior. The calculated heat capacity of MoS3 with SrS3-type structure is about of 39 J/(mol·K).
Galvan, D H
2003-01-01
To get insight into the electronic properties of PrFe4P12 skutterudite, band electronic structure calculations, Total and Projected Density of States, Crystal Orbital Overlap Population and Mulliken Population Analysis were performed. The energy bands yield a semi metallic behavior with a direct gap (at gamma) of 0.02 eV. Total and Projected Density of States provided information of the contribution from each orbital of each atom to the total Density of States. Moreover, the bonding strength between some atoms within the unit cell was obtained. Mulliken Population analysis suggests ionic behavior for this compound.
DEFF Research Database (Denmark)
Enkovaara, J.; Rostgaard, Carsten; Mortensen, Jens Jørgen;
2010-01-01
Electronic structure calculations have become an indispensable tool in many areas of materials science and quantum chemistry. Even though the Kohn-Sham formulation of the density-functional theory (DFT) simplifies the many-body problem significantly, one is still confronted with several numerical...... challenges. In this article we present the projector augmented-wave (PAW) method as implemented in the GPAW program package (https://wiki.fysik.dtu.dk/gpaw) using a uniform real-space grid representation of the electronic wavefunctions. Compared to more traditional plane wave or localized basis set...
Electronic Structure of Organic/Inorganic Interfaces: Insights from First Principles Calculations
Segev, Lior
Electronic devices based on molecules draw a lot of attention in both scientific and industrial activities. Molecules in electronic devices can serve as the heart of the device, featuring versatile physical properties i.e. electronical, optical, magnetic, etc. Molecules can also function as an assist mechanism in which the electronic properties of the underlying material are modified in a predictable fashion according to the molecular monolayer properties. But, the route to applications in both these directions lies in answering fundamental questions related to band offsets between two materials, full electronic structure determination of molecule and substrates, work function modifications, etc. To tackle these questions, we chose to study the interface formed by an alkyl monolayer adsorbed on a Si substrate by utilizing two ab initio methods. First, the density functional theory (DFT) utilizing the local density or the B3LYP approximations for the exchange-correlation potential and, second, the many-body perturbation theory based on the GW approximation. We adapted a "divide and conquer" approach to our system by simulating the infinite counterpart, polyethylene, of our finite alkyl chain to test how the band gap of the two molecules changes when moving from an infinite 1D molecule to a finite length molecule. We find excellent agreement between our GW simulation results for polyethylene and experimental results for the bandstructure, ionization potential and band gap values. From DFT simulations, we analyze the ultra-violet photoelectron spectra (UPS) of odd and even number of carbons alkyl chains and identify the origin of their differences in spectral signature. GW simulations of the full alkyl monolayer/Si(111) system reveal that the projected density of states (DOS) of the upper alkyl chain have an excellent agreement to experimental UPS and inverse-photoemission spectra results. Based on this correspondence, we find the band alignment between the alkyl
First-principles calculations of structural, elastic and electronic properties of Li2B12H12
International Nuclear Information System (INIS)
Highlights: • The fundamental structural parameters and density of states of Li2B12H12 are calculated. • Elastic constants are obtained by the strain energy–strain curves method. • Polycrystalline elastic moduli, Debye temperature and the average elastic wave velocity are determined. • The mechanical stability and elastic anisotropy are analyzed. - Abstract: We investigate the structural, elastic and electronic properties of Li2B12H12 using the first-principles method. Our calculations show that the lowest energy structure of Li2B12H12 is monoclinic C2/m type. We take the monoclinic C2/m Li2B12H12 as a representative to carry out the corresponding theoretical studies. The independent elastic constants are successfully obtained from the strain energy–strain curve calculations. The Shear and Young‘s moduli, as well as Poisson‘s ratio for ideal polycrystalline Li2B12H12 are calculated. The shear anisotropic factors and elastic anisotropy of Li2B12H12 are analyzed. The Debye temperature and the average elastic wave velocity are derived from theoretical elastic constants. According to the obtained results, the monoclinic C2/m Li2B12H12 is found to be mechanically stable and brittle at zero temperature and zero pressure. Furthermore, the density of states and electron charge density distributions are studied. The insulator Li2B12H12 is a technologically interesting indirect hydrogen storage material for further studies
International Nuclear Information System (INIS)
Graphical abstract: UV photoelectron spectrum of the gas-phase thermal decomposition of 5-methyltetrazole (5MTZ), obtained at 195 °C, mechanism and potential energy diagram underlying the thermal dissociation of 5MTZ. Highlights: ► Gas-phase 5-methyltetrazole exists mainly as the 2H-tautomer. ► Thermal decomposition of 5MTZ gives N2, CH3CN and HCN, at 250 °C. ► HCN is formed from secondary reactions. - Abstract: The electronic properties and thermal decomposition of 5-methyltetrazole (5MTZ) are investigated using UV photoelectron spectroscopy (UVPES) and theoretical calculations. Simulated spectra of both 1H- and 2H-5MTZ, based on electron propagator methods, are produced in order to study the relative tautomer population. The thermal decomposition results are rationalized in terms of G2(MP2) results. 5MTZ yields a HOMO ionization energy of 10.82 ± 0.04 eV and the gas-phase 5MTZ assumes predominantly the 2H-form. Its gas-phase thermal decomposition starts at ca. 195 °C and leads to the formation of N2,CH3CN and HCN. N2 is formed from two competing routes, involving 150.2 and 126.2 kJ/mol energy barriers, from 2H- and 1H-5MTZ, respectively. CH3CN is formed also from two competing pathways, requiring activation energies of 218.3 (2H-5MTZ) and 198.6 kJ/mol (1H-5MTZ). Conclusions are also drawn in order to explain the formation of HCN from secondary reactions in the thermal decomposition process.
Institute of Scientific and Technical Information of China (English)
LU Lai-Yu; WEI Dong-Qing; CHEN Xiang-Rong; JI Guang-Fu
2008-01-01
Structures and electronic properties of the pentaerythritol (PE) crystal under volume compression up to 0.85Vo are studied by E - V fitting method using density functional theory (DFT). The compression dependences of the cell volumes, lattice constants, and molecular geometries of solid PE are presented and discussed. It is found that the solid PE presents anisotropy along a- and c-axes, and the c axis is the most compressible. Decreasing anisotropy ratio (c/a) with elevating compression suggests an enhancement of the vdW interaction with increasing compression. The C-C and C-H bonds are significantly reduced under compression, which may be related to the sensitivity. The solid PE has indirect band gap (X - C) in the range of the researched compression and the band gap is decreased with compression.
First-Principles Calculations of Structural, Electronic and Optical Properties of CaTiO3 Crystal
Medeiros, Subênia; Silva, Jusciane; Albuquerque, Eudenilson; Freire, Valder
2013-03-01
The structural, electronic, vibrational, and optical properties of perovskite CaTiO3 in the cubic, orthorhombic, and tetragonal phase are calculated in the framework of density functional theory (DFT) with different exchange-correlation potentials by CASTEP package. The calculated band structure shows an indirect band gap of 1.88 eV at the Γ-R points in the Brillouin zone to the cubic structure, a direct band gap of 2.41 eV at the Γ - Γ points to the orthorhombic structure, and an indirect band gap of 2.31 eV at the M' Γ points to the tetragonal phase. I have concluded that the bonding between Ca and TiO2 is mainly ionic and that the TiO2 entities bond covalently. Unlike some perovskites the CaTiO3 does not exhibit a ferroelectric phase transition down to 4.2 K. It is still known that the CaTiO3 has a static dielectric constant that extrapolates to a value greater than 300 at zero temperature. Our calculated lattice parameters, elastic constants, optical properties, and vibrational frequencies are found to be in good agreement with the available theoretical and experimental values. The results for the effective mass in the electron and hole carriers are also presented in this work.
International Nuclear Information System (INIS)
The structural, elastic, electronic, and thermodynamic properties of ZrxNb1−xC alloys are investigated using the first principles method based on the density functional theory. The results show that the structural properties of ZrxNb1−xC alloys vary continuously with the increase of Zr composition. The alloy possesses both the highest shear modulus (215 GPa) and a higher bulk modulus (294 GPa), with a Zr composition of 0.21. Meanwhile, the Zr0.21 Nb0.79C alloy shows metallic conductivity based on the analysis of the density of states. In addition, the thermodynamic stability of the designed alloys is estimated using the calculated enthalpy of mixing. (condensed matter: electronic structure, electrical, magnetic, and optical properties)
Kong, Bo; Zhang, Yachao
2016-07-01
The electronic structures of the cubic GdH3 are extensively investigated using the ab initio many-body GW calculations treating the Gd 4f electrons either in the core (4f-core) or in the valence states (4f-val). Different degrees of quasiparticle (QP) self-consistent calculations with the different starting points are used to correct the failures of the GGA/GGA + U/HSE03 calculations. In the 4f-core case, GGA + G0W0 calculations give a fundamental band gap of 1.72 eV, while GGA+ GW0 or GGA + GW calculations present a larger band gap. In the 4f-val case, the nonlocal exchange-correlation (xc) functional HSE03 can account much better for the strong localization of the 4f states than the semilocal or Hubbard U corrected xc functional in the Kohn-Sham equation. We show that the fundamental gap of the antiferromagnetic (AFM) or ferromagnetic (FM) GdH3 can be opened up by solving the QP equation with improved starting point of eigenvalues and wave functions given by HSE03. The HSE03 + G0W0 calculations present a fundamental band gap of 2.73 eV in the AFM configuration, and the results of the corresponding GW0 and GW calculations are 2.89 and 3.03 eV, respectively. In general, for the cubic structure, the fundamental gap from G0W0 calculations in the 4f-core case is the closest to the real result. By G0W0 calculations in the 4f-core case, we find that H or Gd defects can strongly affect the band structure, especially the H defects. We explain the mechanism in terms of the possible electron correlation on the hydrogen site. Under compression, the insulator-to-metal transition in the cubic GdH3 occurs around 40 GPa, which might be a satisfied prediction.
International Nuclear Information System (INIS)
The structural, electronic, and optical properties of binary ZnO, ZnSe compounds, and their ternary ZnO1−xSex alloys are computed using the accurate full potential linearized augmented plane wave plus local orbital (FP-LAPW + lo) method in the rocksalt (B1) and zincblende (B3) crystallographic phases. The electronic band structures, fundamental energy band gaps, and densities of states for ZnO1−xSex are evaluated in the range 0 ≤ x ≤ 1 using Wu—Cohen (WC) generalized gradient approximation (GGA) for the exchange—correlation potential. Our calculated results of lattice parameters and bulk modulus reveal a nonlinear variation for pseudo-binary and their ternary alloys in both phases and show a considerable deviation from Vegard's law. It is observed that the predicted lattice parameter and bulk modulus are in good agreement with the available experimental and theoretical data. We establish that the composition dependence of band gap is semi-metallic in B1 phase, while a direct band gap is observed in B3 phase. The calculated density of states is described by taking into account the contribution of Zn 3d, O 2p, and Se 4s, and the optical properties are studied in terms of dielectric functions, refractive index, reflectivity, and energy loss function for the B3 phase and are compared with the available experimental data. (condensed matter: electronic structure, electrical, magnetic, and optical properties)
Baldea, Ioan
2012-01-01
In cases where reorganization is important, present theoretical studies of molecular transport have inherently to resort to models. The Newns-Anderson model is ubiquitous for this purpose but, to author's knowledge, attempts to validate/challenge this model by microscopic calculations are missing in the literature. In this work, results of electronic structure calculations are presented, which demonstrate that the conventional Newns-Anderson model fails to describe redox-active tunneling junctions of recent experimental interest. For the case considered, the ($4, 4^\\prime$)-bipyridine molecule, the failure traces back to the floppy degree of freedom represented by the relative rotation of the two pyridine rings. Expressions that generalize the Newns-Anderson model are deduced, which include significant anharmonicities. These expressions can be straightforwardly utilized as input information in calculations of the partially coherent transport.
Electronic structure of ScN and YN:density-functional theory LDA and GW approximation calculations
Institute of Scientific and Technical Information of China (English)
Lü Tie-Yu; Huang Mei-Chun
2007-01-01
The desirable physical properties of hardness, high temperature stability, and conductivity make the early transition metal nitrides important materials for various technological applications. To learn more about the nature of these materials, the local-density approximation(LDA) and GW approximation i.e. combination of the Green function G and the screened Coulomb interaction W, have been performed. This paper investigates the bulk electronic and physical properties of early transition metal mononitrides, ScN and YN in the rocksalt structure. In this paper, the semicore electrons are regarded as valance electrons. ScN appears to be a semimetal, and YN is semiconductor with band gap of0.142 eV within the LDA, but are in fact semiconductors with indirect band gaps of 1.244 and 0.544 eV respectively, as revealed by calculations performed using GW approximation.
Atomic and electronic structure of hydrogen on ZnO (1bar 100) surface: ab initio hybrid calculations
Usseinov, A. B.; Kotomin, E. A.; Zhukovskii, Yu F.; Purans, J.; Sorokin, A. V.; Akilbekov, A. T.
2013-12-01
Hydrogen atoms unavoidably incorporated into ZnO during growth of bulk samples and thin films considerably affect their electrical conductivity. The results of first principles hybrid LCAO calculations are discussed for hydrogen atoms in the bulk and on the non-polar ZnO (1bar 100) surface. The incorporation energy, the atomic relaxation, the electronic density redistribution and the electronic structure modifications are compared for the surface adsorption and bulk interstitial H positions. It is shown that hydrogen has a strong binding with the surface O ions (2.7 eV) whereas its incorporation into bulk is energetically unfavorable. Surface hydrogen atoms are very shallow donors, thus, contributing to the electronic conductivity.
International Nuclear Information System (INIS)
The focus of this thesis is the study of the electronic and magnetic structure of three representatives of the main Fe-bearing rock-forming silicates: Fe2+2Si2O6, almandine Fe2+3Al2(SiO4)3 and andradite Ca3Fe3+2(SiO4)3. For this purpose the quantum mechanical first principles electronic structure calculations are performed by the most efficient DFT method in the local spin-density approximation for calculating spectroscopic data: the spin-polarized self consistent charge X[alpha] method. These minerals have attracted significant attention due to their abundance in the Earth's crust and mantle, and because crystallised silicates are main components of cosmic dust which is the most abundant raw material in the Universe. The specific feature and strength of these investigations consist in the theoretical characterization of these complex systems based on experimental results. This means that, on one hand, experimental spectroscopic and crystallographic data are being used to judge the reliability of the calculations, whereas, on the other hand, experimental data are interpreted and explained by the theoretical results. This work is divided into seven main parts. Chapter 1 is the introduction to the thesis. Chapter 2 describes the theoretical bases, ideas, approximations and advantages of the SCC- X[alpha] method and basics of the art of cluster construction. Chapter 3 considers physical bases of absorption and Moessbauer spectroscopy, crystal field theory, evaluation of the main spectroscopic values within the frames of the SCC- X[alpha] method and magnetic interaction between atoms. In addition, tetragonally, trigonally and angularly distorted octahedral sites with various degrees of the distortions are calculated and analyzed. The electronic and magnetic structures of orthoferrosilite, almandine and andradite are described in Chapters 4, 5 and 6, respectively. In the case of orthoferrosilite the magnetic interactions between the iron spins within the ribbons and
Long, E. R., Jr.
1979-01-01
The Bethe-Bloch stopping power relations for inelastic collisions were used to determine the absorption of electron and proton energy in cured neat epoxy resin and the absorption of electron energy in a graphite/epoxy composite. Absorption of electron energy due to bremsstrahlung was determined. Electron energies from 0.2 to 4.0 MeV and proton energies from 0.3 to 1.75 MeV were used. Monoenergetic electron energy absorption profiles for models of pure graphite, cured neat epoxy resin, and graphite/epoxy composites are reported. A relation is determined for depth of uniform energy absorption in a composite as a function of fiber volume fraction and initial electron energy. Monoenergetic proton energy absorption profiles are reported for the neat resin model. A relation for total proton penetration in the epoxy resin as a function of initial proton energy is determined. Electron energy absorption in the composite due to bremsstrahlung is reported. Electron and proton energy absorption profiles in cured neat epoxy resin are reported for environments approximating geosynchronous earth orbit.
First-Principles Band Calculations on Electronic Structures of Ag-Doped Rutile and Anatase TiO2
Institute of Scientific and Technical Information of China (English)
HOU Xing-Gang; LIU An-Dong; HUANG Mei-Dong; LIAO Bin; WU Xiao-Ling
2009-01-01
The electronic structures of Ag-doped rutile and anatase TiO2 are studied by first-principles band calculations based on density funetionai theory with the full-potentiai linearized-augraented-plane-wave method.New occupied bands ore found between the band gaps of both Ag-doped rutile and anatase TiO2.The formation of these new bands Capri be explained mainly by their orbitals of Ag 4d states mixed with Ti 3d states and are supposed to contribute to their visible light absorption.
International Nuclear Information System (INIS)
First-principles calculations have been used to study the effect of vacancies on the structural and electronic properties in substoichiometric TiCx and TiNx. The effect of vacancies on equilibrium volumes, bulk moduli, electronic band structures and density of states of the substoichiometric phases was studied using a full-potential linear augmented plane-wave method. A model structure of eight-atom supercells with ordered vacancies within the carbon and nitrogen sublattices is used. We find that the lattice parameters of the studied stoichiometries in both TiCx and TiNx are smaller than that of ideal stoichiometric TiC and TiN. Our results for the variation of the lattice parameters and the bulk moduli for TiCx are found to be in good agreement with experiment. The variation of the energy gaps with the atomic concentration ratio shows that these compounds present the same trends. Results for TiCx are compared to a recent full-potential calculation with relaxed 16-atom supercells
Hamioud, L.; Boumaza, A.; Touam, S.; Meradji, H.; Ghemid, S.; El Haj Hassan, F.; Khenata, R.; Omran, S. Bin
2016-06-01
The present paper aims to study the structural, electronic, optical and thermal properties of the boron nitride (BN) and BAs bulk materials as well as the BNxAs1-x ternary alloys by employing the full-potential-linearised augmented plane wave method within the density functional theory. The structural properties are determined using the Wu-Cohen generalised gradient approximation that is based on the optimisation of the total energy. For band structure calculations, both the Wu-Cohen generalised gradient approximation and the modified Becke-Johnson of the exchange-correlation energy and potential, respectively, are used. We investigated the effect of composition on the lattice constants, bulk modulus and band gap. Deviations of the lattice constants and the bulk modulus from the Vegard's law and the linear concentration dependence, respectively, were observed for the alloys where this result allows us to explain some specific behaviours in the electronic properties of the alloys. For the optical properties, the calculated refractive indices and the optical dielectric constants were found to vary nonlinearly with the N composition. Finally, the thermal effect on some of the macroscopic properties was predicted using the quasi-harmonic Debye model in which the lattice vibrations are taken into account.
Acceleration Schemes for Ab-Initio Molecular Dynamics and Electronic Structure Calculations
Tassone, Francesco; Mauri, Francesco; Car, Roberto
1994-01-01
We study the convergence and the stability of fictitious dynamical methods for electrons. First, we show that a particular damped second-order dynamics has a much faster rate of convergence to the ground-state than first-order steepest descent algorithms while retaining their numerical cost per time step. Our damped dynamics has efficiency comparable to that of conjugate gradient methods in typical electronic minimization problems. Then, we analyse the factors that limit the size of the integ...
Jiang, Hao; Cao, Guanghan; Cao, Chao
2015-01-01
The electronic structure of quasi-one-dimensional superconductor K2Cr3As3 is studied through systematic first-principles calculations. The ground state of K2Cr3As3 is paramagnetic. Close to the Fermi level, the Cr-3dz(2), dxy, and d(x(2)-y(2)) orbitals dominate the electronic states, and three bands cross EF to form one 3D Fermi surface sheet and two quasi-1D sheets. The electronic DOS at EF is less than 1/3 of the experimental value, indicating a large electron renormalization factor around EF. Despite of the relatively small atomic numbers, the antisymmetric spin-orbit coupling splitting is sizable (≈60 meV) on the 3D Fermi surface sheet as well as on one of the quasi-1D sheets. Finally, the imaginary part of bare electron susceptibility shows large peaks at Γ, suggesting the presence of large ferromagnetic spin fluctuation in the compound. PMID:26525099
Electronic structure calculations of positron lifetimes in nuclear materials: SiC and UO2
International Nuclear Information System (INIS)
We present first-principles calculations of positron lifetimes of vacancy-type defects in two nuclear materials: SiC and UO2. We use a self-consistent positron lifetime calculation scheme based on the two-component density functional theory (DFT). Full defect relaxation due to both the creation of the vacancy and the presence of the positron was taken into account. Our results for SiC differ strongly from those published in literature up to now [G. Brauer et al. Phys. Rev. B 54, 2512 (1996)]. This is mostly due to the effect of the relaxation, that was not taken into account before. We also present the first calculated positron lifetimes obtained for UO2 in the DFT+U approach. Results are compared with experimental data. (authors)
Zhang, Yu; Tang, Fu-Ling; Xue, Hong-Tao; Lu, Wen-Jiang; Liu, Jiang-Fei; Huang, Min
2015-02-01
Using first-principles plane-wave calculations within density functional theory, we theoretically studied the atomic structure, bonding energy and electronic properties of the perfect Mo (110)/MoSe2 (100) interface with a lattice mismatch less than 4.2%. Compared with the perfect structure, the interface is somewhat relaxed, and its atomic positions and bond lengths change slightly. The calculated interface bonding energy is about -1.2 J/m2, indicating that this interface is very stable. The MoSe2 layer on the interface has some interface states near the Fermi level, the interface states are mainly caused by Mo 4d orbitals, while the Se atom almost have no contribution. On the interface, Mo-5s and Se-4p orbitals hybridize at about -6.5 to -5.0 eV, and Mo-4d and Se-4p orbitals hybridize at about -5.0 to -1.0 eV. These hybridizations greatly improve the bonding ability of Mo and Se atom in the interface. By Bader charge analysis, we find electron redistribution near the interface which promotes the bonding of the Mo and MoSe2 layer.
Institute of Scientific and Technical Information of China (English)
Muhammad Rashid; Fayyaz Hussain; Muhammad Imran; S A Ahmad; N A Noor; M U Sohaib; S M Alay-e-Abbas
2013-01-01
The structural,electronic,and optical properties of binary ZnO,ZnSe compounds,and their ternary ZnO1-xSex alloys are computed using the accurate full potential linearized augmented plane wave plus local orbital (FP-LAPW + lo) method in the rocksalt (B 1) and zincblende (B3) crystallographic phases.The electronic band structures,fundamental energy band gaps,and densities of states for ZnO1 xSex are evaluated in the range 0 ≤ x ≤ 1 using Wu-Cohen (WC) generalized gradient approximation (GGA) for the exchange-correlation potential.Our calculated results of lattice parameters and bulk modulus reveal a nonlinear variation for pseudo-binary and their ternary alloys in both phases and show a considerable deviation from Vegard's law.It is observed that the predicted lattice parameter and bulk modulus are in good agreement with the available experimental and theoretical data.We establish that the composition dependence of band gap is semi-metallic in B1 phase,while a direct band gap is observed in B3 phase.The calculated density of states is described by taking into account the contribution of Zn 3d,O 2p,and Se 4s,and the optical properties are studied in terms of dielectric functions,refractive index,reflectivity,and energy loss function for the B3 phase and are compared with the available experimental data.
Harish, R. Sugan; Jayalakshmi, D. S.; Viswanathan, E.; Sundareswari, M.
2016-05-01
The mechanical, electronic, thermodynamic properties and structural stability of tetragonal structured CaNi2P2 and CaNi2Sb2 intermetallic compounds has been studied using the FP-LAPW method based on density functional theory. The PBE-GGA exchange correlation has been applied. Using the computed elastic constants, various elastic moduli such as bulk, shear, Young’s modulus, Poisson’s ratio and anisotropy constant are calculated and discussed. Stability of the compounds is confirmed by using their elastic constants. Pugh’s ratio is calculated to analyze the mechanical nature of the compound.
Energy Technology Data Exchange (ETDEWEB)
Tucker, Jon R.; Magyar, Rudolph J.
2012-02-01
High explosives are an important class of energetic materials used in many weapons applications. Even with modern computers, the simulation of the dynamic chemical reactions and energy release is exceedingly challenging. While the scale of the detonation process may be macroscopic, the dynamic bond breaking responsible for the explosive release of energy is fundamentally quantum mechanical. Thus, any method that does not adequately describe bonding is destined to lack predictive capability on some level. Performing quantum mechanics calculations on systems with more than dozens of atoms is a gargantuan task, and severe approximation schemes must be employed in practical calculations. We have developed and tested a divide and conquer (DnC) scheme to obtain total energies, forces, and harmonic frequencies within semi-empirical quantum mechanics. The method is intended as an approximate but faster solution to the full problem and is possible due to the sparsity of the density matrix in many applications. The resulting total energy calculation scales linearly as the number of subsystems, and the method provides a path-forward to quantum mechanical simulations of millions of atoms.
The energetic, electronic and magnetic structures of Fe2−xCoxVSn alloys: Ab-initio calculations
International Nuclear Information System (INIS)
Density Functional Theory (DFT) calculations of a series of the nonstoichiometric Fe2−xCoxVSn full Heusler alloy were carried out utilizing the full potential linearized augmented plane wave (FP-LAPW) method to investigate the electronic, energetic, and magnetic structures of the above systems. Unlike many concentration curves, increasing the cobalt concentration had a crucial effect on the spin polarization as it flattened at 100% at x=1.50, 1.75, and 2.00 where the half- metallic behavior was located with negative formation energy. Moreover, the total magnetic moment of the host material is found to increase with increasing Co concentration. Finally, the half metallic compounds found in some structures of this series might be useful in spintronic devices
DEFF Research Database (Denmark)
Vanin, Marco; Gath, Jesper; Thygesen, Kristian Sommer;
2010-01-01
The stability of graphene nanoribbons in the presence of typical atmospheric molecules is systematically investigated by means of density-functional theory. We calculate the edge formation free energy of five different edge configurations passivated by H, H-2, O, O-2, N-2, CO, CO2, and H2O......, respectively. In addition to the well known hydrogen passivated armchair and zigzag edges, we find the edges saturated by oxygen atoms to be particularly stable under atmospheric conditions. Saturation of the zigzag edge by oxygen leads to the formation of metallic states strictly localized on the oxygen atoms...
First principles total energy calculations of the structural and electronic properties of ScxGa1-xN
International Nuclear Information System (INIS)
Using first principles total energy calculations within the the full-potential linearized augmented plane wave (FP-LAPW) method, we have investigated the structural and electronic properties of ScxGa1-xN, with Sc concentrations varying from 0% up to 100%. In particular we have studied the relative stability of several configurations of ScxGa1-xN in wurtzite-like structures (the ground state configuration of GaN), or in rocksalt-like structures (the ground state configuration of ScN). It is found that for Sc concentrations less than ∼65%, the favored structure is a wurtzite-like one, while for Sc concentrations greater than ∼65%, the favored structure is a NaCl-like structure. It is also found that for the wurtzite-like crystals, the fundamental gap is large and direct. For the rocksalt crystals the fundamental gap is small and indirect, but with an additional larger direct gap. In agreement with the experiments of Little and Kordesch [Appl. Phys. Lett. 78, 2891 (2001)] we found a decrease of the band gap with the increase of the Sc concentration. (Abstract Copyright [2003], Wiley Periodicals, Inc.)
Institute of Scientific and Technical Information of China (English)
Zeng Hui; Zhao Jun; Xiao Xun
2013-01-01
Quantum chemical calculations are performed to investigate the equilibrium C-COOH bond distances and the bond dissociation energies (BDEs) for 15 acids.These compounds are studied by utilizing the hybrid density functional theory (DFT) (B3LYP,B3PW91,B3P86,PBE1PBE) and the complete basis set (CBS-Q) method in conjunction with the 6-31 lG** basis as DFT methods have been found to have low basis sets sensitivity for small and medium molecules in our previous work.Comparisons between the computational results and the experimental values reveal that CBS-Q method,which can produce reasonable BDEs for some systems in our previous work,seems unable to predict accurate BDEs here.However,the B3P86 calculated results accord very well with the experimental values,within an average absolute error of 2.3 kcal/mol.Thus,B3P86 method is suitable for computing the reliable BDEs of C-COOH bond for carboxylic acid compounds.In addition,the energy gaps between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) of studied compounds are estimated,based on which the relative thermal stabilities of the studied acids are also discussed.
International Nuclear Information System (INIS)
Quantum chemical calculations are performed to investigate the equilibrium C—COOH bond distances and the bond dissociation energies (BDEs) for 15 acids. These compounds are studied by utilizing the hybrid density functional theory (DFT) (B3LYP, B3PW91, B3P86, PBE1PBE) and the complete basis set (CBS—Q) method in conjunction with the 6-311G** basis as DFT methods have been found to have low basis sets sensitivity for small and medium molecules in our previous work. Comparisons between the computational results and the experimental values reveal that CBS—Q method, which can produce reasonable BDEs for some systems in our previous work, seems unable to predict accurate BDEs here. However, the B3P86 calculated results accord very well with the experimental values, within an average absolute error of 2.3 kcal/mol. Thus, B3P86 method is suitable for computing the reliable BDEs of C—COOH bond for carboxylic acid compounds. In addition, the energy gaps between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) of studied compounds are estimated, based on which the relative thermal stabilities of the studied acids are also discussed. (atomic and molecular physics)
Electron momentum density and band structure calculations of {alpha}- and {beta}-GeTe
Energy Technology Data Exchange (ETDEWEB)
Vadkhiya, Laxman [Department of Physics, University College of Science, M.L. Sukhadia University, Udaipur 313001, Rajasthan (India); Arora, Gunjan [Department of Physics, Techno India NJR Institute of Technology, Udaipur 313002, Rajasthan (India); Rathor, Ashish [Department of Physics, University College of Science, M.L. Sukhadia University, Udaipur 313001, Rajasthan (India); Ahuja, B.L., E-mail: blahuja@yahoo.com [Department of Physics, University College of Science, M.L. Sukhadia University, Udaipur 313001, Rajasthan (India)
2011-12-15
We have measured isotropic experimental Compton profile of {alpha}-GeTe by employing high energy (662 keV) {gamma}-radiation from a {sup 137}Cs isotope. To compare our experiment, we have also computed energy bands, density of states, electron momentum densities and Compton profiles of {alpha}- and {beta}-phases of GeTe using the linear combination of atomic orbitals method. The electron momentum density is found to play a major role in understanding the topology of bands in the vicinity of the Fermi level. It is seen that the density functional theory (DFT) with generalised gradient approximation is relatively in better agreement with the experiment than the local density approximation and hybrid Hartree-Fock/DFT. - Highlights: > Compton profile of {alpha}-GeTe using a 20 Ci {sup 137}Cs Compton spectrometer. > Compared experimental Compton data with density functional theory. > Reported energy bands and density of states of {alpha}- and {beta}-GeTe. > EVED profiles analysed to check the covalent character.
Nakamura, H.; Hayashi, N.; Nakai, N.; Okumura, M.; Machida, M.
2009-10-01
In order to resolve a discrepancy of the magnetic moment on Fe between the experimental and calculation results, we perform first-principle electronic structure calculations for iron-based superconductors LaFeAsO1-x and LiFeAs also show similar SDW. So far, the first-principle calculations on LaFeAsO actually predicted the SDW state as a ground state. However, the predicted magnetic moment (∼2 μB) per an Fe atom is much larger than the observed one (∼0.35 μB) in experiments [2,4]. The authors suggested that the discrepancy can be resolved by expanding U into a negative U range within LSDA + U framework. In this paper, we revisit the discrepancy and clarify why the negative correction is essential in these compounds. See Ref. [5] for the details of calculation data by LSDA + negative U. In the first-principle calculation on compounds including transition metals, the total energy is frequently corrected by “LSDA + U” approach. The parameter U is theoretically re-expressed as U(≡U-J), where U is the on-site Coulomb repulsion (Hubbard U) and J is the atomic-orbital intra-exchange energy (Hund’s coupling parameter) [6]. The parameter U employed in the electronic structure calculations is usually positive. The positivity promotes the localized character of d-electrons and enhances the magnetic moment in the cases of magnetically ordered compounds. Normally, this positive correction successfully works. In choosing the parameter, one can principally extend the parameter U range to a negative region. The negative case [7] is not popular, but it can occur in the following two cases [8]: (i) the Hubbard U becomes negative and (ii) the intra-exchange J is effectively larger than the Hubbard U. The case (i) has been suggested by many authors based on various theoretical considerations. Here, we note that U should be estimated once screening effects on the long-range Coulomb interaction are taken into account. In fact, small U has been reported [9]. Thus, when the
da Silva, E. Lora; Marinopoulos, A. G.; Vieira, R. B. L.; Vilão, R. C.; Alberto, H. V.; Gil, J. M.; Lichti, R. L.; Mengyan, P. W.; Baker, B. B.
2016-07-01
The electronic structure of hydrogen impurity in Lu2O3 was studied by first-principles calculations and muonium spectroscopy. The computational scheme was based on two methods which are well suited to treat defect calculations in f -electron systems: first, a semilocal functional of conventional density-functional theory (DFT) and secondly a DFT+U approach which accounts for the on-site correlation of the 4 f electrons via an effective Hubbard-type interaction. Three different types of stable configurations were found for hydrogen depending upon its charge state. In its negatively charged and neutral states, hydrogen favors interstitial configurations residing either at the unoccupied sites of the oxygen sublattice or at the empty cube centers surrounded by the lanthanide ions. In contrast, the positively charged state stabilized only as a bond configuration, where hydrogen binds to oxygen ions. Overall, the results between the two methods agree in the ordering of the formation energies of the different impurity configurations, though within DFT+U the charge-transition (electrical) levels are found at Fermi-level positions with higher energies. Both methods predict that hydrogen is an amphoteric defect in Lu2O3 if the lowest-energy configurations are used to obtain the charge-transition, thermodynamic levels. The calculations of hyperfine constants for the neutral interstitial configurations show a predominantly isotropic hyperfine interaction with two distinct values of 926 MHz and 1061 MHz for the Fermi-contact term originating from the two corresponding interstitial positions of hydrogen in the lattice. These high values are consistent with the muonium spectroscopy measurements which also reveal a strongly isotropic hyperfine signature for the neutral muonium fraction with a magnitude slightly larger (1130 MHz) from the ab initio results (after scaling with the magnetic moments of the respective nuclei).
Energy Technology Data Exchange (ETDEWEB)
Li, Ji-Hong [Sichuan Univ., Chengdu (China). Inst. of Atomic and Molecular Physics; Longdong Univ., Qingyang (China). College of Physics and Electronic Engineering; Zhu, Xu-Hui [Sichuan Univ., Chengdu (China). Inst. of Atomic and Molecular Physics; Cheng, Yan [Sichuan Univ., Chengdu (China). Inst. of Atomic and Molecular Physics; Sichuan Univ., Chengdu (China). Key Laboratory of High Energy Density Physics and Technology of Ministry of Education; Ji, Guang-Fu [Chinese Academy of Engineering Physics, Mianyang (China). National Key Laboratory of Shock Wave and Detonation Physics
2015-07-01
Based on the first-principles density functional theory calculations combined with the quasi-harmonic Debye model, the pressure dependencies of the structural, elastic, electronic and thermal properties of Li{sub 2}AgSb were systematically investigated. The calculated lattice parameters and unit cell volume of Li{sub 2}AgSb at the ground state were in good agreement with the available experimental data. The obtained elastic constants, the bulk modulus and the shear modulus revealed that Li{sub 2}AgSb is mechanically stable and behaves in a ductile manner under the applied pressure. The elasticity-relevant properties, the Young's modulus and the Poisson's ratio showed that pressure can enhance the stiffness of Li{sub 2}AgSb and that Li{sub 2}AgSb is mechanically stable up to 20 GPa. The characteristics of the band structure and the partial density of states of Li{sub 2}AgSb were analysed, showing that Li{sub 2}AgSb is a semiconductor with a direct band gap of 217 meV at 0 GPa and that the increasing pressure can make the band structure of Li{sub 2}AgSb become an indirect one. Studies have shown that, unlike temperature, pressure has little effect on the heat capacity and the thermal expansion coefficient of Li{sub 2}AgSb.
Revised self-consistent continuum solvation in electronic-structure calculations
Andreussi, Oliviero; Marzari, Nicola
2011-01-01
The solvation model proposed by Fattebert and Gygi [Journal of Computational Chemistry 23, 662 (2002)] and Scherlis et al. [Journal of Chemical Physics 124, 074103 (2006)] is reformulated, overcoming some of the numerical limitations encountered and extending its range of applicability. We first recast the problem in terms of induced polarization charges that act as a direct mapping of the self-consistent continuum dielectric; this allows to define a functional form for the dielectric that is well behaved both in the high-density region of the nuclear charges and in the low-density region where the electronic wavefunctions decay into the solvent. Second, we outline an iterative procedure to solve the Poisson equation for the quantum fragment embedded in the solvent that does not require multi-grid algorithms, is trivially parallel, and can be applied to any Bravais crystallographic system. Last, we capture some of the non-electrostatic or cavitation terms via a combined use of the quantum volume and quantum s...
Atomic partial charges on CH3NH3PbI3 from first-principles electronic structure calculations
Madjet, Mohamed E.; El-Mellouhi, Fedwa; Carignano, Marcelo A.; Berdiyorov, Golibjon R.
2016-04-01
We calculated the partial charges in methylammonium (MA) lead-iodide perovskite CH3NH3PbI3 in its different crystalline phases using different first-principles electronic charge partitioning approaches, including the Bader, ChelpG, and density-derived electrostatic and chemical (DDEC) schemes. Among the three charge partitioning methods, the DDEC approach provides chemically intuitive and reliable atomic charges for this material, which consists of a mixture of transition metals, halide ions, and organic molecules. The DDEC charges are also found to be robust against the use of hybrid functionals and/or upon inclusion of spin-orbit coupling or dispersive interactions. We calculated explicitly the atomic charges with a special focus on the dipole moment of the MA molecules within the perovskite structure. The value of the dipole moment of the MA is reduced with respect to the isolated molecule due to charge redistribution involving the inorganic cage. DDEC charges and dipole moment of the organic part remain nearly unchanged upon its rotation within the octahedral cavities. Our findings will be of both fundamental and practical importance, as the accurate and consistent determination of the atomic charges is important in order to understand the average equilibrium distribution of the electrons and to help in the development of force fields for larger scale atomistic simulations to describe static, dynamic, and thermodynamic properties of the material.
International Nuclear Information System (INIS)
Graphical abstract: Calculated heats of formation compared to experimental and theoretical data for Ca–Zn system intermetallic compounds. Highlights: ► Ca–Zn system intermetallic compounds have been studied. ► Ca–Zn intermetallic compounds are all conductors. ► Ca–Zn intermetallic compounds are all stable. - Abstract: Structural, elastic and electronic properties, as well as heats of formation, of seven Ca–Zn intermetallic compounds have been studied by using first principles methods. It was found that with increasing Zn concentration, the bulk moduli and shear moduli of Ca–Zn intermetallic compounds increase monotonically. Our results also indicate that Ca3Zn, Ca5Zn3, and CaZn are ductile, while CaZn2, CaZn5, CaZn11, and CaZn13 are brittle. Furthermore, calculations of the electronic properties and heats of formation indicate that seven Ca–Zn intermetallic compounds, considered in this work, are all conductors and thermodynamically stable.
Languages for structural calculations
International Nuclear Information System (INIS)
The differences between human and computing languages are recalled. It is argued that they are to some extent structured in antagonistic ways. Languages in structural calculation, in the past, present, and future, are considered. The contribution of artificial intelligence is stressed
Energy Technology Data Exchange (ETDEWEB)
Kostko, Oleg; Bravaya, Ksenia; Krylov, Anna; Ahmed, Musahid
2009-12-14
We report a combined theoretical and experimental study of ionization of cytosine monomers and dimers. Gas-phase molecules are generated by thermal vaporization of cytosine followed by expansion of the vapor in a continuous supersonic jet seeded in Ar. The resulting species are investigated by single photon ionization with tunable vacuum-ultraviolet (VUV) synchrotron radiation and mass analyzed using reflectron mass spectrometry. Energy onsets for the measured photoionization efficiency (PIE) spectra are 8.60+-0.05 eV and 7.6+-0.1 eV for the monomer and the dimer, respectively, and provide an estimate for the adiabatic ionization energies (AIE). The first AIE and the ten lowest vertical ionization energies (VIEs) for selected isomers of cytosine dimer computed using equation-of-motion coupled-cluster (EOM-IP-CCSD) method are reported. The comparison of the computed VIEs with the derivative of the PIE spectra, suggests that multiple isomers of the cytosine dimer are present in the molecular beam. The calculations reveal that the large red shift (0.7 eV) of the first IE of the lowest-energy cytosine dimer is due to strong inter-fragment electrostatic interactions, i.e., the hole localized on one of the fragments is stabilized by the dipole moment of the other. A sharp rise in the CH+ signal at 9.20+-0.05 eV is ascribed to the formation of protonated cytosine by dissociation of the ionized dimers. The dominant role of this channel is supported by the computed energy thresholds for the CH+ appearance and the barrierless or nearly barrierless ionization-induced proton transfer observed for five isomers of the dimer.
Mitsuoka, Kaoru
2014-11-01
Bacteriorhodopsin (bR) is a light-driven proton pump, which is a membrane protein found in halophilic archeae like Halobacterium salinarum and in eubacteria [1]. When the covalently bound retinal chromophore absorbs the light energy, it changes the conformation from all-trans to 13-cis. This configuration change initiates ion translocation across the cell membrane and a proton moves from inside to outside of the cell. The bR molecules are forming two-dimensional crystals on the membranes of halophilic archeae, and therefore the atomic model of bR was first determined by electron crystallography. The determined structure can be used to determine the pKa values, through which the charge states of ionizable residues in bR determine their pH-dependent properties. The pH-dependent properties are crucial for proton translocation from ionizable residues or to ionizable residues. Detection of the intermediate states of the reaction cycle (photocycle) produced spectroscopic information, which can predict the ionization state of the ionozable residues. In the transition from the L intermediate to the M intermediate, it is known that a proton moves from the Shiff base on the retinal chromophore to Asp85, while a proton is released to the extracellar side from proton-releasing groups including Glu194 and Glu204. Experimentally the pKa value of the proton release is determined to be about 9.7, while the pKa value of Asp85 was measured to change from 2.6 to 7.5 by the proton release from the proton-releasing groups [2]. Here we used the PROPKA program [3] to calculate the pKa values of Asp85 and the proton-releasing groups from the structures at pH 5.5 and at pH 10.0 determined by electron crystallography. The calculation showed that the pKa value of Asp85 changes from 5.3 to 6.1, which qualitatively show the similar changes with the measured difference. The largest change between the structures is the shift of Arg82 by the proton release from the proton-releasing groups
Structural and electronic properties of SrAl2O4:Eu2+ from density functional theory calculations
International Nuclear Information System (INIS)
Highlights: •Persistent phosphor SrAl2O4:Eu2+ was synthesized and studied. •Ab initio calculations of its electronic properties were performed. •Lowest position of the Eu 4f states in the band gap was determined. •Position of the Eu 4f states agrees with the charge transfer transition. -- Abstract: A stoichiometric micro-sized powder SrAl2O4:Eu2+ was synthesized by traditional solid state reaction at 1250 °C. Low-temperature spectroscopic measurements revealed two luminescence bands at 450 nm and 512 nm; their origin was discussed. Theoretical calculations of the structural and optical properties of SrAl2O4:Eu2+ in the framework of the density functional theory (DFT) were carried out; the obtained results were compared with the corresponding experimental data. For the first time, the position of the lowest 4f states of Eu in the host’s band gap was calculated for both available Sr positions to be at about 4.5–5 eV above the top of the valence band. Reliability of this result is confirmed by good agreement with the experimental value of the O(2p)–Eu(4f) charge transfer energy, which is equal to about 4.9 eV
Accurate variational electronic structure calculations with the density matrix renormalization group
Wouters, Sebastian
2014-01-01
During the past 15 years, the density matrix renormalization group (DMRG) has become increasingly important for ab initio quantum chemistry. The underlying matrix product state (MPS) ansatz is a low-rank decomposition of the full configuration interaction tensor. The virtual dimension of the MPS controls the size of the corner of the many-body Hilbert space that can be reached. Whereas the MPS ansatz will only yield an efficient description for noncritical one-dimensional systems, it can still be used as a variational ansatz for other finite-size systems. Rather large virtual dimensions are then required. The two most important aspects to reduce the corresponding computational cost are a proper choice and ordering of the active space orbitals, and the exploitation of the symmetry group of the Hamiltonian. By taking care of both aspects, DMRG becomes an efficient replacement for exact diagonalization in quantum chemistry. DMRG and Hartree-Fock theory have an analogous structure. The former can be interpreted a...
Martin-Samos, Layla; Bussi, Giovanni
2009-08-01
We present here SaX (Self-energies and eXcitations), a plane-waves package aimed at electronic-structure and optical-properties calculations in the GW framework, namely using the GW approximation for quasi-particle properties and the Bethe-Salpeter equation for the excitonic effects. The code is mostly written in FORTRAN90 in a modern style, with extensive use of data abstraction (i.e. objects). SaX employs state of the art techniques and can treat large systems. The package is released with an open source license and can be also download from http://www.sax-project.org/. Program summaryProgram title: SaX (Self-energies and eXcitations) Catalogue identifier: AEDF_v1_0 Program summary URL:http://cpc.cs.qub.ac.uk/summaries/AEDF_v1_0.html Program obtainable from: CPC Program Library, Queen's University, Belfast, N. Ireland Licensing provisions: GNU General Public License No. of lines in distributed program, including test data, etc.: 779 771 No. of bytes in distributed program, including test data, etc.: 4 894 755 Distribution format: tar.gz Programming language: FORTRAN, plus some C utilities Computer: Linux PC, Linux clusters, IBM-SP5 Operating system: Linux, Aix Has the code been vectorised or parallelized?: Yes RAM: depending on the system complexity Classification: 7.3 External routines: Message-Passing Interface (MPI) to perform parallel computations. ESPRESSO ( http://www.quantum-espresso.org) Nature of problem: SaX is designed to calculate the electronic band-structure of semiconductors, including quasi-particle effects and optical properties including excitonic effects. Solution method: The electronic band-structure is calculated using the GW approximation for the self-energy operator. The optical properties are calculated solving the Bethe-Salpeter equation in the GW approximation. The wavefunctions are expanded on a plane-waves basis set, using norm-conserving pseudopotentials. Restrictions: Many objects are non-local matrix represented in plane wave basis
Cao, Jun; Xie, Zhi-Zhong
2016-03-01
The ab initio electronic structure calculations and CASSCF-based nonadiabatic dynamics simulations have been used to investigate the internal conversion and intersystem crossing process of both trans-acrolein and 2-cyclopentenone in the gas phase. Our calculation results show that relaxation from the Franck-Condon region to an S1 minimum is ultrafast and that the S1 state will dominantly undergo intersystem crossing to triplet states due to the existence of significant barriers to access the S1/S0 intersection points and of energetically close-lying triplet states. The S1/T2/T1 three-state intersection is observed in our dynamics simulations to play an important role in the population of the lowest triplet state, which is consistent with previous suggestions. Although the evolution into triplet states involves a similar path and gives rise to a similar triplet quantum yield for these two molecules, the intersystem crossing rate of 2-cyclopentenone is lower owing to the ring constraint that results in a smaller spin-orbital coupling in the singlet-triplet crossing region. The present theoretical study reproduces the experimental results and gives an explanation about the structural factors that govern the excited-state decay of some types of α,β-enones. PMID:26882275
Řezáč, Jan; Huang, Yuanhang; Hobza, Pavel; Beran, Gregory J O
2015-07-14
Many-body noncovalent interactions are increasingly important in large and/or condensed-phase systems, but the current understanding of how well various models predict these interactions is limited. Here, benchmark complete-basis set coupled cluster singles, doubles, and perturbative triples (CCSD(T)) calculations have been performed to generate a new test set for three-body intermolecular interactions. This "3B-69" benchmark set includes three-body interaction energies for 69 total trimer structures, consisting of three structures from each of 23 different molecular crystals. By including structures that exhibit a variety of intermolecular interactions and packing arrangements, this set provides a stringent test for the ability of electronic structure methods to describe the correct physics involved in the interactions. Both MP2.5 (the average of second- and third-order Møller-Plesset perturbation theory) and spin-component-scaled CCSD for noncovalent interactions (SCS-MI-CCSD) perform well. MP2 handles the polarization aspects reasonably well, but it omits three-body dispersion. In contrast, many widely used density functionals corrected with three-body D3 dispersion correction perform comparatively poorly. The primary difficulty stems from the treatment of exchange and polarization in the functionals rather than from the dispersion correction, though the three-body dispersion may also be moderately underestimated by the D3 correction. PMID:26575743
Role of anion doping on electronic structure and magnetism of GdN by first principles calculations
Zhang, Xuejing
2014-01-01
We have investigated the electronic structure and magnetism of anion doped GdN1-yXy (X = B, C, O, F, P, S and As) systems by first-principles calculations based on density functional theory. GdN 1-yXy systems doped by O, C, F, P, and S atoms are more stable than those doped by B and As atoms because of relatively high binding energies. The anion doping and the N defect states modify the density of states at the Fermi level, resulting in a decrease in spin polarization and a slight increase in the magnetic moment at the Gd and N sites. © 2014 The Royal Society of Chemistry.
International Nuclear Information System (INIS)
The electronic structure of the interstitial hydrogen atom in KF, NaCl, KCl, and RbCl cristals has been studied using the self-consistent-field multiple-scattering Xα method. In the present calculation a cluster constituted by the hydrogen atom surrounded by its first anion and cation neighbors in a cubic shell has been used. The optical transition energies and hyperfine contact parameters with the interstitial proton and the first shell nuclei have been evaluated. The agreement obtained with the experimental data and the relative independence of the method under variations of its intrinsic parameters, indicate that this method can be adequate to the study of defects in ionic cristals. (author)
Energy Technology Data Exchange (ETDEWEB)
Hussain, Altaf, E-mail: altafiub@yahoo.com [Department of Physics, Islamia University of Bahawalpur, Punjab 63100 (Pakistan); Aryal, Sitaram; Rulis, Paul [Department of Physics, University of Missouri-Kansas City, MO 64110 (United States); Choudhry, M. Arshad [Department of Physics, Islamia University of Bahawalpur, Punjab 63100 (Pakistan); Chen, Jun [Institute of Applied Physics and Computational Mathematics, Beijing 10088 (China); Ching, W.Y. [Department of Physics, University of Missouri-Kansas City, MO 64110 (United States)
2011-04-28
The electronic structure and optical properties of the Ni{sub 3}Al intermetallic alloy are studied by the first-principles orthogonalized linear combination of atomic orbitals method. Disordered models at different temperatures were constructed using molecular dynamics and the Vienna ab initio simulation package. The average charge transfer from Al to Ni increases steadily with temperature until the liquid phase is reached. The localization index shows the presence of relatively localized states even above the Fermi level in the disordered models. The calculated optical conductivity of the ordered phase is rich in structures and in reasonable agreement with the experimental data. The spectra of the disordered Ni{sub 3}Al models show a single broadened peak at 4.96 eV in the 0 K model which shifts towards 6.62 eV at 1400 K and then down to 5.83 eV in the liquid phase. Other results on the band structure and density of states are also discussed.
Electronic structure and point defect concentrations of C11b MoSi2 by first-principles calculations
International Nuclear Information System (INIS)
Highlights: • The point defects of C11b MoSi2 were studied systematically. • MoSi2 is semimetallic with strong directional covalent bonds. • Some rules of the point defect concentrations were revealed. • Vacancy is a main type of point defect in MoSi2. - Abstract: The electronic structure and point defect concentrations of C11b MoSi2 were studied systematically by the first-principles calculations based on density functional theory. Mo vacancy-induced charge density shows strong directional covalent bonds caused by hybridization of Mo-4d and Si-3p orbitals, which indicates that MoSi2 has low fracture toughness at room temperature. Combining with Wagner–Schottky model, these point defect concentrations of C11b MoSi2 at 2173, 1673, 1223, 773 K as function of composition were also investigated. It is found that the point defect concentrations change drastically for off-stoichiometric compounds. The main structural defects are preferably Mo vacancies or Si anti-structure atoms on the Mo sublattices in Si-rich alloy, and Mo anti-site in Mo-rich alloy, respectively. According to the calculated effective formation enthalpies of point defects, the effective formation enthalpies from big to small in sequence are Mo anti-site, Si anti-site and vacancy (Mo and Si). This result suggests that the vacancy, especially for Si vacancy, is a main type of point defect in C11b MoSi2 system
Joshi, Bhawani Datt; Srivastava, Anubha; Honorato, Sara Braga; Tandon, Poonam; Pessoa, Otília Deusdênia Loiola; Fechine, Pierre Basílio Almeida; Ayala, Alejandro Pedro
2013-09-01
Oncocalyxone A (C17H18O5) is the major secondary metabolite isolated from ethanol extract from the heartwood of Auxemma oncocalyx Taub popularly known as “pau branco”. Oncocalyxone A (Onco A) has many pharmaceutical uses such as: antitumor, analgesic, antioxidant and causative of inhibition of platelet activation. We have performed the optimized geometry, total energy, conformational study, molecular electrostatic potential mapping, frontier orbital energy gap and vibrational frequencies of Onco A employing ab initio Hartree-Fock (HF) and density functional theory (DFT/B3LYP) method with 6-311++G(d, p) basis set. Stability of the molecule arising from hyperconjugative interactions and/or charge delocalization has been analyzed using natural bond orbital (NBO) analysis. UV-vis spectrum of the compound was recorded in DMSO and MeOH solvent. The TD-DFT calculations have been performed to explore the influence of electronic absorption spectra in the gas phase, as well as in solution environment using IEF-PCM and 6-31G basis set. The 13C NMR chemical shifts have been calculated with the B3LYP/6-311++G(d, p) basis set and compared with the experimental values. These methods have been used as tools for structural characterization of Onco A.
Pask, J. E.; Sterne, P. A.
2004-03-01
The finite-element (FE) method is a general approach for the solution of partial differential equations. Like the planewave (PW) method, the FE method is a systematically improvable expansion approach. Unlike the PW method, however, its basis functions are strictly local in real space, which allows for variable resolution in real space and facilitates massively parallel implementation. We discuss the application of the FE method to ab initio electronic-structure calculations.(J.E. Pask, B.M. Klein, C.Y. Fong, and P.A. Sterne, Phys. Rev. B 59), 12352 (1999). In particular, we discuss the use of nonlocal pseudopotentials in bulk calculations, and the handling of long-range interactions in the construction of the Kohn-Sham effective potential and total energy. We show that the total energy converges variationally, and at the optimal theoretical rate consistent with the cubic completeness of the basis. This work was performed under the auspices of the U.S. Department of Energy by University of California, Lawrence Livermore National Laboratory under Contract W-7405-Eng-48.
International Nuclear Information System (INIS)
Structural, electronic and magnetic properties were calculated for the optimized α-U/W(110) thin films (TFs) within the density functional theory. Our optimization for 1U/7W(110) shows that the U–W vertical interlayer spacing (ILS) is expanded by 14.0% compared to our calculated bulk W–W ILS. The spin and orbital magnetic moments (MMs) per U atom were found to be enhanced from zero for the bulk of α-U to 1.4 μB and − 0.4 μB at the interface of the 1U/7W(110), respectively. Inversely, our result for 3U/7W(110) TFs shows that the surface U–U ILS is contracted by 15.7% compared to our obtained bulk U–U spacing. The enhanced spin and orbital MMs in the 1U/7W(110) were then found to be suppressed in 3U/7W(110) to their ignorable bulk values. The calculated density of states (DOS) corroborates the enhancement and suppression of the MMs and shows that the total DOS, in agreement with experiment, is dominated in the vicinity of Fermi level by the 5f U states. Proximity and mismatch effects of the nonmagnetic W(110) substrate were assessed and found to be important for this system. - Highlights: ► We model α-U/W(110) thin films by reconstructing an fcc-like structure for α-U(001). ► The fcc-like α-U(001) has a propensity to wet the fcc-like W(110) substrate. ► In contrast to the freestanding α-U(001), α-U/W(110) is not a magnetic system. ► Interatomic U–U distances are more contracted in α-U/W(110) than α-U(001). ► 5f-U-electrons behave itinerantly and predominate the valence bands in α-U/W(110).
Chae, Kisung; Kim, Hanchul
2013-02-01
We perform spin-polarized density functional theory calculations for a Zn vacancy on the ZnO(10bar 10) surface. Two stable configurations of the surface Zn vacancy are found, and the activation energy barrier is estimated to be ˜0.01 eV. The lower energy configuration has a newly formed surface Zn-O bond to restore the bulk-like structure on the surface. Due to the newly formed bond, the vacancy state in the band gap is characterized by a complicated hybridization of neighboring surface and subsurface atoms and by a more extended electron density. Despite such a hybridization, the surface Zn vacancy is found to have a robust magnetic moment of 1 μ B , implying that surface Zn vacancies may be responsible for the ferromagnetism observed in ZnO thin films and nanoparticles. Simulated scanning tunneling microscope images show that the two structures of the surface Zn vacancy can be distinguished in the filled-state images.
Semidirect algorithms in electron propagator calculations
Energy Technology Data Exchange (ETDEWEB)
Zakrzewski, V.G.; Ortiz, J.V. [Univ. of New Mexico, Albuquerque, NM (United States)
1994-12-31
Electron propagator calculations have been executed with a semi-direct algorithm that generates only a subset of transformed electron repulsion integrals and that takes advantage of Abelian point group symmetry. Diagonal self-energy expressions that are advantageous for large molecules are employed. Illustrative calculations with basis sets in excess of 200 functions include evaluations of the ionization energies of C{sup 2{minus}}{sub 7} and Zn(C{sub 5}H{sub 5}){sub 2}. In the former application, a bound dianion is obtained for a D{sub 3h} structure. In the latter, many final states of the same symmetry are calculated without difficulty.
Energy Technology Data Exchange (ETDEWEB)
Guemou, M., E-mail: guemoumhamed7@gmail.com [Engineering Physics Laboratory, University Ibn Khaldoun of Tiaret, BP 78-Zaaroura, Tiaret 14000 (Algeria); Bouhafs, B. [Modelling and Simulation in Materials Science Laboratory, Physics Department, University of Sidi Bel-Abbes, 22000 Sidi Bel-Abbes (Algeria); Abdiche, A. [Applied Materials Laboratory, Research Center, University of Sidi Bel Abbes, 22000 Sidi Bel Abbes (Algeria); Khenata, R. [Laboratoire de Physique Quantique et de Modelisation Mathematique (LPQ3M), Departement de Technologie, Universite de Mascara, 29000 Mascara (Algeria); Al Douri, Y. [Institute of Nano Electronic Engineering, Universiti Malaysia Perlis, Perlis (Malaysia); Bin Omran, S. [Department of Physics and Astronomy, Faculty of Science, King Saud University, P.O. Box 2455, Riyadh 11451 (Saudi Arabia)
2012-04-15
Density functional calculations are performed to study the structural, electronic and optical properties of technologically important B{sub x}Ga{sub 1-x}As ternary alloys. The calculations are based on the total-energy calculations within the full-potential augmented plane-wave (FP-LAPW) method. For exchange-correlation potential, local density approximation (LDA) and the generalized gradient approximation (GGA) have been used. The structural properties, including lattice constants, bulk modulus and their pressure derivatives, are in very good agreement with the available experimental and theoretical data. The electronic band structure, density of states for the binary compounds and their ternary alloys are given. The dielectric function and the refractive index are also calculated using different models. The obtained results compare very well with previous calculations and experimental measurements.
International Nuclear Information System (INIS)
Density functional calculations are performed to study the structural, electronic and optical properties of technologically important BxGa1-xAs ternary alloys. The calculations are based on the total-energy calculations within the full-potential augmented plane-wave (FP-LAPW) method. For exchange-correlation potential, local density approximation (LDA) and the generalized gradient approximation (GGA) have been used. The structural properties, including lattice constants, bulk modulus and their pressure derivatives, are in very good agreement with the available experimental and theoretical data. The electronic band structure, density of states for the binary compounds and their ternary alloys are given. The dielectric function and the refractive index are also calculated using different models. The obtained results compare very well with previous calculations and experimental measurements.
Institute of Scientific and Technical Information of China (English)
LIU Zhilin; LIN Cheng; LIU Yan; GUO Yanchang
2005-01-01
Based on the phase transformations and strengthening mechanisms during roiling, the strength increments △σb under different strengthening mechanisms are calculated with the covalent electron number nA of the strongest bond in phase cells of alloys and the interface electron density difference △ρ matching the interface stress in alloys. The calculation method of the finishing rolling yield strength is proposed, and it is integrated with the proposed calculation formulas of strength of non quenched-tempered steel. Therefore,the general formulas to simultaneously calculate both the finishing rolling strength and the yield strength of the continuous casting-rolling and non quenched-tempered steel are given. Taken the pipeline steel X70 as an example, the predictions of properties and technological parameters are performed before production or online.
Aryal, Sita Ram
The alumino-silicate solid solution series (Al 4+2xSi2-2 xO10-x) is an important class of ceramics. Except for the end member (x=0), Al2 SiO5 the crystal structures of the other phases, called mullite, have partially occupied sites. Stoichiometric supercell models for the four mullite phases 3Al2O 3 · 2SiO2 · 2Al 2O3 · SiO2, 4 Al2O3· SiO 2, 9Al2O3 · SiO2, and iota-Al2 O3 (iota-alumina) are constructed starting from experimentally reported crystal structures. A large number of models were built for each phase and relaxed using the Vienna ab initio simulation package (VASP) program. The model with the lowest total energy for a given x was chosen as the representative structure for that phase. Electronic structure and mechanical properties of mullite phases were studied via first-principles calculations. Of the various phases of transition alumina, iota-Al 2O3 is the least well known. In addition structural details have not, until now, been available. It is the end member of the aluminosilicate solid solution series with x=1. Based on a high alumina content mullite phase, a structural model for iota- Al2O3 is constructed. The simulated x-ray diffraction (XRD) pattern of this model agrees well with a measured XRD pattern. The iota-Al2 O3 is a highly disordered ultra-low-density phase of alumina with a theoretical density of 2854kg/m3. Using this theoretically constructed model, elastic, thermodynamic, electronic, and spectroscopic properties of iota-Al2 O3 have been calculated and compared it with those of alpha- Al2O3 and gamma- Al2O3. Boron carbide (B4C) undergoes an amorphization under high velocity impacts. The mechanism of amorphization is not clear. Ab initio methods are used to carry out large-scale uniaxial compression simulations on two polytypes of stoichiometric boron carbide (B4C), B 11C-CBC, and B12- CCC where B11C or B12 is the 12-atom icosahedron and CBC or CCC is the three-atom chain. The simulations were performed on large supercells of 180 atoms
Electronics reliability calculation and design
Dummer, Geoffrey W A; Hiller, N
1966-01-01
Electronics Reliability-Calculation and Design provides an introduction to the fundamental concepts of reliability. The increasing complexity of electronic equipment has made problems in designing and manufacturing a reliable product more and more difficult. Specific techniques have been developed that enable designers to integrate reliability into their products, and reliability has become a science in its own right. The book begins with a discussion of basic mathematical and statistical concepts, including arithmetic mean, frequency distribution, median and mode, scatter or dispersion of mea
International Nuclear Information System (INIS)
The efficiency of fast neutron reactors, such as for fusion, breeding and transmutation, depend strongly on the neutron radiation resistance of the materials used in the reactors. The binary Fe-Cr alloy, which has many attractive properties in this regard, is the base for the best steels of today which are, however, still not up to the required standards. Therefore, substantial effort has been devoted to finding new materials that can cope with the demands better. Experimental studies must be complemented with extensive theoretical modelling in order to understand the effects that different alloying elements has on the resistance properties of materials. To this end, the first steps of multi-scale modelling has been taken, starting out with ab initio calculations of the electronic structure of the complete concentration range range of the disordered binary Fe-C alloy. The mixing enthalpy of Fe-Cr has been quantitatively predicted and has, together with data from literature, been used in order to fit two sets of interatomic potentials for the purpose of simulating defect evolution with molecular dynamics and kinetic Monte-Carlo codes. These dedicated Fe-Cr alloy potentials are new and represent important additions to the pure element potentials that can be found in literature
Electronic Structures of S/C-Doped TiO2 Anatase (101 Surface: First-Principles Calculations
Directory of Open Access Journals (Sweden)
Qili Chen
2014-01-01
Full Text Available The electronic structures of sulfur (S or carbon (C-doped TiO2 anatase (101 surfaces have been investigated by density functional theory (DFT plane-wave pseudopotential method. The general gradient approximation (GGA + U (Hubbard coefficient method has been adopted to describe the exchange-correlation effects. All the possible doping situations, including S/C dopants at lattice oxygen (O sites (anion doping, S/C dopants at titanium (Ti sites (cation doping, and the coexisting of anion and cation doping, were studied. By comparing the formation energies, it was found that the complex of anion and cation doping configuration forms easily in the most range of O chemical potential for both S and C doping. The calculated density of states for various S/C doping systems shows that the synergistic effects of S impurities at lattice O and Ti sites lead a sharp band gap narrowing of 1.35 eV for S-doped system comparing with the pure TiO2 system.
Energy Technology Data Exchange (ETDEWEB)
Olsson, Paer
2004-04-01
The efficiency of fast neutron reactors, such as for fusion, breeding and transmutation, depend strongly on the neutron radiation resistance of the materials used in the reactors. The binary Fe-Cr alloy, which has many attractive properties in this regard, is the base for the best steels of today which are, however, still not up to the required standards. Therefore, substantial effort has been devoted to finding new materials that can cope with the demands better. Experimental studies must be complemented with extensive theoretical modelling in order to understand the effects that different alloying elements has on the resistance properties of materials. To this end, the first steps of multi-scale modelling has been taken, starting out with ab initio calculations of the electronic structure of the complete concentration range range of the disordered binary Fe-C alloy. The mixing enthalpy of Fe-Cr has been quantitatively predicted and has, together with data from literature, been used in order to fit two sets of interatomic potentials for the purpose of simulating defect evolution with molecular dynamics and kinetic Monte-Carlo codes. These dedicated Fe-Cr alloy potentials are new and represent important additions to the pure element potentials that can be found in literature.
Electronic structure and rovibrational calculation of the low-lying states of the RbYb molecule
Energy Technology Data Exchange (ETDEWEB)
Tohme, S.N. [Faculty of Science, Beirut Arab University, P.O. Box 11-5020, Riad El Solh, Beirut 1107 2809 (Lebanon); Korek, M., E-mail: fkorek@yahoo.com [Faculty of Science, Beirut Arab University, P.O. Box 11-5020, Riad El Solh, Beirut 1107 2809 (Lebanon)
2013-01-02
Highlights: Black-Right-Pointing-Pointer Potential energy curves of 29 electronic states of YbRb molecule are calculated. Black-Right-Pointing-Pointer We investigated the spectroscopic constants T{sub e}, R{sub e}, {omega}{sub e}, B{sub e}. Black-Right-Pointing-Pointer The rovibrational constants E{sub v}, B{sub v}, D{sub v}, R{sub min}, and R{sub max} have been calculated. Black-Right-Pointing-Pointer We studied 26 new electronic states for the first time. -- Abstract: Complete Active Space Self Consistent Field (CASSCF) method with Multi Reference Configuration Interaction (MRCI) calculations is used to investigate the potential energy curves of the low-lying 29 electronic states in the representation {sup 2s+1}{Lambda}{sup (+/-)} of the RbYb molecule (single and double excitations with Davidson corrections). The harmonic frequency {omega}{sub e}, the internuclear distance R{sub e} and the electronic energy with respect to the ground state T{sub e} have been calculated. The eigenvalues E{sub v}, the rotational constant B{sub v}, and the abscissas of the turning points R{sub min} and R{sub max} have been investigated using the canonical functions approach. The comparison between the values of the present work and those available in the literature for several states shows a very good agreement. Twenty-six new states have been studied here for the first time.
Relativistic calculations for many electron atoms
International Nuclear Information System (INIS)
Many improvements have now been introduced in ab-initio methods for relativistic atomic structure calculations. After a short description of the different methods, we review the various contributions to energy levels and compare the most recent theoretical and experimental results for few electron heavy ions
Energy Technology Data Exchange (ETDEWEB)
Masrour, R., E-mail: rachidmasrour@hotmail.com [Laboratory of Materials, Processes, Environment and Quality, Cady Ayyed University, National School of Applied Sciences, 63 46000 Safi (Morocco); LMPHE (URAC 12), Faculty of Science, Mohammed V-Agdal University, Rabat (Morocco); Hlil, E.K. [Institut Néel, CNRS et Université Joseph Fourier, BP 166, F-38042 Grenoble cedex 9 (France); Hamedoun, M. [Institute of Nanomaterials and Nanotechnologies, MAScIR, Rabat (Morocco); Benyoussef, A. [LMPHE (URAC 12), Faculty of Science, Mohammed V-Agdal University, Rabat (Morocco); Institute of Nanomaterials and Nanotechnologies, MAScIR, Rabat (Morocco); Hassan II Academy of Science and Technology, Rabat (Morocco)
2014-04-01
Self-consistent ab initio calculations, based on Density Functional Theory (DFT) approach and using Full Potential Linear Augmented Plane Wave (FLAPW) method within GGA+U approximation, are performed to investigate both electronic and magnetic properties of the GdS layers. Polarized spin and spin–orbit coupling are included in calculations within the framework of the antiferromagnetic state between two adjacent Gd layers. Magnetic moment considered to lie along (001) axes are computed. Obtained data from ab initio calculations are used as input for the High Temperature Series Expansions (HTSEs) calculations to compute other magnetic parameters. Using the Heisenberg model, the exchange interactions between the magnetic atoms Gd–Gd in the same layer and between the magnetic atoms in the adjacent bilayers are estimated. This estimate is obtained using the antiferromagnetic and ferromagnetic energies computed by abinitio calculations for GdS layers. The High Temperature Series Expansions (HTSEs) of the magnetic susceptibility of GdS with antiferromagnetic moment (m{sub Gd}) is given up to sixth order series versus of (J{sub 11}(Gd–Gd)/k{sub B}T). The Néel temperature T{sub N} is obtained by mean field theory and by HTSEs of the magnetic susceptibility series using the Padé approximant method. The critical exponent γ associated with the magnetic susceptibility is calculated for GdS layers. - Highlights: • Electronic and magnetic properties of GdS are investigated using the ab initio calculations. • Obtained data from abinitio calculations are used as input for HTSEs to compute other magnetic parameters. • Néel temperature and critical exponent are deduced using HTSE method.
Electronic structure and rovibrational calculation of the low-lying states of the RbYb molecule
Tohme, S. N.; Korek, M.
2013-01-01
Complete Active Space Self Consistent Field (CASSCF) method with Multi Reference Configuration Interaction (MRCI) calculations is used to investigate the potential energy curves of the low-lying 29 electronic states in the representation 2s+1Λ(+/-) of the RbYb molecule (single and double excitations with Davidson corrections). The harmonic frequency ωe, the internuclear distance Re and the electronic energy with respect to the ground state Te have been calculated. The eigenvalues Ev, the rotational constant Bv, and the abscissas of the turning points Rmin and Rmax have been investigated using the canonical functions approach. The comparison between the values of the present work and those available in the literature for several states shows a very good agreement. Twenty-six new states have been studied here for the first time.
Energy Technology Data Exchange (ETDEWEB)
Benrekia, A.R., E-mail: benrekia.ahmed@yahoo.com [Faculty of Science and Technology, University of Medea (Algeria); Benkhettou, N. [Laboratoire des Materiaux Magnetiques, Faculte des Sciences, Universite Djillali Liabes de Sidi Bel Abbes (Algeria); Nassour, A. [Laboratoire de Cristallographie, Resonance Magnetique et Modelisations (CRM2, UMR CNRS 7036) Institut Jean Barriol, Nancy Universite BP 239, Boulevard des Aiguillettes, 54506 Vandoeuvre-les-Nancy (France); Driz, M. [Applied Material Laboratory (AML), Electronics Department, University of Sidi bel Abbes (DZ 22000) (Algeria); Sahnoun, M. [Laboratoire de Physique Quantique de la Matiere et Modelisations Mathematique (LPQ3M), Faculty of Science and Technology,University of Mascara (Algeria); Lebegue, S. [Laboratoire de Cristallographie, Resonance Magnetique et Modelisations (CRM2, UMR CNRS 7036) Institut Jean Barriol, Nancy Universite BP 239, Boulevard des Aiguillettes, 54506 Vandoeuvre-les-Nancy (France)
2012-07-01
We present first-principles VASP calculations of the structural, electronic, vibrational, and optical properties of paraelectric SrTiO{sub 3} and KTaO{sub 3}. The ab initio calculations are performed in the framework of density functional theory with different exchange-correlation potentials. Our calculated lattice parameters, elastic constants, and vibrational frequencies are found to be in good agreement with the available experimental values. Then, the bandstructures are calculated with the GW approximation, and the corresponding band gap is used to obtain the optical properties of SrTiO{sub 3} and KTaO{sub 3}.
Institute of Scientific and Technical Information of China (English)
Zhang Ming; Zhang Chuan-Hui; Shen Jiang
2011-01-01
Aiming at developing p-type semiconductors and modulating the band gap for photoelectronic devices and band engineering, we present the ab initio numerical simulation of the effect of codoping ZnO with Al, N and Mg on the crystal lattice and electronic structure. The simulations are based on the Perdew-Burke-Ernzerhof generalised-gradient approximation in density functional theory. Results indicate that electrons close to the Fermi level transfer effectively when Al, Mg, and N replace Zn and O atoms, and the theoretical results were consistent with the experiments. The addition of Mg leads to the variation of crystal lattice, expanse of energy band, and change of band gap. These unusual properties are explained in terms of the computed electronic structure, and the results show promise for the development of next-generation photoconducting devices in optoelectronic information science and technology.
Hornos, Tamas; Gali, Adam; Svensson, Bengt G.
2011-01-01
Large-scale and gap error free calculations of the electronic structure of vacancies in 4H-SiC have been carried out using a hybrid density functional (HSE06) and an accurate charge correction scheme. Based on the results the carbon vacancy is proposed to be responsible for the Z1/2 and EH6/7 DLTS centers
Electronic and magnetic structure of BaCoO2 as obtained from LSDA and LSDA+U calculations
Nazir, Safdar
2011-03-01
Density functional theory is used to study the structural, electronic, and magnetic properties of BaCoO2. Structural relaxation for different collinear magnetic configurations points to a remarkable magneto-elastic coupling in BaCoO2. Although we obtain several stable long range ordered magnetic structures, ferromagnetism is energetically favorable in the case of the LSDA method. In contrast, for the LSDA+U method antiferromagnetic ordering is found to be favorable. © 2011 Elsevier B.V. All rights reserved.
International Nuclear Information System (INIS)
In this work we present valence band studies of LaSb2 using angle-resolved photoelectron spectroscopy with synchrotron radiation and compare these data with band structure calculations. Valence band spectra reveal that Sb 5p states are dominant near the Fermi level and are hybridized with the La 5d states just below. The calculations show a fair agreement with the experimentally determined valence band spectra, allowing an identification of the observed features. We measured some dispersion for kbar, especially for Sb 5p states; no significant dispersion was found for k||. (letter to the editor)
Powell, B J; Bernstein, N; Brake, K; McKenzie, Ross H; Meredith, P; Pederson, M R
2016-01-01
We report first principles density functional calculations for hydroquinone (HQ), indolequinone (IQ) and semiquinone (SQ). These molecules are believed to be the basic building blocks of the eumelanins, a class of bio-macromolecules with important biological functions (including photoprotection) and with potential for certain bioengineering applications. We have used the DeltaSCF (difference of self consistent fields) method to study the energy gap between the highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO), Delta_HL. We show that Delta_HL is similar in IQ and SQ but approximately twice as large in HQ. This may have important implications for our understanding of the observed broad band optical absorption of the eumelanins. The possibility of using this difference in Delta_HL to molecularly engineer the electronic properties of eumelanins is discussed. We calculate the infrared and Raman spectra of the three redox forms from first principles. Each of the molecules ...
Nonlinear stability of E centers in Si_{1-x}Ge_{x}: electronic structure calculations
Chroneos, A.; Bracht, H; Jiang, C; Uberuaga, B. P.
2008-01-01
Electronic structure calculations are used to investigate the binding energies of defect pairs composed of lattice vacancies and phosphorus or arsenic atoms (E centers) in silicon-germanium alloys. To describe the local environment surrounding the E center we have generated special quasirandom structures that represent random silicon-germanium alloys. It is predicted that the stability of E centers does not vary linearly with the composition of the silicon-germanium alloy. Interestingly, we p...
Electronic Structures of S/C-Doped TiO2 Anatase (101) Surface: First-Principles Calculations
2014-01-01
The electronic structures of sulfur (S) or carbon (C)-doped TiO2 anatase (101) surfaces have been investigated by density functional theory (DFT) plane-wave pseudopotential method. The general gradient approximation (GGA) + U (Hubbard coefficient) method has been adopted to describe the exchange-correlation effects. All the possible doping situations, including S/C dopants at lattice oxygen (O) sites (anion doping), S/C dopants at titanium (Ti) sites (cation doping), and the coexisting of an...
Indian Academy of Sciences (India)
N Boukhris; H Meradji; S Amara Korba; S Drablia; S Ghemid; F El Haj Hassan
2014-08-01
The structural, electronic and thermal properties of lead chalcogenides PbS, PbSe and BeTe using full-potential linear augmented plane wave (FP-LAPW) method are investigated. The exchange–correlation energy within the local density approximation (LDA) and the generalized gradient approximation (GGA) are described. The calculated structural parameters are in reasonable agreement with the available experimental and theoretical data. The electronic band structure shows that the fundamental energy gap is direct (L–L) for all the compounds. Thermal effects on some macroscopic properties of these compounds are predicted using the quasi-harmonic Debye model, in which the lattice vibrations are taken into account. The variations of the lattice constant, bulk modulus, heat capacity, volume expansion coefficient and Debye temperature with temperature and pressure are obtained successfully. The effect of spin–orbit interaction is found to be negligible in determining the thermal properties and leads to a richer electronic structure.
Energy Technology Data Exchange (ETDEWEB)
Iwano, Kaoru [Graduate University for Advanced Studies, Institute of Materials Structure Science, High Energy Accelerator Research Organization (KEK), 1-1 Oho, Tsukuba 305-0801 (Japan); Shimoi, Yukihiro, E-mail: y.shimoi@aist.go.j [Nanotechnology Research Institute (NRI), National Institute of Advanced Industrial Science and Technology (AIST), 1-1-1 Umezono, Tsukuba 305-8568 (Japan)
2009-02-01
Density-functional theory (DFT) calculations are performed based on the high-temperature structure of (EDO-TTF){sub 2}PF{sub 6}, a quasi-one-dimensional molecular compound that shows both thermal and photoinduced phase transitions. In this structure, the EDO-TTF molecules are one-dimensionally aligned, accompanied with weak dimerization. Contrary to a common sense, our DFT calculations reveal that the pair having a shorter mutual distance has a weaker intermolecular coupling than the pair with a longer one; the latter is appropriate to be called an electronic dimer. We also estimate the corresponding transfer energies and discuss their relevance to spin correlations and optical excitations.
International Nuclear Information System (INIS)
Many conflicting electron microscopy data for In2O3(ZnO)m indicate that it may have the polymorphous and polytypoid structures. We investigate their stabilities based on four controversial models. The calculated results confirm that the models with the zigzag feature are more stable than the others and it is possible to form different zigzag configurations in the samples as observed in the experiments. The dynamic process of eliminating the dangling bonds and the requirements of maximizing the symmetry and the distances between the In atoms in the slabs can be regarded as the dominant rules to stabilize the system, but the statistical equilibrium processes have the chances to transform it from the ground state structures to the other model structures. The study of the electronic structures based on the plane and zigzag models reveals that their band gaps and effective masses increase monotonically with m. The predicted band gaps are consistent with the experimental results. The anisotropic feature of electron effective mass tensor exhibited in the plane model differs from that of the zigzag one, which is so notable that can be employed to determine which model is more close to the actual structure of a given sample. The calculated results confirm the possibilities of the separation of conduction electrons and defects and the existence of the natural optimized transport channels in the layered structures, which demonstrate its advantage over ZnO to transport electrons and benefit its applications in the optoelectronic devices. - Graphical abstract: The conduction electrons are mainly distributed around the boundaries of the plane or zigzag shape. The optimized transport channels can be formed around the boundaries. - Highlights: • The formation mechanisms for the polytypoid structure of In2O3(ZnO)m are revealed. • The predicted band gaps are consistent with the experimental results. • The natural optimized transport channels in the layered structures are
Auger electron transport calculations in biological matter
International Nuclear Information System (INIS)
The talk briefly discussed physical, biophysical, and biological aspects of Auger emitters. A summary of radiationless transition data available in published literature and databases were presented. Data were presented for electron capture (EC), internal conversions (IC), binding energies of some commonly used radionuclides 123I, 124I, 125I, and 158Gd. For each of these Auger emitting radionuclides some examples of Monte Carlo calculated electron spectra of individual decays were presented. Because most Auger electrons emitted in the decay of radionuclides are short range low energy electrons below 1 keV, a brief discussion was presented on most recent development of physics models for energy loss of electrons in condensed phase and compared with other models and gas phase data. Accuracy of electron spectra calculated in the decay of electron shower by Auger emitting radionuclides depends on availability of accurate physics data. Currently, there are many gaps in physics data as input data to computer codes in need of new evaluation. In addition, comparison should be made between deterministic and Monte Carlo methods to access the accuracy and sensitivity of data to methods and the chosen parameters. It has long been recognized that Auger electron show a high-LET like characteristics when radionuclide is very closely bound to DNA. As most Auger electrons are short range low energy electrons and mostly absorbed with the DNA duplex when in close vicinity to DNA duplex, we believe the physical and biological dosimetry are best achieved by using Monte Carlo track structure simulations able to simulate tracks of low energy electrons below 1keV and in particular sub 100 eV in condensed phas
Energy Technology Data Exchange (ETDEWEB)
Brik, M.G., E-mail: brik@fi.tartu.e [Institute of Physics, University of Tartu, Riia 142, Tartu 51014, Tartu (Estonia)
2011-02-15
Detailed ab initio calculations of the structural, electronic, optical and elastic properties of two crystals - magnesite (MgCO{sub 3}) and calcite (CaCO{sub 3}) - are reported in the present paper. Both compounds are important natural minerals, playing an important role in the carbon dioxide cycling. The optimized crystal structures, band gaps, density of states diagrams, elastic constants, optical absorption spectra and refractive indexes dependence on the wavelength all have been calculated and compared, when available, with literature data. Both crystals are indirect band compounds, with calculated band gaps of 5.08 eV for MgCO{sub 3} and 5.023 eV for CaCO{sub 3}. Both values are underestimated by approximately 1.0 eV with respect to the experimental data. Although both crystals have the same structure, substitution of Mg by Ca ions leads to certain differences, which manifest themselves in noticeable change in the electronic bands profiles and widths, shape of the calculated absorption spectra, and values of the elastic constants. Response of both crystals to the applied hydrostatic pressure was analyzed in the pressure range of phase stability, variations of the lattice parameters and characteristic interionic distances were considered. The obtained dependencies of lattice constants and calculated band gap on pressure can be used for prediction of properties of these two hosts at elevated pressures that occur in the Earth's mantle. -- Research highlights: {yields} Ab initio calculations of physical properties of MgCO{sub 3} and CaCO{sub 3} were performed. {yields} Changes of the calculated properties with replacement of Mg by Ca were followed. {yields} Pressure dependence of the structural and electronic properties was analyzed. {yields} Good agreement with experimental data was demonstrated.
International Nuclear Information System (INIS)
RTAsO (R=rare-earth and T=transition metals) compounds have been studied using density functional theory with the local spin-density approximation (LSDA). In order to take into account the strong on-site Coulomb interaction U present in RTAsO, we have also performed the LSDA+U calculations. We investigated the electronic structure with on-site Coulomb potential for the R-derived 4f and T-derived 3d orbitals to obtain the correct ground state of RTAsO. The structural parameters, band structures and density of states have been given in detail. Overall, the technique developed and tested in this work holds promise in enabling accurate and fully predictive calculations of strongly correlated electron materials. A detailed analysis shows that the LSDA+U method provides the better description of our systems. We note that there is no ab-initio study related to these materials.
Energy Technology Data Exchange (ETDEWEB)
Moussa, M.; Zaoui, A. [Modelling and Simulation in Materials Science Laboratory, Physics Department, University of Sidi Bel-Abbes, 22000 Sidi Bel-Abbes (Algeria); Kacimi, S., E-mail: kacimi200x@yahoo.f [Modelling and Simulation in Materials Science Laboratory, Physics Department, University of Sidi Bel-Abbes, 22000 Sidi Bel-Abbes (Algeria); Boukortt, A. [University of Mostaganem, Abdelhamid Ibn Badis 27000 (Algeria); Bouhafs, B. [Modelling and Simulation in Materials Science Laboratory, Physics Department, University of Sidi Bel-Abbes, 22000 Sidi Bel-Abbes (Algeria)
2010-09-01
RTAsO (R=rare-earth and T=transition metals) compounds have been studied using density functional theory with the local spin-density approximation (LSDA). In order to take into account the strong on-site Coulomb interaction U present in RTAsO, we have also performed the LSDA+U calculations. We investigated the electronic structure with on-site Coulomb potential for the R-derived 4f and T-derived 3d orbitals to obtain the correct ground state of RTAsO. The structural parameters, band structures and density of states have been given in detail. Overall, the technique developed and tested in this work holds promise in enabling accurate and fully predictive calculations of strongly correlated electron materials. A detailed analysis shows that the LSDA+U method provides the better description of our systems. We note that there is no ab-initio study related to these materials.
First-principles calculations of structural, electronic and optical properties of CdxZn1-xS alloys
Noor, Naveed Ahmed
2010-10-01
Structural, electronic and optical properties of ternary alloy system CdxZn1-xS have been studied using first-principles approach based on density functional theory. Electronic structure, density of states and energy band gap values for CdxZn1-xS are estimated in the range 0 ≤ x ≤ 1 using both the standard local density approximation (LDA) as well as the generalized gradient approximations (GGA) of Wu-Cohen (WC) for the exchange-correlation potential. It is observed that the direct band gap EgΓ-Γ of CdxZn1-xS decreases nonlinearly with the compositional parameter x, as observed experimentally. It is also found that Cd s and d, S p and Zn d states play a major role in determining the electronic properties of this alloy system. Furthermore, results for complex dielectric constant ε(ω), refractive index n(ω), normal-incidence reflectivity R(ω), absorption coefficient α(ω) and optical conductivity σ(ω) are also described in a wide range of the incident photon energy and compared with the existing experimental data. © 2010 Elsevier B.V. All rights reserved.
International Nuclear Information System (INIS)
Highlights: • The structures of 2,3,5,6-tetrafluoropyridine for its S0 and S1(π, π∗) states have been calculated. • TFPy is rigidly planar in its ground electronic state, but is quasi-planar and floppy in S1. • The barrier to planarity is 30 cm−1 in the excited state. • The observed vibrational frequencies for both states agree well with the computations. • A ring-bending potential energy function for the S1(π, π∗) state was proposed. - Abstract: Infrared and Raman spectra of 2,3,5,6-tetrafluoropyridine (TFPy) were recorded and vibrational frequencies were assigned for its S0 electronic ground states. Ab initio and density functional theory (DFT) calculations were used to complement the experimental work. The lowest electronic excited state of this molecule was investigated with ultraviolet absorption spectroscopy and theoretical CASSCF calculations. The band origin was found to be at 35,704.6 cm−1 in the ultraviolet absorption spectrum. A slightly puckered structure with a barrier to planarity of 30 cm−1 was predicted by CASSCF calculations for the S1(π, π∗) state. Lower frequencies for the out-of-plane ring bending vibrations for the electronic excited state result from the weaker π bonding within the pyridine ring
Energy Technology Data Exchange (ETDEWEB)
Sheu, Hong-Li; Boopalachandran, Praveenkumar [Department of Chemistry, Texas A& M University, College Station, TX 77843-3255 (United States); Kim, Sunghwan [National Center for Biotechnology Information, National Library of Medicine, National Institutes of Health, Department of Health and Human Services, 8600 Rockville Pike, Bethesda, MD 20894 (United States); Laane, Jaan, E-mail: laane@chem.tamu.edu [Department of Chemistry, Texas A& M University, College Station, TX 77843-3255 (United States)
2015-07-29
Highlights: • The structures of 2,3,5,6-tetrafluoropyridine for its S{sub 0} and S{sub 1}(π, π{sup ∗}) states have been calculated. • TFPy is rigidly planar in its ground electronic state, but is quasi-planar and floppy in S{sub 1}. • The barrier to planarity is 30 cm{sup −1} in the excited state. • The observed vibrational frequencies for both states agree well with the computations. • A ring-bending potential energy function for the S{sub 1}(π, π{sup ∗}) state was proposed. - Abstract: Infrared and Raman spectra of 2,3,5,6-tetrafluoropyridine (TFPy) were recorded and vibrational frequencies were assigned for its S{sub 0} electronic ground states. Ab initio and density functional theory (DFT) calculations were used to complement the experimental work. The lowest electronic excited state of this molecule was investigated with ultraviolet absorption spectroscopy and theoretical CASSCF calculations. The band origin was found to be at 35,704.6 cm{sup −1} in the ultraviolet absorption spectrum. A slightly puckered structure with a barrier to planarity of 30 cm{sup −1} was predicted by CASSCF calculations for the S{sub 1}(π, π{sup ∗}) state. Lower frequencies for the out-of-plane ring bending vibrations for the electronic excited state result from the weaker π bonding within the pyridine ring.
Electronic structure of BaFe2As2 as obtained from DFT/ASW first-principles calculations
Schwingenschlögl, Udo
2010-07-02
We use ab-initio calculations based on the augmented spherical wave method within density functional theory to study the magnetic ordering and Fermi surface of BaFe2As2, the parent compound of the hole-doped iron pnictide superconductors (K,Ba)Fe2As2, for the tetragonal I4/mmm as well as the orthorhombic Fmmm structure. In comparison to full potential linear augmented plane wave calculations, we obtain significantly smaller magnetic energies. This finding is remarkable, since the augmented spherical wave method, in general, is known for a most reliable description of magnetism. © 2010 WILEY-VCH Verlag GmbH & Co. KGaA, Weinheim.
Energy Technology Data Exchange (ETDEWEB)
Khyzhun, O.Y., E-mail: khyzhun@ipms.kiev.ua [Institute for Problems of Materials Science, National Academy of Sciences of Ukraine, 3 Krzhyzhanivsky Street, Kyiv 03142 (Ukraine); Bekenev, V.L.; Denysyuk, N.M. [Institute for Problems of Materials Science, National Academy of Sciences of Ukraine, 3 Krzhyzhanivsky Street, Kyiv 03142 (Ukraine); Parasyuk, O.V. [Department of Inorganic and Physical Chemistry, Eastern European National University, 13 Voli Avenue, Lutsk 43025 (Ukraine); Fedorchuk, A.O. [Department of Inorganic and Organic Chemistry, Lviv National University of Veterinary Medicine and Biotechnologies, Pekarska St., 50, 79010 Lviv (Ukraine)
2014-01-05
Highlights: • Electronic structure of TlPb{sub 2}Cl{sub 5} is calculated by the FP-LAPW method. • The valence band is dominated by contributions of Cl 3p states. • Contributions of Pb 6p{sup *} states dominate at the bottom of the conduction band. • The FP-LAPW data allow concluding that TlPb{sub 2}Cl{sub 5} is an indirect-gap material. • XPS core-level and valence-band spectra of polycrystalline TlPb{sub 2}Cl{sub 5} are measured. -- Abstract: We report on first-principles calculations of total and partial densities of states of atoms constituting TlPb{sub 2}Cl{sub 5} using the full potential linearized augmented plane wave (FP-LAPW) method. The calculations reveal that the valence band of TlPb{sub 2}Cl{sub 5} is dominated by contributions of the Cl 3p-like states, which contribute mainly at the top of the valence band with also significant contributions throughout the whole valence-band region. In addition, the bottom of the conduction band of TlPb{sub 2}Cl{sub 5} is composed mainly of contributions of the unoccupied Pb 6p-like states. Our FP-LAPW data indicate that the TlPb{sub 2}Cl{sub 5} compound is an indirect-gap material with band gap of 3.42 eV. The X-ray photoelectron core-level and valence-band spectra for pristine and Ar{sup +} ion-irradiated surfaces of a TlPb{sub 2}Cl{sub 5} polycrystalline sample were measured. The measurements reveal high chemical stability and confirm experimentally the low hygroscopicity of TlPb{sub 2}Cl{sub 5} surface.