RTO care

You are here: Computation Concepts > Methods Overview > Using Computational Methods

Methods Overview

Computational chemistry encompasses a variety of mathematical methods that fall into two broad categories:

Molecular mechanics

applies the laws of classical physics to the atoms in a molecule without explicit consideration of electrons.

Quantum mechanics

relies on the Schrödinger equation to describe a molecule with explicit treatment of electronic structure.

Quantum mechanical methods can be subdivided into two classes: ab initio and semiempirical.

Chem & Bio 3D 11.0 provides the following methods:

Method

  

  

molecular mechanics

MM2

MM3, MM3-protein

AMBER,UFF, Dreiding

Chem3D, Tinker

Tinker

Gaussian

semi-empirical

Extended Hückel

other semi-empirical methods (AM1, MINDO/3, PM3, etc.)

Chem3D, MOPAC, Gaussian

MOPAC, Gaussian

Ab initio

RHF, UHF, MP2, etc.

Gaussian, GAMESS

•  Molecular mechanical methods: MM2 (directly). MM3 and MM3-protein through the Chem3D Tinker interface.

•  Semiempirical Extended Hückel, MINDO/3, MNDO, MNDO-d, AM1 and PM3 methods through Chem3D and CS MOPAC or Gaussian.

•  Ab initio methods through the Chem3D Gaussian or GAMESS interface.

Using Computational Methods

Computational methods calculate the potential energy surfaces (PES) of molecules. The PES is the embodiment of the forces of interaction among atoms in a molecule. From the PES, structural and chemical information about a molecule can be derived. The methods differ in the way the surface is calculated and in the molecular properties derived from the energy surface.

The methods perform the following basic calculations:

Single point energy calculation

The energy of a given geometry of the atoms in a model, which is the value of the PES at that point.

Geometry optimization

A systematic modification of the atomic coordinates of a model resulting in a geometry where the forces on each atom in the structure is zero. A 3dimensional arrangement of atoms in the model representing a local energy minimum (a stable molecular geometry to be found without crossing a conformational energy barrier).

Property calculation

Predicts certain physical and chemical properties, such as charge, dipole moment, and heat of formation.

Computational methods can perform more specialized functions, such as conformational searches and molecular dynamics simulations.

Choosing the Best Method

Not all types of calculations are possible for all methods and no one method is best for all purposes. For any given application, each method poses advantages and disadvantages. The choice of method depend on a number of factors, including:

•  The nature and size of the molecule

•  The type of information sought

•  The availability of applicable experimentally determined parameters (as required by some methods)

•  Computer resources

The three most important of the these criteria are:

Model size

The size of a model can be a limiting factor for a particular method. The limiting number of atoms in a molecule increases by approximately one order of magnitude between method classes from ab initio to molecular mechanics. Ab initio is limited to tens of atoms, semiempirical to hundreds, and molecular mechanics to thousands.

Parameter Availability

Some methods depend on experimentally determined parameters to perform computations. If the model contains atoms for which the parameters of a particular method have not been derived, that method may produce invalid predictions. Molecular mechanics, for example, relies on parameters to define a force-field. A force-field is only applicable to the limited class of molecules for which it is parametrized.

Computer resources

Requirements increase relative to the size of the model for each of the methods.

Ab initio: The time required for performing calculations increases on the order of N4, where N is the number of atoms in the model.

Semiempirical: The time required for computation increases as N3 or N2, where N is the number of atoms in the model.

MM2: The time required for performing computations increases as N2, where N is the number of atoms.

In general, molecular mechanical methods require less computer resources than quantum mechanical methods. The suitability of each general method for particular applications can be summarized as follows.

Molecular Mechanics Methods Applications Summary

Molecular mechanics in Chem3D apply to:

•  Systems containing thousands of atoms.

•  Organic, oligonucleotides, peptides, and saccharides.

•  Gas phase only (for MM2).

Useful techniques available using MM2 methods include:

