Pseudopotentials and Basis Sets. How to generate and test them
|
|
- Annice Terry
- 5 years ago
- Views:
Transcription
1 Pseudopotentials and Basis Sets How to generate and test them
2 Pseudopotential idea Atomic Si Core electrons highly localized very depth energy are chemically inert 1s 2 2s 2 2p 6 3s 2 3p 2 Valence wave functions must be core valence orthogonal to the core wave functions
3 Cut-off radii Balance between softness and transferability controlled by R c R c SOFTNESS TRANSFERABILITY Larger R c : softer pseudo Si Shorter R c : harder pseudo
4 The atom program pseudopotential generation label xc flavor # valence shells # core orbitals n l Valence configuration Cutoff radii
5 Run the shell script (pg.sh) Procedure (I) Check contents of new directory (Si.tm2) Plot the pseudo-potentials/orbitals
6 Pseudo-wave function Procedure (II) Logarithmic derivative Radial Fourier-T Pseudopotential
7 Procedure (III)
8 Radial charge distribution: Compare all-electron with pseudo-charge
9 Pseudopotential testing (I) The all-electron (ae.sh) and pseudo-test (pt.sh) scripts: run script input pseudopotential to test
10 Pseudopotential testing (II) A transferable pseudo will reproduce the AE energy levels and wave functions in arbitrary environments
11 SUMMER SCHOOL ON COMPUTATIONAL MATERIALS SCIENCE Pseudopotential testing (III) Compute the energy of two different configurations Compute the difference in energy For the pseudopotential to be transferible: 3s2 3p2 (reference) 3s2 3p1 3d1 3s1 3p3 3s1 3p2 3d1 3s0 3p3 3d1
12 Large core-valence overlap Errors due to non-linearity of XC-potential
13 ~ V Non-linear core corrections Standard pseudopotential unscreening: Valence charge only V ps However xc V xc = V xc = V [ ρ ] + v ps scr [ ρ ]( r) V [ ρ ]( r) V [ ρ ]( r) v { V [ ρ + ρ ] V [ ρ ]} xc v c H [ ρ + ρ ]( r) V [ ρ ]( r) + V [ ρ ]( r) v c xc Core-correction v v xc v xc xc v c Keep core charge in pseudopotential generation V ps ps ( r) = V [ ρ + ρ ]( r) V [ ρ ]( r) V [ ρ + ρ ]( r) scr v c H v xc v c
14 New flag PCC input file
15 Pseudo-core & pseudo-valence charge
16 Smooth Fourier Transform The real-space grid required fineness depends on how you define the pseudopotential. The meshcutoff parameter can be determined from the Fourier Transform.
17 Basis generation
18 Key for f linear-scaling scaling: : LOCALITY Large system x2 Divide and Conquer W. Yang, Phys. Rev. Lett. 66, 1438 (1992)
19 Very efficient Atomic Orbitals Lack of systematic for convergence Main features: φ Ilm ( r )=R Il (r I )Y lm ( ˆ r I ) r = r R r I I Numerical Atomic Orbitals (NAOs): Numerical solution of the KS Hamiltonian for the isolated pseudoatom with the same approximations (xc, pseudos) as for the condensed system Size or number of functions Range of localization of these functions Shape or functional form used.
20 Basis Size (I) Depends on the required accuracy and available computational power Quick and dirty calculations Minimal basis set (single-ζ; SZ) Highly converged calculations Complete multiple-ζ + Polarization + Diffuse orbitals
21 Basis Size (II): improving One single radial function per angular momentum shell occupied in the free-atom. Single-ζ (minimal or SZ) Improving the quality? Radial flexibilization: Add more than one radial function within the same angular momentum shell Multiple-ζ Angular flexibilization: Add shells of different angular momentum Polarization
22 Examples
23 Basis Size (III): Polarization Perturbative polarization Apply a small E field to the orbital we want to polarize E Atomic polarization Solve Schrödinger equation for higher angular momentum unbound in the free atom require short cut offs s s+p Si 3d orbitals E. Artacho et al, Phys. Stat. Sol. (b), 215, 809 (1999)
24 Basis Size (IV): Convergence Cohesion curves Bulk Si PW and NAO convergence J. Junquera et. al. PRB 64, (2001).
25 Range (I): How to get sparsity for O(n) Neglecting interactions below a tolerance or beyond some scope of neighbours numerical instablilities for high tolerances. Strictly localized atomic orbitals (zero beyond a given cutoff radius, r c ) Accuracy and computational efficiency depend on the range of the atomic orbitals. Way to define all the cutoff radii in a balanced way.
26 Range (II): Energy Shift Easy approach to define the cutoff radii for the NAOs: 1 2r d dr 2 2 r l( l + 1) + 2 2r + V l ( r) φ l ( r) = ( ε l + δε l ) φl ( r) E. Artacho et al. Phys. Stat. Solidi (b) 215, 809 (1999) A single parameter for all cutoff radii Fireballs O. F. Sankey & D. J. Niklewski, Phys. Rev. B 40, 3979 (1989) BUT, a different cutoff radius for each orbital
27 Range (II): Convergence Bulk Si equal s, p orbitals radii J. Soler et al, J. Phys: Condens. Matter, 14, 2745 (2002)
28 Shape (I) The radial function shape is mainly determined by the pseudopotential. Extra parameters can be introduced to add flexibility: δq : extra charge per atomic specie. Confinement : imposed separately for each angular momentum shell.
29 Shape (II) Soft confinement (J. Junquera et al, Phys. Rev. B 64, (01) ) Shape of the optimal 3s orbital of Mg in MgO for different schemes Corresponding optimal confinement potential Better variational basis sets Removes the discontinuity of the derivative
30 PAO.Basis (I) # shells Add polarization 4s 3d 1 st ζ %block PAO.Basis # Define Basis set Cu 2 # Species label, number of l-shells n=4 0 2 P 1 # n, l, Nzeta, Polarization, NzetaPol Rc n=3 2 2 # n, l, Nzeta Rc nd ζ %endblock PAO.Basis
31 PAO.Basis (II): new generation charge %block PAO.Basis # Define Basis set Cu n=4 0 1 E V c r c n=4 1 1 E Polarization orbital n=3 2 1 E %endblock PAO.Basis
32 Procedure 1. Check the difference in energies involved in your problem 2. For semiquantitative results and general trends use SZ 3. Improve the basis: Automatic DZP (Split Valence & Perturbative Polarization): High quality for most systems Good valence between well converged results & computational cost Standard Rule of thumb in Quantum Chemistry: «a basis should always be doubled before being polarized». 4. Functional optimization of the basis
33 Pseudos & Basis repository Pseudopotentials and basis sets available in the SIESTA web page: Uploaded by users input files to generate them Plots of the radial functions Documentation of the tests done Author s contact information The PAO is pseudopotential-dependent. Check also in the user s mailing list.
Atomic orbitals of finite range as basis sets. Javier Junquera
Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles
More informationBasis sets for SIESTA. Emilio Artacho. Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge
Basis sets for SIESTA Emilio Artacho Nanogune, Ikerbasque & DIPC, San Sebastian, Spain Cavendish Laboratory, University of Cambridge Solving: Basis set Expand in terms of a finite set of basis functions
More informationAtomic orbitals of finite range as basis sets. Javier Junquera
Atomic orbitals of finite range as basis sets Javier Junquera Most important reference followed in this lecture in previous chapters: the many body problem reduced to a problem of independent particles
More informationNorm-conserving pseudopotentials and basis sets in electronic structure calculations. Javier Junquera. Universidad de Cantabria
Norm-conserving pseudopotentials and basis sets in electronic structure calculations Javier Junquera Universidad de Cantabria Outline Pseudopotentials Why pseudopotential approach is useful Orthogonalized
More informationSystematic convergence for realistic projects Fast versus accurate
Systematic convergence for realistic projects Fast versus accurate Daniel Sánchez-Portal Centro de Física de Materiales, Centro Mixto CSIC- UPV/EHU,San Sebastián, Spain Email: sqbsapod@sc.ehu.es Thanks
More informationExercises on basis set generation Full control on the definition of the basis set functions: the PAO.Basis block. Javier Junquera
Exercises on basis set generation Full control on the definition of the basis set functions: the PAO.Basis block Javier Junquera 3.50 0.95 # scalefactor(izeta=1,nzeta) (opt) # scalefactor(izeta=1,nzeta)
More informationThe Siesta program for electronic structure simulations
Université de Liège The Siesta program for electronic structure simulations Javier Junquera Départament de Physique Université de Liège Linear Scaling CPU load ~ N 3 Early 90 s ~ N ~ 100 N (# atoms) G.
More informationThe Plane-Wave Pseudopotential Method
Hands-on Workshop on Density Functional Theory and Beyond: Computational Materials Science for Real Materials Trieste, August 6-15, 2013 The Plane-Wave Pseudopotential Method Ralph Gebauer ICTP, Trieste
More information2. TranSIESTA 1. SIESTA. DFT In a Nutshell. Introduction to SIESTA. Boundary Conditions: Open systems. Greens functions and charge density
1. SIESTA DFT In a Nutshell Introduction to SIESTA Atomic Orbitals Capabilities Resources 2. TranSIESTA Transport in the Nanoscale - motivation Boundary Conditions: Open systems Greens functions and charge
More informationDFT in practice : Part II. Ersen Mete
pseudopotentials Department of Physics Balıkesir University, Balıkesir - Turkey August 13, 2009 - NanoDFT 09, İzmir Institute of Technology, İzmir Outline Pseudopotentials Basic Ideas Norm-conserving pseudopotentials
More informationLDA+U. How to run a simulation within LDA+U framework in SIESTA. Javier Junquera. Daniel Sánchez-Portal
LDA+U How to run a simulation within LDA+U framework in SIESTA Javier Junquera Daniel Sánchez-Portal Transition metal monoxides are prototypes of highly correlated materials Paradigmatic example: NiO Above
More informationThe Linearized Augmented Planewave (LAPW) Method
The Linearized Augmented Planewave (LAPW) Method David J. Singh Oak Ridge National Laboratory E T [ ]=T s [ ]+E ei [ ]+E H [ ]+E xc [ ]+E ii {T s +V ks [,r]} I (r)= i i (r) Need tools that are reliable
More informationWhy use pseudo potentials?
Pseudo potentials Why use pseudo potentials? Reduction of basis set size effective speedup of calculation Reduction of number of electrons reduces the number of degrees of freedom For example in Pt: 10
More informationHow to generate a pseudopotential with non-linear core corrections
How to generate a pseudopotential with non-linear core corrections 14 12 AE core charge AE valence charge PS core charge PS valence charge 10 8 6 4 2 Objectives 0 0 0.5 1 1.5 2 2.5 3 Check whether the
More informationPseudopotential generation and test by the ld1.x atomic code: an introduction
and test by the ld1.x atomic code: an introduction SISSA and DEMOCRITOS Trieste (Italy) Outline 1 2 3 Spherical symmetry - I The Kohn and Sham (KS) equation is (in atomic units): [ 1 ] 2 2 + V ext (r)
More informationOn-the-fly pseudopotential generation in CASTEP
On-the-fly pseudopotential generation in CASTEP Chris J. Pickard School of Physics and Astronomy, University of St Andrews St Andrews, KY16 9SS, United Kingdom September 13, 2006 A quick tutorial Default
More informationConstruction pseudo-potentials for the projector augmentedwave. CSCI699 Assignment 2
Construction pseudo-potentials for the projector augmentedwave Method CSCI699 Assignment 2 Make Your Own PAW Pseudopotentials I. Briefly describe Pseudo-wavefunction (RRKJ2) Local pseudo-potential Non-local
More informationDFT / SIESTA algorithms
DFT / SIESTA algorithms Javier Junquera José M. Soler References http://siesta.icmab.es Documentation Tutorials Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang
More informationTheory of Pseudopotentials. Outline of Talk
Theory of Pseudopotentials David Vanderbilt Rutgers University Outline of Talk Introduction Motivation Basic Idea History and Terminology First-Principles Pseudopotentials Construction Scattering Properties
More informationHow to run SIESTA. Introduction to input & output files
How to run SIESTA Introduction to input & output files Linear-scaling DFT based on Numerical Atomic Orbitals (NAOs) Born-Oppenheimer DFT Pseudopotentials Numerical atomic orbitals relaxations, MD, phonons.
More informationPseudopotentials for hybrid density functionals and SCAN
Pseudopotentials for hybrid density functionals and SCAN Jing Yang, Liang Z. Tan, Julian Gebhardt, and Andrew M. Rappe Department of Chemistry University of Pennsylvania Why do we need pseudopotentials?
More informationThe Projector Augmented Wave method
The Projector Augmented Wave method Advantages of PAW. The theory. Approximations. Convergence. 1 The PAW method is... What is PAW? A technique for doing DFT calculations efficiently and accurately. An
More informationElectron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals
Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals Richard M. Martin UIUC Lecture at Summer School Hands-on introduction to Electronic Structure Materials Computation Center University
More informationReferences. Documentation Manuals Tutorials Publications
References http://siesta.icmab.es Documentation Manuals Tutorials Publications Atomic units e = m e = =1 atomic mass unit = m e atomic length unit = 1 Bohr = 0.5292 Ang atomic energy unit = 1 Hartree =
More informationJune 17, Curtin University of Technology, Perth
U s e r s G u i d e S I E S T A 2.0.2-c11 June 17, 2007 Emilio Artacho Julian D. Gale Alberto García Javier Junquera Richard M. Martin Pablo Ordejón Daniel Sánchez-Portal José M. Soler University of Cambridge
More informationPseudo potential exercises
Pseudo potential exercises Johan M. Carlsson Fritz-Haber-Institut der Max-Planck-Gesellschaft D-14195 Berlin Introduction Castep contains an "on the fly" OTF-pseudo potential generator that can be used
More informationIs Band better than ADF? The core of the issue
Is Band better than ADF? The core of the issue Introduction Introduction What is Band? Introduction What is Band? Slightly different basis set for ADF/Band Introduction What is Band? Slightly different
More informationBasic introduction to SIESTA
Basic introduction to SIESTA Emilio Artacho Department of Earth Sciences University of Cambridge 1 The method Linear-scaling DFT based on NAOs (Numerical Atomic Orbitals) P. Ordejon, E. Artacho & J. M.
More informationLecture on First-principles Computations (14): The Linear Combination of Atomic Orbitals (LCAO) Method
Lecture on First-principles Computations (14): The Linear Combination of Atomic Orbitals (LCAO) Method 任新国 (Xinguo Ren) 中国科学技术大学量子信息重点实验室 Key Laboratory of Quantum Information, USTC Hefei, 2016.11.11 Recall:
More informationPseudopotentials: design, testing, typical errors
Pseudopotentials: design, testing, typical errors Kevin F. Garrity Part 3 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Testing Basis Sets
More informationCHAPTER 3 WIEN2k. Chapter 3 : WIEN2k 50
CHAPTER 3 WIEN2k WIEN2k is one of the fastest and reliable simulation codes among computational methods. All the computational work presented on lanthanide intermetallic compounds has been performed by
More informationATOM User Manual. Version 3.2.8, 29 April Alberto García ICMAB-CSIC, Barcelona vfill
ATOM User Manual Version 3.2.8, 29 April 2008 Alberto García ICMAB-CSIC, Barcelona albertog@icmab.es vfill Contents 1 PREFACE 2 2 A PRIMER ON AB-INITIO PSEUDOPOTENTIALS 2 3 COMPILING THE PROGRAM 3 4 USING
More informationPseudopotential methods for DFT calculations
Pseudopotential methods for DFT calculations Lorenzo Paulatto Scuola Internazionale Superiore di Studi Avanzati and CNR-INFM DEMOCRITOS National Simulation Center Tieste Italy July 9, 2008 Outline pseudopotential
More informationDensity Functional Theory (DFT) modelling of C60 and
ISPUB.COM The Internet Journal of Nanotechnology Volume 3 Number 1 Density Functional Theory (DFT) modelling of C60 and N@C60 N Kuganathan Citation N Kuganathan. Density Functional Theory (DFT) modelling
More informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationALMA: All-scale predictive design of heat management material structures
ALMA: All-scale predictive design of heat management material structures Version Date: 2015.11.13. Last updated 2015.12.02 Purpose of this document: Definition of a data organisation that is applicable
More informationThe Linearized Augmented Planewave (LAPW) Method (WIEN2k, ELK, FLEUR)
The Linearized Augmented Planewave (LAPW) Method (WIEN2k, ELK, FLEUR) David J. Singh Oak Ridge National Laboratory E T [ρ]=t s [ρ]+e ei [ρ]+e H [ρ]+e xc [ρ]+e ii {T s +V ks [ρ,r]}ϕ I (r)=ε i ϕ i (r) Please
More informationBasic introduction to running Siesta. Eduardo Anglada Siesta foundation-uam-nanotec
Basic introduction to running Siesta Eduardo Anglada Siesta foundation-uam-nanotec eduardo.anglada@uam.es eduardo.anglada@nanotec.es Our method Linear-scaling DFT based on NAOs (Numerical Atomic Orbitals)
More informationChapter 3. The (L)APW+lo Method. 3.1 Choosing A Basis Set
Chapter 3 The (L)APW+lo Method 3.1 Choosing A Basis Set The Kohn-Sham equations (Eq. (2.17)) provide a formulation of how to practically find a solution to the Hohenberg-Kohn functional (Eq. (2.15)). Nevertheless
More informationThe LDA+U method: a primer and implementation within SIESTA
The LDA+U method: a primer and implementation within SIESTA Daniel Sánchez-Portal Thanks to Javier Junquera, Sampsa Riikonen and Eduardo Anglada Source of the failure of LDA to describe Mott insulators
More informationKey word: DensityFunctional Theory, Endohedral Energy gap, Electonic properties.
First Principle Studies of Electronic Properties of Nitrogen-doped Endohedral Fullerene 1 M. R. Benam, 2 N. Shahtahmasbi, 1 H. Arabshahi and 1 Z.Zarei 1 Department of Physics, Payame Noor University, P.
More information2.4. Quantum Mechanical description of hydrogen atom
2.4. Quantum Mechanical description of hydrogen atom Atomic units Quantity Atomic unit SI Conversion Ang. mom. h [J s] h = 1, 05459 10 34 Js Mass m e [kg] m e = 9, 1094 10 31 kg Charge e [C] e = 1, 6022
More informationProjector augmented wave Implementation
Projector augmented wave Implementation Peter. E. Blöchl Institute for Theoretical Physics Clausthal University of Technology, Germany http://www.pt.tu-clausthal.de/atp/ 1 = Projector augmented wave +
More informationBasis Sets and Basis Set Notation
Chemistry 46 Fall 215 Dr. Jean M. Standard November 29, 217 Basis Sets and Basis Set Notation Using the LCAO-MO approximation, molecular orbitals can be represented as linear combinations of atomic orbitals,
More informationComputational Modeling Module: Electronic Band Structure Calculations MSE 410, 10/03 & 10/5
Computational Modeling Module: Electronic Band Structure Calculations MSE 410, 10/03 & 10/5 Learning Outcomes and Topics Learning outcomes: Get familiar with density functional theory (DFT) and simulation
More informationData structure of OpenMX and ADPACK
Data structure of OpenMX and ADPACK Structure of ADPACK Solving the 1D Dirac equation Pseudopotentials Localized basis functions Structure of OpenMX Input_std() Total energy Data structure for parallelization
More informationDensity-functional calculations of defect formation energies using the supercell method: Brillouin-zone sampling
Density-functional calculations of defect formation energies using the supercell method: Brillouin-zone sampling Jihye Shim and Eok-Kyun Lee Department of Chemistry and School of Molecular Science (BK21),
More information1 Construction of norm-conserving semi-local pseudopotentials for Si
1 Construction of norm-conserving semi-local pseudopotentials for Si As discussed in class, it is desirable to replace the effective interaction of the valence electrons with the ionic core, i.e. nucleus
More informationarxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 8 Nov 2001
The Siesta method for ab initio order-n materials simulation arxiv:cond-mat/111138v1 [cond-mat.mtrl-sci] 8 Nov 1 José M. Soler, 1 Emilio Artacho, Julian D. Gale, 3 Alberto García, 4 Javier Junquera, 1,5
More informationExplaining the apparent arbitrariness of the LDA-1/2 self-energy. correction method applied to purely covalent systems
Explaining the apparent arbitrariness of the LDA-1/2 self-energy correction method applied to purely covalent systems Kan-Hao Xue, 1,2 Leonardo R. C. Fonseca, 3 and Xiang-Shui Miao 1,2 1 School of Optical
More information1. Hydrogen atom in a box
1. Hydrogen atom in a box Recall H atom problem, V(r) = -1/r e r exact answer solved by expanding in Gaussian basis set, had to solve secular matrix involving matrix elements of basis functions place atom
More informationCP2K: the gaussian plane wave (GPW) method
CP2K: the gaussian plane wave (GPW) method Basis sets and Kohn-Sham energy calculation R. Vuilleumier Département de chimie Ecole normale supérieure Paris Tutorial CPMD-CP2K CPMD and CP2K CPMD CP2K http://www.cpmd.org
More informationECE440 Nanoelectronics. Lecture 07 Atomic Orbitals
ECE44 Nanoelectronics Lecture 7 Atomic Orbitals Atoms and atomic orbitals It is instructive to compare the simple model of a spherically symmetrical potential for r R V ( r) for r R and the simplest hydrogen
More informationHow to compute the projected density of states (PDOS) Javier Junquera
How to compute the projected density of states (PDOS) Javier Junquera Density Of States (DOS) the number of one-electron levels between E and E + de SrTiO 3 bulk Units: (Energy) -1 Projected Density Of
More informationNumerical atomic basis orbitals from H to Kr
PHYSICAL REVIEW B 69, 195113 2004 Numerical atomic basis orbitals from H to Kr T. Ozaki Research Institute for Computational Sciences (RICS), National Institute of Advanced Industrial Science and Technology
More informationElectronic structure calculations with GPAW. Jussi Enkovaara CSC IT Center for Science, Finland
Electronic structure calculations with GPAW Jussi Enkovaara CSC IT Center for Science, Finland Basics of density-functional theory Density-functional theory Many-body Schrödinger equation Can be solved
More informationGaussian Basis Sets for Solid-State Calculations
Gaussian Basis Sets for Solid-State Calculations K. Doll Molpro Quantum Chemistry Software Institute of Theoretical Chemistry, D-70569 Stuttgart, Germany MW-MSSC 2017, Minneapolis, July 10, 2017 Introduction
More informationEfficient projector expansion for the ab initio LCAO method
PHYSICAL REVIEW B 72, 045121 2005 Efficient projector expansion for the ab initio LCAO method T. Ozaki Research Institute for Computational Sciences (RICS), National Institute of Advanced Industrial Science
More informationAb initio study of Mn doped BN nanosheets Tudor Luca Mitran
Ab initio study of Mn doped BN nanosheets Tudor Luca Mitran MDEO Research Center University of Bucharest, Faculty of Physics, Bucharest-Magurele, Romania Oldenburg 20.04.2012 Table of contents 1. Density
More informationFinal Exam Tuesday, May 8, 2012 Starting at 8:30 a.m., Hoyt Hall Duration: 2h 30m
Final Exam Tuesday, May 8, 2012 Starting at 8:30 a.m., Hoyt Hall. ------------------- Duration: 2h 30m Chapter 39 Quantum Mechanics of Atoms Units of Chapter 39 39-1 Quantum-Mechanical View of Atoms 39-2
More informationPseudopotentials: design, testing, typical errors
Pseudopotentials: design, testing, typical errors Kevin F. Garrity Part 1 National Institute of Standards and Technology (NIST) Uncertainty Quantification in Materials Modeling 2015 Parameter free calculations.
More informationS I E S T A 4.0. June 23, CIC-Nanogune and University of Cambridge
U s e r s G u i d e S I E S T A 4.0 June 23, 2016 Emilio Artacho José María Cela Julian D. Gale Alberto García Javier Junquera Richard M. Martin Pablo Ordejón Daniel Sánchez-Portal José M. Soler CIC-Nanogune
More informationAB INITIO STUDY OF NANO STRUCTURED FUNCTIONALIZED GRAPHENE WITH 30C ATOMS
International Journal of Science, Environment and Technology, Vol. 1, No 3, 2012, 108-112 AB INITIO STUDY OF NANO STRUCTURED FUNCTIONALIZED GRAPHENE WITH 30C ATOMS Naveen Kumar* and Jyoti Dhar Sharma Deptt.
More informationIntroduction to Heisenberg model. Javier Junquera
Introduction to Heisenberg model Javier Junquera Most important reference followed in this lecture Magnetism in Condensed Matter Physics Stephen Blundell Oxford Master Series in Condensed Matter Physics
More informationPseudopotential. Meaning and role
Pseudopotential. Meaning and role Jean-Pierre Flament jean-pierre.flament@univ-lille.fr Laboratoire de Physique des Lasers, Atomes et Molécules (PhLAM) Université de Lille-Sciences et technologies MSSC2018
More informationAtomic Structure, Periodic Table, and Other Effects: Chapter 8 of Rex and T. Modern Physics
Atomic Structure, Periodic Table, and Other Effects: Chapter 8 of Rex and T Modern Physics 11/16 and 11/19/2018 1 Introduction In Chapter 7, we studied the hydrogen atom. What about other elements, e.g.,
More informationOpen-Source Pseudopotential Interface/Unification Module (OPIUM): The Basic Ins and Outs of Operation
Open-Source Pseudopotential Interface/Unification Module (OPIUM): The Basic Ins and Outs of Operation Irene K. Metz, Joseph W. Bennett and Sara E. Mason (Dated: May 31, 2018) Learning Objectives 1. Determine
More informationAn Inverse Mass Expansion for Entanglement Entropy. Free Massive Scalar Field Theory
in Free Massive Scalar Field Theory NCSR Demokritos National Technical University of Athens based on arxiv:1711.02618 [hep-th] in collaboration with Dimitris Katsinis March 28 2018 Entanglement and Entanglement
More informationGW quasiparticle energies
Chapter 4 GW quasiparticle energies Density functional theory provides a good description of ground state properties by mapping the problem of interacting electrons onto a KS system of independent particles
More informationAP Chemistry - Problem Drill 10: Atomic Structures and Electron Configuration
AP Chemistry - Problem Drill 10: Atomic Structures and Electron Configuration No. 1 of 10 Instructions: (1) Read the problem statement and answer choices carefully (2) Work the problems on paper as 1.
More informationTight-Binding Model of Electronic Structures
Tight-Binding Model of Electronic Structures Consider a collection of N atoms. The electronic structure of this system refers to its electronic wave function and the description of how it is related to
More informationMulti-reference Density Functional Theory. COLUMBUS Workshop Argonne National Laboratory 15 August 2005
Multi-reference Density Functional Theory COLUMBUS Workshop Argonne National Laboratory 15 August 2005 Capt Eric V. Beck Air Force Institute of Technology Department of Engineering Physics 2950 Hobson
More informationChem 442 Review for Exam 2. Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative (3D) components.
Chem 44 Review for Exam Hydrogenic atoms: The Coulomb energy between two point charges Ze and e: V r Ze r Exact separation of the Hamiltonian of a hydrogenic atom into center-of-mass (3D) and relative
More informationPhysics 606, Quantum Mechanics, Final Exam NAME ( ) ( ) + V ( x). ( ) and p( t) be the corresponding operators in ( ) and x( t) : ( ) / dt =...
Physics 606, Quantum Mechanics, Final Exam NAME Please show all your work. (You are graded on your work, with partial credit where it is deserved.) All problems are, of course, nonrelativistic. 1. Consider
More informationLS coupling. 2 2 n + H s o + H h f + H B. (1) 2m
LS coupling 1 The big picture We start from the Hamiltonian of an atomic system: H = [ ] 2 2 n Ze2 1 + 1 e 2 1 + H s o + H h f + H B. (1) 2m n e 4πɛ 0 r n 2 4πɛ 0 r nm n,m Here n runs pver the electrons,
More informationConquest order N ab initio Electronic Structure simulation code for quantum mechanical modelling in large scale
Fortran Expo: 15 Jun 2012 Conquest order N ab initio Electronic Structure simulation code for quantum mechanical modelling in large scale Lianheng Tong Overview Overview of Conquest project Brief Introduction
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz November 29, 2010 The self interaction error and its correction Perdew-Zunger SIC Average-density approximation Weighted density approximation
More informationDesigned nonlocal pseudopotentials for enhanced transferability
PHYSICAL REVIEW B VOLUME 59, NUMBER 19 15 MAY 1999-I Designed nonlocal pseudopotentials for enhanced transferability Nicholas J. Ramer and Andrew M. Rappe Department of Chemistry and Laboratory for Research
More informationElectrons in Crystals. Chris J. Pickard
Electrons in Crystals Chris J. Pickard Electrons in Crystals The electrons in a crystal experience a potential with the periodicity of the Bravais lattice: U(r + R) = U(r) The scale of the periodicity
More informationElectronic structure, plane waves and pseudopotentials
Electronic structure, plane waves and pseudopotentials P.J. Hasnip Spectroscopy Workshop 2009 We want to be able to predict what electrons and nuclei will do from first principles, without needing to know
More informationLocalized basis methods Theory and implementations. Introduction of OpenMX Implementation of OpenMX
Localized basis methods Theory and implementations Introduction of OpenMX Implementation of OpenMX Total energy Pseudopontials Basis functions Self-consistency Δ-gauge Taisuke Ozaki (ISSP, Univ. of Tokyo)
More informationarxiv:cond-mat/ v1 17 May 1995
Projection of plane-wave calculations into atomic orbitals Daniel Sanchez-Portal, Emilio Artacho, and Jose M. Soler Instituto de Ciencia de Materiales Nicolás Cabrera and Departamento de Física de la Materia
More informationAtomic and molecular physics Revision lecture
Atomic and molecular physics Revision lecture Answer all questions Angular momentum J`2 ` J z j,m = j j+1 j,m j,m =m j,m Allowed values of mgo from j to +jin integer steps If there is no external field,
More informationThe 3 dimensional Schrödinger Equation
Chapter 6 The 3 dimensional Schrödinger Equation 6.1 Angular Momentum To study how angular momentum is represented in quantum mechanics we start by reviewing the classical vector of orbital angular momentum
More informationCOMPUTATIONAL TOOL. Fig. 4.1 Opening screen of w2web
CHAPTER -4 COMPUTATIONAL TOOL Ph.D. Thesis: J. Maibam CHAPTER: 4 4.1 The WIEN2k code In this work, all the calculations presented are performed using the WIEN2k software package (Blaha et al., 2001). The
More informationOPENATOM for GW calculations
OPENATOM for GW calculations by OPENATOM developers 1 Introduction The GW method is one of the most accurate ab initio methods for the prediction of electronic band structures. Despite its power, the GW
More informationarxiv: v1 [cond-mat.mes-hall] 13 Feb 2012
Controlling Band Gap in Silicene Monolayer Using External Electric Field C. Kamal 1 arxiv:1202.2636v1 [cond-mat.mes-hall] 13 Feb 2012 1 Indus Synchrotrons Utilization Division, Raja Ramanna Centre for
More informationAn Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method
Fakultät für Mathematik und Naturwissenschaften - Lehrstuhl für Physikalische Chemie I / Theoretische Chemie An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method Jan-Ole Joswig
More informationTwo implementations of the Projector Augmented Wave (PAW) formalism
Introduction The tools available for detailed first-principles studies of materials have benefited enormously from the development of several international collaborations engaged in developing open source
More informationChapter 5. Atomic spectra
Atomic spectra Sommerfelds relativistic model Sommerfeld succeeded partially in explaining fine structure by extending Bohr Theory i) He allowed the possibility of elliptical orbits for the electrons in
More informationGoal: find Lorentz-violating corrections to the spectrum of hydrogen including nonminimal effects
Goal: find Lorentz-violating corrections to the spectrum of hydrogen including nonminimal effects Method: Rayleigh-Schrödinger Perturbation Theory Step 1: Find the eigenvectors ψ n and eigenvalues ε n
More informationKey concepts in Density Functional Theory (I) Silvana Botti
From the many body problem to the Kohn-Sham scheme European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Temporary Address: Centre
More informationThis is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry.
1 Computational Chemistry (Quantum Chemistry) Primer This is a very succinct primer intended as supplementary material for an undergraduate course in physical chemistry. TABLE OF CONTENTS Methods...1 Basis
More informationSection 3 Electronic Configurations, Term Symbols, and States
Section 3 Electronic Configurations, Term Symbols, and States Introductory Remarks- The Orbital, Configuration, and State Pictures of Electronic Structure One of the goals of quantum chemistry is to allow
More informationPhys 622 Problems Chapter 5
1 Phys 622 Problems Chapter 5 Problem 1 The correct basis set of perturbation theory Consider the relativistic correction to the electron-nucleus interaction H LS = α L S, also known as the spin-orbit
More informationU s e r s G u i d e S I E S T A 3.0-rc2 July 5, 2010 Emilio Artacho José María Cela Julian D. Gale Alberto García Javier Junquera Richard M. Martin Pablo Ordejón Daniel Sánchez-Portal José M. Soler University
More informationHund s rule for monopole harmonics, or why the composite fermion picture works
PERGAMON Solid State Communications 110 (1999) 45 49 Hund s rule for monopole harmonics, or why the composite fermion picture works Arkadiusz Wójs*, John J. Quinn The University of Tennessee, Knoxville,
More informationKondo effect in multi-level and multi-valley quantum dots. Mikio Eto Faculty of Science and Technology, Keio University, Japan
Kondo effect in multi-level and multi-valley quantum dots Mikio Eto Faculty of Science and Technology, Keio University, Japan Outline 1. Introduction: next three slides for quantum dots 2. Kondo effect
More informationLABELING ELECTRONS IN ATOMS
Date: Name: LABELING ELECTRONS IN ATOMS The location of each electron in an atom is determined by a few different factors. Each factor is represented by a QUANTUM NUMBER. Prediction: What do you think
More information8.1 The hydrogen atom solutions
8.1 The hydrogen atom solutions Slides: Video 8.1.1 Separating for the radial equation Text reference: Quantum Mechanics for Scientists and Engineers Section 10.4 (up to Solution of the hydrogen radial
More information