Geometry optimization of solids
|
|
- Brandon Luke Fields
- 6 years ago
- Views:
Transcription
1 The Minnesota Workshop on ab Initio Modeling in Solid State Chemistry with CRYSTAL Minneapolis, MN(U.S.A.) 9-14 July Geometry optimization of solids Bartolomeo Civalleri Dip. di Chimica IFM, Via P.Giuria 5 Università di Torino bartolomeo.civalleri@unito.it Some slides are courtesy of C.M. Zicovich-Wilson, A. Erba, A. Rimola Geometry Optimization for: 2 - Finding equilibrium nuclear configurations; - Finding and characterizing Transition States; - Including the effect of Pressure on structural properties of crystals; - Accounting for nuclear-relaxation effects on tensorial response properties (elastic, piezoelectric, photoelastic, etc.); - 1
2 Chemical structures 3 Molecule/Cluster (0D) Polimer (1D) Carbon nanotube (9,9) τ Surface (2D) Crystal (3D) τ 1 τ 2 (001) τ 1 τ 2 τ 3 Slab model of the (001) MgO surface Urea bulk How can one define a chemical structure? Geometries 4 The same chemical compound (same atoms, same bonds) may have many nuclear configurations. In general these nuclear configurations have different energies, and experimentally they appear as a statistical distribution of nuclear arrangements. Which of such configurations (geometries) can be adopted as a suitable representative of the compound in order to compute its properties and compare them with observations? 2
3 Cation location in proton and alkali-exchanged chabazite Experimental data O1 O4 O3 O2 SIII SIII SII SI Three general cation positions from experiments: 1. SI center of hexagonal prism 2. SII six-membered ring of the prism 3. SIII and SIII eight-membered ring Neutron powder diffraction, MAS-NMR, MQMAS-NMR measurements [1-3] on H, Li, Na and K chabazites show: H - two sittings: O1 (SIII ) and O3 (SII) Li - two sittings: SII and SIII Na - two sittings: SII and SIII K - located at SIII No cations were found at the SI site. [1] Smith, L.J.; Davidson, A.; Cheetham, A.K. Catal. Lett. 1997, 49, 143 [2] Smith, L.J.; Eckert, H.; Cheetham, A.K. J. Am. Chem. Soc. 2000, 122, [3] Smith, L.J.; Eckert, H.; Cheetham, A.K. Chem. Mater. 2001, 13, 385. Instability of site SI in Na-exchanged chabazite Sodium is not stable at site SI and moves towards site SII, A.M. Ferrari, M. Llunell, R. Orlando, M. Mérawa, P. Ugliengo, Chem. Mater.,
4 7 Geometry Optimization Definition (I) Potential Energy Surface (PES) x i : Cartesian coordinates of the N nuclei E=f(x 1,x 2,x 3, x 3N ) PES describes the energy of the system as a function of its geometry ( hypersurface) PES arises naturally from the Born-Oppenheimer approximation PES makes it possible to discuss chemical structures 8 Internal degrees of freedom E=f(x 1,x 2,x 3, x 3N ) Invariant under translations and rotations internal degrees of freedom Molecules: M = 3N 3 3(2) Cartesian translational rotational 3D crystals: M = 3N 3 + ( 9 3 ) lattice 4
5 9 Coordinate systems There are infinite choices to define the M internal coordinates The shape of the PES depends on the adopted coordinate system The best one should be the one that makes easiest the study of the potential energy function (to be discussed later) Let's assume that v = (v 1, v 2,...,v M ) is the nuclear configuration vector in a given coordinate choice PES: relevant information 10 Minimum for product A Transition state A Second order saddle point Minimum for reactant Second order saddle point Transition state B m Valleyridge Inflection point E=f(v 1,v 2,v 3, v M ) Minimum for product B How to describe the features of the PES? Relevant information (invariant under transformations of the coordinate system) is provided by critical points 5
6 PES: relevant information 11 Minimum for product A Transition state A Second order saddle point Minimum for reactant Second order saddle point Transition state B m Valleyridge Inflection point E=f(v 1,v 2,v 3, v M ) Minimum for product B H ij 2 E = vi vj v Relevant information (invariant under transformations of the coordinate system) is provided by critical points 0 12 Geometry Optimization Definition (II) Chemical Structures Equilibrium geometries: are local minima on the PES (e.g. conformers, isomers, polymorphs, reactants, intermediates and products of a chemical reaction) Transition states: correspond to saddle points (firstorder) on the PES (e.g. transition structures of a chemical reaction) 6
7 PES and related quantities (derivatives) 13 E x, a ) Energy differences ( i E E ; Forces (geometry optimisation) x a i 2 E Vibrations (Phonons) x i x j 2 E Elastic constants a i a j E Pressure V 2 E 2 Bulk modulus V x: atomic coordinates a i : lattice constants V: cell volume 14 Which is the meaning of Geometry Optimization? It is the process of finding minima and saddle points on the PES How to do that? 7
8 Geometry Optimization Tools (I) 15 Characterization of the PES finding the critical (or stationary) points First derivatives Gradient (vector) In classical mechanics: F x =- V/ x force acting on the atoms Second derivatives Hessian matrix Condition: E v i v0 = 0 v 0 is a stationary point ij 2 E = vi vj v v i is a nuclear configuration vector in a given coordinate choice H 0 Eigenvalue equation Geometry Optimization Tools (II) 16 Topological analysis of the PES Minima local or global Equilibrium structures forces=0 Hessian matrix positive definite all the eigenvalues must be positive Critical points of order 0 Saddle points forces=0 one or more eigenvalues of the Hessian matrix are negative Critical points of index 1 (one negative eigenvalue; all the others positive) first order saddle points Transition State structures N-order saddle points n negative eigenvalues (N=M maximum) 8
9 Geometry Optimization Methods (I) 17 Finding an equilibrium geometry involves an unconstrained optimization on the PES Algorithms for energy minimization Energy-based methods (e.g. line optimization,...) Simple, widely applicable, very slow convergence Gradient-based methods (e.g. conjugate gradient, quasi-newton,...) Good convergence, expensive Second derivative methods (e.g. Newton, Newton- Raphson,...) Very fast, very expensive Geometry Optimization: Methods (II) 18 Strategy (in Newton and quasi-newton (QN) methods): PES is approximated as a quadratic function (Taylor expansion truncated to second-order local quadratic region) E(v) = E 0 + g 0 t Δv + (1/2) Δv t H 0 Δv g(v) = g 0 +H 0 Δv where Δv = v - v 0 ; g is the gradient and H is the Hessian g(v) = g 0 +H 0 Δv = 0 (at the minimum) Δv = -H -1 g (Newton s step) 9
10 QN methods: some remarks 19 For a quadratic PES, the minimum can be found in one step if the gradient vector and the Hessian matrix are known Few cycles are required to reach the minimum even on non-quadratic PESs The Hessian calculation is computationally expensive, so QN methods are based on the use of gradients (analytically computed) E ( ) gi = 0 = E v 0 < < < 2 E E Hij = vi v vi vj 0 v 0 QN methods: variable metric 20 During the optimization process, the changes in the gradient between steps provide useful information on the surface curvature (second derivatives) An approximation to the Hessian matrix, H, or its inverse, H -1 is built at the beginning and updated during the optimization process (variable metric) by using information from displacements and gradients. Several Hessian updating schemes exist: Fletcher-Powell Murtagh-Sargent (SR1) Powell-Symmetric-Broyden (PSB) Broyden-Fletcher-Goldfarb-Shanno (BFGS) Note: some of them impose a Hessian matrix positive definite 10
11 Flow chart for Quasi-Newton algorithms Choose coordinate system Input starting geometry Obtain initial estimate of Hessian Calculate energy and gradient 21 Update the Hessian Use Hessian and gradient to take a step If necessary, restrict step size Check for convergence on gradient and displacement no Check for maximum cycles no Update geometry yes yes Done Stop 22 Efficiency of a geometry optimization The efficiency (i.e. number of steps) depends on several factors: 1) Initial geometry (experiment, GUI/Molecular Modelling, Molecular Mechanics, lower level QM) 2) Coordinate system (Cartesian, internal (redundant, natural,...)) 3) Hessian initial guess (Identity, Molecular Mechanics, lower level QM, numerical) 4) Hessian updating (Berny, MS, DFP, BFGS,...) 5) Step size control (simple rescaling, trust radius, RFO,...) 11
12 23 Estimating the Hessian matrix A good estimate of the initial Hessian can drastically reduce the optimization process L Identity Matrix [HESSIDEN] No structure information, hence very inefficient K Numerical Hessian [HESSNUM] Close to exact, hence fast convergence. However, high calculation cost (good for difficult cases, e.g. transition state searches) J Model Hessian [HESSMOD1,HESSMOD2] Good and cheap approximation. Based on valence force fields (Lindh, Schlegel). Significant improvement with respect to the Identity matrix (>50% speed-up). Step size control 24 When the step size is too large a simple scaling can be applied (too crude) In non-quadratic functions the Newton step may not be the best choice and should be controlled: Level-shift trust region: A trust radius of an hyper-sphere in which the function is expected to behave quadratically. A parameter, µ, is computed so that the displacement: v i+1 v i = (H i -Iµ) -1 g(v i ) is kept within the trust region. Ensured to be the best direction within the hyper-sphere of radius, τ = v i+1 v i Line search: A scale factor α i is computed as to reach a minimum along the Newton step direction, according to: v i+1 v i = - α i H i -1 g(v i ) 12
13 Choice of coordinate system 25 Good Performance Quadratic behaviour Coordinate transformation: v i = v i (s 1, s 2,..., s M ) Transformation of F(v) results in a new function F': F'(s 1, s 2,..., s M ) = F(v 1 (s 1, s 2,..., s M ), F(v 2 (s 1, s 2,..., s M ),..., F(v M (s 1, s 2,..., s M ) Ideal transformation gives function F' that behaves quadratically close to the critical point Coordinate systems suitable for crystals: Fractionary coordinates + elastic distorsions (symmetry adapted) Internal valence coordinates (redundant) Coordinate systems for crystals (I) 26 Fractionary coordinates + elastic cell distortions Atom coordinates defined in the basis set of the lattice vectors. Cell parameters obtained from elastic distortions (expansion matrix applied on the lattice vectors) J Easy to keep special positions and cell parameters by imposing symmetry K Delocalised by definition. Not useful for chemical reactivity. L In general, non-quadratic L Low dependancy between coordinates and lattice parameters 13
14 Coordinate systems for crystals (II) 27 Redundant internal coordinates All valence coordinates (bond distances, angles and torsions) are considered to form a set of redundant coordinates (N RC > M) Connectivity is usually defined from vicinity criteria (R AB < R vdw,a +R vdw,b ) Cell parameters are embedded into the redundant coordinates Hessian, gradient and geometry displacements are built in terms of the redundant coordinates The redundancies are eliminated to obtain the actual geometry (in Cartesian coordinate space) using numerical approximations (back-transformation) J Quadratic behaviour J Easy choice of geometrical parameters J Easy to constrain 'chemical' degrees of freedom (bond lengths, angles or dihedrals) Good for reactivity studies L The size of redundant space may be very large (improved in CRYSTAL17) L The back-transformation from redundant to non-redundant space is difficult (improved in CRYSTAL17) Comparison between Frac.+Cell and Redundant System Dim Theor. Method Frac+Cell Redundant C-R Nstep Nstep ΔE (µha) α-quartz 3D B3LYP β-quartz 3D LDA All silica Faujasite 3D PBE KNbO 3 tetragonal 3D BLYP BaTiO 3 orthoromb. 3D LDA Bohemite 3D PW Calcite 3D B3PW NaNO 2 3D BP TiO 2 3D PW ZnGeP 2 3D LDA/ECP ZrO 2 3D PBE Formamide 3D B3LYP Ice 3D PW Oxalic acid 3D B3LYP Oxalic acid 3D HF Al 2 O 3 (0001) 12-layer 2D B3LYP EDI(100)/NH 3 2D B3LYP H 2 O infinite chain 1D PBE Li doped PA 1D B3LYP Tot
15 Geometry optimizer in CRYSTAL 30 Quasi-Newton algorithm BFGS Hessian updating Initial model Hessian (Schlegel) Trust Radius (TR) Newton s step only Linear search (if not TR) Symmetry adapted fractional coordinates (symmetrized directions) Symmetry adapted elastic deformations (symmetrized deformations) Analytical forces for Hartree-Fock and DFT 0D, 1D, 2D, 3D systems All-electron and Effective Core Potential (ECP) Atomic positions: K. Doll, V.R. Saunders, N.M. Harrison, IJQC 82 (2001) 1 Cell parameters: K. Doll, R. Dovesi, R. Orlando, TCA 112 (2004) 394 Geometry optimization type Unconstrained optimization of cell and atomic positions Transition state search (atoms, cell) Constraint optimization: constant volume, fragment, Geometry optimization input block 31 E.g.: α-quartz - hexagonal cell RHF/STO-3G a-qua (in. geom. expt.) CRYSTAL OPTGEOM Optional keywords END ENDG Title Dimensionality of the system Crystallographic information (3D only) Space Group (Fm3m 154) Shift of the origin Lattice parameters (hexagonal) Number of non equivalent atoms Atomic number and fractional coordinates Geometry optim. input block Optimization of atomic positions + cell End of geom. opt. input block End of geometry input section Optional keywords are related to: Modification of the algorithm (Initial Hessian, Hessian updating, ) Convergence criteria and tolerances (TOLDEE, TOLDEX, TOLDEG) Constrained geometry optimization See CRYSTAL Tutorials (Geometry Optimization) 15
16 Geometry optimization input block 32 E.g.: α-quartz - hexagonal cell RHF/STO-3G Redundant internal coordinates a-qua (in. geom. expt.) CRYSTAL OPTGEOM INTREDUN Optional keywords END ENDG Full geometry optimization of atomic positions and lattice param.s 33 Geometry Optimization is important for Equilibrium bond lengths and angles Equilibrium cell parameters Discriminating between competing structures Surface relaxation and reconstruction Pressure-driven phase transitions Starting point for many advanced studies 16
17 12/07/17 38 Constrained geometry optimization Optimization in Cartesian/fractionary coordinates (default) 1. Fixing Cartesian coordinates (FRAGMENT) 2. Constant volume geometry optimization (CVOLOPT) 3. Fixing a lattice parameter (CRYDEF FIXDEF) 4. Constant pressure geometry optimization (EXTPRESS) Optimization in redundant internal coordinates (INTREDUN) 1. Fixing a set of internal coordinates (FREEZINT) 2. Adding a new internal coordinate (DEFLNGS, ) 3. Fixing a given internal coordinate (LNGSFROZEN, ) 48 A special acknowledgement Claudio M. Zicovich-Wilson ( ) 17
Reactivity and Organocatalysis. (Adalgisa Sinicropi and Massimo Olivucci)
Reactivity and Organocatalysis (Adalgisa Sinicropi and Massimo Olivucci) The Aldol Reaction - O R 1 O R 1 + - O O OH * * H R 2 R 1 R 2 The (1957) Zimmerman-Traxler Transition State Counterion (e.g. Li
More informationAb initio structure prediction for molecules and solids
Ab initio structure prediction for molecules and solids Klaus Doll Max-Planck-Institute for Solid State Research Stuttgart Chemnitz, June/July 2010 Contents structure prediction: 1) global search on potential
More informationElectronic Supplementary Information
Electronic Supplementary Material (ESI) for CrystEngComm. This journal is The Royal Society of Chemistry 2014 Electronic Supplementary Information Configurational and energetical study of the (100) and
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 informationExploring the energy landscape
Exploring the energy landscape ChE210D Today's lecture: what are general features of the potential energy surface and how can we locate and characterize minima on it Derivatives of the potential energy
More informationInternational Journal of Quantum Chemistry
International Journal of Quantum Chemistry First-principles calculation of second-order elastic constants and equations of state for Lithium Azide, LiN, and Lead Azide, Pb(N ) Journal: International Journal
More informationIntroduction to Geometry Optimization. Computational Chemistry lab 2009
Introduction to Geometry Optimization Computational Chemistry lab 9 Determination of the molecule configuration H H Diatomic molecule determine the interatomic distance H O H Triatomic molecule determine
More informationASSESSMENT OF DFT METHODS FOR SOLIDS
MSSC2009 - Ab Initio Modeling in Solid State Chemistry ASSESSMENT OF DFT METHODS FOR SOLIDS Raffaella Demichelis Università di Torino Dipartimento di Chimica IFM 1 MSSC2009 - September, 10 th 2009 Table
More informationVibrational frequencies in solids: tools and tricks
Vibrational frequencies in solids: tools and tricks Roberto Dovesi Gruppo di Chimica Teorica Università di Torino Torino, 4-9 September 2016 This morning 3 lectures: R. Dovesi Generalities on vibrations
More informationThermodynamics of Solids: Harmonic and Quasi-harmonic Approximations
Thermodynamics of Solids: Harmonic and Quasi-harmonic Approximations, USA, July 9-14, 2017 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it 2017 Outline -
More informationStrain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity
Strain-related Tensorial Properties: Elasticity, Piezoelectricity and Photoelasticity Torino, Italy, September 4-9, 2016 Alessandro Erba Dipartimento di Chimica, Università di Torino (Italy) alessandro.erba@unito.it
More informationSupporting Information. for
Electronic Supplementary Material (ESI) for Physical Chemistry Chemical Physics. This journal is the Owner Societies 2016 Supporting Information for Isoreticular Zirconium-Based Metal-Organic Frameworks:
More informationGeometry optimization
Geometry optimization Trygve Helgaker Centre for Theoretical and Computational Chemistry Department of Chemistry, University of Oslo, Norway European Summer School in Quantum Chemistry (ESQC) 211 Torre
More informationStructural Calculations phase stability, surfaces, interfaces etc
Structural Calculations phase stability, surfaces, interfaces etc Keith Refson STFC Rutherford Appleton Laboratory September 19, 2007 Phase Equilibrium 2 Energy-Volume curves..................................................................
More informationCrystal structure prediction
Crystal structure prediction Klaus Doll Institute for Mathematical Physics, TU Braunschweig Max Planck Institute for Solid State Research, Stuttgart MSSC2011, Turin, September 2011 Motivation structure
More informationExample questions for Molecular modelling (Level 4) Dr. Adrian Mulholland
Example questions for Molecular modelling (Level 4) Dr. Adrian Mulholland 1) Question. Two methods which are widely used for the optimization of molecular geometies are the Steepest descents and Newton-Raphson
More informationThe Interpretation of the Short Range Disorder in the Fluorene- TCNE Crystal Structure
Int. J. Mol. Sci. 2004, 5, 93-100 International Journal of Molecular Sciences ISSN 1422-0067 2004 by MDPI www.mdpi.net/ijms/ The Interpretation of the Short Range Disorder in the Fluorene- TCNE Crystal
More informationTools for QM studies of large systems
Tools for QM studies of large systems Automated, hessian-free saddle point search & characterization QM/MM implementation for zeolites Shaama Mallikarjun Sharada Advisors: Prof. Alexis T Bell, Prof. Martin
More informationAn Introduction to Quantum Chemistry and Potential Energy Surfaces. Benjamin G. Levine
An Introduction to Quantum Chemistry and Potential Energy Surfaces Benjamin G. Levine This Week s Lecture Potential energy surfaces What are they? What are they good for? How do we use them to solve chemical
More informationMaking Electronic Structure Theory work
Making Electronic Structure Theory work Ralf Gehrke Fritz-Haber-Institut der Max-Planck-Gesellschaft Berlin, 23 th June 2009 Hands-on Tutorial on Ab Initio Molecular Simulations: Toward a First-Principles
More informationConvex Optimization CMU-10725
Convex Optimization CMU-10725 Quasi Newton Methods Barnabás Póczos & Ryan Tibshirani Quasi Newton Methods 2 Outline Modified Newton Method Rank one correction of the inverse Rank two correction of the
More informationDFT and beyond: Hands-on Tutorial Workshop Tutorial 1: Basics of Electronic Structure Theory
DFT and beyond: Hands-on Tutorial Workshop 2011 Tutorial 1: Basics of Electronic Structure Theory V. Atalla, O. T. Hofmann, S. V. Levchenko Theory Department, Fritz-Haber-Institut der MPG Berlin July 13,
More informationIntroduction to Tutorial Sessions
Michigan Tech University - Washington State University - University of Torino Ab initio Simulation of Crystalline Systems ASCS2006 September 17-22, 2006 - Spokane, Washington (USA) Introduction to Tutorial
More informationGround-state Structure and Dynamics
Ground-state Structure and Dynamics Jonathan Yates jonathan.yates@materials.ox.ac.uk Materials Modelling Laboratory, Oxford Materials For a given set of atomic positions the ions will experience a force
More informationIntroduction to first-principles modelling and CASTEP
to first-principles modelling and Phil Hasnip to + Atomistic Simulations If we know what the bonding in a material is beforehand, then we can often find good expressions for the forces between atoms, e.g.
More informationOptimization Methods
Optimization Methods Decision making Examples: determining which ingredients and in what quantities to add to a mixture being made so that it will meet specifications on its composition allocating available
More informationOne-Electron Properties of Solids
One-Electron Properties of Solids Alessandro Erba Università di Torino alessandro.erba@unito.it most slides are courtesy of R. Orlando and B. Civalleri Energy vs Wave-function Energy vs Wave-function Density
More informationMethods that avoid calculating the Hessian. Nonlinear Optimization; Steepest Descent, Quasi-Newton. Steepest Descent
Nonlinear Optimization Steepest Descent and Niclas Börlin Department of Computing Science Umeå University niclas.borlin@cs.umu.se A disadvantage with the Newton method is that the Hessian has to be derived
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 informationTransition states and reaction paths
Transition states and reaction paths Lab 4 Theoretical background Transition state A transition structure is the molecular configuration that separates reactants and products. In a system with a single
More informationAB INITIO MODELING OF ALKALI METAL CHALCOGENIDES USING SOGGA THEORY
Int. J. Chem. Sci.: 13(4), 215, 163-1638 ISSN 972-768X www.sadgurupublications.com AB INITIO MODELING OF ALALI METAL CHALCOGENIDES USING SOGGA THEORY HITESH CHANDRA SWARNAR and GUNJAN ARORA a,* Department
More informationFigure 1: Transition State, Saddle Point, Reaction Pathway
Computational Chemistry Workshops West Ridge Research Building-UAF Campus 9:00am-4:00pm, Room 009 Electronic Structure - July 19-21, 2016 Molecular Dynamics - July 26-28, 2016 Potential Energy Surfaces
More informationGeometry Optimization
Geometry Optimization Matt Probert Condensed Matter Dynamics Group Department of Physics, University of York, U.K. http://www-users.york.ac.uk/~mijp1 Overview of lecture n Background Theory n Hellman-Feynman
More informationQuasi-Newton Methods
Newton s Method Pros and Cons Quasi-Newton Methods MA 348 Kurt Bryan Newton s method has some very nice properties: It s extremely fast, at least once it gets near the minimum, and with the simple modifications
More informationPractical calculations using first-principles QM Convergence, convergence, convergence
Practical calculations using first-principles QM Convergence, convergence, convergence Keith Refson STFC Rutherford Appleton Laboratory September 18, 2007 Results of First-Principles Simulations..........................................................
More informationMultivariate Newton Minimanization
Multivariate Newton Minimanization Optymalizacja syntezy biosurfaktantu Rhamnolipid Rhamnolipids are naturally occuring glycolipid produced commercially by the Pseudomonas aeruginosa species of bacteria.
More informationNonlinear Optimization: What s important?
Nonlinear Optimization: What s important? Julian Hall 10th May 2012 Convexity: convex problems A local minimizer is a global minimizer A solution of f (x) = 0 (stationary point) is a minimizer A global
More informationPhonon calculations with SCAN
Workshop on the SCAN density functional: Fundamentals, practices, and extensions Temple university, Philadelphia May 18th, 2017 Hands-on tutorial 3 Phonon calculations with SCAN Yubo Zhang and Jianwei
More informationEAD 115. Numerical Solution of Engineering and Scientific Problems. David M. Rocke Department of Applied Science
EAD 115 Numerical Solution of Engineering and Scientific Problems David M. Rocke Department of Applied Science Multidimensional Unconstrained Optimization Suppose we have a function f() of more than one
More information5 Quasi-Newton Methods
Unconstrained Convex Optimization 26 5 Quasi-Newton Methods If the Hessian is unavailable... Notation: H = Hessian matrix. B is the approximation of H. C is the approximation of H 1. Problem: Solve min
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 informationSupplementary Information for. Universal elastic-hardening-driven mechanical instability in α-quartz and quartz. homeotypes under pressure
Supplementary Information for Universal elastic-hardening-driven mechanical instability in α-quartz and quartz homeotypes under pressure Juncai Dong, Hailiang Zhu, and Dongliang Chen * Beijing Synchrotron
More informationPotentials, periodicity
Potentials, periodicity Lecture 2 1/23/18 1 Survey responses 2 Topic requests DFT (10), Molecular dynamics (7), Monte Carlo (5) Machine Learning (4), High-throughput, Databases (4) NEB, phonons, Non-equilibrium
More informationMolecular Simulation II. Classical Mechanical Treatment
Molecular Simulation II Quantum Chemistry Classical Mechanics E = Ψ H Ψ ΨΨ U = E bond +E angle +E torsion +E non-bond Jeffry D. Madura Department of Chemistry & Biochemistry Center for Computational Sciences
More informationImproving the Convergence of Back-Propogation Learning with Second Order Methods
the of Back-Propogation Learning with Second Order Methods Sue Becker and Yann le Cun, Sept 1988 Kasey Bray, October 2017 Table of Contents 1 with Back-Propagation 2 the of BP 3 A Computationally Feasible
More informationA local MP2 periodic study of crystalline argon
Journal of Physics: Conference Series A local MP2 periodic study of crystalline argon To cite this article: S Casassa et al 2008 J. Phys.: Conf. Ser. 117 012007 Recent citations - Laplace transformed MP2
More informationBefore we start: Important setup of your Computer
Before we start: Important setup of your Computer change directory: cd /afs/ictp/public/shared/smr2475./setup-config.sh logout login again 1 st Tutorial: The Basics of DFT Lydia Nemec and Oliver T. Hofmann
More informationQuasi-Newton methods: Symmetric rank 1 (SR1) Broyden Fletcher Goldfarb Shanno February 6, / 25 (BFG. Limited memory BFGS (L-BFGS)
Quasi-Newton methods: Symmetric rank 1 (SR1) Broyden Fletcher Goldfarb Shanno (BFGS) Limited memory BFGS (L-BFGS) February 6, 2014 Quasi-Newton methods: Symmetric rank 1 (SR1) Broyden Fletcher Goldfarb
More informationC:\Users\Leonardo\Desktop\README_ELMO_NOTES.txt Montag, 7. Juli :48
*********** * Input * *********** Besides the normal GAMESS-UK input directives, to perform an ELMO calculation one needs to specify the elmo keyword and to provide additional instructions that are contained
More informationSelectivity in the initial C-H bond cleavage of n-butane on PdO(101)
Supporting Information for Selectivity in the initial C-H bond cleavage of n-butane on PdO(101) Can Hakanoglu (a), Feng Zhang (a), Abbin Antony (a), Aravind Asthagiri (b) and Jason F. Weaver (a) * (a)
More informationTECHNIQUES TO LOCATE A TRANSITION STATE
32 III. TECHNIQUES TO LOCATE A TRANSITION STATE In addition to the location of minima, the calculation of transition states and reaction pathways is an interesting task for Quantum Chemistry (QC). The
More information22 Path Optimisation Methods
22 Path Optimisation Methods 204 22 Path Optimisation Methods Many interesting chemical and physical processes involve transitions from one state to another. Typical examples are migration paths for defects
More informationOptimization and Root Finding. Kurt Hornik
Optimization and Root Finding Kurt Hornik Basics Root finding and unconstrained smooth optimization are closely related: Solving ƒ () = 0 can be accomplished via minimizing ƒ () 2 Slide 2 Basics Root finding
More informationNonlinearOptimization
1/35 NonlinearOptimization Pavel Kordík Department of Computer Systems Faculty of Information Technology Czech Technical University in Prague Jiří Kašpar, Pavel Tvrdík, 2011 Unconstrained nonlinear optimization,
More informationNANOSTRUCTURED OXIDES: NEW MATERIALS FOR ENERGY AND ENVIRONMENT
NANOSTRUCTURED OXIDES: NEW MATERIALS FOR ENERGY AND ENVIRONMENT Quantum Chemistry Laboratory Dipartimento di Scienza dei Materiali Università Milano-Bicocca http://www.mater.unimib.it/utenti/pacchioni
More informationLine Search Methods for Unconstrained Optimisation
Line Search Methods for Unconstrained Optimisation Lecture 8, Numerical Linear Algebra and Optimisation Oxford University Computing Laboratory, MT 2007 Dr Raphael Hauser (hauser@comlab.ox.ac.uk) The Generic
More informationHigher-Order Methods
Higher-Order Methods Stephen J. Wright 1 2 Computer Sciences Department, University of Wisconsin-Madison. PCMI, July 2016 Stephen Wright (UW-Madison) Higher-Order Methods PCMI, July 2016 1 / 25 Smooth
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 informationNumerical Optimization Professor Horst Cerjak, Horst Bischof, Thomas Pock Mat Vis-Gra SS09
Numerical Optimization 1 Working Horse in Computer Vision Variational Methods Shape Analysis Machine Learning Markov Random Fields Geometry Common denominator: optimization problems 2 Overview of Methods
More informationWANG Bingwu, XU Guangxian & CHEN Zhida
Science in China Ser. B Chemistry 2004 Vol.47 No.2 91 97 91 Study on the covalence of Cu and chemical bonding in an inorganic fullerene-like molecule, [CuCl] 20 [Cp*FeP 5 ] 12 [Cu-(CH 3 CN) 2 Cl ] 5, by
More informationThe rôle of surfaces in prebiotic chemistry. Piero Ugliengo University of Torino Dip. Chimica Via P. Giuria, 7 - Torino
The rôle of surfaces in prebiotic chemistry Piero Ugliengo University of Torino Dip. Chimica Via P. Giuria, 7 - Torino 1 A collaborative work Dip. Chimica - NIS Centre Torino University TO BCN Dep. Quimica
More informationSECTION: CONTINUOUS OPTIMISATION LECTURE 4: QUASI-NEWTON METHODS
SECTION: CONTINUOUS OPTIMISATION LECTURE 4: QUASI-NEWTON METHODS HONOUR SCHOOL OF MATHEMATICS, OXFORD UNIVERSITY HILARY TERM 2005, DR RAPHAEL HAUSER 1. The Quasi-Newton Idea. In this lecture we will discuss
More informationCrystal Relaxation, Elasticity, and Lattice Dynamics
http://exciting-code.org Crystal Relaxation, Elasticity, and Lattice Dynamics Pasquale Pavone Humboldt-Universität zu Berlin http://exciting-code.org PART I: Structure Optimization Pasquale Pavone Humboldt-Universität
More informationComparison of methods for finding saddle points without knowledge of the final states
JOURNAL OF CHEMICAL PHYSICS VOLUME 121, NUMBER 20 22 NOVEMBER 2004 Comparison of methods for finding saddle points without knowledge of the final states R. A. Olsen and G. J. Kroes Leiden Institute of
More informationMATH 4211/6211 Optimization Quasi-Newton Method
MATH 4211/6211 Optimization Quasi-Newton Method Xiaojing Ye Department of Mathematics & Statistics Georgia State University Xiaojing Ye, Math & Stat, Georgia State University 0 Quasi-Newton Method Motivation:
More informationSupplementary material for Electronic Structure of IrO 2 : the Role of the Metal D Orbitals
Supplementary material for Electronic Structure of IrO 2 : the Role of the Metal D Orbitals Yuan Ping 1, Giulia Galli 2 and William A. Goddard III 3 1 Joint Center for Artificial Photosynthesis, Lawrence
More informationAlgorithms for Constrained Optimization
1 / 42 Algorithms for Constrained Optimization ME598/494 Lecture Max Yi Ren Department of Mechanical Engineering, Arizona State University April 19, 2015 2 / 42 Outline 1. Convergence 2. Sequential quadratic
More informationOrbital Density Dependent Functionals
Orbital Density Dependent Functionals S. Kluepfel1, P. Kluepfel1, Hildur Guðmundsdóttir1 and Hannes Jónsson1,2 1. Univ. of Iceland; 2. Aalto University Outline: Problems with GGA approximation (PBE, RPBE,...)
More informationQuasi-Newton methods for minimization
Quasi-Newton methods for minimization Lectures for PHD course on Numerical optimization Enrico Bertolazzi DIMS Universitá di Trento November 21 December 14, 2011 Quasi-Newton methods for minimization 1
More informationApplications: Molecular crystals Graphite MgO(001)/CO MIL-53(Al) 2
Bartolomeo Civalleri Voice: Loredana Valenzano B3LYP augmented with an empirical dispersion term (B3LYP-D*) as applied to solids Università di Torino Dipartimento di Chimica IFM & NIS Torino - MSSC2009-10/09/09
More informationNonlinear Programming
Nonlinear Programming Kees Roos e-mail: C.Roos@ewi.tudelft.nl URL: http://www.isa.ewi.tudelft.nl/ roos LNMB Course De Uithof, Utrecht February 6 - May 8, A.D. 2006 Optimization Group 1 Outline for week
More informationShiqian Ma, MAT-258A: Numerical Optimization 1. Chapter 3. Gradient Method
Shiqian Ma, MAT-258A: Numerical Optimization 1 Chapter 3 Gradient Method Shiqian Ma, MAT-258A: Numerical Optimization 2 3.1. Gradient method Classical gradient method: to minimize a differentiable convex
More informationarxiv:physics/ v1 [physics.comp-ph] 19 Dec 2005
Analytical Hartree-Fock gradients with respect to the cell parameter: systems periodic in one and two dimensions K. Doll Institut für Mathematische Physik, TU Braunschweig, Mendelssohnstraße 3, D-38106
More informationTranslation Symmetry, Space Groups, Bloch functions, Fermi energy
Translation Symmetry, Space Groups, Bloch functions, Fermi energy Roberto Orlando and Silvia Casassa Università degli Studi di Torino July 20, 2015 School Ab initio Modelling of Solids (UniTo) Symmetry
More informationSupplementary Materials
Supplementary Materials Atomistic Origin of Brittle Failure of Boron Carbide from Large Scale Reactive Dynamics Simulations; Suggestions toward Improved Ductility Qi An and William A. Goddard III * Materials
More informationStructure and Dynamics : An Atomic View of Materials
Structure and Dynamics : An Atomic View of Materials MARTIN T. DOVE Department ofearth Sciences University of Cambridge OXFORD UNIVERSITY PRESS Contents 1 Introduction 1 1.1 Observations 1 1.1.1 Microscopic
More informationUnconstrained optimization
Chapter 4 Unconstrained optimization An unconstrained optimization problem takes the form min x Rnf(x) (4.1) for a target functional (also called objective function) f : R n R. In this chapter and throughout
More informationLecture V. Numerical Optimization
Lecture V Numerical Optimization Gianluca Violante New York University Quantitative Macroeconomics G. Violante, Numerical Optimization p. 1 /19 Isomorphism I We describe minimization problems: to maximize
More informationComparative study of Optimization methods for Unconstrained Multivariable Nonlinear Programming Problems
International Journal of Scientific and Research Publications, Volume 3, Issue 10, October 013 1 ISSN 50-3153 Comparative study of Optimization methods for Unconstrained Multivariable Nonlinear Programming
More informationProgramming, numerics and optimization
Programming, numerics and optimization Lecture C-3: Unconstrained optimization II Łukasz Jankowski ljank@ippt.pan.pl Institute of Fundamental Technological Research Room 4.32, Phone +22.8261281 ext. 428
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 informationAndré Schleife Department of Materials Science and Engineering
André Schleife Department of Materials Science and Engineering Yesterday you (should have) learned this: http://upload.wikimedia.org/wikipedia/commons/e/ea/ Simple_Harmonic_Motion_Orbit.gif 1. deterministic
More informationChemistry 334 Part 2: Computational Quantum Chemistry
Chemistry 334 Part 2: Computational Quantum Chemistry 1. Definition Louis Scudiero, Ben Shepler and Kirk Peterson Washington State University January 2006 Computational chemistry is an area of theoretical
More informationNonlinear Optimization for Optimal Control
Nonlinear Optimization for Optimal Control Pieter Abbeel UC Berkeley EECS Many slides and figures adapted from Stephen Boyd [optional] Boyd and Vandenberghe, Convex Optimization, Chapters 9 11 [optional]
More information7/29/2014. Electronic Structure. Electrons in Momentum Space. Electron Density Matrices FKF FKF. Ulrich Wedig
Electron Density Matrices Density matrices Γ, an alternative to the wavefunction Ψ, for the description of a quantum system Electronic Structure The N-particle density matrix Electrons in Momentum Space
More informationSTRUCTURAL DEFECTS IN IMIDATES : AN AB INITIO STUDY
Int. J. Chem. Sci.: 9(4), 2011, 1763-1767 ISSN 0972-768X www.sadgurupublications.com STRUCTURAL DEFECTS IN IMIDATES : AN AB INITIO STUDY M. FATHIMA BEGUM, HEMA TRESA VARGHESE a, Y. SHEENA MARY a, C. YOHANNAN
More informationAb-initio Electronic Structure Calculations β and γ KNO 3 Energetic Materials
ISSN 0974-9373 Vol. 15 No.3 (2011) Journal of International Academy of Physical Sciences pp. 337-344 Ab-initio Electronic Structure Calculations of α, β and γ KNO 3 Energetic Materials Pradeep Jain and
More informationAtom in Molecules a Quantum Theory (AIM)
Atom in Molecules a Quantum Theory (AIM) Richard F. W. Bader* The primary purpose in postulating the existence of atoms in molecules is a consequence of the observation that atoms or groupings of atoms
More informationCHAPTER 2. ELECTRONIC STRUCTURE AND GROUND STATE PROPERTIES OF M 2 O (M: Li, Na, K, Rb)
30 CHAPTER 2 ELECTRONIC STRUCTURE AND GROUND STATE PROPERTIES OF M 2 O (M: Li, Na, K, Rb) 2.1 INTRODUCTION Oxides of alkali metals (M 2 O) (M: Li, Na, K, Rb) play an important role in reducing the work
More informationOptimization for neural networks
0 - : Optimization for neural networks Prof. J.C. Kao, UCLA Optimization for neural networks We previously introduced the principle of gradient descent. Now we will discuss specific modifications we make
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 informationMolecular Modeling of Inorganic Compounds
Peter Comba, Trevor W. Hambley and Bodo Martin Molecular Modeling of Inorganic Compounds Third Completely Revised and Enlarged Edition WILEY- VCH WILEY-VCH Verlag GmbH & Co. KGaA v Contents Preface to
More information2. Quasi-Newton methods
L. Vandenberghe EE236C (Spring 2016) 2. Quasi-Newton methods variable metric methods quasi-newton methods BFGS update limited-memory quasi-newton methods 2-1 Newton method for unconstrained minimization
More informationDFT in practice. Sergey V. Levchenko. Fritz-Haber-Institut der MPG, Berlin, DE
DFT in practice Sergey V. Levchenko Fritz-Haber-Institut der MPG, Berlin, DE Outline From fundamental theory to practical solutions General concepts: - basis sets - integrals and grids, electrostatics,
More informationQuantum-chemical approach to cohesive properties of metallic beryllium
Quantum-chemical approach to cohesive properties of metallic beryllium Elena Voloshina 1, Beate Paulus 2, and Hermann Stoll 3 1 Max-Planck-Institut für Physik komplexer Systeme, Nöthnitzer Straße 38, 01187
More informationOptimization II: Unconstrained Multivariable
Optimization II: Unconstrained Multivariable CS 205A: Mathematical Methods for Robotics, Vision, and Graphics Justin Solomon CS 205A: Mathematical Methods Optimization II: Unconstrained Multivariable 1
More informationNMR and IR spectra & vibrational analysis
Lab 5: NMR and IR spectra & vibrational analysis A brief theoretical background 1 Some of the available chemical quantum methods for calculating NMR chemical shifts are based on the Hartree-Fock self-consistent
More information( ) + ( ) + ( ) = 0.00
Chem 3502/4502 Physical Chemistry II (Quantum Mechanics) 3 Credits Spring Semester 2006 Christopher J. Cramer Lecture 32, April 14, 2006 (Some material in this lecture has been adapted from Cramer, C.
More informationThe Boron Buckyball has an Unexpected T h Symmetry
The Boron Buckyball has an Unexpected T h Symmetry G. Gopakumar, Minh Tho Nguyen, and Arnout Ceulemans* Department of Chemistry and Institute for Nanoscale Physics and Chemistry, University of Leuven,
More informationChapter 4. Unconstrained optimization
Chapter 4. Unconstrained optimization Version: 28-10-2012 Material: (for details see) Chapter 11 in [FKS] (pp.251-276) A reference e.g. L.11.2 refers to the corresponding Lemma in the book [FKS] PDF-file
More information