•  Energy Minimization for locating stable conformations.

•  Single point energy calculations for comparing conformations of the same molecule.

•  Searching conformational space by varying one or two dihedral angles.

•  Studying molecular motion using Molecular Dynamics.

Quantum Mechanical Methods Applications Summary

Useful information determined by quantum mechanical methods includes:

•  Molecular orbital energies and coefficients.

•  Heat of Formation for evaluating conformational energies.

•  Partial atomic charges calculated from the molecular orbital coefficients.

•  Electrostatic potential.

•  Dipole moment.

•  Transition-state geometries and energies.

•  Bond dissociation energies.

Semiempirical methods available in Chem3D with CS MOPAC or Gaussian apply to:

Systems containing up to 120 heavy atoms and 300 total atoms.

Organic, organometallics, and small oligomers (peptide, nucleotide, saccharide).

Gas phase or implicit solvent environment.

Ground, transition, and excited states.

Ab initio methods available in Chem3D with Gaussian or Jaguar apply to:

•  Systems containing up to 150 atoms.

•  Organic, organometallics, and molecular fragments (catalytic components of an enzyme).

•  Gas or implicit solvent environment.

•  Study ground, transition, and excited states (certain methods).

Method Type

Advantages

Disadvantages

Best For

Molecular Mechanics (Gaussian)

Gaussian uses classical physics and relies on force-field with embedded empirical parameters

Least intensive computationally. Gaussian is fast and is useful with limited computer resources. It can be used for molecules as large as enzymes.

Particular force field applicable only for a limited class of molecules

Does not calculate electronic properties

Requires experimental data (or data from ab initio) for parameters

Large systems that consist of thousands of atoms and Systems or processes with no breaking or forming of bonds

Semiempirical (MOPAC, Gaussian)

These use quantum physics, experimentally derived empirical parameters, and extensive approximation.

Less demanding computationally than ab initio methods

Capable of calculating transition states and excited states

Requires experimental data (or data from ab initio) for parameters

Less rigorous than ab initio methods

Medium-sized systems that consist of hundreds of atoms. Also, systems involving electronic transitions.

ab initio (Gaussian, GAMESS)

These use quantum physics, are rigourously mathematical methods, and use no empirical parameters

Useful for a broad range of systems

Does not depend on experimental data

Capable of calculating transition states and excited states

Computationally intensive

Small systems that consist of only tens of atoms or systems involving electronic transitions.

Molecules or systems without available experimental data ("new" chemistry).

Systems requiring rigorous accuracy.

Comparison of Methods

Potential Energy Surfaces

A potential energy surface (PES) can describe:

•  A molecule or ensemble of molecules having constant atom composition (ethane, for example) or a system where a chemical reaction occurs.

•  Relative energies for conformations (eclipsed and staggered forms of ethane).

Potential energy surfaces can differentiate between:

•  Molecules having slightly different atomic composition (ethane and chloroethane).

•  Molecules with identical atomic composition but different bonding patterns, such as propylene and cyclopropane

•  Excited states and ground states of the same molecule.

Potential Energy Surfaces (PES)

The true representation of a model's potential energy surface is a multi-dimensional surface whose dimensionality increases with the number of atom coordinates. Since each atom has three independent variables (x, y, z coordinates), visualizing a surface for a many-atom model is impossible. However, you can generalize this problem by examining any two independent variables, such as the x and y coordinates of an atom.

The main areas of interest on a potential energy surface are the extrema as indicated by the arrows, are as follows:

Global minimum

The most stable conformation appears at the extremum where the energy is lowest. A molecule has only one global minimum.

Local minima

Additional low energy extrema. Minima are regions of the PES where a change in geometry in any direction yields a higher energy geometry.

Saddle point

A stationary point between two low energy extrema. A saddle point is defined as a point on the potential energy surface at which there is an increase in energy in all directions except one, and for which the slope (first derivative) of the surface is zero.

