Density Functional Theory I

Size: px
Start display at page:

Download "Density Functional Theory I"

Transcription

1 Densty Functonal Theoy I cholas M. Hason Depatment of Chemsty Impeal College Lonon & Computatonal Mateals Scence Daesbuy Laboatoy ncholas.hason@c.ac.uk

2 Densty Functonal Theoy I The Many Electon Schönge Equaton Hatee-Fock Theoy Solvng the Schönge Equaton Avong Solvng the Schönge Equaton Densty Matces The basc eas of DFT Kohn-Sham an the non-nteactng system The local ensty appoxmaton LDA Conclusons Densty Functonal Theoy II

3 The Schönge Equaton fo Many Electons Tme nepenent non-elatvstc Bon-Oppenheme A lnea equaton n 3 vaables < j j ext V H E H Ψ Ψ Thee tems knetc enegy extenal potental an electon-electon.

4 The Extenal Potental The nteacton wth the atomc nucle s: V ext at Zα R α α The extenal potental an the numbe of electons completely etemne the Hamltonan.

5 The Vaatonal Pncples o any legal wavefuncton antsymmetc nomalse e enegy s E * [ Ψ] Ψ HΨ Ψ H Ψ - the enegy s a functonal of Ψ E[ Ψ] E 0 - Seach all Ψ to mnmse E > the goun state

6 Hatee Fock Theoy An ansatz fo the stuctue of the wavefuncton [ ] HF... et! 3 Ψ j j j j j j ext HF V E * * * * *

7 The Hatee Fock Equatons non-nteactng electons une the nfluence of a mean fel potental consstng of the classcal Coulomb potental an a nonlocal exchange potental. X ext v V ε j j j X v *

8 Beyon Hatee-Fock Many methos/appoxmatons applcable to small systems. Expensve & scalng wth poblem sze s feocous Eg: MP MP3 MP4 ~ CISD ~ 6 CCSD ~ 6 CCSDT ~ 7

9 Avong Solvng the Schonge Equaton Is t necessay to solve the Schönge equaton an etemne the 3 mensonal wavefuncton n oe to compute the goun state enegy?.. o!

10 The Pa Densty The secon oe ensty matx s efne as P ; 4 3 * Ψ Ψ The agonal elements ae the two patcle ensty matx o pa ensty functon ; P P Ths object n 6 mensons s all we nee to compute the exact total enegy!

11 The Enegy an the Densty Matces The fst oe ensty matx s H only contans one-electon an two-electon opeatos - the enegy can be wtten exactly n tems of P an P P P ; ; R P P Z E at α α α

12 The enegy s a functonal of P : E[P ] Pehaps solve usng by mnmsng E[P ]? Vey ffcult poblem A legal P must be constuctble fom an antsymmetc Ψ. Applyng ths constant n pactce s non tval. Lfe woul be much ease f thee was a way of ong t..!

13 Do we eally nee to know P? o To fn the exact total enegy knowlege of the chage ensty s enough!

14 DFT the theoems Theoem. The extenal potental s unquely etemne by the ensty - - so the total enegy s a unque functonal of the ensty - E[]! Theoem. The ensty whch mnmses the enegy s the goun state ensty an the mnmum enegy s the goun state enegy. [ ] E0 Mn E Hohenbug-Kohn 964 Mel Levy E B Wlson

15 Theoem Wlson s poof The chage ensty has a cusp at the nucleus of any atom such that 0 0 α α α α Z The chage ensty unquely etemnes the Hamltonan thus the wavefuncton an all mateal popetes!!!

16 Densty Functonal Theoy The funamental statement of DFT s δ [ ] E[ ] µ 0 thee s a unvesal functonal E[] whch coul be nsete nto the above equaton an mnmse to obtan the exact goun state ensty an enegy.

17 What s the functonal? Thee ae thee tems n the Hamltonan V ext [] s tval E [ ] T[ ] V ext [ ] V [ ] ee V [ ] V ext ext T an V ee ae vey ffcult to appoxmate!

18 E xc [] - Popetes Does not epen on V ext the specfc system: t s a unvesal functonal. The exact epenence on s unknown E xc can be exactly etemne fo any specfc ensty but the effot s geate than fo usual many-boy calculatons E xc must be appoxmate n applcatons

19 The Homogeneous Electon Gas Fo the non-nteactng gas the knetc an exchange enegy pe patcle can be compute the sngle patcle wavefunctons ae smply plane waves. Pehaps ntegate these enegy enstes fo an nhomogeneous system? 3 T [ ] o chemcal bonng! 3 E x 0.74 E x s OK but what about T? 4 Thomas 97 Fem 98 Dac 930

20 A local functon of the ensty Thomas-Fem-Dac suggests: [ ] Exc ε xc ε ε ε xc x c ε x 3 C J. Slate 95 The functonal s a local functon of the ensty What to o about ε c?

