03964.com

文档资料库 文档搜索专家

文档资料库 文档搜索专家

(a brief introduction)

Molecular Dynamics Simulation

Shaolong Zhang

College of Physics and Electronic Shandong Normal University Jinan Shandong, 250014 China

A Digital Laboratory

“In the real world, this could eventually mean that most chemical experiments are conducted inside the silicon of chips instead of the glassware of laboratories. Turn off that Bunsen burner; it will not be wanted in ten years.”

- The Economist, reporting on the work of the 1998 Chemistry Nobel Prize Awardees

A Different Type of Simulation

?

Many Physically-Based Simulations model easily observable real world phenomena.

Molecular Dynamics Simulations model things too small for us to observe directly.

?

Why Not Quantum Mechanics?

?

Modeling the motion of a complex molecule by solving the wave functions of the various subatomic particles would be accurate…

??

2

? ? ? U ( x, y, z )? ( x, y, z ) ? E? ( x, y, z )

2

2m

?

But it would also be very hard to program and take more computing power than anyone has!

Classical Mechanics

Instead of using Quantum mechanics, we can use classical Newtonian mechanics to model our system. ? This is a simplification of what is actually going on, and is therefore less accurate. ? To alleviate this problem, we use numbers derived from QM for the constants in our classical equations.

?

Molecular Modeling

For each atom in every molecule, we need: ? Position (r) ? Momentum (m + v) ? Charge (q) ? Bond information (which atoms, bond angles, etc.)

From Potential to Movement

To run the simulation, we need the force on each particle.

Fi ? m i a i

We use the gradient of the potential energy function. Now we can find the acceleration.

What is the Potential?

A single atom will be affected by the potential energy functions of every atom in the system:

? ?

Bonded Neighbors Non-Bonded Atoms (either other atoms in the same molecule, or atoms from different molecules)

V ( R ) ? E bonded ? E non ? bonded

Non-Bonded Atoms

There are two potential functions we need to be concerned about between nonbonded atoms:

E non ? bonded ? E van ? der ? Waals ? E electrosta

tic

van der Waals Potential ? Electrostatic Potential

?

E non ? bonded ? E van ? der ? Waals ? E electrosta

tic

The van der Waals Potential

Atoms with no net electrostatic charge will still tend to attract each other at short distances, as long as they don’t get too close. ? Once the atoms are close enough to have overlapping electron clouds, they will repel each other with astounding force

?

The van der Waals Potential

One of the most widely used functions for the van der Waals potential is the Lennard-Jones. It is a compromise between accuracy and computability.

E Lennard

? Jones

?

?

nonbonded pairs

? A ik C ik ? ? 6 ? r 12 rik ? ik

? ? ? ?

The van der Waals Potential

The Constants A and C depend on the atom types, and are derived from experimental data.

?

nonbonded pairs

? A ik C ik ? 12 ? 6 ?r rik ? ik

? ? ? ?

The Electrostatic Potential

Opposite Charges Attract ? Like Charges Repel ? The force of the attraction is inversely proportional to the square of the distance

?

E electrosta

tic

?

?

nonbonded pairs

qiqk Dr ik

Coulomb’s Law

F ?

q1 q 2 4 ?? 0 r

2

The Non-Bonded Potential

Combine the LJ and Electrostatic Potentials

E non ? bonded ? E van ? der ? Waals ? E electrosta

tic

Bonded Atoms

There are three types of interaction between bonded atoms:

Stretching along the bond ? Bending between bonds ? Rotating around bonds

?

E bonded ? E bond ? stretch ? E angle ? bend ? E rotate ? along ? bond

Bond Length Potentials

Both the spring constant and the ideal bond length are dependent on the atoms involved.

E bond ? stretch ?

?

K b (b ? b 0 )

2

1 , 2 pairs

Bond Angle Potentials

The spring constant and the ideal angle are also dependent on the chemical type of the atoms.

E bond ? bend ?

? K ? (?

angles

??0)

2

Torsional Potentials