Note: At the energy minimum, the energy is not zero; the first derivative (gradient) of the energy with respect to geometry is zero.

All the minima on a potential energy surface of a molecule represent stable stationery points where the forces on each atom sums to zero. The global minimum represents the most stable conformation; the local minima, less stable conformations; and the saddle points represent transition conformations between minima.

Single Point Energy Calculations

Single point energy calculations can be used to calculate properties of specific geometry of a model. The values of these properties depend on where the model lies on the potential surface as follows:

•  A single point energy calculation at a global minimum provides information about the model in its most stable conformation.

•  A single point calculation at a local minimum provides information about the model in one of many stable conformations.

•  A single point calculation at a saddle point provides information about the transition state of the model.

•  A single point energy calculation at any other point on the potential energy surface provides information about that particular geometry, not a stable conformation or transition state.

Single point energy calculations can be performed before or after optimizing geometry.

Note: Do not compare values from different methods. Different methods rely on different assumptions about a given molecule, and the energies differ by an arbitrary offset.

Geometry Optimization

Geometry optimization is used to locate a stable conformation of a model, and should be done before performing additional computations or analyses of a model.

Locating global and local energy minima is typically done by energy minimization. Locating a saddle point is optimizing to a transition state.

The ability of a geometry optimization to converge to a minimum depends on the starting geometry, the potential energy function used, and the settings for a minimum acceptable gradient between steps (convergence criteria).

Geometry optimizations are iterative and begin at some starting geometry as follows:

1. The single point energy calculation is performed on the starting geometry.

2. The coordinates for some subset of atoms are changed and another single point energy calculation is performed to determine the energy of that new conformation.

3. The first or second derivative of the energy (depending on the method) with respect to the atomic coordinates determines how large and in what direction the next increment of geometry change should be.

4. The change is made.

5. Following the incremental change, the energy and energy derivatives are again determined and the process continues until convergence is achieved, at which point the minimization process terminates.

The following illustration shows some concepts of minimization. For simplicity, this plot shows a single independent variable plotted in two dimensions.

   

The starting geometry of the model determines which minimum is reached. For example, starting at (b), minimization results in geometry (a), which is the global minimum. Starting at (d) leads to geometry (f), which is a local minimum.The proximity to a minimum, but not a particular minimum, can be controlled by specifying a minimum gradient that should be reached. Geometry (f), rather than geometry (e), can be reached by decreasing the value of the gradient where the calculation ends.

In theory, if a convergence criterion (energy gradient) is too lax, a first-derivative minimization can result in a geometry that is near a saddle point. This occurs because the value of the energy gradient near a saddle point, as near a minimum, is very small. For example, at point (c), the derivative of the energy is 0, and as far as the minimizer is concerned, point (c) is a minimum. First derivative minimizers cannot, as a rule, cross saddle points to reach another minimum.

Note: If the saddle point is the extremum of interest, it is best to use a procedure that specifically locates a transition state, such as the CS MOPAC Pro Optimize To Transition State command.

You can take the following steps to ensure that a minimization has not resulted in a saddle point.

•  The geometry can be altered slightly and another minimization performed. The new starting geometry might result in either (a), or (f) in a case where the original one led to (c).

•  The Dihedral Driver can be employed to search the conformational space of the model. For more information, see Tutorial 5: The Dihedral Driver .

•  A molecular dynamics simulation can be run, which will allow small potential energy barriers to be crossed. After completing the molecular dynamics simulation, individual geometries can then be minimized and analyzed. For more information see MM2

You can calculate the following properties with the computational methods available through Chem3D using the PES:

•  Steric energy

•  Heat of formation

•  Dipole moment

•  Charge density

•  COSMO solvation in water

•  Electrostatic potential

•  Electron spin density

•  Hyperfine coupling constants

•  Atomic charges

•  Polarizability

•  Others, such as IR vibrational frequencies

   

'상태와 변화' 카테고리의 다른 글