21 E[] The Kohn Sham Appoach Wte the ensty n tems of a set of non-nteactng obtals The non nteactng knetc enegy an the classcal Coulomb nteacton s T ] [ ] [ V H Allow us to ecast the enegy functonal as: ] [ ] [ ] [ ] [ ] [ xc H ext s E V V T E ] [ ] [ ] [ ] [ ] [ H ee s xc V V T T E Whee we have ntouce

22 Vaaton Theoem > Kohn Sham Equatons Vay the enegy wth espect to the obtals an. o appoxmatons So If we knew E xc [] we coul solve fo the exact goun state enegy an ensty! Cost 3.. In pncple. ] [ xc xc E V ext V xc V ε

23 The on-nteactng system Thee exsts an effectve mean fel potental whch when apple to a system of non-nteactng femons wll geneate the exact goun state enegy an chage ensty!!! Ψ j... V xc E[]

24 Hatee Fock s a ensty functonal theoy Hatee-Fock theoy s a ensty functonal theoy! - but wth a non-local potental. j j j xc V V

25 Quantum Monte Calo Smulatons Fo the Homogeneous electon gas the exact epenence of ε xc can be compute. ε xc Cepeley an Ale 980

26 The Local Densty Appoxmaton - LDA LDA Exc [ ] ε xc ε xc Pctue coutesy of Aneas Savn

27 Conclusons I Fo the goun state enegy an ensty thee s an exact mappng between the many boy system an a fcttous nonnteactng system. - DFT-people stuy the fcttous system! The fcttous system s subject to an unknown potental eve fom the exchange-coelaton functonal The enegy functonal may be appoxmate as a local functon of the ensty!

28 Densty Functonal Theoy II Why oes the LDA wok? The exchange coelaton hole Compason wth exact exchange an coelaton enegy enstes Genealse gaent appoxmatons GGA s Sem-local nteactons: Meta-GGA s Hyb-exchange functonals Pefomance n molecules an sols

Slide 1. Quantum Mechanics: the Practice

Slide 1. Quantum Mechanics: the Practice Slde Quantum Mecancs: te Pactce Slde Remnde: Electons As Waves Wavelengt momentum = Planck? λ p = = 6.6 x 0-34 J s Te wave s an exctaton a vbaton: We need to know te ampltude of te exctaton at evey pont

More information

Electron density: Properties of electron density (non-negative): => exchange-correlation functionals should respect these conditions.

Electron density: Properties of electron density (non-negative): => exchange-correlation functionals should respect these conditions. lecton densty: ρ ( =... Ψ(,,..., ds d... d Pobablty of fndng one electon of abtay spn wthn a volume element d (othe electons may be anywhee. s Popetes of electon densty (non-negatve:.. 3. ρ ( d = ρ( =

More information

THE MANY-BODY PROBLEM

THE MANY-BODY PROBLEM 3.30: Lectue 5 Feb 5 005 THE MANY-BODY PROBLEM Feb 5 005 3.30 Atomstc Modelng of Mateals -- Geband Cede and Ncola Maza When s a patcle lke a wave? Wavelength momentum = Planck λ p = h h = 6.6 x 0-34 J

More information

E For K > 0. s s s s Fall Physical Chemistry (II) by M. Lim. singlet. triplet

E For K > 0. s s s s Fall Physical Chemistry (II) by M. Lim. singlet. triplet Eneges of He electonc ψ E Fo K > 0 ψ = snglet ( )( ) s s+ ss αβ E βα snglet = ε + ε + J s + Ks Etplet = ε + ε + J s Ks αα ψ tplet = ( s s ss ) ββ ( αβ + βα ) s s s s s s s s ψ G = ss( αβ βα ) E = ε + ε

More information

19 The Born-Oppenheimer Approximation

19 The Born-Oppenheimer Approximation 9 The Bon-Oppenheme Appoxmaton The full nonelatvstc Hamltonan fo a molecule s gven by (n a.u.) Ĥ = A M A A A, Z A + A + >j j (883) Lets ewte the Hamltonan to emphasze the goal as Ĥ = + A A A, >j j M A

More information

Supersymmetry in Disorder and Chaos (Random matrices, physics of compound nuclei, mathematics of random processes)

Supersymmetry in Disorder and Chaos (Random matrices, physics of compound nuclei, mathematics of random processes) Supesymmety n Dsoe an Chaos Ranom matces physcs of compoun nucle mathematcs of anom pocesses Lteatue: K.B. Efetov Supesymmety n Dsoe an Chaos Cambge Unvesty Pess 997999 Supesymmety an Tace Fomulae I.V.

More information

Introduction to Density Functional Theory. Jeremie Zaffran 2 nd year-msc. (Nanochemistry)

Introduction to Density Functional Theory. Jeremie Zaffran 2 nd year-msc. (Nanochemistry) Introducton to Densty Functonal Theory Jereme Zaffran nd year-msc. (anochemstry) A- Hartree appromatons Born- Oppenhemer appromaton H H H e The goal of computatonal chemstry H e??? Let s remnd H e T e

More information

DENSITY-FUNCTIONAL THEORY, AND DENSITY-FUNCTIONAL PRACTICE

DENSITY-FUNCTIONAL THEORY, AND DENSITY-FUNCTIONAL PRACTICE 3.320: Lectue 7 Feb 24 2005 DENSITY-FUNCTIONAL THEORY, AND DENSITY-FUNCTIONAL PRACTICE Feb 24 2005 3.320 Atomistic Modeling of Mateials -- Geband Cede and Nicola Mazai Hatee-Fock Equations ψ,,..., = 1

More information

Complex atoms and the Periodic System of the elements

Complex atoms and the Periodic System of the elements Complex atoms and the Peodc System of the elements Non-cental foces due to electon epulson Cental feld appoxmaton electonc obtals lft degeneacy of l E n l = R( hc) ( n δ ) l Aufbau pncple Lectue Notes

More information

Energy in Closed Systems

Energy in Closed Systems Enegy n Closed Systems Anamta Palt palt.anamta@gmal.com Abstact The wtng ndcates a beakdown of the classcal laws. We consde consevaton of enegy wth a many body system n elaton to the nvese squae law and

More information

8 Baire Category Theorem and Uniform Boundedness

8 Baire Category Theorem and Uniform Boundedness 8 Bae Categoy Theoem and Unfom Boundedness Pncple 8.1 Bae s Categoy Theoem Valdty of many esults n analyss depends on the completeness popety. Ths popety addesses the nadequacy of the system of atonal

More information

Part V: Velocity and Acceleration Analysis of Mechanisms

Part V: Velocity and Acceleration Analysis of Mechanisms Pat V: Velocty an Acceleaton Analyss of Mechansms Ths secton wll evew the most common an cuently pactce methos fo completng the knematcs analyss of mechansms; escbng moton though velocty an acceleaton.

More information

PHYS 705: Classical Mechanics. Derivation of Lagrange Equations from D Alembert s Principle

PHYS 705: Classical Mechanics. Derivation of Lagrange Equations from D Alembert s Principle 1 PHYS 705: Classcal Mechancs Devaton of Lagange Equatons fom D Alembet s Pncple 2 D Alembet s Pncple Followng a smla agument fo the vtual dsplacement to be consstent wth constants,.e, (no vtual wok fo

More information

Teoría del Funcional de la Densidad (Density Functional Theory)

Teoría del Funcional de la Densidad (Density Functional Theory) Teoía del Funcional de la Densidad (Density Functional Theoy) Motivation: limitations of the standad appoach based on the wave function. The electonic density n() as the key vaiable: Functionals & Thomas-Femi

More information

Density functional theory

Density functional theory Densty functonal theoy Tadtonal ab nto: fndng the N-electon wavefunton Ψ(1,,,N) dependng on 4N cood. DFT: fndng the total electon spn-denstes dependng on 8 coodnates Hohenbeg & Kohn: E0[ ρ0] = ρ0( ) VNed

More information

Large scale magnetic field generation by accelerated particles in galactic medium

Large scale magnetic field generation by accelerated particles in galactic medium Lage scale magnetc feld geneaton by acceleated patcles n galactc medum I.N.Toptygn Sant Petesbug State Polytechncal Unvesty, depatment of Theoetcal Physcs, Sant Petesbug, Russa 2.Reason explonatons The

More information

Nuclear Chart. Takashi NAKATSUKASA Theoretical Nuclear Physics Laboratory RIKEN Nishina Center. Real-space, real-time approaches ) Few-body model

Nuclear Chart. Takashi NAKATSUKASA Theoretical Nuclear Physics Laboratory RIKEN Nishina Center. Real-space, real-time approaches ) Few-body model Takash NAKATSUKASA Theoetcal Nuclea Physcs Laboatoy RIKEN Nshna Cente 009.3.5-6 Mn-WS: Real-space, eal-tme appoaches DFT, TDDFT ((Q)RPA ) Few-body model (CDCC ) Nuclea Chat 70 Los Alamos Natonal Laboatoy's

More information

Set of square-integrable function 2 L : function space F

Set of square-integrable function 2 L : function space F Set of squae-ntegable functon L : functon space F Motvaton: In ou pevous dscussons we have seen that fo fee patcles wave equatons (Helmholt o Schödnge) can be expessed n tems of egenvalue equatons. H E,

More information

UNIT10 PLANE OF REGRESSION

UNIT10 PLANE OF REGRESSION UIT0 PLAE OF REGRESSIO Plane of Regesson Stuctue 0. Intoducton Ojectves 0. Yule s otaton 0. Plane of Regesson fo thee Vaales 0.4 Popetes of Resduals 0.5 Vaance of the Resduals 0.6 Summay 0.7 Solutons /

More information

Physics 2A Chapter 11 - Universal Gravitation Fall 2017

Physics 2A Chapter 11 - Universal Gravitation Fall 2017 Physcs A Chapte - Unvesal Gavtaton Fall 07 hese notes ae ve pages. A quck summay: he text boxes n the notes contan the esults that wll compse the toolbox o Chapte. hee ae thee sectons: the law o gavtaton,

More information

2 dependence in the electrostatic force means that it is also

2 dependence in the electrostatic force means that it is also lectc Potental negy an lectc Potental A scala el, nvolvng magntues only, s oten ease to wo wth when compae to a vecto el. Fo electc els not havng to begn wth vecto ssues woul be nce. To aange ths a scala

More information

THE TIME-DEPENDENT CLOSE-COUPLING METHOD FOR ELECTRON-IMPACT DIFFERENTIAL IONIZATION CROSS SECTIONS FOR ATOMS AND MOLECULES

THE TIME-DEPENDENT CLOSE-COUPLING METHOD FOR ELECTRON-IMPACT DIFFERENTIAL IONIZATION CROSS SECTIONS FOR ATOMS AND MOLECULES Intenatonal The Tme-Dependent cence Pess Close-Couplng IN: 9-59 Method fo Electon-Impact Dffeental Ionzaton Coss ectons fo Atoms... REVIEW ARTICE THE TIME-DEPENDENT COE-COUPING METHOD FOR EECTRON-IMPACT

More information

4.4 Continuum Thermomechanics

4.4 Continuum Thermomechanics 4.4 Contnuum Themomechancs The classcal themodynamcs s now extended to the themomechancs of a contnuum. The state aables ae allowed to ay thoughout a mateal and pocesses ae allowed to be eesble and moe

More information

Machine Learning. Spectral Clustering. Lecture 23, April 14, Reading: Eric Xing 1

Machine Learning. Spectral Clustering. Lecture 23, April 14, Reading: Eric Xing 1 Machne Leanng -7/5 7/5-78, 78, Spng 8 Spectal Clusteng Ec Xng Lectue 3, pl 4, 8 Readng: Ec Xng Data Clusteng wo dffeent ctea Compactness, e.g., k-means, mxtue models Connectvty, e.g., spectal clusteng

More information

Chapter 23: Electric Potential

Chapter 23: Electric Potential Chapte 23: Electc Potental Electc Potental Enegy It tuns out (won t show ths) that the tostatc foce, qq 1 2 F ˆ = k, s consevatve. 2 Recall, fo any consevatve foce, t s always possble to wte the wok done

More information

9/12/2013. Microelectronics Circuit Analysis and Design. Modes of Operation. Cross Section of Integrated Circuit npn Transistor

9/12/2013. Microelectronics Circuit Analysis and Design. Modes of Operation. Cross Section of Integrated Circuit npn Transistor Mcoelectoncs Ccut Analyss and Desgn Donald A. Neamen Chapte 5 The pola Juncton Tanssto In ths chapte, we wll: Dscuss the physcal stuctue and opeaton of the bpola juncton tanssto. Undestand the dc analyss

More information

A NOTE ON ELASTICITY ESTIMATION OF CENSORED DEMAND

A NOTE ON ELASTICITY ESTIMATION OF CENSORED DEMAND Octobe 003 B 003-09 A NOT ON ASTICITY STIATION OF CNSOD DAND Dansheng Dong an Hay. Kase Conell nvesty Depatment of Apple conomcs an anagement College of Agcultue an fe Scences Conell nvesty Ithaca New

More information

3.320: Lecture 6 (Feb ) Photos of Hartree, Fock, Hohenberg, Kohn, and Sham removed for copyright reasons.

3.320: Lecture 6 (Feb ) Photos of Hartree, Fock, Hohenberg, Kohn, and Sham removed for copyright reasons. 3.30: Lectue 6 Feb 7 005 HARTREE-FOCK AND DFT Photos of Hatee, Fock, Hohenbeg, Kohn, and Sham emoved fo copyight easons. Feb 7 005 3.30 Atomistic Modeling of Mateials -- Geband Cede and Nicola Mazai Matix

More information

Khintchine-Type Inequalities and Their Applications in Optimization

Khintchine-Type Inequalities and Their Applications in Optimization Khntchne-Type Inequaltes and The Applcatons n Optmzaton Anthony Man-Cho So Depatment of Systems Engneeng & Engneeng Management The Chnese Unvesty of Hong Kong ISDS-Kolloquum Unvestaet Wen 29 June 2009

More information

A Lattice Energy Calculation for LiH

A Lattice Energy Calculation for LiH A Lattice Enegy Calculation fo LiH Fank Riou Lithium hyie is a white cystalline soli with the face-centee cubic cystal stuctue (see lattice shown below). The moel fo LiH(s) popose in this stuy constists

More information

DFT functionals in WIEN2k

DFT functionals in WIEN2k d E xc DFT functionals in WIEN2k ] [ 2 1 ] [ xc ext o E d d d V T E 1-electon equations Kohn Sham } 2 1 { 2 V V V i i i xc C ext E F i i 2 vay de/d=0 -Z/ E xc = E x + E c : exchange-coelation enegy V xc

More information

Physics 11b Lecture #2. Electric Field Electric Flux Gauss s Law

Physics 11b Lecture #2. Electric Field Electric Flux Gauss s Law Physcs 11b Lectue # Electc Feld Electc Flux Gauss s Law What We Dd Last Tme Electc chage = How object esponds to electc foce Comes n postve and negatve flavos Conseved Electc foce Coulomb s Law F Same

More information

Multistage Median Ranked Set Sampling for Estimating the Population Median

Multistage Median Ranked Set Sampling for Estimating the Population Median Jounal of Mathematcs and Statstcs 3 (: 58-64 007 ISSN 549-3644 007 Scence Publcatons Multstage Medan Ranked Set Samplng fo Estmatng the Populaton Medan Abdul Azz Jeman Ame Al-Oma and Kamaulzaman Ibahm

More information

Integral Vector Operations and Related Theorems Applications in Mechanics and E&M

Integral Vector Operations and Related Theorems Applications in Mechanics and E&M Dola Bagayoko (0) Integal Vecto Opeatons and elated Theoems Applcatons n Mechancs and E&M Ι Basc Defnton Please efe to you calculus evewed below. Ι, ΙΙ, andιιι notes and textbooks fo detals on the concepts

More information

Stellar Astrophysics. dt dr. GM r. The current model for treating convection in stellar interiors is called mixing length theory:

Stellar Astrophysics. dt dr. GM r. The current model for treating convection in stellar interiors is called mixing length theory: Stella Astophyscs Ovevew of last lectue: We connected the mean molecula weght to the mass factons X, Y and Z: 1 1 1 = X + Y + μ 1 4 n 1 (1 + 1) = X μ 1 1 A n Z (1 + ) + Y + 4 1+ z A Z We ntoduced the pessue

More information

Lecture 14: Forces and Stresses

Lecture 14: Forces and Stresses The Nuts and Bolts of Frst-Prncples Smulaton Lecture 14: Forces and Stresses Durham, 6th-13th December 2001 CASTEP Developers Group wth support from the ESF ψ k Network Overvew of Lecture Why bother? Theoretcal

More information

Photodisintegration of light nuclei

Photodisintegration of light nuclei Photodsntegaton of lght nucle N Banea The Hebew Unvesty Jeusalem Guseppna Olandn Tento (Italy) Wnfed Ledemann Tento (Italy) Sona Bacca Damstadt (Gemany) Doon Gazt Jeusalem (Isael) Yael Ronen Jeusalem (Isael)

More information

4. Compare the electric force holding the electron in orbit ( r = 0.53

4. Compare the electric force holding the electron in orbit ( r = 0.53 Electostatics WS Electic Foce an Fiel. Calculate the magnitue of the foce between two 3.60-µ C point chages 9.3 cm apat.. How many electons make up a chage of 30.0 µ C? 3. Two chage ust paticles exet a

More information

V. Principles of Irreversible Thermodynamics. s = S - S 0 (7.3) s = = - g i, k. "Flux": = da i. "Force": = -Â g a ik k = X i. Â J i X i (7.

V. Principles of Irreversible Thermodynamics. s = S - S 0 (7.3) s = = - g i, k. Flux: = da i. Force: = -Â g a ik k = X i. Â J i X i (7. Themodynamcs and Knetcs of Solds 71 V. Pncples of Ievesble Themodynamcs 5. Onsage s Teatment s = S - S 0 = s( a 1, a 2,...) a n = A g - A n (7.6) Equlbum themodynamcs detemnes the paametes of an equlbum

More information

Test 1 phy What mass of a material with density ρ is required to make a hollow spherical shell having inner radius r i and outer radius r o?

Test 1 phy What mass of a material with density ρ is required to make a hollow spherical shell having inner radius r i and outer radius r o? Test 1 phy 0 1. a) What s the pupose of measuement? b) Wte all fou condtons, whch must be satsfed by a scala poduct. (Use dffeent symbols to dstngush opeatons on ectos fom opeatons on numbes.) c) What

More information

PHYS Week 5. Reading Journals today from tables. WebAssign due Wed nite

PHYS Week 5. Reading Journals today from tables. WebAssign due Wed nite PHYS 015 -- Week 5 Readng Jounals today fom tables WebAssgn due Wed nte Fo exclusve use n PHYS 015. Not fo e-dstbuton. Some mateals Copyght Unvesty of Coloado, Cengage,, Peason J. Maps. Fundamental Tools

More information

Chapter 8. Linear Momentum, Impulse, and Collisions

Chapter 8. Linear Momentum, Impulse, and Collisions Chapte 8 Lnea oentu, Ipulse, and Collsons 8. Lnea oentu and Ipulse The lnea oentu p of a patcle of ass ovng wth velocty v s defned as: p " v ote that p s a vecto that ponts n the sae decton as the velocty

More information

Tian Zheng Department of Statistics Columbia University

Tian Zheng Department of Statistics Columbia University Haplotype Tansmsson Assocaton (HTA) An "Impotance" Measue fo Selectng Genetc Makes Tan Zheng Depatment of Statstcs Columba Unvesty Ths s a jont wok wth Pofesso Shaw-Hwa Lo n the Depatment of Statstcs at

More information

COMPLEMENTARY ENERGY METHOD FOR CURVED COMPOSITE BEAMS

COMPLEMENTARY ENERGY METHOD FOR CURVED COMPOSITE BEAMS ultscence - XXX. mcocd Intenatonal ultdscplnay Scentfc Confeence Unvesty of skolc Hungay - pl 06 ISBN 978-963-358-3- COPLEENTRY ENERGY ETHOD FOR CURVED COPOSITE BES Ákos József Lengyel István Ecsed ssstant

More information

A Scheme for Calculating Atomic Structures beyond the Spherical Approximation

A Scheme for Calculating Atomic Structures beyond the Spherical Approximation Jounal of Moden Physcs,,, 4-4 do:.46/jmp..55 Publshed Onlne May (http://www.crp.og/jounal/jmp) A cheme fo Calculatng Atomc tuctues beyond the phecal Appoxmaton Abstact Mtyasu Myasta, Katsuhko Hguch, Masahko

More information

Chapter Fifiteen. Surfaces Revisited

Chapter Fifiteen. Surfaces Revisited Chapte Ffteen ufaces Revsted 15.1 Vecto Descpton of ufaces We look now at the vey specal case of functons : D R 3, whee D R s a nce subset of the plane. We suppose s a nce functon. As the pont ( s, t)

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 5 Lecture 3 Contnuous Systems an Fels (Chapter 3) Where Are We Now? We ve fnshe all the essentals Fnal wll cover Lectures through Last two lectures: Classcal Fel Theory Start wth wave equatons

More information

Arbitrary Choice of Basic Variables in Density Functional Theory. I. Formalism

Arbitrary Choice of Basic Variables in Density Functional Theory. I. Formalism Abtay Choce of Basc Vaables n Densty Functonal Theoy. I. Foals Masaho Hguch * Insttute of Sol State an Mateals Reseach P.O. Box 76 D-7 Desen Geany Katsuho Hguch Deatent of lectcal ngneeng Hosha Unvesty

More information

Plasmonics and non-local interactions from TDDFT: graphene and metal surfaces

Plasmonics and non-local interactions from TDDFT: graphene and metal surfaces Plasmonics and non-local inteactions fom TDDFT: gaphene and metal sufaces Thomas Olsen Cente fo Atomic-scale Mateials Design CAMD Depatment of Physics Technical Univesity of Denmak Outline Linea esponse

More information

EE 5337 Computational Electromagnetics (CEM)

EE 5337 Computational Electromagnetics (CEM) 7//28 Instucto D. Raymond Rumpf (95) 747 6958 cumpf@utep.edu EE 5337 Computatonal Electomagnetcs (CEM) Lectue #6 TMM Extas Lectue 6These notes may contan copyghted mateal obtaned unde fa use ules. Dstbuton

More information

A. Thicknesses and Densities

A. Thicknesses and Densities 10 Lab0 The Eath s Shells A. Thcknesses and Denstes Any theoy of the nteo of the Eath must be consstent wth the fact that ts aggegate densty s 5.5 g/cm (ecall we calculated ths densty last tme). In othe

More information

Molecular Dynamics and Density Functional Theory

Molecular Dynamics and Density Functional Theory Molecular Dynamcs and Densty Functonal Theory What do we need? An account n pemfc cluster: Host name: pemfc.chem.sfu.ca I wll take care of that. Ths can be usually a common account for all of you but please

More information

Select/Special Topics in Atomic Physics

Select/Special Topics in Atomic Physics Select/Specal Topcs n Atomc Physcs P. C. Deshmukh Depatment of Physcs Indan Insttute of Technology Madas Chenna 600036 pcd@physcs.tm.ac.n Unt 4() Lectue 9 Many-Electon Atoms: Hatee-Fock Self-Consstent

More information

Thermodynamics of solids 4. Statistical thermodynamics and the 3 rd law. Kwangheon Park Kyung Hee University Department of Nuclear Engineering

Thermodynamics of solids 4. Statistical thermodynamics and the 3 rd law. Kwangheon Park Kyung Hee University Department of Nuclear Engineering Themodynamcs of solds 4. Statstcal themodynamcs and the 3 d law Kwangheon Pak Kyung Hee Unvesty Depatment of Nuclea Engneeng 4.1. Intoducton to statstcal themodynamcs Classcal themodynamcs Statstcal themodynamcs

More information

Event Shape Update. T. Doyle S. Hanlon I. Skillicorn. A. Everett A. Savin. Event Shapes, A. Everett, U. Wisconsin ZEUS Meeting, October 15,

Event Shape Update. T. Doyle S. Hanlon I. Skillicorn. A. Everett A. Savin. Event Shapes, A. Everett, U. Wisconsin ZEUS Meeting, October 15, Event Shape Update A. Eveett A. Savn T. Doyle S. Hanlon I. Skllcon Event Shapes, A. Eveett, U. Wsconsn ZEUS Meetng, Octobe 15, 2003-1 Outlne Pogess of Event Shapes n DIS Smla to publshed pape: Powe Coecton

More information

Contact, information, consultations

Contact, information, consultations ontact, nfomaton, consultatons hemsty A Bldg; oom 07 phone: 058-347-769 cellula: 664 66 97 E-mal: wojtek_c@pg.gda.pl Offce hous: Fday, 9-0 a.m. A quote of the week (o camel of the week): hee s no expedence

More information

Many Electron Atoms. Electrons can be put into approximate orbitals and the properties of the many electron systems can be catalogued

Many Electron Atoms. Electrons can be put into approximate orbitals and the properties of the many electron systems can be catalogued Many Electon Atoms The many body poblem cannot be solved analytically. We content ouselves with developing appoximate methods that can yield quite accuate esults (but usually equie a compute). The electons

More information

A capacitor is simply two pieces of metal near each other, separated by an insulator or air. A capacitor is used to store charge and energy.

A capacitor is simply two pieces of metal near each other, separated by an insulator or air. A capacitor is used to store charge and energy. -1 apactors A capactor s smply two peces of metal near each other, separate by an nsulator or ar. A capactor s use to store charge an energy. A parallel-plate capactor conssts of two parallel plates separate

More information

CHAPTER 4 EVALUATION OF FORCE-CONSTANT MATRIX

CHAPTER 4 EVALUATION OF FORCE-CONSTANT MATRIX CHAPTER 4 EVALUATION OF FORCE-CONSTANT MATRIX 4.- AIM OF THE WORK As antcpated n the ntodcton the am of the pesent ok s to obtan the nmecal vale of the foce-constant matx fo tantalm. In a fst step expesson

More information

Solving the Dirac Equation: Using Fourier Transform

Solving the Dirac Equation: Using Fourier Transform McNa Schola Reeach Jounal Volume Atcle Solvng the ac quaton: Ung oue Tanfom Vncent P. Bell mby-rddle Aeonautcal Unvety, Vncent.Bell@my.eau.edu ollow th and addtonal wok at: http://common.eau.edu/na Recommended

More information

CSJM University Class: B.Sc.-II Sub:Physics Paper-II Title: Electromagnetics Unit-1: Electrostatics Lecture: 1 to 4

CSJM University Class: B.Sc.-II Sub:Physics Paper-II Title: Electromagnetics Unit-1: Electrostatics Lecture: 1 to 4 CSJM Unvesty Class: B.Sc.-II Sub:Physcs Pape-II Ttle: Electomagnetcs Unt-: Electostatcs Lectue: to 4 Electostatcs: It deals the study of behavo of statc o statonay Chages. Electc Chage: It s popety by

More information

Mechanics Physics 151

Mechanics Physics 151 Mechancs Physcs 151 Lectue 18 Hamltonan Equatons of Moton (Chapte 8) What s Ahead We ae statng Hamltonan fomalsm Hamltonan equaton Today and 11/6 Canoncal tansfomaton 1/3, 1/5, 1/10 Close lnk to non-elatvstc

More information

Analytical classical dynamics

Analytical classical dynamics Analytcal classcal ynamcs by Youun Hu Insttute of plasma physcs, Chnese Acaemy of Scences Emal: yhu@pp.cas.cn Abstract These notes were ntally wrtten when I rea tzpatrck s book[] an were later revse to

More information

Quantum Mechanics I - Session 5

Quantum Mechanics I - Session 5 Quantum Mechanics I - Session 5 Apil 7, 015 1 Commuting opeatos - an example Remine: You saw in class that Â, ˆB ae commuting opeatos iff they have a complete set of commuting obsevables. In aition you

More information

Efficiency of the principal component Liu-type estimator in logistic

Efficiency of the principal component Liu-type estimator in logistic Effcency of the pncpal component Lu-type estmato n logstc egesson model Jbo Wu and Yasn Asa 2 School of Mathematcs and Fnance, Chongqng Unvesty of Ats and Scences, Chongqng, Chna 2 Depatment of Mathematcs-Compute

More information

Chapter 13 - Universal Gravitation

Chapter 13 - Universal Gravitation Chapte 3 - Unesal Gataton In Chapte 5 we studed Newton s thee laws of moton. In addton to these laws, Newton fomulated the law of unesal gataton. Ths law states that two masses ae attacted by a foce gen

More information

Lecture 10: Condensed matter systems

Lecture 10: Condensed matter systems Lectue 0: Codesed matte systems Itoducg matte ts codesed state.! Ams: " Idstgushable patcles ad the quatum atue of matte: # Cosequeces # Revew of deal gas etopy # Femos ad Bosos " Quatum statstcs. # Occupato

More information

Correspondence Analysis & Related Methods

Correspondence Analysis & Related Methods Coespondence Analyss & Related Methods Ineta contbutons n weghted PCA PCA s a method of data vsualzaton whch epesents the tue postons of ponts n a map whch comes closest to all the ponts, closest n sense

More information

Applied Statistical Mechanics Lecture Note - 13 Molecular Dynamics Simulation

Applied Statistical Mechanics Lecture Note - 13 Molecular Dynamics Simulation Appled Statstcal Mechancs Lectue Note - 3 Molecula Dynamcs Smulaton 고려대학교화공생명공학과강정원 Contents I. Basc Molecula Dynamcs Smulaton Method II. Popetes Calculatons n MD III. MD n Othe Ensembles I. Basc MD Smulaton

More information

Atomic units The atomic units have been chosen such that the fundamental electron properties are all equal to one atomic unit.

Atomic units The atomic units have been chosen such that the fundamental electron properties are all equal to one atomic unit. tomc uts The atomc uts have bee chose such that the fudametal electo popetes ae all equal to oe atomc ut. m e, e, h/, a o, ad the potetal eegy the hydoge atom e /a o. D3.33564 0-30 Cm The use of atomc

More information

CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122

CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 19122 CLIMBING THE LADDER OF DENSITY FUNCTIONAL APPROXIMATIONS JOHN P. PERDEW DEPARTMENT OF PHYSICS TEMPLE UNIVERSITY PHILADELPHIA, PA 191 THANKS TO MANY COLLABORATORS, INCLUDING SY VOSKO DAVID LANGRETH ALEX

More information

Matrix Elements of Many-Electron Wavefunctions. noninteger principal quantum number. solutions to Schröd. Eq. outside sphere of radius r

Matrix Elements of Many-Electron Wavefunctions. noninteger principal quantum number. solutions to Schröd. Eq. outside sphere of radius r 30 - Matx Eements of Many-Eecton Wavefunctons Last tme: ν = R En, f ( ν, ) g ( ν, ) need both f and g to satsfy bounday condton fo E < 0 as ν = n µ πµ s phase shft of f ( ν, ) nonntege pncpa quantum numbe

More information

Physical & Interfacial Electrochemistry 2013

Physical & Interfacial Electrochemistry 2013 Physcal & Intefacal Electochemsty 013 Lectue 3. Ion-on nteactons n electolyte solutons. Module JS CH3304 MoleculaThemodynamcs and Knetcs Ion-Ion Inteactons The themodynamc popetes of electolyte solutons

More information

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations

Robert Eisberg Second edition CH 09 Multielectron atoms ground states and x-ray excitations Quantum Physcs 量 理 Robert Esberg Second edton CH 09 Multelectron atoms ground states and x-ray exctatons 9-01 By gong through the procedure ndcated n the text, develop the tme-ndependent Schroednger equaton

More information

Nuclear size corrections to the energy levels of single-electron atoms

Nuclear size corrections to the energy levels of single-electron atoms Nuclea size coections to the enegy levels of single-electon atoms Babak Nadii Nii a eseach Institute fo Astonomy and Astophysics of Maagha (IAAM IAN P. O. Box: 554-44. Abstact A study is made of nuclea

More information

If there are k binding constraints at x then re-label these constraints so that they are the first k constraints.

If there are k binding constraints at x then re-label these constraints so that they are the first k constraints. Mathematcal Foundatons -1- Constaned Optmzaton Constaned Optmzaton Ma{ f ( ) X} whee X {, h ( ), 1,, m} Necessay condtons fo to be a soluton to ths mamzaton poblem Mathematcally, f ag Ma{ f ( ) X}, then

More information

4 SingularValue Decomposition (SVD)

4 SingularValue Decomposition (SVD) /6/00 Z:\ jeh\self\boo Kannan\Jan-5-00\4 SVD 4 SngulaValue Decomposton (SVD) Chapte 4 Pat SVD he sngula value decomposton of a matx s the factozaton of nto the poduct of thee matces = UDV whee the columns

More information

CSU ATS601 Fall Other reading: Vallis 2.1, 2.2; Marshall and Plumb Ch. 6; Holton Ch. 2; Schubert Ch r or v i = v r + r (3.

CSU ATS601 Fall Other reading: Vallis 2.1, 2.2; Marshall and Plumb Ch. 6; Holton Ch. 2; Schubert Ch r or v i = v r + r (3. 3 Eath s Rotaton 3.1 Rotatng Famewok Othe eadng: Valls 2.1, 2.2; Mashall and Plumb Ch. 6; Holton Ch. 2; Schubet Ch. 3 Consde the poston vecto (the same as C n the fgue above) otatng at angula velocty.

More information

24-2: Electric Potential Energy. 24-1: What is physics

24-2: Electric Potential Energy. 24-1: What is physics D. Iyad SAADEDDIN Chapte 4: Electc Potental Electc potental Enegy and Electc potental Calculatng the E-potental fom E-feld fo dffeent chage dstbutons Calculatng the E-feld fom E-potental Potental of a

More information

Let us next consider how one can calculate one- electron elements of the long- range nuclear- electron potential:

Let us next consider how one can calculate one- electron elements of the long- range nuclear- electron potential: Hatee- Fock Long Range In the following assume the patition of 1/ into shot ange and long ance components 1 efc( / a) ef ( / a) + Hee a is a paamete detemining the patitioning between shot ange and long

More information

A Hartree-Fock Example Using Helium

A Hartree-Fock Example Using Helium Univesity of Connecticut DigitalCommons@UConn Chemisty Education Mateials Depatment of Chemisty June 6 A Hatee-Fock Example Using Helium Cal W. David Univesity of Connecticut, Cal.David@uconn.edu Follow

More information

Smoothed Particle Hydrodynamics

Smoothed Particle Hydrodynamics Smooted Patcle Hydodynamcs Applcaton Example Alan Hec Noembe 9, 010 Fluds Nae-Stokes equatons Smooted Patcle Hydodynamcs Smooted Patcle Hydodynamcs Noembe 010 Fluds Lquds, e.g. wate Gasses, e.g. a Plasmas

More information

Distinct 8-QAM+ Perfect Arrays Fanxin Zeng 1, a, Zhenyu Zhang 2,1, b, Linjie Qian 1, c

Distinct 8-QAM+ Perfect Arrays Fanxin Zeng 1, a, Zhenyu Zhang 2,1, b, Linjie Qian 1, c nd Intenatonal Confeence on Electcal Compute Engneeng and Electoncs (ICECEE 15) Dstnct 8-QAM+ Pefect Aays Fanxn Zeng 1 a Zhenyu Zhang 1 b Lnje Qan 1 c 1 Chongqng Key Laboatoy of Emegency Communcaton Chongqng

More information

Chapter 2. A Brief Review of Electron Diffraction Theory

Chapter 2. A Brief Review of Electron Diffraction Theory Chapte. A Bef Revew of Electon Dffacton Theoy 8 Chapte. A Bef Revew of Electon Dffacton Theoy The theoy of gas phase electon dffacton s hadly a new topc. t s well establshed fo decades and has been thooughly

More information

Professor Wei Zhu. 1. Sampling from the Normal Population

Professor Wei Zhu. 1. Sampling from the Normal Population AMS570 Pofesso We Zhu. Samplg fom the Nomal Populato *Example: We wsh to estmate the dstbuto of heghts of adult US male. It s beleved that the heght of adult US male follows a omal dstbuto N(, ) Def. Smple

More information

Lecture 7: Angular Momentum, Hydrogen Atom

Lecture 7: Angular Momentum, Hydrogen Atom Lectue 7: Angula Momentum, Hydogen Atom Vecto Quantization of Angula Momentum and Nomalization of 3D Rigid Roto wavefunctions Conside l, so L 2 2 2. Thus, we have L 2. Thee ae thee possibilities fo L z

More information

HYBRID EXCHANGE CORRELATION FUNCTIONALS AND POTENTIALS: CONCEPT ELABORATION INTRODUCTION. A. V. Arbuznikov UDC

HYBRID EXCHANGE CORRELATION FUNCTIONALS AND POTENTIALS: CONCEPT ELABORATION INTRODUCTION. A. V. Arbuznikov UDC Jounal of Stuctual Chemsty, Vol. 48, Supplement, pp. S1-S31, 007 Ognal Russan Tet Copyght 007 by A. V. Abuznkov HYBRID EXCHANGE CORRELATION FUNCTIONALS AND POTENTIALS: CONCEPT ELABORATION A. V. Abuznkov

More information

Hamiltonian multivector fields and Poisson forms in multisymplectic field theory

Hamiltonian multivector fields and Poisson forms in multisymplectic field theory JOURNAL OF MATHEMATICAL PHYSICS 46, 12005 Hamltonan multvecto felds and Posson foms n multsymplectc feld theoy Mchael Foge a Depatamento de Matemátca Aplcada, Insttuto de Matemátca e Estatístca, Unvesdade

More information

An Approach to Inverse Fuzzy Arithmetic

An Approach to Inverse Fuzzy Arithmetic An Appoach to Invese Fuzzy Athmetc Mchael Hanss Insttute A of Mechancs, Unvesty of Stuttgat Stuttgat, Gemany mhanss@mechaun-stuttgatde Abstact A novel appoach of nvese fuzzy athmetc s ntoduced to successfully

More information

Doublet structure of Alkali spectra:

Doublet structure of Alkali spectra: Doublet stuctue of : Caeful examination of the specta of alkali metals shows that each membe of some of the seies ae closed doublets. Fo example, sodium yellow line, coesponding to 3p 3s tansition, is

More information

9.1 The multiplicative group of a finite field. Theorem 9.1. The multiplicative group F of a finite field is cyclic.

9.1 The multiplicative group of a finite field. Theorem 9.1. The multiplicative group F of a finite field is cyclic. Chapte 9 Pimitive Roots 9.1 The multiplicative goup of a finite fld Theoem 9.1. The multiplicative goup F of a finite fld is cyclic. Remak: In paticula, if p is a pime then (Z/p) is cyclic. In fact, this

More information

A Crash Course in (2 2) Matrices

A Crash Course in (2 2) Matrices A Cash Couse in ( ) Matices Seveal weeks woth of matix algeba in an hou (Relax, we will only stuy the simplest case, that of matices) Review topics: What is a matix (pl matices)? A matix is a ectangula

More information

APPLICATIONS OF SEMIGENERALIZED -CLOSED SETS

APPLICATIONS OF SEMIGENERALIZED -CLOSED SETS Intenatonal Jounal of Mathematcal Engneeng Scence ISSN : 22776982 Volume Issue 4 (Apl 202) http://www.mes.com/ https://stes.google.com/ste/mesounal/ APPLICATIONS OF SEMIGENERALIZED CLOSED SETS G.SHANMUGAM,

More information

P 365. r r r )...(1 365

P 365. r r r )...(1 365 SCIENCE WORLD JOURNAL VOL (NO4) 008 www.scecncewoldounal.og ISSN 597-64 SHORT COMMUNICATION ANALYSING THE APPROXIMATION MODEL TO BIRTHDAY PROBLEM *CHOJI, D.N. & DEME, A.C. Depatment of Mathematcs Unvesty

More information

GENERALIZATION OF AN IDENTITY INVOLVING THE GENERALIZED FIBONACCI NUMBERS AND ITS APPLICATIONS

GENERALIZATION OF AN IDENTITY INVOLVING THE GENERALIZED FIBONACCI NUMBERS AND ITS APPLICATIONS #A39 INTEGERS 9 (009), 497-513 GENERALIZATION OF AN IDENTITY INVOLVING THE GENERALIZED FIBONACCI NUMBERS AND ITS APPLICATIONS Mohaad Faokh D. G. Depatent of Matheatcs, Fedows Unvesty of Mashhad, Mashhad,

More information

Groupoid and Topological Quotient Group

Groupoid and Topological Quotient Group lobal Jounal of Pue and Appled Mathematcs SSN 0973-768 Volume 3 Numbe 7 07 pp 373-39 Reseach nda Publcatons http://wwwpublcatoncom oupod and Topolocal Quotent oup Mohammad Qasm Manna Depatment of Mathematcs

More information

This chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density.

This chapter illustrates the idea that all properties of the homogeneous electron gas (HEG) can be calculated from electron density. 1 Unform Electron Gas Ths chapter llustrates the dea that all propertes of the homogeneous electron gas (HEG) can be calculated from electron densty. Intutve Representaton of Densty Electron densty n s

More information

Multipole Radiation. March 17, 2014

Multipole Radiation. March 17, 2014 Multpole Radaton Mach 7, 04 Zones We wll see that the poblem of hamonc adaton dvdes nto thee appoxmate egons, dependng on the elatve magntudes of the dstance of the obsevaton pont,, and the wavelength,

More information

Three-dimensional systems with spherical symmetry

Three-dimensional systems with spherical symmetry Thee-dimensiona systems with spheica symmety Thee-dimensiona systems with spheica symmety 006 Quantum Mechanics Pof. Y. F. Chen Thee-dimensiona systems with spheica symmety We conside a patice moving in

More information

Introduction to density-functional theory. Emmanuel Fromager

Introduction to density-functional theory. Emmanuel Fromager Institut de Chimie, Strasbourg, France Page 1 Emmanuel Fromager Institut de Chimie de Strasbourg - Laboratoire de Chimie Quantique - Université de Strasbourg /CNRS M2 lecture, Strasbourg, France. Institut

More information