First-principles calculations of insulators in a. finite electric field
|
|
- Edwin Tyler
- 5 years ago
- Views:
Transcription
1 Université de Liège First-principles calculations of insulators in a finite electric field M. Veithen, I. Souza, J. Íñiguez, D. Vanderbilt, K. M. Rabe and h. Ghosez Supported by: FNRS Belgium, VW Stiftung, Communauté Française & Rutgers University Many thanks: Donald R. Hamann
2 Outline Introduction: Difficulties related to the electric field perturbation 1 st part: Structural response to macroscopic electric fields Formalism Berry phase calculation of and the ddk Implementation In practice: how to optimize the structure for a fixed polarization Applications nd 2 part: First-principles approach to insulators in finite electric fields Formalism Definition of an electric field dependent energy functional For < c : minimization of the energy functional Computation of forces and stresses Implementation in ABINIT In practice: how to do a finite electric field calculation Applications Conclusions
3 Nonperiodic lectric field perturbation: difficulties The perturbation - r is nonperiodic and unbound from below Methods based on Bloch's theorem and eigenstates Unbound from below nk do not apply nergy can always be lowered by transferring charge from valence states in one region to conduction states in a distant region c v e - L g lectric field bends the energy bands Dielectric breakdown by interband (Zener) tunneling An infinite crystal in the presence of an electric field does not have a ground-state
4 1 st Method: Structural response to macroscopic electric fields lectric field derivatives of arbitrary order can be computed from DFT R. W. Nunes and X. Gonze, RB 63, (2001) Low-order Taylor expansion of the energy with respect to the electric field Dependence of the energy on the structural degrees of freedom is preserved to all orders Two step approach: 1) Map out the energy as a function of the polarization 2) Use this energy surface & field-coupling term to compute the ground state structure in the presence of the electric field More informations Na Sai, K. M. Rabe and D. Vanderbilt, RB 66, (2002) H. Fu and R.. Cohen, Nature 403, 281 (2000)
5 " #! Constrained polarization approach: formalism Structural response at constant electric field F min R, F R,, olarization: thermodynamic conjugate of : Legendre transformation ( ): R eq eq Strain and atomic positions R at fixed electric field 1 F F R, F, min R, min F R, F R,,, Strain and atomic positions at fixed polarization: R eq & eq Inverse Legendre transformation F min F
6 $ & % A 4 < 0 = 1, + ) * ; / > 2 B 3 9 ' J L H F G K I ( Q U M T O N Constrained polarization approach: formalism Taylor expansion of F(R, ) around = 0 F R,, F R,, 0 -. F R,, , 7 8 : 2 F R,?, 0 CD Truncation of the Taylor expansion at the lowest order F 1 R,, F R,, 0 R,, 0 Approximation is supposed to be valid in systems where the ionic contribution to the polarization dominates (f. ex. ferroelectrics) Resulting expression F min, R, F R,, 0 RS R,, 0
7 V ^ Z _ [ ] ` \ f d j e i o p W Y X ~ ƒ u } v s q r t Ž ˆ š Constrained polarization approach: formalism System of equations to be solved F R,, 0 R F R,, 0 fhg R,, 0 R R,, 0 khl 0 acb 0 mcn R,, 0 In practice: Taylor expansion of F(R,,0) and (R,,0) F R, R,, 0, 0 F R 0, R 0, 0, 0 0, 0 f R 1 ŠZ Œ wcx R y{z e v1 2 K { R 2 2 C { 2 R { R 0, 0 : trial guess of the initial coordinates and strains ( f, : Hellmann-Feynmann forces and stresses K : interatomic force constants C : rigid atom elastic constants : coupling parameters between R and Z * : Born effective charges e : clamped ion piezoelectric tensor R = R R 0 & = - 0 )
8 Ÿ ª ž œ µ ¹ º ½ ¼» Á À ¾ Constrained polarization approach: formalism Linear systems of equations 1 K Z 1 C e 0 Z e R «f ± ² ³ = difference between initial and target values of Iterative solution Choose R 0 & 0 Compute f, (Hellmann-Feynmann theorem) (Berry phase) K, C,, e & Z * (linear response) Solve linear system of equations to obtain new atomic positions and strains R 0 + R R ÂÄà repeat until converged
9 Ñ Ð Å Ò Æ Ö Ô Ø Ù Ó Ç ã á Ý ß à ä â È É å è ë ê æ é ì Ê Ì Berry phase calculation of and the ddk Definition of the Born effective charges Z Computation of : string averaged Berry phase ÍÏÎ R k-b k k+b 2 2 e Õ3 Computation of Z * (linear response, non-stationarry expression) Z ÚÜÛ 2 i e ROBLM: If the ddk is computed from linear response, the relation between Z * and is only satisfied in a the limit of a dense k-point mesh Finite k-point grid: finite difference formula of the ddk 2 dk Þ3 dk ln occ m det S k, k R k b Comments: Ë No violoation of the charge neutrality k ç1 2b b b b b Same arguments in case of the piezoelectric tensor
10 í î ï ñ ò ó ô õ ö ø Implementation ABINIT: new Berry phase routine (berryphase_new.f) asier to use olarization in cartesian coordinates ð MI parallelization over k-points ANADDB: Structure of the DDB Header Total energy 1 st -order energy derivatives 2 nd -order energy derivatives New routine relaxpol.f Called from anaddb.f Compute in cartesian coordinates Solve linear system of equations Compute residual forces and stresses Update atomic positions and lattice constants
11 ù ú û ü ý þ ÿ In practice (v.4.3.x) 1 st step: Ground state calculation of the forces and the stress tensor nd 2 step: Berryphase calculation of and the ddk kptopt = 3 berryopt = -3: use berryphase_new.f routine rfdir = 1 1 1: compute projection of & ddk along x, y and z nband = number of occupied bands th 3 step: Linear response calculation th 4 step: Use ANADDB to compute new atomic positions and lattice constants polflag = 1 targetpol = target value of the polarization (cartesian coordinates & C/m 2 ) relaxat, relaxstr natfix, iatfix, istrfix specify which degrees of freedom are allowed to relax Repeat steps 1, 2, 3 & 4 until convergence is reached
12 Application: double well in BaTiO3
13 Application: c/a of btio 3 Motivation: btio 3 thin films grown on a SrTiO 3 substrate as a function of SrTiO 3 btio 3 xperiment: c/a decreases with decreasing film thickness Theory (effective hamiltonian): depolarizing electric field reduction of the polarization reduction of c/a Result a has been fixed at the theoretical lattice constant of SrTiO 3 C. Lichtensteiger et al., cond-mat/
14 Application: non-linear dielectric response of tetragonal btio 3 z as a function of an electric field along z Ionic dielectric constant: = (linear response = 17.37) Tunability: d zz d z pm V
15 2 nd method: First-principles approach to insulators in finite electric fields For a supercell of size L = N k a (a = unit cell size, N k = number of k- points), use only fields smaller than For < c, we can minimize the energy functional c g N k a F F KS a 3 FKS = Kohn-Sham energy at zero electric field = macroscopic polarization that is computed as a Berry phase of field polarized Bloch functions Long-lived metastable state For > c, F has no minimum More informations I. Souza, J. ĺñiguez and D. Vanderbilt, hys. Rev. Lett. 89, (2002)
16 ; 8 9 : = < # ' & % $! " 5. ( 7 3, - 4 / 6 Minimization of the energy functional F Strategy: use the preconditione cg minimization implemented in ABINIT to minimiez F KS Gradient of the energy functional F G nk F > G nk H KS k W nk usual zero-field term electric field term lectric field term k-b k k+b W nk )+* C m 1 S mn k, k b u mk b S mn k, k b u mk b Transform G nk into a preconditioned cg search direction D nk
17 @ A? I R Q N M L J O K C B H G Minimization of the energy functional F Update of a state new cos old sin D nk Minimize F( ) D = 0: analytic formula for min M. C. ayne et al., Rev. Mod. hys. 64, 1045 (1992) 0: no analytic formula for performed numerically F, line minimization must be min
18 S Y WX UV Z ^ [ ] \ a c _` b T Computation of forces and stresses Field polarized Bloch functions are stationary points of F, the Hellmann- Feynmann theorem yields for the force on atom j f j df dr j F r j (no implicit dependence on r j via the wavefunctions) F F KS a 3 el ion no explicit dependence on r j f j F KS r j ez j Similar arguments show that the expression used to compute the stress under zero electric field remains valid
19 d e f g h gstate.f initberry.f scfcv.f Implementation in ABINIT (4.3.x) Initialization of Berry phase calculation: olarization, ddk & electric field store informations in the efield_type structured datatype berryphase_new.f Initialize polarization pel_cg(1:3) & Initialize overlap matrices Start SCF optimization of wavefunctions vtorho.f vtowfk.f cgwf.f Non-SCF line minimizations. For each iteration update overlap matrices add the electric field contribution to the gradient perform line minimization numerically (linemin.f) Update polarization: add change in Berry phase to pel_cg(1:3) nd SCF optimization of wavefunctions berryphase_new.f Recompute polarization, check that it is consistent with the value obtained from the update of the wavefuntions
20 i j k l m n o p q s u In practice (4.3.x) 1 st step: erform a calculation under zero electric field nd 2 step: lectric fiel calculation (use wavefunctions computed during the first step to initialize the calculation) berryopt = 4 efield(1:3) = cartesian coordinates of the electric field in atomic units nsym = 1 kptopt = 3 nband = number of valence bands COMMNTS: insulators only r MI parallelization not yet implemented no spin polarization (spin polarization probably available in ABINITv4.4.x) t You should increase the amplitude of the electric field in small steps Suggestion: use multiple datasets. For each dataset, use a slightly larger value of the electric field and take wavefunctions of the previous dataset to initialize the SCF cycle.
21 ~ y v x } w Applications: AlAs, electric field along x Force on Al along the x direction Born effective charges: Z f z { Finit electric field calculation: Z * = Linear response calculation: Z * =
22 Š ˆ Applications: AlAs, electric field along x Applications: AlAs, electric field along x Optical dielectric constant F ƒ Ž ŒŽ Ÿž œ š Ÿž œ š «Ÿž ª Ž
23 ± Applications: AlAs, electric field along (1,1,1) olarziation as a function of = *(1,1,1) Non-linear optical susceptibilities Finite electric field calculation: (2) = pm/v Non-linear response calculation: (2) = pm/v
24 I. Souza et al., RL 89, (2002) Application to cubic semiconductors
25 ²» Conclusions First-principles calculations of insulators in finite electric fields 1 st method: structural response to macroscopic electric fields ³ Derivatives of the energy with respect to the electric field can be accurately computed from DFT Low order Taylor expansion of the energy with respect to the electric field µ Map out the energy as a function of the polarization Iterative optimization of the structure under the constrained of a fixed polarization nd 2 method: first-principles approach to insulators in finite electric fields ractical scheme for computing the electronic structure of insulators under a finite bias ¹ For < c: minimization of a modified energy functional º Accurate computation of the total energy, the polarization, the forces and stresses Finite difference calculation of electric field derivatives
LA PRISE DE CALAIS. çoys, çoys, har - dis. çoys, dis. tons, mantz, tons, Gas. c est. à ce. C est à ce. coup, c est à ce
> ƒ? @ Z [ \ _ ' µ `. l 1 2 3 z Æ Ñ 6 = Ð l sl (~131 1606) rn % & +, l r s s, r 7 nr ss r r s s s, r s, r! " # $ s s ( ) r * s, / 0 s, r 4 r r 9;: < 10 r mnz, rz, r ns, 1 s ; j;k ns, q r s { } ~ l r mnz,
More informationFramework for functional tree simulation applied to 'golden delicious' apple trees
Purdue University Purdue e-pubs Open Access Theses Theses and Dissertations Spring 2015 Framework for functional tree simulation applied to 'golden delicious' apple trees Marek Fiser Purdue University
More informationOptimal Control of PDEs
Optimal Control of PDEs Suzanne Lenhart University of Tennessee, Knoville Department of Mathematics Lecture1 p.1/36 Outline 1. Idea of diffusion PDE 2. Motivating Eample 3. Big picture of optimal control
More informationAn Example file... log.txt
# ' ' Start of fie & %$ " 1 - : 5? ;., B - ( * * B - ( * * F I / 0. )- +, * ( ) 8 8 7 /. 6 )- +, 5 5 3 2( 7 7 +, 6 6 9( 3 5( ) 7-0 +, => - +< ( ) )- +, 7 / +, 5 9 (. 6 )- 0 * D>. C )- +, (A :, C 0 )- +,
More informationTheory of structural response to macroscopic electric fields in ferroelectric systems
PHYSICAL REVIEW B 66, 104108 2002 Theory of structural response to macroscopic electric fields in ferroelectric systems Na Sai, Karin M. Rabe, and David Vanderbilt Department of Physics and Astronomy,
More informationAn Introduction to Optimal Control Applied to Disease Models
An Introduction to Optimal Control Applied to Disease Models Suzanne Lenhart University of Tennessee, Knoxville Departments of Mathematics Lecture1 p.1/37 Example Number of cancer cells at time (exponential
More informationOC330C. Wiring Diagram. Recommended PKH- P35 / P50 GALH PKA- RP35 / RP50. Remarks (Drawing No.) No. Parts No. Parts Name Specifications
G G " # $ % & " ' ( ) $ * " # $ % & " ( + ) $ * " # C % " ' ( ) $ * C " # C % " ( + ) $ * C D ; E @ F @ 9 = H I J ; @ = : @ A > B ; : K 9 L 9 M N O D K P D N O Q P D R S > T ; U V > = : W X Y J > E ; Z
More information! " # $! % & '! , ) ( + - (. ) ( ) * + / 0 1 2 3 0 / 4 5 / 6 0 ; 8 7 < = 7 > 8 7 8 9 : Œ Š ž P P h ˆ Š ˆ Œ ˆ Š ˆ Ž Ž Ý Ü Ý Ü Ý Ž Ý ê ç è ± ¹ ¼ ¹ ä ± ¹ w ç ¹ è ¼ è Œ ¹ ± ¹ è ¹ è ä ç w ¹ ã ¼ ¹ ä ¹ ¼ ¹ ±
More informationHomogeneous Electric and Magnetic Fields in Periodic Systems
Electric and Magnetic Fields in Periodic Systems Josef W. idepartment of Chemistry and Institute for Research in Materials Dalhousie University Halifax, Nova Scotia June 2012 1/24 Acknowledgments NSERC,
More informationJ. Phys. Chem. Solids 61, 147 (2000). Berry-phase theory of proper piezoelectric response David Vanderbilt Department of Physics and Astronomy, Rutgers University, Piscataway, New Jersey 08855-0849 (July
More informationObtain from theory/experiment a value for the absorbtion cross-section. Calculate the number of atoms/ions/molecules giving rise to the line.
Line Emission Observations of spectral lines offer the possibility of determining a wealth of information about an astronomical source. A usual plan of attack is :- Identify the line. Obtain from theory/experiment
More informationOptimization of a parallel 3d-FFT with non-blocking collective operations
Optimization of a parallel 3d-FFT with non-blocking collective operations Chair of Computer Architecture Technical University of Chemnitz Département de Physique Théorique et Appliquée Commissariat à l
More informationA Functional Quantum Programming Language
A Functional Quantum Programming Language Thorsten Altenkirch University of Nottingham based on joint work with Jonathan Grattage and discussions with V.P. Belavkin supported by EPSRC grant GR/S30818/01
More informationGeometry Optimisation
Geometry Optimisation Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, UK http://www.cmt.york.ac.uk/cmd http://www.castep.org Motivation Overview of Talk Background
More information3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 REDCLIFF MUNICIPAL PLANNING COMMISSION FOR COMMENT/DISCUSSION DATE: TOPIC: April 27 th, 2018 Bylaw 1860/2018, proposed amendments to the Land Use Bylaw regarding cannabis
More informationproperties Michele Catti Dipartimento di Scienza dei Materiali Università di Milano Bicocca, Italy
Elastic and piezoelectric tensorial properties Michele Catti Dipartimento di Scienza dei Materiali Università di Milano Bicocca, Italy (catti@mater.unimib.it) 1 Tensorial physical properties of crystals
More informationChapter 134 : First-Principles Modeling of Ferroelectric Oxides Nanostructures
HANDBOOK OF THEORETICAL AND COMPUTATIONAL NANOTECHNOLOGY Chapter 134 : First-Principles Modeling of Ferroelectric Oxides Nanostructures by Philippe GHOSEZ (Université de Liège, Belgium), and Javier JUNQUERA
More informationMax. Input Power (W) Input Current (Arms) Dimming. Enclosure
Product Overview XI025100V036NM1M Input Voltage (Vac) Output Power (W) Output Voltage Range (V) Output urrent (A) Efficiency@ Max Load and 70 ase Max ase Temp. ( ) Input urrent (Arms) Max. Input Power
More informationSystematic treatment of displacements, strains, and electric fields in density-functional perturbation theory
Systematic treatment of displacements, strains, and electric fields in density-functional perturbation theory Xifan Wu, 1 David Vanderbilt, 1 and D. R. Hamann 2 1 Department of Physics and Astronomy, Rutgers
More informationFundamentals 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 informationElectric displacement as the fundamental variable in electronic-structure calculations
Electric displacement as the fundamental variable in electronic-structure calculations CECAM - Centre Européen de Calcul Atomique et Moléculaire EPF Lausanne, Switzerland Conference UC Davis, 6/23/2009
More informationNon-linear optics, k p perturbation theory, and the Sternheimer equation
Non-linear optics, k p perturbation theory, and the Sternheimer equation David A. Strubbe Department of Materials Science and Engineering Massachusetts Institute of Technology, Cambridge, MA Formerly Department
More informationPhonon Dispersion, Interatomic Force Constants Thermodynamic Quantities
Phonon Dispersion, Interatomic Force Constants Thermodynamic Quantities Umesh V. Waghmare Theoretical Sciences Unit J N C A S R Bangalore ICMR OUTLINE Vibrations and interatomic force constants (IFC) Extended
More informationALTER TABLE Employee ADD ( Mname VARCHAR2(20), Birthday DATE );
!! "# $ % & '( ) # * +, - # $ "# $ % & '( ) # *.! / 0 "# "1 "& # 2 3 & 4 4 "# $ % & '( ) # *!! "# $ % & # * 1 3 - "# 1 * #! ) & 3 / 5 6 7 8 9 ALTER ; ? @ A B C D E F G H I A = @ A J > K L ; ?
More informationPlanning for Reactive Behaviors in Hide and Seek
University of Pennsylvania ScholarlyCommons Center for Human Modeling and Simulation Department of Computer & Information Science May 1995 Planning for Reactive Behaviors in Hide and Seek Michael B. Moore
More informationQUESTIONS ON QUARKONIUM PRODUCTION IN NUCLEAR COLLISIONS
International Workshop Quarkonium Working Group QUESTIONS ON QUARKONIUM PRODUCTION IN NUCLEAR COLLISIONS ALBERTO POLLERI TU München and ECT* Trento CERN - November 2002 Outline What do we know for sure?
More informationVectors. Teaching Learning Point. Ç, where OP. l m n
Vectors 9 Teaching Learning Point l A quantity that has magnitude as well as direction is called is called a vector. l A directed line segment represents a vector and is denoted y AB Å or a Æ. l Position
More information" #$ P UTS W U X [ZY \ Z _ `a \ dfe ih j mlk n p q sr t u s q e ps s t x q s y i_z { U U z W } y ~ y x t i e l US T { d ƒ ƒ ƒ j s q e uˆ ps i ˆ p q y
" #$ +. 0. + 4 6 4 : + 4 ; 6 4 < = =@ = = =@ = =@ " #$ P UTS W U X [ZY \ Z _ `a \ dfe ih j mlk n p q sr t u s q e ps s t x q s y i_z { U U z W } y ~ y x t i e l US T { d ƒ ƒ ƒ j s q e uˆ ps i ˆ p q y h
More informationConnection equations with stream variables are generated in a model when using the # $ % () operator or the & ' %
7 9 9 7 The two basic variable types in a connector potential (or across) variable and flow (or through) variable are not sufficient to describe in a numerically sound way the bi-directional flow of matter
More informationETIKA V PROFESII PSYCHOLÓGA
P r a ž s k á v y s o k á š k o l a p s y c h o s o c i á l n í c h s t u d i í ETIKA V PROFESII PSYCHOLÓGA N a t á l i a S l o b o d n í k o v á v e d ú c i p r á c e : P h D r. M a r t i n S t r o u
More informationFirst-Principles Wannier Functions of Silicon and Gallium. Arsenide arxiv:cond-mat/ v1 [cond-mat.mtrl-sci] 22 Nov 1996.
First-Principles Wannier Functions of Silicon and Gallium Arsenide arxiv:cond-mat/9611176v1 [cond-mat.mtrl-sci] 22 Nov 1996 Pablo Fernández 1, Andrea Dal Corso 1, Francesco Mauri 2, and Alfonso Baldereschi
More informationA Beamforming Method for Blind Calibration of Time-Interleaved A/D Converters
A Beamforming Method for Blind Calibration of Time-nterleaved A/D Converters Bernard C. Levy University of California, Davis Joint wor with Steve Huang September 30, 200 Outline 1 Motivation Problem Formulation
More informationMatrices and Determinants
Matrices and Determinants Teaching-Learning Points A matri is an ordered rectanguar arra (arrangement) of numbers and encosed b capita bracket [ ]. These numbers are caed eements of the matri. Matri is
More informationVector analysis. 1 Scalars and vectors. Fields. Coordinate systems 1. 2 The operator The gradient, divergence, curl, and Laplacian...
Vector analysis Abstract These notes present some background material on vector analysis. Except for the material related to proving vector identities (including Einstein s summation convention and the
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 informationQuantum Molecular Dynamics Basics
Quantum Molecular Dynamics Basics Aiichiro Nakano Collaboratory for Advanced Computing & Simulations Depts. of Computer Science, Physics & Astronomy, Chemical Engineering & Materials Science, and Biological
More informationJuan Juan Salon. EH National Bank. Sandwich Shop Nail Design. OSKA Beverly. Chase Bank. Marina Rinaldi. Orogold. Mariposa.
( ) X é X é Q Ó / 8 ( ) Q / ( ) ( ) : ( ) : 44-3-8999 433 4 z 78-19 941, #115 Z 385-194 77-51 76-51 74-7777, 75-5 47-55 74-8141 74-5115 78-3344 73-3 14 81-4 86-784 78-33 551-888 j 48-4 61-35 z/ zz / 138
More informationPrincipal Secretary to Government Haryana, Town & Country Planning Department, Haryana, Chandigarh.
1 From To Principal Secretary to Government Haryana, Town & Country Planning Department, Haryana, Chandigarh. The Director General, Town & Country Planning Department, Haryana, Chandigarh. Memo No. Misc-2339
More informationGeneral Neoclassical Closure Theory: Diagonalizing the Drift Kinetic Operator
General Neoclassical Closure Theory: Diagonalizing the Drift Kinetic Operator E. D. Held eheld@cc.usu.edu Utah State University General Neoclassical Closure Theory:Diagonalizing the Drift Kinetic Operator
More informationAb initio study of the nonlinear optics of III-V semiconductors in the terahertz regime
Ab initio study of the nonlinear optics of III-V semiconductors in the terahertz regime Eric Roman, 1 Jonathan R. Yates, 1,2 Marek Veithen, 3 David Vanderbilt, 4 and Ivo Souza 1,2 1 Department of Physics,
More informationAbinit electron phonon interaction calculations for geniuses
Abinit electron phonon interaction calculations for geniuses Matthieu Verstraete, Bin Xu, and Momar Diakhate Unité Nanomat, Dept of Physics, Université de Liège, Belgium 1 I. INTRODUCTION This document
More informationTowards a High Level Quantum Programming Language
MGS Xmas 04 p.1/?? Towards a High Level Quantum Programming Language Thorsten Altenkirch University of Nottingham based on joint work with Jonathan Grattage and discussions with V.P. Belavkin supported
More informationarxiv: v1 [cond-mat.mtrl-sci] 22 Apr 2010
First-principles study of competing ferroelectric and antiferroelectric instabilities in BaTiO 3 /BaO superlattices Eric Bousquet 1,2, Javier Junquera 3, and Philippe Ghosez 1 1 Institut de physique (B5),
More informationGlasgow eprints Service
Kalna, K. and Asenov, A. (1) Multiple delta doping in aggressively scaled PHEMTs. In, Ryssel, H. and Wachutka, G. and Grunbacher, H., Eds. Solid-State Device Research Conference, 11-13 September 1, pages
More informationTowards a numerical solution of the "1/2 vs. 3/2" puzzle
Towards a numerical solution of the /2 s 3/2 puzzle Benoît Blossier DESY Zeuthen Motiations Hea Quark Effectie Theor Theoretical expectations Experimental measurements umerical computation of Conclusion
More informationThe University of Bath School of Management is one of the oldest established management schools in Britain. It enjoys an international reputation for
The University of Bath School of Management is one of the oldest established management schools in Britain. It enjoys an international reputation for the quality of its teaching and research. Its mission
More informationExamination paper for TFY4240 Electromagnetic theory
Department of Physics Examination paper for TFY4240 Electromagnetic theory Academic contact during examination: Associate Professor John Ove Fjærestad Phone: 97 94 00 36 Examination date: 16 December 2015
More informationBerry Phases and Curvatures in Electronic-Structure Theory. David Vanderbilt Rutgers University
Berry Phases and Curvatures in Electronic-Structure Theory David Vanderbilt Rutgers University Rahman prize for: Theory of polarization (King-Smith & Vanderbilt) Ultrasoft pseudopotentials Three quick
More informationStabilization of highly polar BiFeO 3 like structure: a new interface design route for enhanced ferroelectricity in artificial perovskite superlattice
Stabilization of highly polar BiFeO 3 like structure: a new interface design route for enhanced ferroelectricity in artificial perovskite superlattice Speaker: Xifan Wu Temple University Outline Part I.
More informationOptimizing Stresses for Testing DRAM Cell Defects Using Electrical Simulation
Optimizing Stresses for Testing RAM Cell efects Using Electrical Simulation Zaid Al-Ars Ad J. van de Goor Faculty of Information Technology and Systems Section of Computer Engineering elft University of
More informationB œ c " " ã B œ c 8 8. such that substituting these values for the B 3 's will make all the equations true
System of Linear Equations variables Ð unknowns Ñ B" ß B# ß ÞÞÞ ß B8 Æ Æ Æ + B + B ÞÞÞ + B œ, "" " "# # "8 8 " + B + B ÞÞÞ + B œ, #" " ## # #8 8 # ã + B + B ÞÞÞ + B œ, 3" " 3# # 38 8 3 ã + 7" B" + 7# B#
More informationAN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS
AN IDENTIFICATION ALGORITHM FOR ARMAX SYSTEMS First the X, then the AR, finally the MA Jan C. Willems, K.U. Leuven Workshop on Observation and Estimation Ben Gurion University, July 3, 2004 p./2 Joint
More informationComplex Analysis. PH 503 Course TM. Charudatt Kadolkar Indian Institute of Technology, Guwahati
Complex Analysis PH 503 Course TM Charudatt Kadolkar Indian Institute of Technology, Guwahati ii Copyright 2000 by Charudatt Kadolkar Preface Preface Head These notes were prepared during the lectures
More informationCurie-Weiss law in thin-film ferroelectrics
JOURNAL OF APPLIED PHYSICS 100, 044114 006 Curie-Weiss law in thin-film ferroelectrics Biao Wang School of Physics and Engineering, Sun Yat-sen University, Guangzhou, China C. H. Woo a Department of Electronic
More informationF O R SOCI AL WORK RESE ARCH
7 TH EUROPE AN CONFERENCE F O R SOCI AL WORK RESE ARCH C h a l l e n g e s i n s o c i a l w o r k r e s e a r c h c o n f l i c t s, b a r r i e r s a n d p o s s i b i l i t i e s i n r e l a t i o n
More informationFirst principles studies of the Born effective charges and electronic dielectric tensors for the relaxor PMN (PbMg 1/3 Nb 2/3 O 3 )
First principles studies of the Born effective charges and electronic dielectric tensors for the relaxor PMN (PbMg 1/3 Nb 2/3 O 3 ) Narayani Choudhury 1, R. E. Cohen 2 and Eric J. Walter 3 1 Solid State
More informationSelf Consistent Cycle
Self Consistent Cycle Step 0 : defining your system namelist SYSTEM How to specify the System All periodic systems can be specified by a Bravais Lattice and and atomic basis How to specify the Bravais
More informationMustafa Uludogan 1, Tahir Cagin, William A. Goddard, III Materials and Process Simulation Center, Caltech, Pasadena, CA 91125, U.S.A.
Ab Initio Studies On Phase Behavior of Barium Titanate Mustafa Uludogan 1, Tahir Cagin, William A. Goddard, III Materials and Process Simulation Center, Caltech, Pasadena, CA 91125, U.S.A. 1 Physics Department,
More informationTime-Dependent Density-Functional Theory
Summer School on First Principles Calculations for Condensed Matter and Nanoscience August 21 September 3, 2005 Santa Barbara, California Time-Dependent Density-Functional Theory X. Gonze, Université Catholique
More informationRedoing the Foundations of Decision Theory
Redoing the Foundations of Decision Theory Joe Halpern Cornell University Joint work with Larry Blume and David Easley Economics Cornell Redoing the Foundations of Decision Theory p. 1/21 Decision Making:
More informationThermodynamic theory of epitaxial ferroelectric thin films with dense domain structures
Thermodynamic theory of epitaxial ferroelectric thin films with dense domain structures V. G. Koukhar, N. A. ertsev,,, and R. Waser, A.F. Ioffe hysico-technical Institute, Russian Academy of ciences, 9
More informationLattice Vibrations. Chris J. Pickard. ω (cm -1 ) 200 W L Γ X W K K W
Lattice Vibrations Chris J. Pickard 500 400 300 ω (cm -1 ) 200 100 L K W X 0 W L Γ X W K The Breakdown of the Static Lattice Model The free electron model was refined by introducing a crystalline external
More informationFrom 180º stripe domains to more exotic patterns of polarization in ferroelectric nanostructures. A first principles view
From 180º stripe domains to more exotic patterns of polarization in ferroelectric nanostructures. A first principles view Pablo Aguado-Puente Javier Junquera Ferroelectricity: Basic definitions Existence
More informationTHE ANALYSIS OF THE CONVECTIVE-CONDUCTIVE HEAT TRANSFER IN THE BUILDING CONSTRUCTIONS. Zbynek Svoboda
THE NLSIS OF THE CONECTIE-CONDUCTIE HET TRNSFER IN THE BUILDING CONSTRUCTIONS Zbynek Soboda Czech Technical Uniersity in Prague Faculty of Ciil Engineering 166 29 Prague 6 - Czech Republic BSTRCT The numerical
More informationNumerical construction of Wannier functions
July 12, 2017 Internship Tutor: A. Levitt School Tutor: É. Cancès 1/27 Introduction Context: Describe electrical properties of crystals (insulator, conductor, semi-conductor). Applications in electronics
More informationVASP: Dielectric response Perturbation theory, linear response, and finite electric fields
VASP: Dielectric response Perturbation theory, linear response, and finite electric fields University of Vienna, Faculty of Physics and Center for Computational Materials Science, Vienna, Austria Outline
More informationDept 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 informationConstructive Decision Theory
Constructive Decision Theory Joe Halpern Cornell University Joint work with Larry Blume and David Easley Economics Cornell Constructive Decision Theory p. 1/2 Savage s Approach Savage s approach to decision
More informationContinuing Education and Workforce Training Course Schedule Winter / Spring 2010
Continuing Education and Workforce Training Course Schedule Winter / Spring 2010 www.southeastmn.edu/training Creating the Solution that Fits! Information Information Registration Registration via our
More informationElectronic Structure of Crystalline Solids
Electronic Structure of Crystalline Solids Computing the electronic structure of electrons in solid materials (insulators, conductors, semiconductors, superconductors) is in general a very difficult problem
More informationThe Fourier series are applicable to periodic signals. They were discovered by the
3.1 Fourier Series The Fourier series are applicable to periodic signals. They were discovered by the famous French mathematician Joseph Fourier in 1807. By using the Fourier series, a periodic signal
More informationCalculating GW corrections. with ABINIT
Calculating GW corrections S = igw with ABINIT Valerio Olevano (1), Rex Godby (2), Lucia Reining (1), Giovanni Onida (3), Marc Torrent (4) and Gian-Marco Rignanese (5) 1) Laboratoire des Solides Irradiés,
More informationSOLAR MORTEC INDUSTRIES.
SOLAR MORTEC INDUSTRIES www.mortecindustries.com.au Wind Regions Region A Callytharra Springs Gascoyne Junction Green Head Kununurra Lord Howe Island Morawa Toowoomba Wittanoom Bourke Region B Adelaide
More informationAutomating the Finite Element Method
Automating the Finite Element Method Anders Logg logg@tti-c.org Toyota Technological Institute at Chicago Johan Hoffman, Johan Jansson, Claes Johnson, Robert C. Kirby, Matthew Knepley, Ridgway Scott Simula
More informationDamping Ring Requirements for 3 TeV CLIC
Damping Ring Requirements for 3 TeV CLIC Quantity Symbol Value Bunch population N b 9 4.1 10 No. of bunches/train k bt 154 Repetition frequency f r 100 Hz Horizontal emittance γε x 7 4.5 10 m Vertical
More information4.3 Laplace Transform in Linear System Analysis
4.3 Laplace Transform in Linear System Analysis The main goal in analysis of any dynamic system is to find its response to a given input. The system response in general has two components: zero-state response
More informationBand Structure Calculations; Electronic and Optical Properties
; Electronic and Optical Properties Stewart Clark University of Outline Introduction to band structures Calculating band structures using Castep Calculating optical properties Examples results Some applications
More informationTable of Contents. Table of Contents Computing the piezoelectric tensor for AlN. Introduction. Computing the piezoelectric tensor
Table of Contents Table of Contents Computing the piezoelectric tensor for AlN Introduction The piezoelectric tensor Computing the piezoelectric tensor Analysis of the e 33 e33 coefficient Alternative
More informationSome emission processes are intrinsically polarised e.g. synchrotron radiation.
Polarisation Some emission processes are intrinsically polarised e.g. synchrotron radiation. B e linearly polarised emission circularly polarised emission Scattering processes can either increase or decrease
More informationPlasma diagnostics and abundance determinations for planetary nebulae current status. Xiaowei Liu Department of Astronomy, Peking University
3 ' 4 2 - '. %! A?@ @ 3 +.?; = ; %,9 ' ü / :, %4 ' 9 8 +! (76 5 +. *". 34 %, /, /10 *$+,,' -+!. *$+ () ', ' & "$# Plasma diagnostics and abundance determinations for planetary nebulae current status
More information1. Allstate Group Critical Illness Claim Form: When filing Critical Illness claim, please be sure to include the following:
Dear Policyholder, Per your recent request, enclosed you will find the following forms: 1. Allstate Group Critical Illness Claim Form: When filing Critical Illness claim, please be sure to include the
More informationA General Procedure to Design Good Codes at a Target BER
A General Procedure to Design Good odes at a Target BER Speaker: Xiao Ma 1 maxiao@mail.sysu.edu.cn Joint work with: hulong Liang 1, Qiutao Zhuang 1, and Baoming Bai 2 1 Dept. Electronics and omm. Eng.,
More informationNew BaBar Results on Rare Leptonic B Decays
New BaBar Results on Rare Leptonic B Decays Valerie Halyo Stanford Linear Accelerator Center (SLAC) FPCP 22 valerieh@slac.stanford.edu 1 Table of Contents Motivation to look for and Analysis Strategy Measurement
More informationPharmacological and genomic profiling identifies NF-κB targeted treatment strategies for mantle cell lymphoma
CORRECTION NOTICE Nat. Med. 0, 87 9 (014) Pharmacoogica and genomic profiing identifies NF-κB targeted treatment strategies for mante ce ymphoma Rami Raha, Mareie Fric, Rodrigo Romero, Joshua M Korn, Robert
More informationRarefied Gas FlowThroughan Orifice at Finite Pressure Ratio
Rarefied Gas FlowThroughan Orifice at Finite Pressure Ratio Felix Sharipov Departamento de Física, Universidade Federal do Paraná Caixa Postal 90, 85-990 Curitiba, Brazil Email: sharipov@fisica.ufpr.br;
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 informationALMA: All-scale predictive design of heat management material structures
ALMA: All-scale predictive design of heat management material structures Version Date: 2015.11.13. Last updated 2015.12.02 Purpose of this document: Definition of a data organisation that is applicable
More informationElectrons in a periodic potential
Chapter 3 Electrons in a periodic potential 3.1 Bloch s theorem. We consider in this chapter electrons under the influence of a static, periodic potential V (x), i.e. such that it fulfills V (x) = V (x
More information236 Chapter 4 Applications of Derivatives
26 Chapter Applications of Derivatives Î$ &Î$ Î$ 5 Î$ 0 "Î$ 5( 2) $È 26. (a) g() œ ( 5) œ 5 Ê g () œ œ Ê critical points at œ 2 and œ 0 Ê g œ ± )(, increasing on ( _ß 2) and (!ß _), decreasing on ( 2 ß!)!
More informationLecture 4: Basic elements of band theory
Phys 769 Selected Topics in Condensed Matter Physics Summer 010 Lecture 4: Basic elements of band theory Lecturer: Anthony J. Leggett TA: Bill Coish 1 Introduction Most matter, in particular most insulating
More informationComputational Methods. Chem 561
Computational Methods Chem 561 Lecture Outline 1. Ab initio methods a) HF SCF b) Post-HF methods 2. Density Functional Theory 3. Semiempirical methods 4. Molecular Mechanics Computational Chemistry " Computational
More informationMARKET AREA UPDATE Report as of: 1Q 2Q 3Q 4Q
Year: 2013 Market Area (Cit, State): Greater Portland, OR MARKET AREA UPDATE Report as of: 1Q 2Q 3Q 4Q Provided b (Compan / Companies): Coldwell Banker Seal, Archibald Relocation, Prudential Northwest
More informationDensity Functional Theory
Density Functional Theory Iain Bethune EPCC ibethune@epcc.ed.ac.uk Overview Background Classical Atomistic Simulation Essential Quantum Mechanics DFT: Approximations and Theory DFT: Implementation using
More informationMott metal-insulator transition on compressible lattices
Mott metal-insulator transition on compressible lattices Markus Garst Universität zu Köln T p in collaboration with : Mario Zacharias (Köln) Lorenz Bartosch (Frankfurt) T c Mott insulator p c T metal pressure
More informationSpontaneous polarization of ferroelectric BaTiO3
Spontaneous polarization of ferroelectric BaTiO3 Tutorial on how to use ATK for polarization calculations Version 2014.0 Spontaneous polarization of ferroelectric BaTiO3: Tutorial on how to use ATK for
More informationUNIQUE FJORDS AND THE ROYAL CAPITALS UNIQUE FJORDS & THE NORTH CAPE & UNIQUE NORTHERN CAPITALS
Q J j,. Y j, q.. Q J & j,. & x x. Q x q. ø. 2019 :. q - j Q J & 11 Y j,.. j,, q j q. : 10 x. 3 x - 1..,,. 1-10 ( ). / 2-10. : 02-06.19-12.06.19 23.06.19-03.07.19 30.06.19-10.07.19 07.07.19-17.07.19 14.07.19-24.07.19
More informationSymbols and dingbats. A 41 Α a 61 α À K cb ➋ à esc. Á g e7 á esc. Â e e5 â. Ã L cc ➌ ã esc ~ Ä esc : ä esc : Å esc * å esc *
Note: Although every effort ws tken to get complete nd ccurte tble, the uhtor cn not be held responsible for ny errors. Vrious sources hd to be consulted nd MIF hd to be exmined to get s much informtion
More informationLecture 4: Hartree-Fock Theory
Lecture 4: Hartree-Fock Theory One determinant to rule them all, One determinant to find them, One determinant to bring them all and in the darkness bind them Second quantization rehearsal The formalism
More informationOptimal data-independent point locations for RBF interpolation
Optimal dataindependent point locations for RF interpolation S De Marchi, R Schaback and H Wendland Università di Verona (Italy), Universität Göttingen (Germany) Metodi di Approssimazione: lezione dell
More informationFront-end. Organization of a Modern Compiler. Middle1. Middle2. Back-end. converted to control flow) Representation
register allocation instruction selection Code Low-level intermediate Representation Back-end Assembly array references converted into low level operations, loops converted to control flow Middle2 Low-level
More information