Fireball.

Size: px
Start display at page:

Download "Fireball."

Transcription

1 Fireball

2 SIESTA Fireball atomic-like atomic-like Numerical Numerical basis basis set set Precomputed 2 and 3 center integrals Precomputed 2 and 3 center integrals atomic-like atomic-like Numerical Numerical basis basis set set Density and Potential on grid Density and Potential on grid Fast Fast Fourier Fourier Transform Transform to to solve solve Poisson's Poisson's equation equation saved saved to to disk disk before before computation computation Kohn-Sham or McWEDA SCF-loop Kohn-Sham or McWEDA SCF-loop ρρ VVhhaarrttrreeee Kohn-Sham Kohn-Sham is is like like SIESTA SIESTA (on (on grid) grid) more more precise precise calculations calculations McWEDA McWEDA use use atomic atomic charges charges faster, faster, approximative approximative

3 Pseudopotentials "begin" Basis Fuctions pseudo-atomic orbitals solutions with pseudopotential in finite box "create" Interaction integrals Vμνij (r) = < μ (R) V[ ρij ] ν(r+r) > "fireball-prog" System of interest

4 McWEDA Multi-center Weighted Exchangecorrelation Density Approximation 4-centre integrals Vμνij (r,r',r'') r',r'' = < μ (0) V [ ρij ( r', r'') ] ν(r) > r' 0 3-centre integrals Vμνs (r,r') r' = < μ (0) V'[ Charge s (r') ] ν(r) > r'' r'' r i,j... index of atomic orbitals i> iatom, n, L, m > 0 effective charge r s... index of atomic shell Charge(s) Charge( iatom, n, L)... Lowdin, Mulliken...

5 Fireball2Smeagol

6 Fireball SCF Fireball2Smeagol SCF Real Real space space Hamiltonian Hamiltonian Real Real space space Hamiltonian Hamiltonian k-space k-space Hamiltonian(k) Hamiltonian(k) k-space k-space Hamiltonian(k) Hamiltonian(k) Diagonalization Diagonalization Eigenstates Eigenstates (k) (k) Smeagol Smeagol contour contour integration integration Density Density Matrix Matrix (k) (k) Real Real space space Density Density Matrix Matrix Kohn-Sham Kohn-Sham Grid Grid density density Grid Grid potantial potantial Density Density Matrix Matrix (k) (k) Real Real space space Density Density Matrix Matrix McWEDA McWEDA or or CHARGES CHARGES Kohn-Sham Kohn-Sham Grid Grid density density Grid Grid potantial potantial or or McWEDA McWEDA CHARGES CHARGES

7 Basic scheme of Fireball2Smeagol Run PRINCIPAL LAYER directory SCF SCF CHARGES export export leads leads ELECTRODE.left ELECTRODE.right CHARGES.left CHARGES.right EXTENDED MOLECULE directory SCF SCF H.scf Smeagol Smeagol run run

8 BULK calculation fireball.in geometry ansewer.bas element fireball.in x y z cel.lvs lattice input.kpts k-points in x,y,z (for bulk) SCF export leads Output MOLECULE.kpts k-points in x,y (exported leads) Projection H, S, ρ ELECTRODE H_k (μ,ν) S_k (μ,ν) rho_k (μ,ν)

9 ELECTRODE.left μ ν Re H0 or.right Im H0 Re H1 Im H1 Re S0 Im S0 Re S1 Im S1 Re ρ0 Im ρ0 H0 H1 Principal layer Re ρ1 Im ρ1

10 ansewer.bas Extended molecule calculation smeagol.optional LEADS must be in the same order as in BULK calculation cel.lvs input.kpts Output smeagol.cur... current... I(V) curve smeagol.trc... Transmission... T(E)

11 HL ΣL decimation ΣL HL M decimation HM GM Effective hamiltonian: Green's function: Spectral function: HMR HR ΣR Conductance: ΣR

12 Example 1 H22 quantum dot in non-selfconsistent potential

13 smeagol.optional LUMO 12.0 Å 0.7 Å 15.0 Å 5.0 ev HOMO

14 LUMO 12.0 Å 0.7 Å 15.0 Å 5.0 ev HOMO Why you don't see this state? V bias

15 HOMO LUMO 1.0 Å 11.0 Å 12.0 Å 12.5 Å 13.5 Å 0.7 Å G/G Å V bias

16 H2 quantum dot (exlanation)... HOMO bonding... LUMO antibonding T(E) V = 0 ev V = +5 ev tu nn el lin g forbidden by symmetry V = -5 ev tic t12 ~ <s3 H s1 - s2> =0 ba llis t12 ~ <s3 H s1 + s2> 0

17 antibonding 1.0 Å bonding 11.0 Å 12.0 Å 13.5 Å 15.0 Å 0.7 Å G/G Å V bias

18 antibonding 1.0 Å bonding 11.0 Å 12.0 Å 12.5 Å 13.5 Å 0.7 Å G/G Å V bias

19 Example 2 Breaking Gold Chain with extended hopping

20 extended hopping Original basis function Extended basis function 7Å no overlap interaction.optional This is replaced to H and S matrix in apex region before before construction of Green's function

21 T(E) for 1D Gold chain by breaking 3.0 Å G/G0 d d s p

22 function used for smearing f(r) = 1/ (1+ exp (r - r_start / r_scale r_start = 4.0 r_scale = 0.25 t12 = f. H12 + (1-f). H'12 original extended S-matrix as well ))

23 Example 3 Conductance of molecules

24 Au AuT1 AuT3 AuSH S: AuNH2 H Au :S: H AuB1 AuB3 AuBSH AuBNH2 Au :N H

25 Total energy by pulling of molecules weaker bonds covalent bonds

26

27 I-z curve ev

28 Suplementary

29 Instalation of Smeagol into Fireball We will provide makefile to incorporate Smeagol to Fireball new SCF loop for non-equlibirum calculation module_smeagol.f90... global variables Directories of smeagol solve H for new density

Introduction to DFTB. Marcus Elstner. July 28, 2006

Introduction to DFTB. Marcus Elstner. July 28, 2006 Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD

More information

Electronic Structure Models

Electronic Structure Models Electronic Structure Models Hückel Model (1933) Basic Assumptions: (a) One orbital per atom contributes to the basis set; all orbitals "equal" (b) The relevant integrals involving the Hamiltonian are α

More information

Quantum Chemistry. NC State University. Lecture 5. The electronic structure of molecules Absorption spectroscopy Fluorescence spectroscopy

Quantum Chemistry. NC State University. Lecture 5. The electronic structure of molecules Absorption spectroscopy Fluorescence spectroscopy Quantum Chemistry Lecture 5 The electronic structure of molecules Absorption spectroscopy Fluorescence spectroscopy NC State University 3.5 Selective absorption and emission by atmospheric gases (source:

More information

ELEMENTARY BAND THEORY

ELEMENTARY BAND THEORY ELEMENTARY BAND THEORY PHYSICIST Solid state band Valence band, VB Conduction band, CB Fermi energy, E F Bloch orbital, delocalized n-doping p-doping Band gap, E g Direct band gap Indirect band gap Phonon

More information

Support information: (a) + + (b)

Support information: (a) + + (b) 1 upport information: Mechanism of the toddart-heath bistable rotaxane molecular switch Abstract The supporting information is organized as follows: ection 1 describes the real device and our strategy

More information

Density Functional Theory (DFT) modelling of C60 and

Density 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 information

Tight-Binding Model of Electronic Structures

Tight-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 information

Atomic orbitals of finite range as basis sets. Javier Junquera

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 information

Title. I-V curve? e-e interactions? Conductance? Electrical Transport Through Single Molecules. Vibrations? Devices?

Title. I-V curve? e-e interactions? Conductance? Electrical Transport Through Single Molecules. Vibrations? Devices? Electrical Transport Through Single Molecules Harold U. Baranger, Duke University Title with Rui Liu, San-Huang Ke, and Weitao Yang Thanks to S. Getty, M. Fuhrer and L. Sita, U. Maryland Conductance? I-V

More information

2. TranSIESTA 1. SIESTA. DFT In a Nutshell. Introduction to SIESTA. Boundary Conditions: Open systems. Greens functions and charge density

2. 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 information

Momentum filtering effect in molecular wires

Momentum filtering effect in molecular wires PHYSICAL REVIEW B 70, 195309 (2004) Momentum filtering effect in molecular wires Chao-Cheng Kaun, 1, * Hong Guo, 1 Peter Grütter, 1 and R. Bruce Lennox 1,2 1 Center for the Physics of Materials and Department

More information

How to run SIESTA. Introduction to input & output files

How 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 information

Speed-up of ATK compared to

Speed-up of ATK compared to What s new @ Speed-up of ATK 2008.10 compared to 2008.02 System Speed-up Memory reduction Azafulleroid (molecule, 97 atoms) 1.1 15% 6x6x6 MgO (bulk, 432 atoms, Gamma point) 3.5 38% 6x6x6 MgO (k-point sampling

More information

Pseudopotentials for hybrid density functionals and SCAN

Pseudopotentials 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 information

Time-dependent density functional theory (TDDFT)

Time-dependent density functional theory (TDDFT) Advanced Workshop on High-Performance & High-Throughput Materials Simulations using Quantum ESPRESSO ICTP, Trieste, Italy, January 16 to 27, 2017 Time-dependent density functional theory (TDDFT) Ralph

More information

DFT / SIESTA algorithms

DFT / 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 information

Intro to ab initio methods

Intro to ab initio methods Lecture 2 Part A Intro to ab initio methods Recommended reading: Leach, Chapters 2 & 3 for QM methods For more QM methods: Essentials of Computational Chemistry by C.J. Cramer, Wiley (2002) 1 ab initio

More information

Comparison of various abinitio codes used in periodic calculations

Comparison of various abinitio codes used in periodic calculations Comparison of various abinitio codes used in periodic calculations 1 Prof.P. Ravindran, Department of Physics, Central University of Tamil Nadu, India & Center for Materials Science and Nanotechnology,

More information

3.1 Electron tunneling theory

3.1 Electron tunneling theory Scanning Tunneling Microscope (STM) was invented in the 80s by two physicists: G. Binnig and H. Rorher. They got the Nobel Prize a few years later. This invention paved the way for new possibilities in

More information

References. Documentation Manuals Tutorials Publications

References.   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 information

Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 28, December 08, 2014

Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi. Lecture 28, December 08, 2014 Chem 4502 Introduction to Quantum Mechanics and Spectroscopy 3 Credits Fall Semester 2014 Laura Gagliardi Lecture 28, December 08, 2014 Solved Homework Water, H 2 O, involves 2 hydrogen atoms and an oxygen

More information

LUMO + 1 LUMO. Tómas Arnar Guðmundsson Report 2 Reikniefnafræði G

LUMO + 1 LUMO. Tómas Arnar Guðmundsson Report 2 Reikniefnafræði G Q1: Display all the MOs for N2 in your report and classify each one of them as bonding, antibonding or non-bonding, and say whether the symmetry of the orbital is σ or π. Sketch a molecular orbital diagram

More information

Spin-crossover molecules: puzzling systems for electronic structure methods. Andrea Droghetti School of Physics and CRANN, Trinity College Dublin

Spin-crossover molecules: puzzling systems for electronic structure methods. Andrea Droghetti School of Physics and CRANN, Trinity College Dublin Spin-crossover molecules: puzzling systems for electronic structure methods Andrea Droghetti School of Physics and CRANN, Trinity College Dublin Introduction Introduction Can we read a molecule spin? Can

More information

Toward an O(N) method for the non-equilibrium Green's function calculation

Toward an O(N) method for the non-equilibrium Green's function calculation Toward an O(N) method for the non-equilibrium Green's function calculation Masahiro FUKUDA Institute for Solid State Physics, University of Tokyo, Japan Table of Content Introduction Non equilibrium Green

More information

Time-dependent density functional theory (TDDFT)

Time-dependent density functional theory (TDDFT) 04/05/16 Hands-on workshop and Humboldt-Kolleg: Density-Functional Theory and Beyond - Basic Principles and Modern Insights Isfahan University of Technology, Isfahan, Iran, May 2 to 13, 2016 Time-dependent

More information

arxiv:cond-mat/ v2 [cond-mat.mtrl-sci] 14 Aug 2002

arxiv:cond-mat/ v2 [cond-mat.mtrl-sci] 14 Aug 2002 Density functional method for nonequilibrium electron transport arxiv:cond-mat/65v [cond-mat.mtrl-sci] 4 Aug Mads Brandbyge,, José-Luis Mozos, Pablo Ordejón, Jeremy Taylor, and Kurt Stokbro Mikroelektronik

More information

Dept of Mechanical Engineering MIT Nanoengineering group

Dept of Mechanical Engineering MIT Nanoengineering group 1 Dept of Mechanical Engineering MIT Nanoengineering group » Recap of HK theorems and KS equations» The physical meaning of the XC energy» Solution of a one-particle Schroedinger equation» Pseudo Potentials»

More information

General Physical Chemistry II

General Physical Chemistry II General Physical Chemistry II Lecture 13 Aleksey Kocherzhenko October 16, 2014" Last time " The Hückel method" Ø Used to study π systems of conjugated molecules" Ø π orbitals are treated separately from

More information

Electronic communication through molecular bridges Supporting Information

Electronic communication through molecular bridges Supporting Information Electronic communication through molecular bridges Supporting Information Carmen Herrmann and Jan Elmisz Institute of Inorganic and Applied Chemistry, University of Hamburg, Martin-Luther-King-Platz 6,

More information

SCF calculation on HeH +

SCF calculation on HeH + SCF calculation on HeH + Markus Meuwly Department of Chemistry, University of Basel, Basel, Switzerland Abstract This document describes the main steps involved in carrying out a SCF calculation on the

More information

IFM Chemistry Computational Chemistry 2010, 7.5 hp LAB2. Computer laboratory exercise 1 (LAB2): Quantum chemical calculations

IFM Chemistry Computational Chemistry 2010, 7.5 hp LAB2. Computer laboratory exercise 1 (LAB2): Quantum chemical calculations Computer laboratory exercise 1 (LAB2): Quantum chemical calculations Introduction: The objective of the second computer laboratory exercise is to get acquainted with a program for performing quantum chemical

More information

Molecular Orbital Theory. WX AP Chemistry Chapter 9 Adapted from: Luis Bonilla Abel Perez University of Texas at El Paso

Molecular Orbital Theory. WX AP Chemistry Chapter 9 Adapted from: Luis Bonilla Abel Perez University of Texas at El Paso Molecular Orbital Theory WX AP Chemistry Chapter 9 Adapted from: Luis Bonilla Abel Perez University of Texas at El Paso Molecular Orbital Theory The goal of molecular orbital theory is to describe molecules

More information

Molecular Simulation I

Molecular Simulation I Molecular Simulation I Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences

More information

Strong Correlation Effects in Fullerene Molecules and Solids

Strong Correlation Effects in Fullerene Molecules and Solids Strong Correlation Effects in Fullerene Molecules and Solids Fei Lin Physics Department, Virginia Tech, Blacksburg, VA 2461 Fei Lin (Virginia Tech) Correlations in Fullerene SESAPS 211, Roanoke, VA 1 /

More information

Lecture B6 Molecular Orbital Theory. Sometimes it's good to be alone.

Lecture B6 Molecular Orbital Theory. Sometimes it's good to be alone. Lecture B6 Molecular Orbital Theory Sometimes it's good to be alone. Covalent Bond Theories 1. VSEPR (valence shell electron pair repulsion model). A set of empirical rules for predicting a molecular geometry

More information

Bonding and Physical Properties The Molecular Orbital Theory

Bonding and Physical Properties The Molecular Orbital Theory Bonding and Physical Properties The Molecular Orbital Theory Ø Developed by F. Hund and R. S. Mulliken in 1932 Ø Diagram of molecular energy levels Ø Magnetic and spectral properties Paramagnetic vs. Diamagnetic

More information

Basis 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 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 information

VALENCE Hilary Term 2018

VALENCE Hilary Term 2018 VALENCE Hilary Term 2018 8 Lectures Prof M. Brouard Valence is the theory of the chemical bond Outline plan 1. The Born-Oppenheimer approximation 2. Bonding in H + 2 the LCAO approximation 3. Many electron

More information

Lecture 6. Tight-binding model

Lecture 6. Tight-binding model Lecture 6 Tight-binding model In Lecture 3 we discussed the Krönig-Penny model and we have seen that, depending on the strength of the potential barrier inside the unit cell, the electrons can behave like

More information

Topological insulator with time-reversal symmetry

Topological insulator with time-reversal symmetry Phys620.nb 101 7 Topological insulator with time-reversal symmetry Q: Can we get a topological insulator that preserves the time-reversal symmetry? A: Yes, with the help of the spin degree of freedom.

More information

Time-Dependent Electron Localization Function! (TD-ELF)! GOAL! Time-resolved visualization of the breaking and formation of chemical bonds.!

Time-Dependent Electron Localization Function! (TD-ELF)! GOAL! Time-resolved visualization of the breaking and formation of chemical bonds.! Time-Dependent Electron Localization Function! (TD-ELF)! GOAL! Time-resolved visualization of the breaking and formation of chemical bonds.! How can one give a rigorous mathematical meaning to chemical

More information

Winmostar tutorial Quantum ESPRESSO Spin Polarization V X-Ability Co,. Ltd. 2017/8/8

Winmostar tutorial Quantum ESPRESSO Spin Polarization V X-Ability Co,. Ltd. 2017/8/8 Winmostar tutorial Quantum ESPRESSO Spin Polarization V7.025 X-Ability Co,. Ltd. question@winmostar.com 2017/8/8 Contents I. SCF calculation II. Bands calculation III. Fermi surface 2 Environment setting

More information

Theoretical Chemistry - Level II - Practical Class Molecular Orbitals in Diatomics

Theoretical Chemistry - Level II - Practical Class Molecular Orbitals in Diatomics Theoretical Chemistry - Level II - Practical Class Molecular Orbitals in Diatomics Problem 1 Draw molecular orbital diagrams for O 2 and O 2 +. E / ev dioxygen molecule, O 2 dioxygenyl cation, O 2 + 25

More information

Supplementary Figure S1. STM image of monolayer graphene grown on Rh (111). The lattice

Supplementary Figure S1. STM image of monolayer graphene grown on Rh (111). The lattice Supplementary Figure S1. STM image of monolayer graphene grown on Rh (111). The lattice mismatch between graphene (0.246 nm) and Rh (111) (0.269 nm) leads to hexagonal moiré superstructures with the expected

More information

Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering

Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering Fundamentals and applications of Density Functional Theory Astrid Marthinsen PhD candidate, Department of Materials Science and Engineering Outline PART 1: Fundamentals of Density functional theory (DFT)

More information

DFT EXERCISES. FELIPE CERVANTES SODI January 2006

DFT EXERCISES. FELIPE CERVANTES SODI January 2006 DFT EXERCISES FELIPE CERVANTES SODI January 2006 http://www.csanyi.net/wiki/space/dftexercises Dr. Gábor Csányi 1 Hydrogen atom Place a single H atom in the middle of a largish unit cell (start with a

More information

... surface Green function. Figure 1:

... surface Green function. Figure 1: 1 1 Introduction The program calculates the transmission in the L-C or C-R interface. The system the program can calculate is written in the top part of Fig.1. The effects of semi-infinite electrodes,

More information

Conductance Calculations for Small Molecules

Conductance Calculations for Small Molecules Conductance Calculations for Small Molecules Karsten W. Jacobsen CAMP, Dept. of Physics Technical University of Denmark Outline Framework Density Functional Theory Conductance calculations methodology

More information

Table of Contents. Table of Contents Spin-orbit splitting of semiconductor band structures

Table of Contents. Table of Contents Spin-orbit splitting of semiconductor band structures Table of Contents Table of Contents Spin-orbit splitting of semiconductor band structures Relavistic effects in Kohn-Sham DFT Silicon band splitting with ATK-DFT LSDA initial guess for the ground state

More information

Efficient implementation of the nonequilibrium Green function method for electronic transport calculations

Efficient implementation of the nonequilibrium Green function method for electronic transport calculations PHYSICAL REVIEW B 81, 35116 21 Efficient implementation of the nonequilibrium Green function method for electronic transport calculations Taisuke Ozaki Research Center for Integrated Science (RCIS), Japan

More information

Pseudopotentials and Basis Sets. How to generate and test them

Pseudopotentials and Basis Sets. How to generate and test them Pseudopotentials and Basis Sets How to generate and test them 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

More information

Linking electronic and molecular structure: Insight into aqueous chloride solvation. Supplementary Information

Linking electronic and molecular structure: Insight into aqueous chloride solvation. Supplementary Information Linking electronic and molecular structure: Insight into aqueous chloride solvation Ling Ge, Leonardo Bernasconi, and Patricia Hunt Department of Chemistry, Imperial College London, London SW7 2AZ, United

More information

Apparent reversal of molecular orbitals reveals entanglement

Apparent reversal of molecular orbitals reveals entanglement Apparent reversal of molecular orbitals reveals entanglement Andrea Donarini P.Yu, N. Kocic, B.Siegert, J.Repp University of Regensburg and Shanghai Tech University Entangled ground state Spectroscopy

More information

Key concepts in Density Functional Theory (II)

Key concepts in Density Functional Theory (II) Kohn-Sham scheme and band structures European Theoretical Spectroscopy Facility (ETSF) CNRS - Laboratoire des Solides Irradiés Ecole Polytechnique, Palaiseau - France Present Address: LPMCN Université

More information

Chemistry 543--Final Exam--Keiderling May 5, pm SES

Chemistry 543--Final Exam--Keiderling May 5, pm SES Chemistry 543--Final Exam--Keiderling May 5,1992 -- 1-5pm -- 174 SES Please answer all questions in the answer book provided. Make sure your name is clearly indicated and that the answers are clearly numbered,

More information

Excitation energies of dissociating H 2 : A problematic case for the adiabatic approximation of time-dependent density functional theory

Excitation energies of dissociating H 2 : A problematic case for the adiabatic approximation of time-dependent density functional theory JOURNAL OF CHEMICAL PHYSICS VOLUME 113, NUMBER 19 15 NOVEMBER 2000 Eitation energies of dissociating H 2 : A problematic case for the adiabatic approximation of time-dependent density functional theory

More information

An Approximate DFT Method: The Density-Functional Tight-Binding (DFTB) Method

An 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 information

1. Motivation 1. Motivation 1. Motivation L d R - V + 2. Bulk Transport E( k ) = 2 k 2 2m * ε µ) 2. Bulk Transport 2. Bulk Transport k d = -eeτ 3. Some basic concepts τ l m =vτ l ϕ λ = 2π/k 3. Some basic

More information

Practical Guide to Density Functional Theory (DFT)

Practical Guide to Density Functional Theory (DFT) Practical Guide to Density Functional Theory (DFT) Brad Malone, Sadas Shankar Quick recap of where we left off last time BD Malone, S Shankar Therefore there is a direct one-to-one correspondence between

More information

Electron bands in crystals Pseudopotentials, Plane Waves, Local Orbitals

Electron 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 information

Recent advances in development of single-point orbital-free kinetic energy functionals

Recent advances in development of single-point orbital-free kinetic energy functionals PacifiChem 2010 p. 1/29 Recent advances in development of single-point orbital-free kinetic energy functionals Valentin V. Karasiev vkarasev@qtp.ufl.edu Quantum Theory Project, Departments of Physics and

More information

Vibronic Coupling in Quantum Wires: Applications to Polydiacetylene

Vibronic Coupling in Quantum Wires: Applications to Polydiacetylene Vibronic Coupling in Quantum Wires: Applications to Polydiacetylene An Exhaustively Researched Report by Will Bassett and Cole Johnson Overall Goal In order to elucidate the absorbance spectra of different

More information

Band calculations: Theory and Applications

Band calculations: Theory and Applications Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation

More information

Wavelets for density functional calculations: Four families and three. applications

Wavelets for density functional calculations: Four families and three. applications Wavelets for density functional calculations: Four families and three Haar wavelets Daubechies wavelets: BigDFT code applications Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Interpolating

More information

2.1 Introduction: The many-body problem

2.1 Introduction: The many-body problem Chapter 2 Smeagol: Density Functional Theory and NEGF s 2.1 Introduction: The many-body problem In solid state physics one is interested in systems comprising many atoms, and consequently many electrons.

More information

CP2K: the gaussian plane wave (GPW) method

CP2K: 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 information

6: Plane waves, unit cells, k- points and all that

6: Plane waves, unit cells, k- points and all that The Nuts and Bolts of First-Principles Simulation 6: Plane waves, unit cells, k- points and all that Durham, 6th- 13th December 2001 CASTEP Developers Group with support from the ESF ψ k Network Overview

More information

Valence bond theory accounts, at least qualitatively, for the stability of the covalent bond in terms of overlapping atomic orbitals.

Valence bond theory accounts, at least qualitatively, for the stability of the covalent bond in terms of overlapping atomic orbitals. Molecular Orbital Theory Valence bond theory accounts, at least qualitatively, for the stability of the covalent bond in terms of overlapping atomic orbitals. Using the concept of hybridization, valence

More information

EXAM INFORMATION. Radial Distribution Function: B is the normalization constant. d dx. p 2 Operator: Heisenberg Uncertainty Principle:

EXAM INFORMATION. Radial Distribution Function: B is the normalization constant. d dx. p 2 Operator: Heisenberg Uncertainty Principle: EXAM INFORMATION Radial Distribution Function: P() r RDF() r Br R() r B is the normalization constant., p Operator: p ^ d dx Heisenberg Uncertainty Principle: n ax n! Integrals: xe dx n1 a x p Particle

More information

We have arrived to the question: how do molecular bonds determine the band gap? We have discussed that the silicon atom has four outer electrons.

We have arrived to the question: how do molecular bonds determine the band gap? We have discussed that the silicon atom has four outer electrons. ET3034Tux - 2.2.2 - Band Gap 2 - Electrons in Molecular Bonds We have arrived to the question: how do molecular bonds determine the band gap? We have discussed that the silicon atom has four outer electrons.

More information

Functionalized Carbon Nanotubes a key to nanotechnology?

Functionalized Carbon Nanotubes a key to nanotechnology? 1 27th Max Born Symposium Multiscale Modeling of Real Materials Wroclaw, Sep 19, 2010 Functionalized Carbon Nanotubes a key to nanotechnology? Karolina Milowska, Magda Birowska & Jacek A. Majewski Faculty

More information

Lecture. Ref. Ihn Ch. 3, Yu&Cardona Ch. 2

Lecture. Ref. Ihn Ch. 3, Yu&Cardona Ch. 2 Lecture Review of quantum mechanics, statistical physics, and solid state Band structure of materials Semiconductor band structure Semiconductor nanostructures Ref. Ihn Ch. 3, Yu&Cardona Ch. 2 Reminder

More information

Ab-initio Molecular Electronics: are we reaching the limit? Yamila.García Applied Physics Department, Alicante University, Spain

Ab-initio Molecular Electronics: are we reaching the limit? Yamila.García Applied Physics Department, Alicante University, Spain Ab-initio Molecular Electronics: are we reaching the limit? Yamila.García Applied Physics Department, Alicante University, Spain To be, or not to be: that is the question: Whether tis nobler in mind to

More information

PAPER No. 7: Inorganic chemistry II MODULE No. 5: Molecular Orbital Theory

PAPER No. 7: Inorganic chemistry II MODULE No. 5: Molecular Orbital Theory Subject Chemistry Paper No and Title Module No and Title Module Tag 7, Inorganic chemistry II 5, Molecular Orbital Theory CHE_P7_M5 TABLE OF CONTENTS 1. Learning Outcomes 2. Introduction to Ligand Field

More information

LoProp: local property calculations with quantum chemical methods

LoProp: local property calculations with quantum chemical methods LoProp: local property calculations with quantum chemical methods Laura Gagliardi Dipartimento di Chimica Fisica F. Accascina Università degli Studi di Palermo Italy Gunnar Karlström, Jesper Krogh, and

More information

I 2 Vapor Absorption Experiment and Determination of Bond Dissociation Energy.

I 2 Vapor Absorption Experiment and Determination of Bond Dissociation Energy. I 2 Vapor Absorption Experiment and Determination of Bond Dissociation Energy. What determines the UV-Vis (i.e., electronic transitions) band appearance? Usually described by HOMO LUMO electron jump LUMO

More information

Computational Nanoscience

Computational Nanoscience Computational Nanoscience Applications for Molecules, Clusters, and Solids KALMÄN VARGA AND JOSEPH A. DRISCOLL Vanderbilt University, Tennessee Щ CAMBRIDGE HP UNIVERSITY PRESS Preface Part I One-dimensional

More information

Introduction to Hartree-Fock Molecular Orbital Theory

Introduction to Hartree-Fock Molecular Orbital Theory Introduction to Hartree-Fock Molecular Orbital Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Origins of Mathematical Modeling in Chemistry Plato (ca. 428-347

More information

Ab initio tight-binding LMTO method for nonequilibrium electron transport in nanosystems

Ab initio tight-binding LMTO method for nonequilibrium electron transport in nanosystems PHYSICA REVIEW B 71, 195422 2005 Ab initio tight-binding MTO method for nonequilibrium electron transport in nanosystems Sergey V. Faleev, 1, * François éonard, 1 Derek A. Stewart, 1, and Mark van Schilfgaarde

More information

List of Figures Page Figure No. Figure Caption No. Figure 1.1.

List of Figures Page Figure No. Figure Caption No. Figure 1.1. List of Figures Figure No. Figure Caption Page No. Figure 1.1. Cation- interactions and their modulations. 4 Figure 1.2. Three conformations of benzene dimer, S is not a minimum on the potential energy

More information

GaAs -- MLWF. Special thanks to Elias Assmann (TU Graz) for the generous help in preparation of this tutorial

GaAs -- MLWF. Special thanks to Elias Assmann (TU Graz) for the generous help in preparation of this tutorial GaAs -- MLWF + + Special thanks to Elias Assmann (TU Graz) for the generous help in preparation of this tutorial YouTube video: https://youtu.be/r4c1yhdh3ge 1. Wien2k SCF Create a tutorial directory, e.g.

More information

CHEM6085: Density Functional Theory Lecture 10

CHEM6085: Density Functional Theory Lecture 10 CHEM6085: Density Functional Theory Lecture 10 1) Spin-polarised calculations 2) Geometry optimisation C.-K. Skylaris 1 Unpaired electrons So far we have developed Kohn-Sham DFT for the case of paired

More information

Supporting information for: Effect of Structural Fluctuations on. Elastic Lifetimes of Adsorbate States. Isonicotinic Acid on Rutile(110)

Supporting information for: Effect of Structural Fluctuations on. Elastic Lifetimes of Adsorbate States. Isonicotinic Acid on Rutile(110) Supporting information for: Effect of Structural Fluctuations on Elastic Lifetimes of Adsorbate States Isonicotinic Acid on Rutile(110) Moritz Müller,,,, Daniel Sánchez-Portal,,, He Lin,,# Gian Paolo Brivio,

More information

Theoretical Concepts of Spin-Orbit Splitting

Theoretical Concepts of Spin-Orbit Splitting Chapter 9 Theoretical Concepts of Spin-Orbit Splitting 9.1 Free-electron model In order to understand the basic origin of spin-orbit coupling at the surface of a crystal, it is a natural starting point

More information

Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory.

Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn was awarded with the Nobel Prize in Chemistry in 1998 for his development of the density functional theory. Walter Kohn receiving his Nobel Prize from His Majesty the King at the Stockholm

More information

Chapter 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 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 information

The Plane-Wave Pseudopotential Method

The 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 information

Simulating charge transfer in (organic) solar cells

Simulating charge transfer in (organic) solar cells Simulating charge transfer in (organic) solar cells July 24, 2013 Organic solar cells 1 c Fraunhofer ISE Organic solar cells Step 1: Microscopic structure Organic solar cells 1 c Fraunhofer ISE Organic

More information

MOLECULAR STRUCTURE. Molecular Structure - B. Molecular Structure - B. Molecular Structure - B. Molecular Structure - B. Molecular Structure - B

MOLECULAR STRUCTURE. Molecular Structure - B. Molecular Structure - B. Molecular Structure - B. Molecular Structure - B. Molecular Structure - B MOLECULAR STRUCTURE Molecular Orbital all orbitals of the appropriate symmetry contribute to a molecular orbital. Bundet Boekfa Chem Div, Faculty Lib Arts & Sci Kasetsart University Kamphaeng Saen Campus

More information

Introduction to Density Functional Theory (DFT)

Introduction to Density Functional Theory (DFT) Introduction to Density Functional Theory (DFT) Brad Malone, Sadas Shankar The Problem: What's the big deal? Materials we are often interested in contain a macroscopically large number of particles (~1023

More information

Introduction to Hartree-Fock calculations in Spartan

Introduction to Hartree-Fock calculations in Spartan EE5 in 2008 Hannes Jónsson Introduction to Hartree-Fock calculations in Spartan In this exercise, you will get to use state of the art software for carrying out calculations of wavefunctions for molecues,

More information

Systematic convergence for realistic projects Fast versus accurate

Systematic 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 information

Problem with Kohn-Sham equations

Problem with Kohn-Sham equations Problem with Kohn-Sham equations (So much time consuming) H s Ψ = E el Ψ ( T + V [ n] + V [ n] + V [ n]) ϕ = Eϕ i = 1, 2,.., N s e e ext XC i i N nr ( ) = ϕi i= 1 2 The one-particle Kohn-Sham equations

More information

Electronic structure calculations: fundamentals George C. Schatz Northwestern University

Electronic structure calculations: fundamentals George C. Schatz Northwestern University Electronic structure calculations: fundamentals George C. Schatz Northwestern University Electronic Structure (often called Quantum Chemistry) calculations use quantum mechanics to determine the wavefunctions

More information

Luigi Paolasini

Luigi Paolasini Luigi Paolasini paolasini@esrf.fr LECTURE 4: MAGNETIC INTERACTIONS - Dipole vs exchange magnetic interactions. - Direct and indirect exchange interactions. - Anisotropic exchange interactions. - Interplay

More information

NH 3 H 2 O N 2. Why do they make chemical bonds? Molecular Orbitals

NH 3 H 2 O N 2. Why do they make chemical bonds? Molecular Orbitals N 2 NH 3 H 2 O Why do they make chemical bonds? 5 Molecular Orbitals Why do they make chemical bonds? Stabilization Bond energy Types of Chemical Bonds Metallic Bond Ionic Bond Covalent Bond Covalent Bond

More information

arxiv: v1 [cond-mat.mes-hall] 28 Aug 2009

arxiv: v1 [cond-mat.mes-hall] 28 Aug 2009 Efficient implementation of the nonequilibrium Green function method for electronic transport calculations Taisuke Ozaki Research Center for Integrated Science (RCIS), Japan Advanced Institute of Science

More information

Formation and characterization of a. molecule-metal-molecule bridge in real space SUPPORTING INFORMATION

Formation and characterization of a. molecule-metal-molecule bridge in real space SUPPORTING INFORMATION Formation and characterization of a molecule-metal-molecule bridge in real space SUPPORTING INFORMATION Florian Albrecht,, Mathias Neu, Christina Quest, Ingmar Swart,, and Jascha Repp Institute of Experimental

More information

Lecture 4: Band theory

Lecture 4: Band theory Lecture 4: Band theory Very short introduction to modern computational solid state chemistry Band theory of solids Molecules vs. solids Band structures Analysis of chemical bonding in Reciprocal space

More information

Electron States of Diatomic Molecules

Electron States of Diatomic Molecules IISER Pune March 2018 Hamiltonian for a Diatomic Molecule The hamiltonian for a diatomic molecule can be considered to be made up of three terms Ĥ = ˆT N + ˆT el + ˆV where ˆT N is the kinetic energy operator

More information