17.11.2014 Views

ChemOffice.Com - CambridgeSoft

ChemOffice.Com - CambridgeSoft

ChemOffice.Com - CambridgeSoft

SHOW MORE
SHOW LESS

You also want an ePaper? Increase the reach of your titles

YUMPU automatically turns print PDFs into web optimized ePapers that Google loves.

Administrator<br />

6. The modified values of K sres and BL res are<br />

used in the molecular mechanics portion of the<br />

MM2 computation to further refine the<br />

molecule.<br />

Stretch-Bend Cross Terms<br />

Stretch-bend cross terms are used when a coupling<br />

occurs between bond stretching and angle bending.<br />

For example, when an angle is compressed, the<br />

MM2 force field uses the stretch-bend force<br />

constants to lengthen the bonds from the central<br />

atom in the angle to the other two atoms in the<br />

angle.<br />

1<br />

Ε= ∑<br />

2 K sb( r − r o<br />

)θ − θ o<br />

Stretch / Bend<br />

( )<br />

The force constant (K sb ) differs for different atom<br />

combinations.<br />

The seven different atom combinations where<br />

force constants are available for describing the<br />

situation follow:<br />

• X-B, C, N, O-Y<br />

• B-B, C, N, O-H<br />

• X-Al, S-Y<br />

• X-Al, S-H<br />

• X-Si, P-Y<br />

• X-Si, P-H<br />

• X-Ga, Ge, As, Se-Y, P-Y<br />

where X and Y are any non-hydrogen atom.<br />

User-Imposed Constraints<br />

Additional terms are included in the force field<br />

when constraints are applied to torsional angles and<br />

non-bonded distances by the Optimal field in the<br />

Measurements table. These terms use a harmonic<br />

potential function, where the force constant has<br />

been set to a large value (4 for torsional constraints<br />

and 10 6 for non-bonded distances) in order to<br />

enforce the constraint.<br />

For torsional constraints the additional term and<br />

force constant is described by:<br />

Ε= ∑ 4(θ − θ o<br />

) 2<br />

Torsions<br />

For non-bonded distance constraints the additional<br />

term and force constant is:<br />

Ε= ∑ 10 6<br />

(r − r o<br />

) 2<br />

Distance<br />

Molecular Dynamics<br />

Simulation<br />

In its broadest sense, molecular dynamics is<br />

concerned with simulating molecular motion.<br />

Motion is inherent to all chemical processes. Simple<br />

vibrations, like bond stretching and angle bending,<br />

give rise to IR spectra. Chemical reactions,<br />

hormone-receptor binding, and other complex<br />

processes are associated with many kinds of<br />

intramolecular and intermolecular motions. The<br />

MM2 method of molecular dynamics simulation<br />

uses Newton’s equations of motion to simulate the<br />

movement of atoms.<br />

Conformational transitions and local vibrations are<br />

the usual subjects of molecular dynamics studies.<br />

Molecular dynamics alters the values of the<br />

intramolecular degrees of freedom in a stepwise<br />

fashion. The steps in a molecular dynamics<br />

simulation represent the changes in atom position<br />

over time, for a given amount of kinetic energy.<br />

The driving force for chemical processes is<br />

described by thermodynamics. The mechanism by<br />

which chemical processes occur is described by<br />

kinetics. Thermodynamics describes the energetic<br />

relationships between different chemical states,<br />

whereas the sequence or rate of events that occur as<br />

molecules transform between their various possible<br />

states is described by kinetics.<br />

138•<strong>Com</strong>putation Concepts<br />

<strong>CambridgeSoft</strong><br />

Molecular Mechanics Theory in Brief

Hooray! Your file is uploaded and ready to be published.

Saved successfully!

Ooh no, something went wrong!