Molecular structure and properties calculations


Disadvantages of MM methods



Download 2.13 Mb.
Page2/3
Date09.06.2018
Size2.13 Mb.
#53901
1   2   3

Disadvantages of MM methods:

  • each force field provides good results for a limited class of molecules (for wich it was parametrized). No force field is general!

  • since empirical methods do not consider the electrons, they can not describe bond formation and breaking

  • many molecular properties depending on subtle electronic effects are not reproducible by MM methods.


See AMBER parameter set!


b) molecular dynamics (MD methods)



    • consists in following the temporal evolution of a system of interacting atoms

      • by integrating their equations of motions: Fi=miai

      • being given a set of positions and velocities => evolution in time is completely determined

  • forces acting on atoms are derived from the potential:


The simplest choice for V is to write it as a sum of pairwise interactions:  

The potential functions representing the non-bonded energy are formulated as a sum over interactions between the particles of the system. The simplest choice is the "pair potential", in which the total potential energy can be calculated from the sum of energy contributions between pairs of atoms.



Lennard-Jones potential:

The intramolecular bonding interactions must be also considered:



      • stretching energy:

      • bending energy:

      • torsion energy:

Force field are also nedeed in this case (CHARMM)

Advantages of MD methods

  • can treat very large molecular systems (thousands of atoms)

  • used in biophysics (proteins, viruses)


Disadvantages of MD methods

  • not able to model bond forming/breaking since electronic structure methods does not enter these models


2. QUANTUM MECHANICAL APROACHES

  1. semiempirical methods (AM1, PM3, PPP, INDO, MINDO, ...)

  • approximate methods in which some quantities (electron repulsion integrals) are taken from experiment, some small quantities are neglected and some quantities are estimated by fitting to experimental data; empirical parameters and functions are used to compensate for errors introduced by removing the integrals

  • model only the valence electrons

  • limited to hundred of atoms

  • can be used to study ground and excited molecular states

  • parametrized to reproduce experimental reference data or accurate high-level theoretical predictions

Ab Initio wave function



  1. non empirical methods

DFT electronic density

(quantum mechanics methods)

(electronic structure methods)


  • do not require empirical parameters

  • can be used for any molecular system

  • limited to tens of atoms

  • can be used to study ground and excited molecular states

Quantum mechanics methods are based on the following principles:



  • nuclei and electrons are distinguished from each other

  • electron-electron and electron-nuclear interactions are considered explicitly.

  • interactions are governed by nuclear and electron charges (potential energy) and electron motions

  • interactions determine the spatial distribution of nuclei and electrons and their energies


What level of theory?

1. The choice of level of theory to use is usually determined by the property being calculated and the cost of the calculation in terms of computer usage.

2. In general the higher the level of theory the greater the computer cost. For example Hartree-Fock based ab-initio methods are more costly than semiempirical methods. Correlated methods, depending on the sophistication, can add additional cost to the Hartree-Fock method.

3. For the Hartree-Fock and Density Functional methods an additional factor is the quality (size) of the basis sets used. In general the larger the number of functions used in the basis set the more flexible (better) the basis set. Each additional function adds extra cost to the calculation.



T
ype of calculations



  1. Single point calculations

energy, wave-function and other requested properties at a single fixed geometry

  • usually done first at the beginning of a study on a new molecule to check out the nature of the wave-function

  • it is also frequently carried out after a geometry optimization, with a larger basis set or more superior method than is possible with the basis set and method used to optimize geometry


These calculations are performed:

  • to obtain basic information about a given molcule

  • to check a molecular geometry to be used as a starting geometry in an optimization process

  • to compute more accurate energies and other properties for an optimized geometry at a lower theoretical level.





  1. Geometry optimization calculations

►calculation of the wave-function and energy at a starting geometry and then searching for the geometry which will give the lowest energy (local or global energy minimum)
Potential energy surface (PES)

  • specify the way in which the energy of a molecular system varies with small changes in its structure

  • mathematical relationship linking the molecular structure and the resultant energy

local minimum – the lowest energy point in a limited region of PES

global minimum – the lowest energy point anywhere on the PES
saddle points  maximum in one direction and minimum in all other directions

correspond to transition state structures



all minima + saddle points = stationary points ()
Geometry optimizations usually locate the stationary points closest to the starting geometry!!!

Download 2.13 Mb.

Share with your friends:
1   2   3




The database is protected by copyright ©ininet.org 2024
send message

    Main page