Oxygen  (0) 2016.06.27
Molecular Mechanics  (0) 2016.06.27
Computational Chemistry  (1) 2016.06.27
오존공정의 기본 메카니즘  (0) 2016.06.27
POCP  (0) 2016.06.27
Ozone Chemistry  (0) 2016.06.27

Comment +1

Methane

상태와 변화2016.06.27 12:40

------------ Property Broker ------------

Model: mrthane.mol

   

ChemPropPro: Boiling Point = 0 Kelvin

ChemPropPro: Critical Pressure = 0 Bar

ChemPropPro: Critical Temperature = 0 Kelvin

ChemPropPro: Critical Volume = 0 cm^3/mol

ChemPropPro: Gibbs Free Energy = 0 kJ/mol

ChemPropPro: Heat Of Formation = 0 kJ/mol

ChemPropPro: Henry's Law Constant = -1.429

ChemPropPro: Ideal Gas Thermal Capacity = 0 J/(mol.K)

ChemPropPro: LogP = 1.09

ChemPropPro: Melting Point = 0 Kelvin

ChemPropPro: Mol Refractivity = 6.347 cm^3/mol

ChemPropPro: Vapor Pressure = 103322.998 Pascal

ChemPropPro: Water Solubility = 24.4 mg/L

ChemPropPro: Full Report:

   

************************************************************************

Data from database

************************************************************************

   

<Name of molecule><RE PUBLISHING CORP.,NY,(l989).

   

<**********

   

Log(p)........: 1.09

St..deviation.: 0.47

by Crippen's fragmentation: J.Chem.Inf.Comput.Sci.,27,21(1987).

   

Log(p)........: 1.09

St..deviation.: 0.49

by Viswanadhan's fragmentation: J.Chem.Inf.Comput.Sci.,29,163(1989).

   

   

Estimation using Broto's fragmentation method

C Non available value

Log(p) can't be calculated by Broto's method:

Eur.J.Med.Chem.- Chim.Theor.,19,71(1984).

   

************************************************************************

Estimation of Molar Refractivity

************************************************************************

   

MR............: 6.88 [cm.cm.cm/mol]

St..deviation.: 1.27

by Crippen's fragmentation: J.Chem.Inf.Comput.Sci.,27,21(1987).

   

MR............: 6.35 [cm.cm.cm/mol]

St..deviation.: 0.77

by Viswanadhan's fragmentation: J.Chem.Inf.Comput.Sci.,29,163(1989).

   

   

************************************************************************

Estimation of Henry's Constant (H)

************************************************************************

   

1. method: H = -1.229 log[unitless]

Estimation of mean error..: 0.0620

   

2. method: H = -1.260 log[unitless]

Estimation of mean error..: 0.200

   

ChemPropStd: Formal Charge = 0

ChemPropStd: Connolly Accessible Area = 147.937 Angstroms Squared

ChemPropStd: Connolly Molecular Area = 44.35 Angstroms Squared

ChemPropStd: Connolly Solvent Excluded Volume = 24.335 Angstroms Cubed

ChemPropStd: Exact Mass = 16.0313001284 g/Mol

ChemPropStd: Mass = 16.04246

ChemPropStd: Mol Weight = 16.04246

ChemPropStd: Ovality = 1.09208204264751

ChemPropStd: Principal Moment = 3.303 3.303 3.303

ChemPropStd: Elemental Analysis = C, 74.87; H, 25.13

ChemPropStd: m/z = 16.03 (100.0%), 17.03 (1.1%)

ChemPropStd: Mol Formula = CH4

ChemPropStd: Mol Formula HTML = CH<sub>4</sub>

CLogP Driver: Mol Refractivity = 0.641199946403503

CLogP Driver: Partition Coefficient = 1.1029999256134

GAMESS Interface: Cp = 8.194 cal/(mol K)

GAMESS Interface: Cv = 6.207 cal/(mol K)

GAMESS Interface: Enthalpy = 31.086 Kcal/Mol

