Ab initio molecular dynamics : BO
|
|
- Melina Farmer
- 5 years ago
- Views:
Transcription
1 School on First Principles Simulations, JNCASR, 2010 Ab initio molecular dynamics : BO Vardha Srinivasan IISER Bhopal
2 Why ab initio MD Free of parametrization and only fundamental constants required. Bond forming and breaking processes can be treated accurately. Transferable unlike many classical MD methods. On the fly computation of potential energy surface, so avoids dimensionality bottleneck.
3 Quantum dynamics i Ψ(r, R; t) =ĤΨ(r, R; t) t r electrons R (RI) ions Ĥ = 2 2M I I 2 2 2m e e 2 + ˆV e e + ˆV I I + ˆV e I electron-electron repulsion Ψ(r, R; t) Complicated Full quantum dynamics quite messy and not needed in most cases. ion-ion repulsion electron-ion attraction
4 Decoupling electrons and ions Ĥ e = 2 2m e 2 e + ˆV e e + ˆV I I + ˆV e I First solve for clamped nuclei Ĥ e ψ k (r; R) =E k (R)ψ k (r; R) Instantaneous representation for full solution Ψ(r, R; t) = l ψ l (r; R)χ l (R; t) Integrate out r dependence dr ψ k (r; R) ĤΨ(r, R; t)
5 Decoupling electrons and ions i t χ(r; t) = 2 2 I + E k (R) 2M I Exact non-adiabatic coupling C kl = drψk 2 2 I 2M I C kl = C k δ kl C kl =0 Adiabatic Approx. ψ k + 1 M I χ(r; t)+ l C kl χ l (R; t) drψk [ i I ] [ i I ] leads to the Born-Oppenheimer approximation
6 Born-Oppenheimer Dynamics i t χ(r; t) = 2 2 I + E k (R) 2M I χ(r; t) Effective Kinetic energy potential energy Dynamics of ions does not change the state of the electronic sub-system. Thus, Ψ(r, R; t) =χ(r; t)ψ(r; R) Ek(R) defines an energy eigenvalue for every configuration R - Potential Energy Surface (of the k th electronic eigenstate) Classical ions would move on the ground-state surface.
7 Classical Dynamics for ions Let s define a classical Lagrangean for B-O dynamics L BO = 1 2 M IṘ2 I E BO (R) Equations of motion are obtained by d dt ṘL BO R L BO =0 Thus we have M I R = I E BO (R) E BO = ψ 0 Ĥe ψ 0 =minψ Ĥe ψ ψ The potential energy at every instant is obtained by solving the electronic problem variationally.
8 DFT Energy and Forces E BO (R) =E DFT (R)+E II (R) =E tot Where the first term is just the DFT energy computed from the ground-state density n(r) by methods discussed earlier in this school. Hellmann-Feynman forces : F i = I E BO = drn(r) I V ext (r) I E II (R) External potential on electrons
9
10
11 B-O Molecular Dynamics At every time step we need to perform a minimization to reach self consistency. This could be time-consuming (especially around 1985 where iterative diagonalization schemes were not employed for first-principles calculations). The accuracy of the simulation critically depends on the accuracy of the SCF minimization. (Energy drifts) The performance depends on the algorithms used to extrapolate the wavefunctions from the previous steps.
12 B-O Molecular Dynamics BO MD of butadiene molecule by two different routes to guess initial wavefunction at each instant P. Pulay and G. Fogarasi, Chem. Phys. Lett. 386, 272 (2004). Energy drifts are a concern in BO MD but can be reduced by choosing good extrapolation schemes and other clever ways of propagating the wavefunction.
13 Choosing parameters for MD A large time step means lesser overall computational cost. But need to sample fastest ionic motion. Roughly, dt ~ dtmax where dtmax=1/ωmax=period of fastest phonon. SCF convergence at each step should be very good to conserve energy, avoid systematic drifts and ensure accurate forces. Error on DFT energy is quadratic in SCF error of charge density whereas error on forces is linear.
14 Some extensions
15
16 &control calculation='md' restart_mode='from_scratch', pseudo_dir = '$PSEUDO_DIR/', outdir='$tmp_dir/', dt=20, nstep=100, 1 a.u. = 2.4 * s disk_io='high' / &system ibrav= 2, celldm(1)=10.18, nat= 2, ntyp= 1, ecutwfc = 8.0, nosym=.true. / &electrons conv_thr = 1.0d-8 mixing_beta = 0.7 / &ions pot_extrapolation='second-order' wfc_extrapolation='second-order' / ATOMIC_SPECIES Si Si.pz-vbc.UPF ATOMIC_POSITIONS Si Si K_POINTS {automatic}
17 iteration # 2 ecut= 8.00 Ry beta=0.70 Davidson diagonalization with overlap ethr = 4.32E-10, avg # of iterations = 2.0 total cpu time spent up to now is 0.10 secs End of self-consistent calculation k = ( 113 PWs) bands (ev): ! total energy = Ry Harris-Foulkes estimate = Ry estimated scf accuracy < 2.8E-09 Ry convergence has been achieved in 2 iterations Forces acting on atoms (Ry/au): atom 1 type 1 force = atom 2 type 1 force = Total force = Total SCF correction = Entering Dynamics: iteration = 3 time = pico-seconds ATOMIC_POSITIONS (alat) Si Si kinetic energy (Ekin) = Ry temperature = K Ekin + Etot (const) = Ry Linear momentum : Writing output data file pwscf.save second order wave-functions extrapolation
18 References Most of the equations have been taken from D. Marx and J. Hutter, Ab initio molecular dynamics, Cambridge University Press, 2009 Part of the material has been adapted from P. Giannozzi s tutorial on QE The original CP paper is R. Car and M. Parrinello, Phys. Rev. Lett. 55, 2471 (1985).
Ab initio molecular dynamics. Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy. Bangalore, 04 September 2014
Ab initio molecular dynamics Simone Piccinin CNR-IOM DEMOCRITOS Trieste, Italy Bangalore, 04 September 2014 What is MD? 1) Liquid 4) Dye/TiO2/electrolyte 2) Liquids 3) Solvated protein 5) Solid to liquid
More informationGWL tutorial. Paolo Umari, Università degli Studi di Padova, Italy Democritos, Trieste
GWL tutorial Paolo Umari, Università degli Studi di Padova, Italy Democritos, Trieste GW calculation with QE and GWL Benfits: Optimal basis (reduced) for representing polarizabilty operators Full convergence
More informationQuantum ESPRESSO. Input and Output description
Quantum ESPRESSO Input and Output description Where can I find useful information about Quantum ESPRESSO? Where can I find useful information about Quantum ESPRESSO? prompt > cd $espresso_dir/doc; ls *.html
More informationExamples of Defects in Solids from ESPRESSO
Examples of Defects in Solids from ESPRESSO Alessandra Satta SLACS-INFMCNR Sardinian LAboratory for Computational Materials Science Hands-on Tutorial on the Quantum-ESPRESSO Package Cagliari, September
More informationMetals Magnetic systems
Metals Magnetic systems Ralph Gebauer Cape Town, July 2008 Metallic systems: k-points and smearing: Let us consider iron, in the bcc phase. In the first part of this exercise, we neglect magnetism and
More informationPWSCF First examples
PWSCF First examples (much more in espresso 3.1.1/examples directory!) Guido Fratesi (Università di Milano) Urbana, August 2006 A pw.x input file &CONTROL / &SYSTEM / &ELECTRONS title calculation restart_mode
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 informationStructural optimizations
Structural optimizations Guido Fratesi (Università di Milano) Urbana, August 2006 Acetylene molecule ( ) We want to use pw.x to find the optimized geometry of acetylene. Let us suppose we have the following
More informationThe quasi-harmonic approximation (QHA)
The quasi-harmonic approximation (QHA) M. Palumbo 19/01/2017 Trieste, Italy Limitations of the harmonic approximation E tot(r I, u I )=E tot(r I )+ I,α E tot u u Iα + 1 2 E tot Iα 2 I,α u u Iα u Iα u Jβ
More informationDFT calculation of pressure induced phase transition in silicon
DFT calculation of pressure induced phase transition in silicon Michael Scherbela scherbela@student.tugraz.at 2015-12-07 Contents 1 Introduction 2 2 Structure of the input file 2 3 Calculation of phase
More informationEnergy and Forces in DFT
Energy and Forces in DFT Total Energy as a function of nuclear positions {R} E tot ({R}) = E DF T ({R}) + E II ({R}) (1) where E DF T ({R}) = DFT energy calculated for the ground-state density charge-density
More informationAb initio molecular dynamics: Basic Theory and Advanced Met
Ab initio molecular dynamics: Basic Theory and Advanced Methods Uni Mainz October 30, 2016 Bio-inspired catalyst for hydrogen production Ab-initio MD simulations are used to learn how the active site
More information23 The Born-Oppenheimer approximation, the Many Electron Hamiltonian and the molecular Schrödinger Equation M I
23 The Born-Oppenheimer approximation, the Many Electron Hamiltonian and the molecular Schrödinger Equation 1. Now we will write down the Hamiltonian for a molecular system comprising N nuclei and n electrons.
More informationModule 2: Quantum Espresso Walkthrough
Module 2: Quantum Espresso Walkthrough Energy and Geometry Optimization of the H 2 Molecule We will be using the PWSCF code for quantum mechanical calculations of extended systems. The PWSCF program is
More informationLab 3: Handout Quantum-ESPRESSO: a first principles code, part 2.
1 Lab 3: Handout Quantum-ESPRESSO: a first principles code, part 2. In this lab, we will be using Quantum-ESPRESSO as our first-principles code again. In problem 1, we will compare energy between allotropes
More informationThe Energetics of the Hydrogenation of a Single-Walled Carbon Nanotube. Janet Ryu Nicola Marzari May 13, J
The Energetics of the Hydrogenation of a Single-Walled Carbon Nanotube Janet Ryu Nicola Marzari May 13, 2005 3.021J Janet Ryu May 13, 2005 Final Paper 3.021J Carbon Nanotubes The Energetics of the Hydrogenation
More informationAb initio Molecular Dynamics Born Oppenheimer and beyond
Ab initio Molecular Dynamics Born Oppenheimer and beyond Reminder, reliability of MD MD trajectories are chaotic (exponential divergence with respect to initial conditions), BUT... With a good integrator
More informationCar-Parrinello Molecular Dynamics
Car-Parrinello Molecular Dynamics Eric J. Bylaska HPCC Group Moral: A man dreams of a miracle and wakes up with loaves of bread Erich Maria Remarque Molecular Dynamics Loop (1) Compute Forces on atoms,
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 informationConical Intersections. Spiridoula Matsika
Conical Intersections Spiridoula Matsika The Born-Oppenheimer approximation Energy TS Nuclear coordinate R ν The study of chemical systems is based on the separation of nuclear and electronic motion The
More informationSolid 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 informationIntroduction to the QUANTUM ESPRESSO package and its application to computational catalysis Input/Output description
School on computational materials modeling in catalysis Bangalore, 1-5 September 2014 Introduction to the QUANTUM ESPRESSO package and its application to computational catalysis Input/Output description
More informationProblem Set 2: First-Principles Energy Methods
Problem Set 2: First-Principles Energy Methods Problem 1 (10 points): Convergence of absolute energies with respect to cutoff energies. A Using the Quantum ESPRESSO PWscf package, calculate the energy
More informationComputing NMR parameters using the GIPAW method
Computing NMR parameters using the GIPAW method DFT and NMR with Quantum Espresso (QE) thibault.charpentier@cea.fr & Ari.P.Seitsonen@iki.fi Welcome to the hands-on session on the GIPAW method. The idea
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 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 informationThe Overhauser Instability
The Overhauser Instability Zoltán Radnai and Richard Needs TCM Group ESDG Talk 14th February 2007 Typeset by FoilTEX Introduction Hartree-Fock theory and Homogeneous Electron Gas Noncollinear spins and
More informationSome Review & Introduction to Solar PV
1.021, 3.021, 10.333, 22.00 : Introduction to Modeling and Simulation : Spring 2012 Part II Quantum Mechanical Methods : Lecture 9 Some Review & Introduction to Solar PV Jeffrey C. Grossman Department
More informationCHEM3023: Spins, Atoms and Molecules
CHEM3023: Spins, Atoms and Molecules Lecture 3 The Born-Oppenheimer approximation C.-K. Skylaris Learning outcomes Separate molecular Hamiltonians to electronic and nuclear parts according to the Born-Oppenheimer
More informationLecture 9: Molecular Orbital theory for hydrogen molecule ion
Lecture 9: Molecular Orbital theory for hydrogen molecule ion Molecular Orbital Theory for Hydrogen Molecule Ion We have seen that the Schrödinger equation cannot be solved for many electron systems. The
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 informationDiatomic Molecules. 7th May Hydrogen Molecule: Born-Oppenheimer Approximation
Diatomic Molecules 7th May 2009 1 Hydrogen Molecule: Born-Oppenheimer Approximation In this discussion, we consider the formulation of the Schrodinger equation for diatomic molecules; this can be extended
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 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 informationThe electronic structure of materials 1
Quantum mechanics 2 - Lecture 9 December 18, 2013 1 An overview 2 Literature Contents 1 An overview 2 Literature Electronic ground state Ground state cohesive energy equilibrium crystal structure phase
More informationAccelerated Quantum Molecular Dynamics
Accelerated Quantum Molecular Dynamics Enrique Martinez, Christian Negre, Marc J. Cawkwell, Danny Perez, Arthur F. Voter and Anders M. N. Niklasson Outline Quantum MD Current approaches Challenges Extended
More informationDensity Functional Theory for Electrons in Materials
Density Functional Theory for Electrons in Materials Richard M. Martin Department of Physics and Materials Research Laboratory University of Illinois at Urbana-Champaign 1 Density Functional Theory for
More informationIV. Classical Molecular Dynamics
IV. Classical Molecular Dynamics Basic Assumptions: 1. Born-Oppenheimer Approximation 2. Classical mechanical nuclear motion Unavoidable Additional Approximations: 1. Approximate potential energy surface
More informationAb initio molecular dynamics
Ab initio molecular dynamics Molecular dynamics Why? allows realistic simulation of equilibrium and transport properties in Nature ensemble averages can be used for statistical mechanics time evolution
More informationAb initio molecular dynamics and nuclear quantum effects
Ab initio molecular dynamics and nuclear quantum effects Luca M. Ghiringhelli Fritz Haber Institute Hands on workshop density functional theory and beyond: First principles simulations of molecules and
More informationCourse SS18. Organization
0-0 Course SS18 Stefan.Goedecker@unibas.ch Organization Lecture with script: http://comphys.unibas.ch/teaching.htm Simple exercises: Traditional analytic problems and numerical problems (on your own laptop?)
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 informationFirst-principles modeling: The evolution of the field from Walter Kohn s seminal work to today s computer-aided materials design
First-principles modeling: The evolution of the field from Walter Kohn s seminal work to today s computer-aided materials design Peter Kratzer 5/2/2018 Peter Kratzer Abeokuta School 5/2/2018 1 / 34 Outline
More informationab initio Electronic Structure Calculations
ab initio Electronic Structure Calculations New scalability frontiers using the BG/L Supercomputer C. Bekas, A. Curioni and W. Andreoni IBM, Zurich Research Laboratory Rueschlikon 8803, Switzerland ab
More informationTime reversible Born Oppenheimer molecular dynamics
Time reversible Born Oppenheimer molecular dynamics Jianfeng Lu Mathematics Department Department of Physics Duke University jianfeng@math.duke.edu KI-Net Conference, CSCAMM, University of Maryland, May
More informationDFT+U practical session
DFT+U practical session Matteo Cococcioni GGA and GGA+U calculations in FeO Calculation of U for bulk Fe Calculation of U for NiO Exercise I: evaluating U for Cu 2 O Exercise II: evaluating U for FePO
More informationQuantum ESPRESSO. PWSCF: first steps
Quantum ESPRESSO PWSCF: first steps What can I learn in this tutorial? What can I learn in this tutorial? How to run PWscf (pw.x) in self-consistent mode for Silicon How to get the band structure of Silicon
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 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 informationQuick reference guide on PLUMED with Quantum ESPRESSO
Quick reference guide on PLUMED with Quantum ESPRESSO Changru Ma SISSA, Trieste March 30, 2011 Contents 1 Introduction 2 1.1 Overview................................... 2 1.2 Collective variables.............................
More informationFast Eigenvalue Solutions
Fast Eigenvalue Solutions Techniques! Steepest Descent/Conjugate Gradient! Davidson/Lanczos! Carr-Parrinello PDF Files will be available Where HF/DFT calculations spend time Guess ρ Form H Diagonalize
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 informationFirst-principles Molecular Dynamics Simulations
First-principles Molecular Dynamics Simulations François Gygi University of California, Davis fgygi@ucdavis.edu http://eslab.ucdavis.edu http://www.quantum-simulation.org MICCoM Computational School, Jul
More informationRecent advances in development of single-point orbital-free kinetic energy functionals
PacifiChem 2010 p. 1/29 Recent advances in development of single-point orbital-free kinetic energy functionals Valentin V. Karasiev vkarasev@qtp.ufl.edu Quantum Theory Project, Departments of Physics and
More informationSet the initial conditions r i. Update neighborlist. Get new forces F i
v Set the initial conditions r i ( t 0 ), v i ( t 0 ) Update neighborlist Quantum mechanical models Get new forces F i ( r i ) Solve the equations of motion numerically over time step Δt : r i ( t n )
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 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 informationBorn-Oppenheimer Approximation
Born-Oppenheimer Approximation Adiabatic Assumption: Nuclei move so much more slowly than electron that the electrons that the electrons are assumed to be obtained if the nuclear kinetic energy is ignored,
More informationMolecular Dynamics. Park City June 2005 Tully
Molecular Dynamics John Lance Natasa Vinod Xiaosong Dufie Priya Sharani Hongzhi Group: August, 2004 Prelude: Classical Mechanics Newton s equations: F = ma = mq = p Force is the gradient of the potential:
More information(1/2) M α 2 α, ˆTe = i. 1 r i r j, ˆV NN = α>β
Chemistry 26 Spectroscopy Week # The Born-Oppenheimer Approximation, H + 2. Born-Oppenheimer approximation As for atoms, all information about a molecule is contained in the wave function Ψ, which is the
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 informationK-Points and Metals. Ralph Gebauer ICTP, Trieste. (slides courtesy of Shobhana Narasimhan)
Joint ICTP-TWAS Caribbean School on Electronic Structure Fundamentals and Methodologies (an Ab-initio Perspective) Cartagena Colombia, 27.08. to 21.09.2012 K-Points and Metals Ralph Gebauer ICTP, Trieste
More informationAb Initio Molecular Dynamics: Theory and Implementation
John von Neumann nstitute for Computing Ab nitio Molecular Dynamics: Theory and mplementation Dominik Marx and Jürg Hutter published in Modern Methods and Algorithms of Quantum Chemistry, Proceedings,
More informationv(r i r j ) = h(r i )+ 1 N
Chapter 1 Hartree-Fock Theory 1.1 Formalism For N electrons in an external potential V ext (r), the many-electron Hamiltonian can be written as follows: N H = [ p i i=1 m +V ext(r i )]+ 1 N N v(r i r j
More informationIntroduction to Hartree-Fock Molecular Orbital Theory
Introduction to Hartree-Fock Molecular Orbital Theory C. David Sherrill School of Chemistry and Biochemistry Georgia Institute of Technology Origins of Mathematical Modeling in Chemistry Plato (ca. 428-347
More informationExchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride. Dimer. Philip Straughn
Exchange Correlation Functional Investigation of RT-TDDFT on a Sodium Chloride Dimer Philip Straughn Abstract Charge transfer between Na and Cl ions is an important problem in physical chemistry. However,
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 informationIntro to ab initio methods
Lecture 2 Part A Intro to ab initio methods Recommended reading: Leach, Chapters 2 & 3 for QM methods For more QM methods: Essentials of Computational Chemistry by C.J. Cramer, Wiley (2002) 1 ab initio
More informationAb-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University
Ab-initio molecular dynamics: from the basics up to quantum effects Roberto Car Princeton University Hands-on Tutorial Workshop on Ab-Initio Molecular Simulations, Fritz- Haber-Institut, Berlin, July 12-21,
More informationDensity Functional Theory: from theory to Applications
Density Functional Theory: from theory to Applications Uni Mainz May 27, 2012 Large barrier-activated processes time-dependent bias potential extended Lagrangian formalism Basic idea: during the MD dynamics
More informationDay 1 : Introduction to AIMD and PIMD
Day 1 : Introduction to AIMD and PIMD Aug 2018 In today s exercise, we perform ab-initio molecular dynamics (AIMD) and path integral molecular dynamics (PIMD) using CP2K[?]. We will use the Zundel s cation
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 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 informationModels for Time-Dependent Phenomena
Models for Time-Dependent Phenomena I. Phenomena in laser-matter interaction: atoms II. Phenomena in laser-matter interaction: molecules III. Model systems and TDDFT Manfred Lein p.1 Outline Phenomena
More informationAb initio molecular dynamics: Propagating the density matrix with Gaussian orbitals
JOURNAL OF CHEMICAL PHYSICS VOLUME 114, NUMBER 8 JUNE 001 Ab initio molecular dynamics: Propagating the density matrix with Gaussian orbitals H. Bernhard Schlegel and John M. Millam Department of Chemistry,
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 informationBasic tutorial to CPMD calculations
Car and Parinnello Molecular Dynamics http://www.cpmd.org/ Basic tutorial to CPMD calculations Sébastien LE ROUX sebastien.leroux@ipcms.unistra.fr INSTITUT DE PHYSIQUE ET DE CHIMIE DES MATÉRIAUX DE STRASBOURG,
More informationIs there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich
Is there a future for quantum chemistry on supercomputers? Jürg Hutter Physical-Chemistry Institute, University of Zurich Chemistry Chemistry is the science of atomic matter, especially its chemical reactions,
More informationAb initio molecular dynamics
Ab initio molecular dynamics Kari Laasonen, Physical Chemistry, Aalto University, Espoo, Finland (Atte Sillanpää, Jaakko Saukkoriipi, Giorgio Lanzani, University of Oulu) Computational chemistry is a field
More informationLarge Scale Electronic Structure Calculations
Large Scale Electronic Structure Calculations Jürg Hutter University of Zurich 8. September, 2008 / Speedup08 CP2K Program System GNU General Public License Community Developers Platform on "Berlios" (cp2k.berlios.de)
More informationMulti-Scale Modeling from First Principles
m mm Multi-Scale Modeling from First Principles μm nm m mm μm nm space space Predictive modeling and simulations must address all time and Continuum Equations, densityfunctional space scales Rate Equations
More informationInteratomic Potentials. The electronic-structure problem
Interatomic Potentials Before we can start a simulation, we need the model! Interaction between atoms and molecules is determined by quantum mechanics: Schrödinger Equation + Born-Oppenheimer approximation
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 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 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 informationFast and accurate Coulomb calculation with Gaussian functions
Fast and accurate Coulomb calculation with Gaussian functions László Füsti-Molnár and Jing Kong Q-CHEM Inc., Pittsburgh, Pennysylvania 15213 THE JOURNAL OF CHEMICAL PHYSICS 122, 074108 2005 Received 8
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 informationQuantum Final Project by Anustup Poddar and Cody Tripp 12/10/2013
Quantum Final Project by Anustup Poddar and Cody Tripp 12102013 Introduction The Hamiltonian in the Schrӧdinger equation is the sum of a kinetic and potential energy operator. The Fourier grid Hamiltonian
More informationLecture 10. Central potential
Lecture 10 Central potential 89 90 LECTURE 10. CENTRAL POTENTIAL 10.1 Introduction We are now ready to study a generic class of three-dimensional physical systems. They are the systems that have a central
More informationLecture 8: Introduction to Density Functional Theory
Lecture 8: Introduction to Density Functional Theory Marie Curie Tutorial Series: Modeling Biomolecules December 6-11, 2004 Mark Tuckerman Dept. of Chemistry and Courant Institute of Mathematical Science
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 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 informationLecture 6. Tight-binding model
Lecture 6 Tight-binding model In Lecture 3 we discussed the Krönig-Penny model and we have seen that, depending on the strength of the potential barrier inside the unit cell, the electrons can behave like
More informationTime-dependent density functional theory (TDDFT)
Advanced Workshop on High-Performance & High-Throughput Materials Simulations using Quantum ESPRESSO ICTP, Trieste, Italy, January 16 to 27, 2017 Time-dependent density functional theory (TDDFT) Ralph
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 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 informationFaddeev Random Phase Approximation (FRPA) Application to Molecules
Faddeev Random Phase Approximation (FRPA) Application to Molecules Matthias Degroote Center for Molecular Modeling (CMM) Ghent University INT 2011 Spring Program Fermions from Cold Atoms to Neutron Stars:
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 informationIntroduction to Density Functional Theory
Introduction to Density Functional Theory S. Sharma Institut für Physik Karl-Franzens-Universität Graz, Austria 19th October 2005 Synopsis Motivation 1 Motivation : where can one use DFT 2 : 1 Elementary
More informationAb Ini'o Molecular Dynamics (MD) Simula?ons
Ab Ini'o Molecular Dynamics (MD) Simula?ons Rick Remsing ICMS, CCDM, Temple University, Philadelphia, PA What are Molecular Dynamics (MD) Simulations? Technique to compute statistical and transport properties
More information