On the adaptive finite element analysis of the Kohn-Sham equations
|
|
- Edwin Harvey
- 5 years ago
- Views:
Transcription
1 On the adaptive finite element analysis of the Kohn-Sham equations Denis Davydov, Toby Young, Paul Steinmann Denis Davydov, LTM, Erlangen, Germany August 2015
2 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 2 /29 Outline 1. Density Functional Theory and its discretisation with FEM 2. Mesh motion approach 3. A posteriori mesh refinement 4. Numerical examples and discussion 5. Conclusions and future work
3 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 3 /29 Different Scales time scale d deal.ii (a finite element Differential Equations Analysis Library) atomistic simulation quantum physics continuum mechanics Our goal is to solve quantum mechanics using h-fem. length scale Bangerth, W., Hartmann, R., & Kanschat, G. (2007). deal.ii---a general-purpose object-oriented finite element library. ACM Transactions on Mathematical Software, 33(4), 24 es. doi: /
4 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 4 /29 Kohn-Sham theorem 1. The ground state property of many-electron system is uniquely determined by an electron density (x) x 2 R 3 2. There exists an energy functional of the system. In the ground state electron density minimises this functional. Born-Oppenheimer approximation: nuclei are fixed in space, electrons are moving in a static external potential generated by nuclei. O R J R I Electronic state is described by a wavefunction (x)
5 Density Functional Theory Total energy of the system E 0 = E kin + E Hartree + E ion + E xc + E zz The electron density: (x) := X f (x) 2 The kinetic energy of non-interacting electrons: E kin := X partial occupancy orthonormal electronic Z apple 1 f number, X such that wavefunctions 2 r2 dx f N e The exchange-correlation Z energy which represents quantum many-body interactions: E xc := (x)" xc ( (x)) dx The Hartree energy E Hartree := 1 Z (electrostatic Z interaction between electrons): (x) (x 0 ) 2 x x 0 dx 0 dx. Electrostatic interaction between electrons Z Z " and the external # potential of nuclei: X Z I E ion := V ion (x) (x)dx =: (x)dx x R I Repulsive nuclei-nuclei E zz := 1 X electrostatic interaction energy: Z I Z J 2 R I R J I,J6=I I Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 5 /29
6 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 6 /29 Kohn-Sham eigenvalue problem stationary conditions " # 1 2 r2 + V e (x; ) (x) = (x) V e (x) :=V ion (x)+v Hartree (x)+v xc (x) X I Z I x R I + Z (x 0 ) x x 0 dx0 + E xc Explicit Approach: compute Hartree potential more efficiently via solving the Poisson equation r 2 V Hartree (x) =4 (x) V ion need to evaluate for every quadrature point. Hence, bad scaling w.r.t. number of atoms. Implicit Approach: recall that 1/r is Green s functions of the Laplace operator " # X r 2 [V Hartree (x)+v ion (x)] = 4 (x) Z I (x R I ) solution is not in H 1! I
7 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 7 /29 Kohn-Sham eigenvalue problem (cont.) stationary conditions " # 1 2 r2 + V e (x; ) (x) = (x) V e (x) :=V S ion(x)+ V L ion(x)+v Hartree (x) + V xc (x) As a result evaluation of V ion scales linearly w.r.t. the number of atoms. Gaussian Approach: split Coulomb potential into fully-local short-range and smooth long-range parts Vion(x) L erf ( x R I / ) := Z I x R I X I This potential corresponds to the Gaussian charge distribution, i.e. r 2 Vion(x) L 4 X Z I x exp RI 2 3/2 3 2 I The remaining short-range part is used directly during assembly of the eigenvalue problem Vion(x) S Z I := V L 1 erf ( x R I / ) x R I ion(x) = Z I x R I X I X I whereas the Hartree potential and the long-range parts are evaluated by solving the Poisson problem r 2 V Hartree (x)+vion(x) L " X Z I x =4 (x) exp RI 2 # 3/2 3 2 I
8 Discretisation in space Introduce a finite (x) = X element basis for the wave-functions and the potential fields in i (x) '(x) = X i i ' i N ' i (x) (x) := X f (x) 2 The generalised eigenvalue problem [K ij + V ij ] j = M ij j solve iteratively until convergence i.e. split into a sequence of linear problems The generalised Poisson problem L ij ' j = R H j + R V j the potential Zmatrix V ij := N i (x)v e (x; ', )N j (x)dx contribution to the RHS from the electron density Z Rj H := 4 N ' i (x) (x)dx the mass (or overlap) matrix Z M ij := N i (x) N j (x)dx the kinetic matrix K ij := 1 Z 2 rn i (x) rn j (x)dx possibly non-zero contribution to the RHS from nuclei density - R V j the Laplace matrix Z L ij := rn ' i (x) rn ' j (x)dx Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 8 /29
9 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 9 /29 Outline 1. Density Functional Theory and its discretisation with FEM 2. Mesh motion approach 3. A posteriori mesh refinement 4. Numerical examples and discussion 5. Conclusions and future work
10 Mesh motion approach In order to increase the accuracy of numerical solution and avoid singularities of the Coulomb potential, we need to obtain a mesh such that ions are located at vertexes. O V I R I I u(v I ) - find the closest vertex - prescribe it s motion to nucleus position - move other nodes as to Hence we face a problem of determining the mesh motion field u(x) according to the specified point-wise constraints. Approach 1: Laplace equation [c(x)ru(x)] r = 0 on, u(x) =0 such that u(v I )=R I V I, initial position of the vertex closest to nuclei I. c(x) = 1 min I x R I. Approach 1I: small-strain linear elasticity [c(x)c : r sym u(x)] r = 0 on, u(x) =0 such that u(v I )=R I V I. fourth order isotropic elastic tensor with unit bulk and shear moduli. Jasak, H. & Tuković, Z (2006). Automatic Mesh Motion for the Unstructured Finite Volume Method Transactions of Famena, 30, Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 10/29
11 Mesh motion (scaled Jacobian element quality) Naive displacement of the closest vertex to the position of nuclei lead to a unacceptably destroyed elements. Laplace approach was found to lead to a better element quality near nuclei. Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 11/29
12 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 12/29 Outline 1. Density Functional Theory and its discretisation with FEM 2. Mesh motion approach 3. A posteriori mesh refinement 4. Numerical examples and discussion 5. Conclusions and future work
13 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 13/29 A posteriori mesh refinement Approach 1: Residual based error-estimator for each eigenpair: Z apple 2 Z 1 e, 2 := h 2 e 2 r2 + V e dx + h e [r n] 2 e e 2 e X 2 e, Approach 1I: Kelly error estimator applied to electron density or the electrostatic potential '. 2 e( ) :=h e Z@ e [[r( ) n] da Marking strategy: marks a minimum subset of a triangulation M T squared error is more than a given fraction of the squared total error " # X apple X e 2 e2t e2m 2 e whose
14 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 14/29 Outline 1. Density Functional Theory and its discretisation with FEM 2. Mesh motion approach 3. A posteriori mesh refinement 4. Numerical examples and discussion 5. Conclusions and future work
15 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 15/29 Treatment of the ionic potential Consider hydrogen atom centered at the origin. V e V ion H V ion (x) V L ion(x) The 1/r singularity can not be represented exactly by non-enhanced FE trial spaces. For the Gaussian charge approach a convergence to the analytical solution is evident.
16 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 16/29 Treatment of the ionic potential (cont.) Consider hydrogen and helium atoms that are centered at the origin. H He Each approach converges to the expected solution with subsequent mesh refinement. The split of the potential into fully local short-range and smooth long-range part is beneficial from the computational points of view. Otherwise scales badly w.r.t. the number of atoms.
17 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 17/29 Quadrature order For explicit treatment of ionic potential nonlinearities mainly come from the exchange-correlation potential and the Coulomb (1/r) potential. Consider helium and carbon atoms with different quadrature orders used to integrate Hamiltonian matrix and RHS of the Laplace problem. He C For both atoms energy convergence w.r.t. quadrature order is observed
18 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 18/29 Quadrature order (cont.) He C A posteriori refinement has a bigger influence on the resulting energy in comparison to the quadrature order used in calculations. Insufficient quadrature rule may result in energies lower than those obtained form the radial solution (i.e. violation of the variational principle).
19 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 19/29 Choice of the polynomial degree in FE spaces 1) which FE basis is the best (DoFs per element vs interpolation properties)? 2) can we use the same FE basis for Poisson and eigenvalue problem? O Linear FEs are too costly. The FE basis in the Poisson problem should be twice the polynomial degree of the eigenvalue problem in order to maintain variational character of the solution (the difference between energies becomes negative for 2/2 combination). The efficiency of the basis to reach the chemical accuracy (1e-3) depends on the problem.
20 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 20/29 Error estimators Consider neon (Ne) atom. Ne
21 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 21/29 Error estimators (cont.) Compare different error estimators: Ne Residual error estimator and Kelly estimator applied to eigenvectors result in a good convergence w.r.t. the expected total energy. The error estimator based on jumps in density (Kelly) does not lead to convergence. The Kelly estimator applied to the solution of the Poisson problem shows suboptimal convergence.
22 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 22/29 Towards optimisation of nucleus positions Consider a hydrogen molecule ion. Atoms are placed on the x axis symmetrically about the origin with varying interatomic distance. The same structured mesh is used as an input. H H
23 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 23/29 Towards optimisation of nucleus positions Consider carbon monoxide. Atoms are placed on the x axis symmetrically about the origin with varying interatomic distance. The same structured mesh is used as an input. O C a triple bond (+) :O C : ( ) 7 (x) 2
24 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 24/29 Outline 1. Density Functional Theory and its discretisation with FEM 2. Mesh motion approach 3. A posteriori mesh refinement 4. Numerical examples and discussion 5. Conclusions and future work
25 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 25/29 Summary Conclusions: Solution of the quantum many-body problem within the context of DFT using a realspace method based on h-adaptive FE discretisation of the space has been developed Different treatments of the ionic potential produce adequate results given fine enough meshes and are well behaved with respect to refinement. The quadrature rule is not the most influential factor so long as it is accurate enough to integrate the mass matrix. In order to maintain variational nature of the solution, the polynomial degree of the Poisson FE basis should be twice of that used for the eigenvalue problem. The effectiveness of the polynomial degree in the FE basis depends on the problem considered. A residual-based error estimator or Kelly-error estimator applied to eigenvectors result in adequate convergence. Error estimators based on the density should be avoided. The proposed mesh motion approach to move vertices in such a way that each atom has an associated vertex at its position is well suited to be applied to the structural optimisation of molecules.
26 Denis Davydov, LTM, Erlangen, Germany College Station (August 2015) 26/29 Thank you for your attention References: Davydov, D.; Young, T. & Steinmann, P. (2015) On the adaptive finite element analysis of the Kohn-Sham equations: Methods, algorithms, and implementation. Journal for Numerical Methods in Engineering. Advanced Grant MOCOPOLY
Solid State Theory: Band Structure Methods
Solid State Theory: Band Structure Methods Lilia Boeri Wed., 11:15-12:45 HS P3 (PH02112) http://itp.tugraz.at/lv/boeri/ele/ Plan of the Lecture: DFT1+2: Hohenberg-Kohn Theorem and Kohn and Sham equations.
More informationCHEM6085: Density Functional Theory
Lecture 5 CHEM6085: Density Functional Theory Orbital-free (or pure ) DFT C.-K. Skylaris 1 Consists of three terms The electronic Hamiltonian operator Electronic kinetic energy operator Electron-Electron
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 informationDensity Functional Theory
Density Functional Theory March 26, 2009 ? DENSITY FUNCTIONAL THEORY is a method to successfully describe the behavior of atomic and molecular systems and is used for instance for: structural prediction
More informationTeoría del Funcional de la Densidad (Density Functional Theory)
Teoría del Funcional de la Densidad (Density Functional Theory) Motivation: limitations of the standard approach based on the wave function. The electronic density n(r) as the key variable: Functionals
More informationElectron Correlation
Electron Correlation Levels of QM Theory HΨ=EΨ Born-Oppenheimer approximation Nuclear equation: H n Ψ n =E n Ψ n Electronic equation: H e Ψ e =E e Ψ e Single determinant SCF Semi-empirical methods Correlation
More informationThe Plane Stress Problem
The Plane Stress Problem Martin Kronbichler Applied Scientific Computing (Tillämpad beräkningsvetenskap) February 2, 2010 Martin Kronbichler (TDB) The Plane Stress Problem February 2, 2010 1 / 24 Outline
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 informationDensity Functional Theory. Martin Lüders Daresbury Laboratory
Density Functional Theory Martin Lüders Daresbury Laboratory Ab initio Calculations Hamiltonian: (without external fields, non-relativistic) impossible to solve exactly!! Electrons Nuclei Electron-Nuclei
More informationQuantum Mechanical Simulations
Quantum Mechanical Simulations Prof. Yan Wang Woodruff School of Mechanical Engineering Georgia Institute of Technology Atlanta, GA 30332, U.S.A. yan.wang@me.gatech.edu Topics Quantum Monte Carlo Hartree-Fock
More informationOVERVIEW OF QUANTUM CHEMISTRY METHODS
OVERVIEW OF QUANTUM CHEMISTRY METHODS Outline I Generalities Correlation, basis sets Spin II Wavefunction methods Hartree-Fock Configuration interaction Coupled cluster Perturbative methods III Density
More informationSolving Many-Body Schrödinger Equation Using Density Functional Theory and Finite Elements
Solving Many-Body Schrödinger Equation Using Density Functional Theory and Finite Elements Institute of Physics, Academy of Sciences of the Czech Republic June 21, 2008 Introduction Contens Density Functional
More informationGEM4 Summer School OpenCourseWare
GEM4 Summer School OpenCourseWare http://gem4.educommons.net/ http://www.gem4.org/ Lecture: Molecular Mechanics by Ju Li. Given August 9, 2006 during the GEM4 session at MIT in Cambridge, MA. Please use
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 informationScientific Computing I
Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Neckel Winter 2013/2014 Module 8: An Introduction to Finite Element Methods, Winter 2013/2014 1 Part I: Introduction to
More informationElectrochemistry project, Chemistry Department, November Ab-initio Molecular Dynamics Simulation
Electrochemistry project, Chemistry Department, November 2006 Ab-initio Molecular Dynamics Simulation Outline Introduction Ab-initio concepts Total energy concepts Adsorption energy calculation Project
More informationLehrstuhl Informatik V. Lehrstuhl Informatik V. 1. solve weak form of PDE to reduce regularity properties. Lehrstuhl Informatik V
Part I: Introduction to Finite Element Methods Scientific Computing I Module 8: An Introduction to Finite Element Methods Tobias Necel Winter 4/5 The Model Problem FEM Main Ingredients Wea Forms and Wea
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 informationFinite Element Method in Geotechnical Engineering
Finite Element Method in Geotechnical Engineering Short Course on + Dynamics Boulder, Colorado January 5-8, 2004 Stein Sture Professor of Civil Engineering University of Colorado at Boulder Contents Steps
More informationCHEMISTRY Topic #1: Bonding What Holds Atoms Together? Spring 2012 Dr. Susan Lait
CHEMISTRY 2000 Topic #1: Bonding What Holds Atoms Together? Spring 2012 Dr. Susan Lait Why Do Bonds Form? An energy diagram shows that a bond forms between two atoms if the overall energy of the system
More informationLarge-scale real-space electronic structure calculations
Large-scale real-space electronic structure calculations YIP: Quasi-continuum reduction of field theories: A route to seamlessly bridge quantum and atomistic length-scales with continuum Grant no: FA9550-13-1-0113
More informationMeasurement of deformation. Measurement of elastic force. Constitutive law. Finite element method
Deformable Bodies Deformation x p(x) Given a rest shape x and its deformed configuration p(x), how large is the internal restoring force f(p)? To answer this question, we need a way to measure deformation
More informationElectronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory
Electronic structure theory: Fundamentals to frontiers. 1. Hartree-Fock theory MARTIN HEAD-GORDON, Department of Chemistry, University of California, and Chemical Sciences Division, Lawrence Berkeley National
More informationQuantum Continuum Mechanics for Many-Body Systems
Workshop on High Performance and Parallel Computing Methods and Algorithms for Materials Defects, 9-13 February 2015 Quantum Continuum Mechanics for Many-Body Systems J. Tao 1,2, X. Gao 1,3, G. Vignale
More informationIndex. C 2 ( ), 447 C k [a,b], 37 C0 ( ), 618 ( ), 447 CD 2 CN 2
Index advection equation, 29 in three dimensions, 446 advection-diffusion equation, 31 aluminum, 200 angle between two vectors, 58 area integral, 439 automatic step control, 119 back substitution, 604
More informationCHEM6085: Density Functional Theory Lecture 10
CHEM6085: Density Functional Theory Lecture 10 1) Spin-polarised calculations 2) Geometry optimisation C.-K. Skylaris 1 Unpaired electrons So far we have developed Kohn-Sham DFT for the case of paired
More informationIntroduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić
Introduction to Density Functional Theory with Applications to Graphene Branislav K. Nikolić Department of Physics and Astronomy, University of Delaware, Newark, DE 19716, U.S.A. http://wiki.physics.udel.edu/phys824
More informationMolecular Bonding. Molecular Schrödinger equation. r - nuclei s - electrons. M j = mass of j th nucleus m 0 = mass of electron
Molecular onding Molecular Schrödinger equation r - nuclei s - electrons 1 1 W V r s j i j1 M j m i1 M j = mass of j th nucleus m = mass of electron j i Laplace operator for nuclei Laplace operator for
More informationKarhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques
Institut für Numerische Mathematik und Optimierung Karhunen-Loève Approximation of Random Fields Using Hierarchical Matrix Techniques Oliver Ernst Computational Methods with Applications Harrachov, CR,
More informationVALENCE Hilary Term 2018
VALENCE Hilary Term 2018 8 Lectures Prof M. Brouard Valence is the theory of the chemical bond Outline plan 1. The Born-Oppenheimer approximation 2. Bonding in H + 2 the LCAO approximation 3. Many electron
More informationIntroduction to Computational Chemistry
Introduction to Computational Chemistry Vesa Hänninen Laboratory of Physical Chemistry room B430, Chemicum 4th floor vesa.hanninen@helsinki.fi September 3, 2013 Introduction and theoretical backround September
More informationMolecular Mechanics: The Ab Initio Foundation
Molecular Mechanics: The Ab Initio Foundation Ju Li GEM4 Summer School 2006 Cell and Molecular Mechanics in BioMedicine August 7 18, 2006, MIT, Cambridge, MA, USA 2 Outline Why are electrons quantum? Born-Oppenheimer
More informationDFT calculations of NMR indirect spin spin coupling constants
DFT calculations of NMR indirect spin spin coupling constants Dalton program system Program capabilities Density functional theory Kohn Sham theory LDA, GGA and hybrid theories Indirect NMR spin spin coupling
More informationLecture 9 Approximations of Laplace s Equation, Finite Element Method. Mathématiques appliquées (MATH0504-1) B. Dewals, C.
Lecture 9 Approximations of Laplace s Equation, Finite Element Method Mathématiques appliquées (MATH54-1) B. Dewals, C. Geuzaine V1.2 23/11/218 1 Learning objectives of this lecture Apply the finite difference
More informationModule 6 1. Density functional theory
Module 6 1. Density functional theory Updated May 12, 2016 B A DDFT C K A bird s-eye view of density-functional theory Authors: Klaus Capelle G http://arxiv.org/abs/cond-mat/0211443 R https://trac.cc.jyu.fi/projects/toolbox/wiki/dft
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 informationBehind the "exciting" curtain: The (L)APW+lo method
Behind the "exciting" curtain: The (L)APW+lo method Aug 7, 2016 Andris Gulans Humboldt-Universität zu Berlin Kohn-Sham equation Potential due to nuclei Exchange-correlation potential Potential due to electron
More informationInverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros
Inverse Design (and a lightweight introduction to the Finite Element Method) Stelian Coros Computational Design Forward design: direct manipulation of design parameters Level of abstraction Exploration
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University Variational Problems of the Dirichlet BVP of the Poisson Equation 1 For the homogeneous
More informationThe Plane-Wave Pseudopotential Method
Hands-on Workshop on Density Functional Theory and Beyond: Computational Materials Science for Real Materials Trieste, August 6-15, 2013 The Plane-Wave Pseudopotential Method Ralph Gebauer ICTP, Trieste
More informationDept of Mechanical Engineering MIT Nanoengineering group
1 Dept of Mechanical Engineering MIT Nanoengineering group » To calculate all the properties of a molecule or crystalline system knowing its atomic information: Atomic species Their coordinates The Symmetry
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 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 informationFVM for Fluid-Structure Interaction with Large Structural Displacements
FVM for Fluid-Structure Interaction with Large Structural Displacements Željko Tuković and Hrvoje Jasak Zeljko.Tukovic@fsb.hr, h.jasak@wikki.co.uk Faculty of Mechanical Engineering and Naval Architecture
More informationOslo node. Highly accurate calculations benchmarking and extrapolations
Oslo node Highly accurate calculations benchmarking and extrapolations Torgeir Ruden, with A. Halkier, P. Jørgensen, J. Olsen, W. Klopper, J. Gauss, P. Taylor Explicitly correlated methods Pål Dahle, collaboration
More informationAdvanced Quantum Mechanics
Advanced Quantum Mechanics University of York Lecturer: Rex Godby Notes by Victor Naden Robinson Lecture 1: TDSE Lecture 2: TDSE Lecture 3: FMG Lecture 4: FMG Lecture 5: Ehrenfest s Theorem and the Classical
More informationChapter 1: The Finite Element Method
Chapter 1: The Finite Element Method Michael Hanke Read: Strang, p 428 436 A Model Problem Mathematical Models, Analysis and Simulation, Part Applications: u = fx), < x < 1 u) = u1) = D) axial deformation
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 5 The Hartree-Fock method C.-K. Skylaris Learning outcomes Be able to use the variational principle in quantum calculations Be able to construct Fock operators
More information3: Density Functional Theory
The Nuts and Bolts of First-Principles Simulation 3: Density Functional Theory CASTEP Developers Group with support from the ESF ψ k Network Density functional theory Mike Gillan, University College London
More informationIntroduction to DFTB. Marcus Elstner. July 28, 2006
Introduction to DFTB Marcus Elstner July 28, 2006 I. Non-selfconsistent solution of the KS equations DFT can treat up to 100 atoms in routine applications, sometimes even more and about several ps in MD
More informationPoisson Solver, Pseudopotentials, Atomic Forces in the BigDFT code
CECAM Tutorial on Wavelets in DFT, CECAM - LYON,, in the BigDFT code Kernel Luigi Genovese L_Sim - CEA Grenoble 28 November 2007 Outline, Kernel 1 The with Interpolating Scaling Functions in DFT for Interpolating
More informationMO Calculation for a Diatomic Molecule. /4 0 ) i=1 j>i (1/r ij )
MO Calculation for a Diatomic Molecule Introduction The properties of any molecular system can in principle be found by looking at the solutions to the corresponding time independent Schrodinger equation
More informationInstructor background for the discussion points of Section 2
Supplementary Information for: Orbitals Some fiction and some facts Jochen Autschbach Department of Chemistry State University of New York at Buffalo Buffalo, NY 14260 3000, USA Instructor background for
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 4 Molecular orbitals C.-K. Skylaris Learning outcomes Be able to manipulate expressions involving spin orbitals and molecular orbitals Be able to write down
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 informationDensity Functional Theory - II part
Density Functional Theory - II part antonino.polimeno@unipd.it Overview From theory to practice Implementation Functionals Local functionals Gradient Others From theory to practice From now on, if not
More informationVibration of Thin Beams by PIM and RPIM methods. *B. Kanber¹, and O. M. Tufik 1
APCOM & ISCM -4 th December, 23, Singapore Vibration of Thin Beams by PIM and RPIM methods *B. Kanber¹, and O. M. Tufik Mechanical Engineering Department, University of Gaziantep, Turkey. *Corresponding
More informationContinuum mechanics of beam-like structures using one-dimensional finite element based on Serendipity Lagrange cross-sectional discretisation, Mayank Patni, Prof. Paul Weaver, Dr Alberto Pirrera Bristol
More informationGeneralized Finite Element Methods for Three Dimensional Structural Mechanics Problems. C. A. Duarte. I. Babuška and J. T. Oden
Generalized Finite Element Methods for Three Dimensional Structural Mechanics Problems C. A. Duarte COMCO, Inc., 7800 Shoal Creek Blvd. Suite 290E Austin, Texas, 78757, USA I. Babuška and J. T. Oden TICAM,
More informationHandbook of Computational Quantum Chemistry. DAVID B. COOK The Department of Chemistry, University of Sheffield
Handbook of Computational Quantum Chemistry DAVID B. COOK The Department of Chemistry, University of Sheffield Oxford New York Tokyo OXFORD UNIVERSITY PRESS 1998 CONTENTS 1 Mechanics and molecules 1 1.1
More informationIntroduction to Computational Materials Science
Introduction to Computational Materials Science Fundamentals to Applications RICHARD LESAR lowa State University.CAMBRIDGE ::: UNIVERSITY PRESS CONTENTS Preface 1 Introduction to materials modeling and
More informationBAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS
Journal of Computational and Applied Mechanics, Vol.., No. 1., (2005), pp. 83 94 BAR ELEMENT WITH VARIATION OF CROSS-SECTION FOR GEOMETRIC NON-LINEAR ANALYSIS Vladimír Kutiš and Justín Murín Department
More informationEstimating Convergence Probability for the Hartree-Fock Algorithm for Calculating Electronic Wavefunctions for Molecules
Estimating Convergence Probability for the Hartree-Fock Algorithm for Calculating Electronic Wavefunctions for Molecules Sofia Izmailov, Fang Liu December 14, 2012 1 Introduction In quantum mechanics,
More informationA One-Slide Summary of Quantum Mechanics
A One-Slide Summary of Quantum Mechanics Fundamental Postulate: O! = a! What is!?! is an oracle! operator wave function (scalar) observable Where does! come from?! is refined Variational Process H! = E!
More informationMSE 102, Fall 2014 Midterm #2. Write your name here [10 points]:
MSE 102, Fall 2014 Midterm #2 Write your name here [10 points]: Instructions: Answer all questions to the best of your abilities. Be sure to write legibly and state your answers clearly. The point values
More informationIntroduction 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 informationElectron States of Diatomic Molecules
IISER Pune March 2018 Hamiltonian for a Diatomic Molecule The hamiltonian for a diatomic molecule can be considered to be made up of three terms Ĥ = ˆT N + ˆT el + ˆV where ˆT N is the kinetic energy operator
More informationComputational Physics. J. M. Thijssen
Computational Physics J. M. Thijssen Delft University of Technology CAMBRIDGE UNIVERSITY PRESS Contents Preface xi 1 Introduction 1 1.1 Physics and computational physics 1 1.2 Classical mechanics and statistical
More informationNumerical Solutions to Partial Differential Equations
Numerical Solutions to Partial Differential Equations Zhiping Li LMAM and School of Mathematical Sciences Peking University The Residual and Error of Finite Element Solutions Mixed BVP of Poisson Equation
More informationAlgorithms and Computational Aspects of DFT Calculations
Algorithms and Computational Aspects of DFT Calculations Part I Juan Meza and Chao Yang High Performance Computing Research Lawrence Berkeley National Laboratory IMA Tutorial Mathematical and Computational
More informationBlock Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems
Mitglied der Helmholtz-Gemeinschaft Block Iterative Eigensolvers for Sequences of Dense Correlated Eigenvalue Problems Birkbeck University, London, June the 29th 2012 Edoardo Di Napoli Motivation and Goals
More informationConstitutive models. Constitutive model: determines P in terms of deformation
Constitutive models Constitutive model: determines P in terms of deformation Elastic material: P depends only on current F Hyperelastic material: work is independent of path strain energy density function
More informationQuantum Chemical Simulations and Descriptors. Dr. Antonio Chana, Dr. Mosè Casalegno
Quantum Chemical Simulations and Descriptors Dr. Antonio Chana, Dr. Mosè Casalegno Classical Mechanics: basics It models real-world objects as point particles, objects with negligible size. The motion
More informationQuantum Chemistry. NC State University. Lecture 5. The electronic structure of molecules Absorption spectroscopy Fluorescence spectroscopy
Quantum Chemistry Lecture 5 The electronic structure of molecules Absorption spectroscopy Fluorescence spectroscopy NC State University 3.5 Selective absorption and emission by atmospheric gases (source:
More informationDISPENSA FEM in MSC. Nastran
DISPENSA FEM in MSC. Nastran preprocessing: mesh generation material definitions definition of loads and boundary conditions solving: solving the (linear) set of equations components postprocessing: visualisation
More informationSame idea for polyatomics, keep track of identical atom e.g. NH 3 consider only valence electrons F(2s,2p) H(1s)
XIII 63 Polyatomic bonding -09 -mod, Notes (13) Engel 16-17 Balance: nuclear repulsion, positive e-n attraction, neg. united atom AO ε i applies to all bonding, just more nuclei repulsion biggest at low
More informationChapter 3 Variational Formulation & the Galerkin Method
Institute of Structural Engineering Page 1 Chapter 3 Variational Formulation & the Galerkin Method Institute of Structural Engineering Page 2 Today s Lecture Contents: Introduction Differential formulation
More informationSession 1. Introduction to Computational Chemistry. Computational (chemistry education) and/or (Computational chemistry) education
Session 1 Introduction to Computational Chemistry 1 Introduction to Computational Chemistry Computational (chemistry education) and/or (Computational chemistry) education First one: Use computational tools
More informationA Hybrid Method for the Wave Equation. beilina
A Hybrid Method for the Wave Equation http://www.math.unibas.ch/ beilina 1 The mathematical model The model problem is the wave equation 2 u t 2 = (a 2 u) + f, x Ω R 3, t > 0, (1) u(x, 0) = 0, x Ω, (2)
More informationStructure of diatomic molecules
Structure of diatomic molecules January 8, 00 1 Nature of molecules; energies of molecular motions Molecules are of course atoms that are held together by shared valence electrons. That is, most of each
More informationREVIEW : INTRODUCTION TO THE MOLECULAR ORIGINS OF MECHANICAL PROPERTIES QUANTITATIVE TREATMENT OF INTERATOMIC BONDING : THE LENNARD-JONES POTENTIAL
LECTURE #19 : 3.11 MECANICS OF MATERIALS F3 INSTRUCTOR : Professor Christine Ortiz OFFICE : 13-422 PONE : 452-384 WWW : http://web.mit.edu/cortiz/www REVIEW : INTRODUCTION TO TE MOLECULAR ORIGINS OF MECANICAL
More informationConvergence and optimality of an adaptive FEM for controlling L 2 errors
Convergence and optimality of an adaptive FEM for controlling L 2 errors Alan Demlow (University of Kentucky) joint work with Rob Stevenson (University of Amsterdam) Partially supported by NSF DMS-0713770.
More informationYingwei Wang Computational Quantum Chemistry 1 Hartree energy 2. 2 Many-body system 2. 3 Born-Oppenheimer approximation 2
Purdue University CHM 67300 Computational Quantum Chemistry REVIEW Yingwei Wang October 10, 2013 Review: Prof Slipchenko s class, Fall 2013 Contents 1 Hartree energy 2 2 Many-body system 2 3 Born-Oppenheimer
More information1. Hydrogen atom in a box
1. Hydrogen atom in a box Recall H atom problem, V(r) = -1/r e r exact answer solved by expanding in Gaussian basis set, had to solve secular matrix involving matrix elements of basis functions place atom
More informationThe calculation of the universal density functional by Lieb maximization
The calculation of the universal density functional by Lieb maximization Trygve Helgaker, Andy Teale, and Sonia Coriani Centre for Theoretical and Computational Chemistry (CTCC), Department of Chemistry,
More informationPerformance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis
Performance Evaluation of Various Smoothed Finite Element Methods with Tetrahedral Elements in Large Deformation Dynamic Analysis Ryoya IIDA, Yuki ONISHI, Kenji AMAYA Tokyo Institute of Technology, Japan
More informationGame Physics. Game and Media Technology Master Program - Utrecht University. Dr. Nicolas Pronost
Game and Media Technology Master Program - Utrecht University Dr. Nicolas Pronost Soft body physics Soft bodies In reality, objects are not purely rigid for some it is a good approximation but if you hit
More informationBand calculations: Theory and Applications
Band calculations: Theory and Applications Lecture 2: Different approximations for the exchange-correlation correlation functional in DFT Local density approximation () Generalized gradient approximation
More information3: Many electrons. Orbital symmetries. l =2 1. m l
3: Many electrons Orbital symmetries Atomic orbitals are labelled according to the principal quantum number, n, and the orbital angular momentum quantum number, l. Electrons in a diatomic molecule experience
More informationAlkali metals show splitting of spectral lines in absence of magnetic field. s lines not split p, d lines split
Electron Spin Electron spin hypothesis Solution to H atom problem gave three quantum numbers, n,, m. These apply to all atoms. Experiments show not complete description. Something missing. Alkali metals
More informationVariational Methods for Electronic Structure
Variational Methods for Electronic Structure The hydrogen atom is a two-body system consisting of a proton and an electron. If spin and relativistic effects are ignored, then the Schrödinger equation for
More informationThe successful wavefunction can be written as a determinant: # 1 (2) # 2 (2) Electrons. This can be generalized to our 2N-electron wavefunction:
T2. CNDO to AM1: The Semiempirical Molecular Orbital Models The discussion in sections T2.1 T2.3 applies also to ab initio molecular orbital calculations. T2.1 Slater Determinants Consider the general
More informationNon-Born-Oppenheimer Effects Between Electrons and Protons
Non-Born-Oppenheimer Effects Between Electrons and Protons Kurt Brorsen Department of Chemistry University of Illinois at Urbana-Champaign PI: Sharon Hammes-Schiffer Funding: NSF, AFOSR Computer time:
More informationCoupling atomistic and continuum modelling of magnetism
Coupling atomistic and continuum modelling of magnetism M. Poluektov 1,2 G. Kreiss 2 O. Eriksson 3 1 University of Warwick WMG International Institute for Nanocomposites Manufacturing 2 Uppsala University
More informationContinuum Mechanics and the Finite Element Method
Continuum Mechanics and the Finite Element Method 1 Assignment 2 Due on March 2 nd @ midnight 2 Suppose you want to simulate this The familiar mass-spring system l 0 l y i X y i x Spring length before/after
More informationEXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD
Journal of KONES Powertrain and Transport, Vol. 7, No. EXPERIMENTAL IDENTIFICATION OF HYPERELASTIC MATERIAL PARAMETERS FOR CALCULATIONS BY THE FINITE ELEMENT METHOD Robert Czabanowski Wroclaw University
More information3. Numerical integration
3. Numerical integration... 3. One-dimensional quadratures... 3. Two- and three-dimensional quadratures... 3.3 Exact Integrals for Straight Sided Triangles... 5 3.4 Reduced and Selected Integration...
More informationArbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact
Arbitrary Normal and Tangential Loading Sequences for Circular Hertzian Contact Philip P. Garland 1 and Robert J. Rogers 2 1 School of Biomedical Engineering, Dalhousie University, Canada 2 Department
More informationPDEs, part 1: Introduction and elliptic PDEs
PDEs, part 1: Introduction and elliptic PDEs Anna-Karin Tornberg Mathematical Models, Analysis and Simulation Fall semester, 2013 Partial di erential equations The solution depends on several variables,
More information1 Density functional theory (DFT)
1 Density functional theory (DFT) 1.1 Introduction Density functional theory is an alternative to ab initio methods for solving the nonrelativistic, time-independent Schrödinger equation H Φ = E Φ. The
More information