Force Fields and Simulations Flashcards

1
Q

For a diatomic molecule the potential energy only depends on

A

one variable - the bond length

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
2
Q

Plotting oscillator potential to approximate potential energy represents

A

how often energy varies as a function of the distortion of the bond length from its equilibrium value

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
3
Q

1 d potential energy (E) =

A

f(l) = 1/2 * k (l-l0)^2 for bond length l, min bond length l0 and bond force constant k

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
4
Q

If potential energy is a function of more than one variable we have a…

A

potential energy surface

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
5
Q

potential energy wells correspond to

A

lowest energy states for reactants and products

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
6
Q

Saddle point on PES correspond to

A

transition state for reaction

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
7
Q

hypersurface depends on

A

many variables

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
8
Q

molecules spend most of the time in

A

potential wells on PES

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
9
Q

energy minima on PES are

A

STABLE CONFORMATIONS

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
10
Q

Anti conformation (180) is more stable than gauche (60) because

A

reduced steric hindrance

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
11
Q

PES calculated by considering

A

rotation about backbone dihedral angles

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
12
Q

To calculate potential energy of molecular system write

A

as a sum of all the different contributions from bond stretching/bending, dihedral angle distortions and nonbonded interaction (eg steric repulsion and van der Waals) and electrostatic interaction

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
13
Q

Can use computational algorithms to minimise

A

steric energy and find stable conformations corresponding to minima on the PES

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
14
Q

High force constant for bond corresponds to

A

stiff bonds where distortion from equilibrium value corresponds to large energies

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
15
Q

Energy in torsional motion of dihedral angle =

A

sum up to number of terms (n) in series expansion of: (force constants / 2)*(1+cos(nangle))

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
16
Q

2 neutral atoms approaching at long distances

A

attract each other through dispersion interactions and repel at close distances

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
17
Q

Lennard jones potential

A

epsilon is well depth and sigma is distance where attractive and repulsive energies cancel equation on cheat sheet

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
18
Q

Energy due to electrostatic interactions can be given by

A

a charge-charge interaction equation on cheat sheet

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
19
Q

in some force fields electrostatic interactions can be represents by

A

bond dipoles. equation on cheat sheet

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
20
Q

Steric energy should not be compared between

A

molecules. Compared between different conformation of a molecule instead.

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
21
Q

Individual force fields are not normally needed because

A

so similarin different molecules so can use the same parameters

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
22
Q

making force fields:

A

using computational fitting procedures. using experimental data, using other calculations eg monte carlo (details on pg 10-11)

How well did you know this?
1
Not at all
2
3
4
5
Perfectly
23
Q

Structure with the lowest energy is the most stable and therefore the

A

most abundant structure.

24
Q

IF first derivative of energy is 0 and second derivative is > 0

25
Formal derivation of minimisation algorithm
start at x. x(min) = x + delta x. differential of x(min) = 0. expand as taylor series. equation. Find x(min) as roughly x - E'/E''.
26
For higher dimensions replace E' by
vector (force vector)
27
For higher dimensions replace E'' by
hessian matrix, giving a vector of small changes for each iteration.
28
Boltzmann Law gives
relative population of states with different energies.
29
If more than 2 states in system we can use boltzman distribution to find
probability of each state
30
partition function is the
sum over boltzmann factors for each energy level and energy of lowest energy state is subtracted from each energy level
31
To find probability of different conformations:
Find energy relative to lowest energy state for each conformation. Work our Boltzmann factor for each conformation. Sum up boltzmann factor to give partition function. Work out probability of each level using partition. Us relative energy to get correct partition function.
32
Uses of Molecular mechanics (3)
finding global energy minima, conformational searching and modelling of interactions between molecules.
33
Advantages of MM over QM: (2):
MM faster, possible for large molecules, simple organic molecule stuctures are better for MM as poor representation of dispresion in Qm
34
Disadvantages of MM over QM: (3)
No account of electrons, not good where electronic effects dominate over steric effects, and polarisability not taken into account in most force fields.
35
Limitations of MM(4):
in large molecules, the number of accessible conformations is great, the number of molecular arrangements grows exponentially as the number of molecules increases. if only considering energy minima then ignore entropy and free energy, molecules are moving all the time so sampling diff energy states.
36
Average energy =
average over representation number of states. equation on cheat sheet.
37
Velocity Verlet algorithm
breaks time down into a sequence of small discrete time-steps (delta t). method on page 18. delta t is very small and should be 1/20 of the fastest motion in system.
38
Leapfrog algorithm
velocities are half a time step behind the positions but then leap ahead by half a time step. Algorithm of pg 18
39
Monte carlo simulation
uses random numbers to make small changes from current configuration, calculating Blotzmann factor for energy change and acceptingrejecting the change. diagram on pg 19
40
MC vs MD (5):
Both generate energy states according to their Bolztmann weights. Only coordinates enter MC, md has velocities and forces. Md good for big molecules. More good.free programs available for MD.
41
LARGE molecules hard to simulate via MC because
hard to generate MC mobes that sample conformational space efficiently
42
Mc has many advantages over MD such as
use of special MC moves and use of thermodynamic ensembles for simulation. Even possible to have MC moves whree molecules appear/disappear. Move can be used to speed up simulation of mixing or demixing in a binary mixture of 2 liquids.
43
When simulating bulk systems with many molecules avoid
serious problem caused by fact we can only simulate a few hundred molecules (at most a few thousand).
44
For a small system surface molecules dominate over
bulk molecules
45
Periodic boundary conditions
employed alongside minimum image convention
46
when simulating in the bulk it is usually necessary to
conserve certain quantities that define the thermodynamic ensemble used. important ensembles on page 21
47
Thermodyanmic properties can be obtained using
MC and MD techniques. examples on page 22
48
increasing the temperature in MC will
change the moves that are accepted with additional moves to higher energy configurations being accepted
49
For MD temperature can be obtained from velocities in a simulation by the
equipartition principle equation on page 22
50
dihedral angle distributions plot the
frequency of dihedral angles occuring during a simulation.
51
Can take natural log of the dihedral angle distribution to obtain the
torsional free energy
52
Radial distribution function g(r) shows
the average structure in a liquid. 1st peak corresponds to solvation shell, 2nd solvation shell is less well defined. g(r) decays to 1 at long distances as no long range structure in a liquid
53
Bulk diffusion coefficient can be estimated in MD simulation from the
mean squared displacement at long times or the velocity autocorrection function equation on page 23
54
Molecular mechanics is for single molecules and
solids and 'clusters'
55
MC and MD for
single molecules and fluids