Convergence of DFT eigenvalues with cell volume and vacuum level
|
|
- Beverly Anderson
- 5 years ago
- Views:
Transcription
1 Converence of DFT eienvalues with cell volume and vacuum level Sohrab Ismail-Beii October 4, 2013 Computin work functions or absolute DFT eienvalues (e.. ionization potentials) requires some care. Obviously, we as discussin systems with surfaces in contact with vacuum computin the absolute eienvalue for the interior of an infinite crystal is not possible since the alinment with vacuum depends on the dipoles on the surface which are not bein modeled. Therefore, what we really have in mind are 2D (slabs, sheets), 1D (polymers, nanotubes, nanowires) and 0D (molecules, atoms) where there are some periodic directions and some directions oin into the vacuum. et the lattice vector alon the periodic direction be a and the lattice vector alon the vacuum direction be. We will be investiatin the limit below. In eneral, the DFT Kohn-Sham Hamiltonian can be written as H = 2 2m + φ(r) + V xc(r) where φ is the total electrostatic potential from all chares, ions and electrons, φ(r) = d 3 r ρ(r ) r r = d 3 q 4π ρ(q)e iq r (2π) 3 q 2. Here, ρ(r) is the total (ion+electron) chare density and ρ(q) is its Fourier transform. Most of the problems with determinin absolute eneries and the vacuum level are because the Coulomb interaction is lon raned and so computin the potential in a periodic system of chares can be ill-defined. In could say it is about handlin the q = 0 diverence in periodic systems. Before continuin forward, there are two types of problems in determinin absolute eienvalues and potentials in a DFT calculation, and they are of different oriin. First, some DFT codes will set the averae potential to zero over the unit cell. This not only includes the electrostatic potential but also the exchane-correlation potential V xc. This means that the averae value of V xc over the volume of the unit cell Ω, 1 Ω Ω d 3 r V xc (r), is bein added or subtracted from the output potential. As the volume of the unit cell is chaned by increasin or decreasin the vacuum, this averae chanes as well. Whether a 1
2 iven code does this or not can only be determined by lookin inside. If one is usin a local or semi-local functional (e.. DA or GGA) and has enouh vacuum for the electron density to be localized on the system, then the interal is just some number and the averae will scale like 1/Ω. What this means is that the value of the potential in the vacuum reion could have such a contribution that depends on Ω and must be corrected for or extrapolated away. Second, all DFT codes run on periodic systems will have to do somethin with the diverence at q = 0 and thus adjust the averae potential in some way. This is what we are discussin below. Before we bein, we will need some eneral orientation and one mathematical result. It has to do with the fact that for a function of compact support with well-defined Fourier transform, interatin over the Fourier space or summin over the discrete reciprocal (Fourier) lattice ive exactly the same thin. We will do this in one dimension to keep the math simple. Consider two functions F (x) and G(x) whose product F ()G(x) is of compact support: it is non-zero only for 0 < x < a. We want the interal of their product (inner product) which can be done in real or reciprocal space I = where the Fourier convention is dx F (x) G(x) = F (x) = dq 2π F (q)e iqx. dq 2π F (q) G(q) Instead of interatin over all x, let s take a unit cell of lenth > a. Then the interal of F (x)g(x) over the finite rane 0 < x < is also I. Once we are workin on a finite interval of lenth, we can make F (x) and G(x) periodic without chanin I: F p (x) = n= F (x n) = 1 F ()e ix and similarly for G p (x). Here the reciprocal lattice vectors are = 2πn/ for inteers n. I is unchaned if we use the periodic functions What we found is that I = I = 0 dx F p (x) G p (x) = 1 F () G(). dq 2π F (q) 1 G(q) = F () G(). The discrete sum and the continuous interal are exactly the same: for a function of compact support, one can sample discrete or continuously and no information is lost (Nyquist theorem). We can write this in a slihtly different way which is more useful below. We can choose to write q = k + where π/ < k < π/ is in the first Brillouin zone (BZ) and a reciprocal lattice vector. Therefore I = π/ dk π/ 2π F (k + ) G(k + ). 2
3 Since the BZ has lenth 2π/ we chane replace the interal by the averae over the BZ and write I = 1 F (k + ) G(k + ) BZ where h BZ π/ π/ dk 2π/ h(k). Thus another way to write out main findin is I = 1 F () 1 G() = F (k + ) G(k + ) BZ. Now we bein by focusin on the slab case. We have to periodic directions in the xy plane and one direction oin into the vacuum. So if we let G xy denote the reciprocal lattice vectors in the xy lattice, all functions of interest have periodicity in the xy plane and have a Fourier series in the xy plane: for example, the chare density ρ(r) has Fourier representation ρ(x, y, z) = G xy dq z 2π ρ(g xy, q z ) e igxy rxy+iqzz. A The vector r xy = (x, y) is the position in the plane. A is the area of the periodic cell in xy. Instead of lookin at the potential φ(r) itself, consider the eienvalue E n for some bound stae ψ n (r): [ ] E n = d 3 r ψ n (r) 2 2m + φ(r) + V xc(r) ψ n (r). The part of interest is the electrostatic contribution Φ n = d 3 r φ(r) ψ n (r) 2. To see this, let s define the Fourier transform of ψ n (r) 2 to be p n (q). Then workin in Fourier space Φ n = dqz 4π ρ(g xy, q z ) p n (G xy, q z ) 2π A(G 2 G xy + q 2. xy z) Now, in an actual plane wave calculation, we have a periodic cell alon z of lenth. Both the bound state ψ n and the chare density ρ will decay exponentially into the vacuum outside the slab. Thus for practical purposes, they have compact support. So by choosin a lare enouh, the periodic imaes of these densities will not have any sinificant overlap so that the kinetic and exchane-correlation contributions to E n will be very well convered for any reasonable. It is aain the electrostatic part that ives us trouble. What the plane wave code actually calculates for Φ n uses the discrete reciprocal lattice vectors G z alon z instead of a continuous q z. The computed result from the periodic calculation is Φ p n and is Φ p n = G xy 4π ρ(g xy, G z ) p n (G xy, G z ) A(G 2 G xy + G 2 z z) where G xy = 0, G z = 0 is excluded from the sum. 3
4 (Most plane wave codes will write Ω = A the cell volume instead of A in the denominator.) We have to exclude the 0 wave vector contribution since the summand is very ill defined: the denominator oes to zero; while the numerator ρ(q) also oes to zero for q 0 for a neutral system, it may only o to zero linearly for a system with a dipole moment so the result is poorly defined. This missin term and what to do with it is the basic problem. Usin our mathematical result, for any G xy 0, all the interands or summands are quite smooth so we can use the equivalence of the continuous interal and discrete sum alon z to find dq z 4π ρ(g xy, q z ) p n (G xy, q z ) 2π A(G 2 xy + qz) 2 = 4π ρ(g xy, G z ) p n (G xy, G z ) A(G 2 G xy + G 2. z z) However, for G xy = 0, we can t just blindly apply the mathematical equivalence since the interand is actually diverent for G z = 0 while the interal is enerally well defined for physical systems (see below). Thus far, we have shown that all G xy 0 terms are identical between Φ n and Φ p n so that Φ n = Φ p dq z 4π ρ(q z ) p n (q z ) n + 2π Aqz 2 4π ρ(g z ) p n (G z ) AG 2. G z z 0 To clarify the difference further, we can write an arbitrary q z as the sum of the closest G z plus some k z that is in the first BZ q z = k z + G z as we did in our mathematical result above. We then have the equivalent expression Φ n = Φ p n + 4π ρ(kz ) p n (k z ) A BZ + 4π A [ ρ(kz + G z ) p n (k z + G z ) (k z + G z ) 2 ρ(g ] z) p n (G z ) BZ G 2. z G z 0 To proceed, we need to look at some of the physical characteristics of ρ and p n. We expand both in series ρ(k) = αk + βk 2 + O(k 3 ), p n (k) = 1 + α n k + β n k 2 + O(k 3 ). The expansion for ρ starts at dipolar order because our system is net neutral. The p n expansion starts at 1 since ψ n (r) 2 is normalized. The product will behave like This product divided by k 2 is then ρ(k) p n (k) = αk + Qk 2 + O(k 3 ) + O(k 4 ) + ρ(k) p n (k) k 2 = α k + Q + O(k) + O(k2 ) + This looks badly diverent, but the reciprocal lattice vectors and the BZ have inversion symmetry: for any valid G z, G z is also valid and for any k z in the BZ, so is k z. Thus the BZ averae as well as the sum over G z in the expression for Φ n are both samplin the diverent α/k term equally alon positive and neative values so its contribution actually vanishes around the oriin. Thus the effective summands are actually quite smooth here. What this means is that the sum over G z 0 terms above is very small since we are subtractin the averae of a smooth function over the BZ from the value at the center of the BZ; anyways, the 4
5 difference converes to zero as ets larer and the averain is over a narrower reion. Therefore, any type of discrepancy between Φ n and Φ p is really comin from the very first term: Φ n = Φ p n + 4π ρ(kz ) p n (k z ) A Usin the Taylor series in k from above, the desired averae is 4π ρ(kz ) p n (k z ) A BZ + terms converin very rapidly in. = 4π α BZ A k + Q + Rk + Sk2 + Only the even terms survive, and the averae of k 2 is proportional to (π/) 2. So we et 4π ρ(kz ) p n (k z ) A = C BZ + D 3 + E 5 + for some constants in the numerators. Therefore, we expect for the eienvalue E n () = E n ( ) + C 1 + D 3 + By performin a set of computations at various, we can fit E n () to the above form and extrapolate to = to et E n ( ). (This is an example of Richardson extrapolation.) This concludes the slab eometry. BZ. Considerin the other eometries (1D wire and 0D molecule) is mainly takin the above derivation and chanin some discrete sums to continuous ones. For example, for a wire eometry with periodic axis alon z, the x and y directions are those oin into the vacuum. Thus we will have a discrete sum over G z and interals over q x and q y ; when we put in a periodic box we et a discretized sum over G x and G y instead; the volume factor A is replaced by a 2 where a is the periodic lattice lenth alon z; the G z 0 contributions will be well behaved while G z = 0 will be problematic; for G z = 0, the contributions for (G x, G y ) 0 and those from averain over a BZ centered around (G x, G y ) will be very close; so the difference between Φ n and Φ p n will aain be due to the averain in the BZ around q = 0. The main difference is that now we have vector k so where α is a vector, Q is a matrix, etc., and ρ(k) p n (k) = α T k + k T Qk + O(k 3 ) ρ(k) p n (k) k 2 = αt ˆk k + ˆk T Qˆk + O(ˆkˆkk) + O(ˆkˆkkk) + where ˆk = k/ k is the unit vector. The averain proceeds as before notin that the odd power cancel by symmetry and the even terms ive for the averae of k 2 somethin of order (π/) 2. The overall volume prefactor is now 1/(a 2 ) so our expansion will look like E n () = E n ( ) + C 2 + D 4 + The 0D point eometry has vacuum in all directions so the q interals are continuous alon all three directions; the discretization is to replace q by a three vector G; the volume factor is Ω = 3 ; and we proceed as before. So the series starts with 3 and oes in powers of two from there. 5
6 Our final summary is thus that E n () E n ( ) = constant C 1 + D subleadin terms = E 2 + F 4 + Ω G 3 + H 5 + 2D slab case 1D wire case 0D point case. for constants C, D,... Therefore, we expect the eienvalue to convere to the limit as E n () = E n ( ) + C m + D m+2 + where m = 1 for 2D, m = 2 for 1D, and m = 3 for 0D. By performin a set of computations at various, we can fit E n () to the above form and extrapolate to = to et E n ( ). (This is an example of Richardson extrapolation.) 6
Another possibility is a rotation or reflection, represented by a matrix M.
1 Chapter 25: Planar defects Planar defects: orientation and types Crystalline films often contain internal, 2-D interfaces separatin two reions transformed with respect to one another, but with, otherwise,
More informationJournal of Inequalities in Pure and Applied Mathematics
Journal of Inequalities in Pure and Applied Mathematics L HOSPITAL TYPE RULES FOR OSCILLATION, WITH APPLICATIONS IOSIF PINELIS Department of Mathematical Sciences, Michian Technoloical University, Houhton,
More informationConsider a square wave potential with lattice constant a. In the first unit cell, this means:
Consider a square wave potential with lattice constant a. In the first unit cell, this means: = for - a/4 < x < a/4, =0 for -a/ < x < - a/4 and a/4 < x < a/, This is repeated throuhout the crystal. 1.
More informationWhat types of isometric transformations are we talking about here? One common case is a translation by a displacement vector R.
1. Planar Defects Planar defects: orientation and types Crystalline films often contain internal, -D interfaces separatin two reions transformed with respect to one another, but with, otherwise, essentially,
More informationRenormalization Group Theory
Chapter 16 Renormalization Group Theory In the previous chapter a procedure was developed where hiher order 2 n cycles were related to lower order cycles throuh a functional composition and rescalin procedure.
More informationEfficient method for obtaining parameters of stable pulse in grating compensated dispersion-managed communication systems
3 Conference on Information Sciences and Systems, The Johns Hopkins University, March 12 14, 3 Efficient method for obtainin parameters of stable pulse in ratin compensated dispersion-manaed communication
More informationIntegration in General Relativity
arxiv:physics/9802027v1 [math-ph] 14 Feb 1998 Interation in General Relativity Andrew DeBenedictis Dec. 03, 1995 Abstract This paper presents a brief but comprehensive introduction to certain mathematical
More informationCh 14: Feedback Control systems
Ch 4: Feedback Control systems Part IV A is concerned with sinle loop control The followin topics are covered in chapter 4: The concept of feedback control Block diaram development Classical feedback controllers
More informationPhysics 20 Homework 1 SIMS 2016
Physics 20 Homework 1 SIMS 2016 Due: Wednesday, Auust 17 th Problem 1 The idea of this problem is to et some practice in approachin a situation where you miht not initially know how to proceed, and need
More informationExam 2A Solution. 1. A baseball is thrown vertically upward and feels no air resistance. As it is rising
Exam 2A Solution 1. A baseball is thrown vertically upward and feels no air resistance. As it is risin Solution: Possible answers: A) both its momentum and its mechanical enery are conserved - incorrect.
More informationDynamical Diffraction
Dynamical versus Kinematical Di raction kinematical theory valid only for very thin crystals Dynamical Diffraction excitation error s 0 primary beam intensity I 0 intensity of other beams consider di racted
More information11. Excitons in Nanowires and Nanotubes
Excitons in Nanowires and Nanotubes We have seen that confinement in quantum wells leads to enhanced excitonic effects in the optical response of semiconductors The bindin enery of the stronest bound excitons
More informationStochastic simulations of genetic switch systems
Stochastic simulations of enetic switch systems Adiel Loiner, 1 Azi Lipshtat, 2 Nathalie Q. Balaban, 1 and Ofer Biham 1 1 Racah Institute of Physics, The Hebrew University, Jerusalem 91904, Israel 2 Department
More informationReciprocal Symmetric Boltzmann Function and Unified Boson-Fermion Statistics
Reciprocal Symmetric Boltzmann Function and Unified Boson-Fermion Statistics Mushfiq Ahmad Department of Physics, Rajshahi University, Rajshahi, Banladesh E-mail: mushfiqahmad@ru.ac.bd Muhammad O. G. Talukder
More information2.3. PBL Equations for Mean Flow and Their Applications
.3. PBL Equations for Mean Flow and Their Applications Read Holton Section 5.3!.3.1. The PBL Momentum Equations We have derived the Reynolds averaed equations in the previous section, and they describe
More informationBoundary value problems - an introduction
Boundary value problems - an introduction Up till now, we ve been worryin about initial value problems in DEs. Also quite common, and enerally much more difficult, are boundary value problems. One oriin
More informationIntroduction to density functional perturbation theory for lattice dynamics
Introduction to density functional perturbation theory for lattice dynamics SISSA and DEMOCRITOS Trieste (Italy) Outline 1 Lattice dynamic of a solid: phonons Description of a solid Equations of motion
More informationV DD. M 1 M 2 V i2. V o2 R 1 R 2 C C
UNVERSTY OF CALFORNA Collee of Enineerin Department of Electrical Enineerin and Computer Sciences E. Alon Homework #3 Solutions EECS 40 P. Nuzzo Use the EECS40 90nm CMOS process in all home works and projects
More informationC. Non-linear Difference and Differential Equations: Linearization and Phase Diagram Technique
C. Non-linear Difference and Differential Equations: Linearization and Phase Diaram Technique So far we have discussed methods of solvin linear difference and differential equations. Let us now discuss
More informationNANO 703-Notes. Chapter 12-Reciprocal space
1 Chapter 1-Reciprocal space Conical dark-field imain We primarily use DF imain to control imae contrast, thouh STEM-DF can also ive very hih resolution, in some cases. If we have sinle crystal, a -DF
More information5 Shallow water Q-G theory.
5 Shallow water Q-G theory. So far we have discussed the fact that lare scale motions in the extra-tropical atmosphere are close to eostrophic balance i.e. the Rossby number is small. We have examined
More informationCourse , General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 8: Lorenz Energy Cycle
Course.8, General Circulation of the Earth's Atmosphere Prof. Peter Stone Section 8: Lorenz Enery Cycle Enery Forms: As we saw in our discussion of the heat budet, the enery content of the atmosphere per
More informationLinearized optimal power flow
Linearized optimal power flow. Some introductory comments The advantae of the economic dispatch formulation to obtain minimum cost allocation of demand to the eneration units is that it is computationally
More information1 CHAPTER 7 PROJECTILES. 7.1 No Air Resistance
CHAPTER 7 PROJECTILES 7 No Air Resistance We suppose that a particle is projected from a point O at the oriin of a coordinate system, the y-axis bein vertical and the x-axis directed alon the round The
More information6 Mole Concept. g mol. g mol. g mol ) + 1( g : mol ratios are the units of molar mass. It does not matter which unit is on the
What is a e? 6 Mole Concept The nature of chemistry is to chane one ecule into one or more new ecules in order to create new substances such as paints, fertilizers, food additives, medicines, etc. When
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 informationExperiment 3 The Simple Pendulum
PHY191 Fall003 Experiment 3: The Simple Pendulum 10/7/004 Pae 1 Suested Readin for this lab Experiment 3 The Simple Pendulum Read Taylor chapter 5. (You can skip section 5.6.IV if you aren't comfortable
More informationMASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics. Physics 8.01x Fall Term 2001 EXAM 1 SOLUTIONS
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics Physics 8.01x Fall Term 2001 EXAM 1 SOLUTIONS Problem 1: We define a vertical coordinate system with positive upwards. The only forces actin
More informationNotes: Most of the material presented in this chapter is taken from Jackson, Chap. 2, 3, and 4, and Di Bartolo, Chap. 2. 2π nx i a. ( ) = G n.
Chapter. Electrostatic II Notes: Most of the material presented in this chapter is taken from Jackson, Chap.,, and 4, and Di Bartolo, Chap... Mathematical Considerations.. The Fourier series and the Fourier
More informationFirst- and Second Order Phase Transitions in the Holstein- Hubbard Model
Europhysics Letters PREPRINT First- and Second Order Phase Transitions in the Holstein- Hubbard Model W. Koller 1, D. Meyer 1,Y.Ōno 2 and A. C. Hewson 1 1 Department of Mathematics, Imperial Collee, London
More informationProblem Set 2 Solutions
UNIVERSITY OF ALABAMA Department of Physics and Astronomy PH 125 / LeClair Sprin 2009 Problem Set 2 Solutions The followin three problems are due 20 January 2009 at the beinnin of class. 1. (H,R,&W 4.39)
More informationSTOCHASTICALLY GENERATED MULTIGROUP DIFFUSION COEFFICIENTS
STOCHASTICALLY GENERATED MULTIGROUP DIFFUSION COEFFICIENTS A Thesis Presented to The Academic Faculty by Justin M. Pounders In Partial Fulfillment of the Requirements for the Deree Master of Science in
More informationAdvanced Methods Development for Equilibrium Cycle Calculations of the RBWR. Andrew Hall 11/7/2013
Advanced Methods Development for Equilibrium Cycle Calculations of the RBWR Andrew Hall 11/7/2013 Outline RBWR Motivation and Desin Why use Serpent Cross Sections? Modelin the RBWR Axial Discontinuity
More informationPhysicsAndMathsTutor.com
. A particle of mass m is projected vertically upwards, at time t =, with speed. The particle is mv subject to air resistance of manitude, where v is the speed of the particle at time t and is a positive
More informationPhysics 18 Spring 2011 Homework 2 - Solutions Wednesday January 26, 2011
Physics 18 Sprin 011 Homework - s Wednesday January 6, 011 Make sure your name is on your homework, and please box your final answer. Because we will be ivin partial credit, be sure to attempt all the
More informationAVERAGE AND STANDARD DEVIATION
AVERAGE AND STANDARD DEVIATION Link to: physicspaes home pae. To leave a comment or report an error, please use the auxiliary blo. Post date: 2 Jun 25. Reference: Griffiths, David J. (25), Introduction
More information9. v > 7.3 mi/h x < 2.5 or x > x between 1350 and 5650 hot dogs
.5 Etra Practice. no solution. (, 0) and ( 9, ). (, ) and (, ). (, 0) and (, 0) 5. no solution. ( + 5 5 + 5, ) and ( 5 5 5, ) 7. (0, ) and (, 0). (, ) and (, 0) 9. (, 0) 0. no solution. (, 5). a. Sample
More informationModeling for control of a three degrees-of-freedom Magnetic. Levitation System
Modelin for control of a three derees-of-freedom Manetic evitation System Rafael Becerril-Arreola Dept. of Electrical and Computer En. University of Toronto Manfredi Maiore Dept. of Electrical and Computer
More information4.3. Solving Friction Problems. Static Friction Problems. Tutorial 1 Static Friction Acting on Several Objects. Sample Problem 1.
Solvin Friction Problems Sometimes friction is desirable and we want to increase the coefficient of friction to help keep objects at rest. For example, a runnin shoe is typically desined to have a lare
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 informationExperiment 4: Electrostatic Force
MASSACHUSETTS INSTITUTE OF TECHNOLOGY Department of Physics 8. Sprin 3 OBJECTIVE Experiment 4: Electrostatic Force To measure ε, the permittivity of free space. INTRODUCTION Electrostatic force plays a
More informationStrong Interference and Spectrum Warfare
Stron Interference and Spectrum Warfare Otilia opescu and Christopher Rose WILAB Ruters University 73 Brett Rd., iscataway, J 8854-86 Email: {otilia,crose}@winlab.ruters.edu Dimitrie C. opescu Department
More information#A52 INTEGERS 15 (2015) AN IMPROVED UPPER BOUND FOR RAMANUJAN PRIMES
#A5 INTEGERS 15 (015) AN IMPROVED UPPER BOUND FOR RAMANUJAN PRIMES Anitha Srinivasan Dept. of Mathematics, Saint Louis University Madrid Campus, Madrid, Spain asriniv@slu.edu John W. Nicholson Decatur,
More informationCmpt 250 Unsigned Numbers January 11, 2008
Cmt 25 Unsined Numbers Januar, 28 These notes serve two uroses in the contet of Cmt 25: as we develo the basic desin of an inteer ALU, we ll review the basics of number reresentation and combinational
More informationLP Rounding and Combinatorial Algorithms for Minimizing Active and Busy Time
LP Roundin and Combinatorial Alorithms for Minimizin Active and Busy Time Jessica Chan, Samir Khuller, and Koyel Mukherjee University of Maryland, Collee Park {jschan,samir,koyelm}@cs.umd.edu Abstract.
More informationWhat Causes Image Intensity Changes?
Ede Detection l Why Detect Edes? u Information reduction l Replace imae by a cartoon in which objects and surface markins are outlined create line drawin description l These are the most informative parts
More informationMatrix multiplication: a group-theoretic approach
CSG399: Gems of Theoretical Computer Science. Lec. 21-23. Mar. 27-Apr. 3, 2009. Instructor: Emanuele Viola Scribe: Ravi Sundaram Matrix multiplication: a roup-theoretic approach Given two n n matrices
More informationCalculation of Reflection and Transmission Coefficients in scuff-transmission
Calculation of Reflection and Transmission Coefficients in scuff-transmission Homer Reid May 9, 2015 Contents 1 The Setup 2 2 Scattering coefficients from surface currents 4 2.1 Computation of b(q).........................
More informationULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor
ULOF Accident Analysis for 300 MWt Pb-Bi Coolled MOX Fuelled SPINNOR Reactor Ade afar Abdullah Electrical Enineerin Department, Faculty of Technoloy and Vocational Education Indonesia University of Education
More informationFig.1: Non-stationary temperature distribution in a circular plate.
LECTURE Problem 1 13. 5 : Fi.1: Non-stationary temperature distribution in a circular plate. The non-stationary, radial symmetric temperature distribution ur, t in a circular plate of radius c is determined
More informationCentral Limit Theorems and Proofs
Math/Stat 394, Winter 0 F.W. Scholz Central Limit Theorems and Proos The ollowin ives a sel-contained treatment o the central limit theorem (CLT). It is based on Lindeber s (9) method. To state the CLT
More informationRESISTANCE STRAIN GAGES FILLAMENTS EFFECT
RESISTANCE STRAIN GAGES FILLAMENTS EFFECT Nashwan T. Younis, Younis@enr.ipfw.edu Department of Mechanical Enineerin, Indiana University-Purdue University Fort Wayne, USA Bonsu Kan, kan@enr.ipfw.edu Department
More informationGeopotential tendency and vertical motion
Geopotential tendency and vertical motion Recall PV inversion Knowin the PV, we can estimate everythin else! (Temperature, wind, eopotential ) In QG, since the flow is eostrophic, we can obtain the wind
More informationParametric Equations
Parametric Equations Suppose a cricket jumps off of the round with an initial velocity v 0 at an anle θ. If we take his initial position as the oriin, his horizontal and vertical positions follow the equations:
More informationAccording to the Ratio of Height to Least Lateral Dimension: (ACI Code and ACI Code ) < 22, when member is braced against sideways.
CHAPTER 5 COLUMN 5.1 INTRODUCTION Columns are structural elements used primarily to support compressive loads. All practical columns are members subjected not only to axial load but also to bendin moments,
More informationMechanics Cycle 3 Chapter 12++ Chapter 12++ Revisit Circular Motion
Chapter 12++ Revisit Circular Motion Revisit: Anular variables Second laws for radial and tanential acceleration Circular motion CM 2 nd aw with F net To-Do: Vertical circular motion in ravity Complete
More informationAltitude measurement for model rocketry
Altitude measurement for model rocketry David A. Cauhey Sibley School of Mechanical Aerospace Enineerin, Cornell University, Ithaca, New York 14853 I. INTRODUCTION In his book, Rocket Boys, 1 Homer Hickam
More information2.2 Differentiation and Integration of Vector-Valued Functions
.. DIFFERENTIATION AND INTEGRATION OF VECTOR-VALUED FUNCTIONS133. Differentiation and Interation of Vector-Valued Functions Simply put, we differentiate and interate vector functions by differentiatin
More informationCHEM6085: Density Functional Theory
Lecture 11 CHEM6085: Density Functional Theory DFT for periodic crystalline solids C.-K. Skylaris 1 Electron in a one-dimensional periodic box (in atomic units) Schrödinger equation Energy eigenvalues
More informationPractical 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 informationf 1. (8.1.1) This means that SI unit for frequency is going to be s 1 also known as Hertz d1hz
ecture 8-1 Oscillations 1. Oscillations Simple Harmonic Motion So far we have considered two basic types of motion: translational motion and rotational motion. But these are not the only types of motion
More informationDurham Research Online
Durham Research Online Deposited in DRO: 15 March 2011 Version of attached file: Published Version Peer-review status of attached file: Peer-reviewed Citation for published item: Brand, S. and Kaliteevski,
More informationThis relationship is known as the ideal gas law and is mathematically described with the formula below:
Chemistry 20 Ideal as law If we combine all the information contained in Boyle s, Charles and Avoadro s laws, we can derive an expression that describes the temperature, pressure and volume of a as. This
More informationSupplementary Information for: Exciton Radiative Lifetimes in. Layered Transition Metal Dichalcogenides
Supplementary Information for: Exciton Radiative Lifetimes in Layered Transition Metal Dichalcogenides Maurizia Palummo,, Marco Bernardi,, and Jeffrey C. Grossman, Dipartimento di Fisica, Università di
More informationHomework # 2. SOLUTION - We start writing Newton s second law for x and y components: F x = 0, (1) F y = mg (2) x (t) = 0 v x (t) = v 0x (3)
Physics 411 Homework # Due:..18 Mechanics I 1. A projectile is fired from the oriin of a coordinate system, in the x-y plane (x is the horizontal displacement; y, the vertical with initial velocity v =
More informationTheoretische Physik 2: Elektrodynamik (Prof. A-S. Smith) Home assignment 11
WiSe 22..23 Prof. Dr. A-S. Smith Dipl.-Phys. Matthias Saba am Lehrstuhl für Theoretische Physik I Department für Physik Friedrich-Alexander-Universität Erlangen-Nürnberg Problem. Theoretische Physik 2:
More informationJ.L. Kirtley Jr. September 4, 2010
Massachusetts Institute of Technoloy Department of Electrical Enineerin and Computer Science 6.007 Electromanetic Enery: From Motors to Lasers Supplemental Class Notes Manetic Circuit Analo to Electric
More informationA-Level Mathematics. MM05 Mechanics 5 Final Mark scheme June Version/Stage: v1.0
A-Level Mathematics MM0 Mechanics Final Mark scheme 6360 June 07 Version/Stae: v.0 Mark schemes are prepared by the Lead Assessment Writer and considered, toether with the relevant questions, by a panel
More informationCode_Aster. Element CABLE_GAINE
Titre : Élément CABLE_GAINE Date : 28/07/2015 Pae : 1/12 Element CABLE_GAINE Summary: The element CABLE_GAINE presented in this document is to model cables of prestressin bein able to rub or slip into
More informationMATHCHEM 1.0. By Madhavan Narayanan Graduate Student, Department of Chemistry, Temple University
MATHCHEM.0 By Madhavan Narayanan Graduate Student, Department of Chemistry, Temple University Preface I dedicate this document to my beloved parents and my teachers without whom I would not have had the
More informationMassachusetts Institute of Technology Department of Electrical Engineering and Computer Science
Massachusetts nstitute of Technoloy Department of Electrical Enineerin and Computer Science 6.685 Electric Machines Class Notes 2 Manetic Circuit Basics September 5, 2005 c 2003 James L. Kirtley Jr. 1
More information7. QCD. Particle and Nuclear Physics. Dr. Tina Potter. Dr. Tina Potter 7. QCD 1
7. QCD Particle and Nuclear Physics Dr. Tina Potter Dr. Tina Potter 7. QCD 1 In this section... The stron vertex Colour, luons and self-interactions QCD potential, confinement Hadronisation, jets Runnin
More information1. According to the U.S. DOE, each year humans are. of CO through the combustion of fossil fuels (oil,
1. Accordin to the U.S. DOE, each year humans are producin 8 billion metric tons (1 metric ton = 00 lbs) of CO throuh the combustion of fossil fuels (oil, coal and natural as). . Since the late 1950's
More informationIsotropic diffeomorphisms: solutions to a differential system for a deformed random fields study
Isotropic diffeomorphisms: solutions to a differential system for a deformed random fields study Marc Briant, Julie Fournier To cite this version: Marc Briant, Julie Fournier. Isotropic diffeomorphisms:
More informationMODULI SPACES OF COHERENT SYSTEMS OF SMALL SLOPE ON ALGEBRAIC CURVES
MODULI SPACES OF COHERENT SYSTEMS OF SMALL SLOPE ON ALGEBRAIC CURVES S. B. BRADLOW, O. GARCÍA-PRADA, V. MERCAT, V. MUÑOZ, AND P. E. NEWSTEAD Abstract. Let C be an alebraic curve of enus 2. A coherent system
More informationUsing Quasi-Newton Methods to Find Optimal Solutions to Problematic Kriging Systems
Usin Quasi-Newton Methos to Fin Optimal Solutions to Problematic Kriin Systems Steven Lyster Centre for Computational Geostatistics Department of Civil & Environmental Enineerin University of Alberta Solvin
More informationarxiv: v2 [astro-ph.co] 7 Feb 2017
Nonzero Density-Velocity Consistency Relations for Lare Scale Structures Luca Alberto Rizzo 1,2,3, David F. Mota 3, Patrick Valaeas 1,2 1 Institut de Physique Théorique, CEA, IPhT, F-91191 Gif-sur-Yvette,
More informationAccuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation
Accuracy benchmarking of DFT results, domain libraries for electrostatics, hybrid functional and solvation Stefan Goedecker Stefan.Goedecker@unibas.ch http://comphys.unibas.ch/ Multi-wavelets: High accuracy
More informationSolution sheet 9. λ < g u = 0, λ = g α 0 : u = αλ. C(u, λ) = gλ gλ = 0
Advanced Finite Elements MA5337 - WS17/18 Solution sheet 9 On this final exercise sheet about variational inequalities we deal with nonlinear complementarity functions as a way to rewrite a variational
More informationKernel B Splines and Interpolation
Kernel B Splines and Interpolation M. Bozzini, L. Lenarduzzi and R. Schaback February 6, 5 Abstract This paper applies divided differences to conditionally positive definite kernels in order to generate
More informationPhysics 11 Fall 2012 Practice Problems 2 - Solutions
Physics 11 Fall 01 Practice Problems - s 1. True or false (inore any effects due to air resistance): (a) When a projectile is fired horizontally, it takes the same amount of time to reach the round as
More informationv( t) g 2 v 0 sin θ ( ) ( ) g t ( ) = 0
PROJECTILE MOTION Velocity We seek to explore the velocity of the projectile, includin its final value as it hits the round, or a taret above the round. The anle made by the velocity vector with the local
More informationWater Formation on Pt and Pt-based Alloys: A TheoreticalDescription of a Catalytic Reaction
DOI: 10.1002/cphc.200500613 Water Formation on Pt and Pt-based Alloys: A TheoreticalDescription of a Catalytic Reaction Timo Jacob* [a] and William A. Goddard III [b] In the past, the modelinof catalytic
More informationAssignment 6. Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question.
Assinment 6 Name MULTIPLE CHOICE. Choose the one alternative that best completes the statement or answers the question. Round your answer, if appropriate. 1) A man 6 ft tall walks at a rate of 3 ft/sec
More informationATOMIC BONDING Atomic Bonding
ATOMIC BONDING Atomic Bonding Primary Bonds Secondary Bonds Ionic Covalent Metallic van der Waals 1. IONIC BONDING q 11 Na & 17 Cl These two ions are attracted to eachother by the electrostatic force developed
More informationTopic 9-1: Bloch Theorem and the Central Equation Kittel Pages:
Topic 9-1: Bloch Theorem and the Central Equation Kittel Pages: 167-174 Summary: We begin here by postulating Bloch s theorems which develop the form of the wavefunction in a periodic solid. We then show
More information(a) 1m s -2 (b) 2 m s -2 (c) zero (d) -1 m s -2
11 th Physics - Unit 2 Kinematics Solutions for the Textbook Problems One Marks 1. Which one of the followin Cartesian coordinate system is not followed in physics? 5. If a particle has neative velocity
More informationGLY Geomorphology Notes
GLY 5705 - Geomorpholoy Notes Dr. Peter N. Adams Revised: Oct. 2012 24 Hillslope Processes Associated Readins: Anderson and Anderson (2010), pp. 313-328 To examine processes operatin on hillslopes, we
More informationChem/Biochem 471 Exam 3 12/18/08 Page 1 of 7 Name:
Che/Bioche 47 Exa /8/08 Pae of 7 Please leave the exa paes stapled toether. The forulas are on a separate sheet. This exa has 5 questions. You ust answer at least 4 of the questions. You ay answer ore
More informationbefore the collision and v 1 f and v 2 f after the collision. Since conservation of the linear momentum
Lecture 7 Collisions Durin the preious lecture we stared our discussion of collisions As it was stated last tie a collision is an isolated eent in which two or ore odies (the collidin odies) exert relatiely
More informationExperiment 1: Simple Pendulum
COMSATS Institute of Information Technoloy, Islamabad Campus PHY-108 : Physics Lab 1 (Mechanics of Particles) Experiment 1: Simple Pendulum A simple pendulum consists of a small object (known as the bob)
More informationStat260: Bayesian Modeling and Inference Lecture Date: March 10, 2010
Stat60: Bayesian Modelin and Inference Lecture Date: March 10, 010 Bayes Factors, -priors, and Model Selection for Reression Lecturer: Michael I. Jordan Scribe: Tamara Broderick The readin for this lecture
More informationWire antenna model of the vertical grounding electrode
Boundary Elements and Other Mesh Reduction Methods XXXV 13 Wire antenna model of the vertical roundin electrode D. Poljak & S. Sesnic University of Split, FESB, Split, Croatia Abstract A straiht wire antenna
More informationChapter K. Oscillatory Motion. Blinn College - Physics Terry Honan. Interactive Figure
K. - Simple Harmonic Motion Chapter K Oscillatory Motion Blinn Collee - Physics 2425 - Terry Honan The Mass-Sprin System Interactive Fiure Consider a mass slidin without friction on a horizontal surface.
More informationDisclaimer: This lab write-up is not
Disclaimer: This lab write-up is not to be copied, in whole or in part, unless a proper reference is made as to the source. (It is stronly recommended that you use this document only to enerate ideas,
More informationThe Bionic Lightweight Design of the Mid-rail Box Girder Based on the Bamboo Structure
Wenmin Chen 1, Weian Fu 1, Zeqian Zhan 1, Min Zhan 2 Research Institute of Mechanical Enineerin, Southwest Jiaoton University, Chendu, Sichuan, China (1), Department of Industrial and System Enineerin,
More informationINFLUENCE OF TUBE BUNDLE GEOMETRY ON HEAT TRANSFER TO FOAM FLOW
HEFAT7 5 th International Conference on Heat Transfer, Fluid Mechanics and Thermodynamics Sun City, South Africa Paper number: GJ1 INFLUENCE OF TUBE BUNDLE GEOMETRY ON HEAT TRANSFER TO FOAM FLOW Gylys
More informationFoundation on Compressible Fluid Flow
Chapter Foundation on Compressible Fluid Flow. COMRESSIBLE FLUIDS In everyday life one reconizes three states of matter, namely, solid, liquid and aseous. Solids, liquids and ases are all comprised of
More informationShort Sample Solutions to the Sample Exam for (3rd Year Course 6) Hilary Term 2011
Short Sample Solutions to the Sample Exam for (3rd Year Course 6) Hilary Term 0. [4] A Lattice is an infinite set of points in space where the environment of any given point is identical to the enivironment
More informationIntroduction to Parallelism in CASTEP
to ism in CASTEP Stewart Clark Band University of Durham 21 September 2012 Solve for all the bands/electrons (Band-) Band CASTEP solves the Kohn-Sham equations for electrons in a periodic array of nuclei:
More information