Described by a dihedral angle and coefficient of symmetry (n=1,2,3), around the middle bond.

E ro ta te ? a lo n g ? b o n d ?

?

1, 4 p a irs

K ? (1 ? co s( n? ))

Integration Algorithms

Forces like the LJ potential have powers of 12, which would make Euler horribly unstable (even worse than usual) ? RK and Midpoint algorithms would seem to help ? However, force calculations are extremely expensive, and RK and Midpoint require multiple force calculations per timestep

?

Integration Algorithms

Using RK is justifiable in other circumstances, because it allows you to take larger timesteps (two force calculations allowing you to more than double the timestep) ? This is normally not achievable in MD simulations, because the forces are very rapidly changing non-linear functions. ? We need an algorithm with the stability benefits of RK without the cost of extra force calculations!

?

Verlet Algorithm

?

First, take a third-order Taylor step:

r (t ? ? t ) ? r (t ) ? v (t ) ? t ? 1 2 a (t ) ? t ?

2

1 3!

?

r (t ) ? t ? O ( ? t )

3 4

?

Now take a step backward:

r (t ? ? t ) ? r (t ) ? v (t ) ? t ? 1 2 a (t ) ? t ?

2

1 3!

?

r (t ) ? t ? O ( ? t )

3 4

Verlet Algorithm

?

When adding the two formulas, the first and third derivatives cancel out:

r (t ? ? t ) ? r (t ? ? t ) ? 2 r (t ) ? a (t ) ? t ? O ( ? t )

2 4

?

And we can express the next timestep in terms of the previous position and the current acceleration:

r (t ? ? t ) ? 2 r (t ) ? r (t ? ? t ) ? a (t ) ? t ? O ( ? t )

2 4

Verlet Algorithm

Pros: ? Simple & Effective ? Low Memory & CPU Requirements (don’t need to store velocities or perform multiple force calculations) ? Time Reversible ? Very stable even with large numbers of interacting particles Cons: ? Not as accurate as RK ? We never calculate velocities! (when would we need them?)

Obtaining Velocities

?

We can estimate the velocities using a finite difference:

v (t ) ?

? ?

1 2?t

[ r ( t ? ? t ) ? r ( t ? ? t )] ? O ( ? t )

2

This has a second order error, while our algorithm has a fourth order error There are variations of the Verlet algorithm, such as the leapfrog algorithm, which seek to improve velocity estimations.

Molecules in Solution

?

?

?

In real situations, a molecule is rarely isolated. In biological systems, proteins, RNA, and DNA are immersed in a sea of water molecules To accurately portray the effect of the solvent molecules on a system, the solvent molecules must be free flowing How do we establish computational boundaries while keeping a realistic solvent simulation?

Periodic Boundary Conditions

Simulate a segment of molecules in a larger solution by having repeatable regions ? Potential calculations are run only on each atom’s closest counterpart in the 27 cubes ? When an atom moves off the edge, it reappears on the other side (like in asteroids)

?

Cutoff Methods

Ideally, every atom should interact with every other atom ? This creates a force calculation algorithm of quadratic order ? We may be able to ignore atoms at large distances from each other without suffering too much loss of accuracy

?

Cutoff Methods

Truncation – cuts off calculation at a predefined distance ? Shift – alters the entire function as to be zero at the cutoff distance ? Switch – begins tapering to zero as the function approaches the cutoff distance

?

Live Demo

Using NAMD (Not Another Molecular Dynamics Simulation) and VMD (Visual Molecular Dynamics) from the University of Illinois at Urbana Champaign

Resources

Books ? Tamar Schlick Molecular Modeling and Simulation: An Interdisciplinary Guide 2002 ? Alan Hinchliffe Molecular Modelling for Beginners 2003 ? D. C. Rapaport The Art of Molecular Dynamics Simulation. 2004 ? Daan Frenkel, B. Smit Understanding Molecular Simulation 2001 Websites

? ?

Molecular Dynamics Tutorial at EMBnet