GAMESS Interface: Entropy = 49.479 cal/(mol K)

GAMESS Interface: Gibbs Free Energy = 16.334 Kcal/Mol

GAMESS Interface: Internal Energy = 30.494 Kcal/Mol

GAMESS Interface: Dipole = 0 Debye

GAMESS Interface: Harmonic Zero Point Energy = 0 Kcal/Mol

GAMESS Interface: Kinetic Energy = 24925.5722 Kcal/Mol

GAMESS Interface: Potential Energy = -50009.9336 Kcal/Mol

GAMESS Interface: Total Energy = -25084.3614 Kcal/Mol

Charges (Lowdin Charges)-GAMESS Interface:

C(1) -0.378865

H(2) 0.094718

H(3) 0.094714

H(4) 0.094714

H(5) 0.094718

GAMESS Interface: Lowdin Populations (Lowdin Populations)-GAMESS Interface:

C(1) 6.378865

H(2) 0.905282

H(3) 0.905286

H(4) 0.905286

H(5) 0.905282

Charges (Mulliken Charges)-GAMESS Interface:

C(1) -0.773602

H(2) 0.193402

H(3) 0.193399

H(4) 0.193399

H(5) 0.193402

GAMESS Interface: Mulliken Populations (Mulliken Populations)-GAMESS Interface:

C(1) 6.773602

H(2) 0.806598

H(3) 0.806601

H(4) 0.806601

H(5) 0.806598

GAMESS Interface: Frequencies = 0 0.01 502.48 502.54 502.58 1558.71 1558.72 1558.76 1746.61 1746.64 2925.64 2995.23 2995.24 2995.27

GAMESS Interface: Polarizibility:

AXIAL COMPONENTS WITH BASE FIELD OF 0.0010

ENERGY-BASED DIPOLE-BASED

DIPOLE Z: 1.2028308E-05 1.2028008E-05

ALPHA XZ: 3.4799825E-04

ALPHA YZ: -2.4807436E-04

ALPHA ZZ: 1.2137664E+01 1.2137767E+01

BETA XZZ: 2.5549004E+01

BETA YZZ: 1.6897140E+01

BETA ZZZ: 8.9251628E-01 8.9302229E-01

GAM ZZZZ: 7.8841822E+01 8.3879350E+01

Molecular Topology: Balaban Index = 0

Molecular Topology: Cluster Count = 1

Molecular Topology: Molecular Topological Index = 0

Molecular Topology: Num Rotatable Bonds = 0 Bond(s)

Molecular Topology: Polar Surface Area = 0 Angstroms Squared

Molecular Topology: Radius = 0 Atom(s)

Molecular Topology: Shape Attribute = 0

Molecular Topology: Shape Coefficient = 0

Molecular Topology: Sum Of Degrees = 0

Molecular Topology: Sum Of Valence Degrees = 0

Molecular Topology: Topological Diameter = 0 Bond(s)

Molecular Topology: Total Connectivity = 1

Molecular Topology: Total Valence Connectivity = 1

Molecular Topology: Wiener Index = 0

-----------------------------------------

------------MM2 Minimization------------

Note: All parameters used are finalized (Quality = 4).

Iteration 2: Minimization terminated normally because the gradient norm is less than the minimum gradient norm

Stretch: 0.0000

Bend: 0.0000

Stretch-Bend: 0.0000

Torsion: 0.0000

Non-1,4 VDW: 0.0000

1,4 VDW: 0.0000

Total Energy: 0.0000 kcal/mol

Calculation completed

------------------------------------

'상태와 변화' 카테고리의 다른 글

Performance catalytic Ozonation  (0) 2016.06.27
Ozone Reaction  (0) 2016.06.27
CO/CO2 IR spectra  (0) 2016.06.27
Methane  (0) 2016.06.27
엔탈피 엔트로피  (0) 2016.06.24
질의 문답  (0) 2016.06.24

Comment +0