http://www.ch.embnet.org/MD_tutorial/index.html http://www.ks.uiuc.edu/

Theoretical and Computational Biophysics Group at UIUC (home of VMD and NAMD)

赞助商链接

相关文章:

- Typhoon 8600激光扫描成像仪
- Typhoon 8600 激光扫描成像仪 产品简介: Typhoon 8600 是
*molecular**dynamics*公司产品,由仪器、仪 器控制软件和其他多种附件组成,配有苹果 Power Mac G4 电脑...

- AutoDock&ADT教程_图文
- (GA) 遗传算法 Tabu Search 禁忌搜索
*Molecular**dynamics*(MD) 分子动力学 Energy minimization 能量最小化*Molecular**dynamics*(MD) 分子模拟 Flexible-Protein ...

- NAMD入门教程(一)
- 不仅如此,利用 NAMD 还可以进行极具特色的 IMD(Interactive
*Molecular**Dynamics*, 交互式分子动力学模拟)和 SMD(Steered*Molecular**Dynamics*,可控式分子动力学模拟) 。...

- 入射粒子能量对沉积铝膜形状影响的二维截面模拟
- 【KEY WORDS】 :
*molecular**dynamics*; deposition; groove 一、引言铝膜的淀积主要有两种方法, 真空蒸发法和溅射法。 其中溅射法由于较高的入射粒子能 量而应用...

- A COURSE ON THEORY OF CHARGED POLYMERS 课堂讲义 Ref...
- Strongly charged flexible polyelectrolytes in poor solvents:
*Molecular**dynamics*simulations with explicit solvent. Journal of Chemical Physics 118, 6634-6647 (...

- 多尺度方法综述
- BD法(2004,Xiao S P and Belytschko T 2004 A bridging domain method for coupling continua with
*molecular**dynamics*Comput. Methods Appl. Mech. Eng. 193 ...

- AutoDock中文教程_图文
- (GA) 遗传算法 Tabu Search 禁忌搜索
*Molecular**dynamics*(MD) 分子动力学 Energy minimization 能量最小化*Molecular**dynamics*(MD) 分子模拟 Flexible-Protein ...

- NAMD入门教程(一)_图文
- 分子动力学模拟概论 分子动力学模拟(
*Molecular**Dynamics*Simulation)是指利用计算机软件,根据牛 顿力学的基本原理, 模拟大分子的相互作用和运动变化的研究方法。 生命...

- 分子动力学模拟及其在材料科学中的研究进展
- 基于QM的分子动力学称之为现代 分子动力学(ModernMD), 也称之为从头分子动力学(Abinitio
*moleculardynamics*,AIMD)。密 度泛函分子动力学(DFMD)和第一原理分子动力...

- 甘油水溶液氢键特性的分子动力学模拟3
- 分子动力学模拟(
*molecular**dynamics*simulation, MDS)方法已经被广 泛用来研究物质和体系的微观结构及其氢键特性。甘油分子小且柔性强,易于被 玻璃化,被认为是最有...

更多相关标签:

- Introduction+to+Molecular+Dynamics+Simulation
- Ultra-large Scale Molecular Dynamics Simulation for
- Molecular dynamics simulation of a thin water layer evaporation and evaporation coefficient
- The ReaxFF Polarizable Reactive Force Fields for Molecular Dynamics Simulation of Ferroelec
- Molecular Dynamics Simulation of Solvent-Polymer Interdiffusion. I. Fickian diffusion
- Molecular Dynamics Simulation of the Structure I Empty Gas Hydrate
- Molecular dynamics simulation of polymer insertion into lipid bilayers
- Molecular dynamics simulation of multiphase flow on micro-scale using CUDA
- Combining Molecular Dynamics with Lattice-Boltzmann A Hybrid Method for the Simulation of (
- Molecular dynamics simulation high entropy films
- Molecular Dynamics simulation
- Molecular Dynamics Simulations
- Introduction to Molecular Dynamics Simulation-NIC
- Molecular dynamics simulations_05
- Molecular dynamics simulations_04